d1 <- scan("dat_acc_TCGH16_EW2.csv",skip=3) d2 <- scan("dat_acc_TCGH16_NS2.csv",skip=3) d3 <- scan("dat_acc_TCGH16_UD2.csv",skip=3) nn <- length(d1) tt <- seq(1,nn,by=1) tt <- tt*0.01 pltdat<-function(tt,dd,str){ plot( tt,dd, type="l", xlab="Time (sec)", ylab="Acceleration (gal)", main=str, xlim=c(0,300), ylim=c(-max(abs(dd)),max(abs(dd))), axes=FALSE ) c1 <- seq(0,300,by=30) c2 <- c(round(-max(abs(dd))),0,round(max(abs(dd)))) axis(1,c1) axis(2,c2) } par(mfrow=c(3,1)) dd <-d1; pltdat(tt,dd,"TCGH16_EW2") dd <-d2; pltdat(tt,dd,"TCGH16_NS2") dd <-d3; pltdat(tt,dd,"TCGH16_UD2")