source("http://www.uvm.edu/~rsingle/stat5210/data/scripts-5210.R")

#concentration (ug/dl) of a chemical recorded at 6 different times

 dat <- otherdata("lackfit.dat")
 
 mod1 <- lm(CONC ~ TIME, data=dat)
 anova(mod1)
 
 plot( CONC ~ TIME, data=dat) 
 
 pure.error(dat$CONC, dat$TIME) #SYNTAX: pure.error(data$y.var, data$x.var)
 
 multifig(2,2)
 plot( CONC ~ TIME, data=dat)
 plot((CONC)^(1/2) ~ TIME, data=dat)
 plot((CONC)^(1/4) ~ TIME, data=dat)
 plot(   log(CONC) ~ TIME, data=dat)
 multifig(1,1)

#Tukey's Rule of the Bulge --> transform Y down