Return to etf.S CVS log | Up to [ELWIX - Embedded LightWeight unIX -] / embedaddon / ntp / scripts / stats |
1.1 ! misho 1: options(digits=4) ! 2: file2 <- "etf_summary" ! 3: etf <- scan(file1, list(day=0, sec=0, offset=0, stab=0)) ! 4: r <- lsfit(etf$sec, etf$offset) ! 5: count<-length(etf$sec) ! 6: mean<-r$coef[[1]] ! 7: std<-sqrt(var(r$residuals)) ! 8: slope<-r$coef[[2]] * 1000 ! 9: cat("\n", file=file2 , append=TRUE, fill=FALSE, sep="") ! 10: cat(file1, "\n", file=file2, append=TRUE, fill=FALSE, sep="") ! 11: cat("etf1 ", count, ", T ", mean, " ns, R ", slope, " ps/s, std ", std, " us\n", file=file2, append=TRUE, fill=FALSE, sep="") ! 12: str <- paste("eps/", file1, ".eps", sep="") ! 13: postscript(str, , , , 5, pointsize=18) ! 14: par(mgp=c(1, 0, 0), tck=0.03, mar=c(2, 2, 1, 1)) ! 15: plot(etf$sec, etf$offset, type="l", xlab=paste("MJD", etf$day, "Time (s)"), ylab="External Offset (ns)", ylim=c(-400, 400))