1options(digits=4) 2loop <- scan(file1, list(day=0, sec=0, offset=0, freq=0, tc=0)) 3loop$offset <- loop$offset * 1e6 4str <- paste("eps/", file1, ".eps", sep="") 5postscript(str, , , , 5, pointsize=18) 6par(mgp=c(1, 0, 0), tck=0.03, mar=c(2, 2, 1, 1)) 7plot(loop$sec, loop$offset, type="l", xlab=paste("MJD", loop$day, "Time (s)"), ylab="PLL Offset (us)", ylim=c(-400, 400)) 8