# # Written by: # -- # John L. Weatherwax 2009-04-21 # # email: wax@alum.mit.edu # # Please send comments and especially bug reports to the # above email address. # #----- # Fort Collins snowfall #ftcollinssnow <- read.table("data/ftcollinssnow.txt", header=T) if(is.null(version$language) == FALSE) data(ftcollinssnow) attach(ftcollinssnow) # Duplicates Fig 1.6 "../EPS-figs/ftcollins1.eps" postscript("../../WriteUp/Graphics/Chapter2/ftcollins1.eps", onefile=FALSE, horizontal=FALSE) plot(Early,Late,xlim=c(0,60),ylim=c(0,60),bty="l") abline(h=mean(Late)) abline(lm(Late~Early),lty=2) dev.off() m1 <- lm(Late~Early) summary(m1) detach("ftcollinssnow")