pdf("kuehl_02_02.pdf")

trt <- c(rep(1,9),rep(2,6),rep(3,8))
tempdiff <- c(2.2,1.6,0.8,1.8,1.4,0.4,0.6,1.5,0.5,
0.3,0.0,0.6,0.0,-0.3,0.2,
0.1,0.1,0.2,-0.4,0.3,0.1,0.1,-0.5)

rabbits <- data.frame(trt,tempdiff)

attach(rabbits)

trt <- factor(trt)

rabbits

tapply(tempdiff,trt,mean)
tapply(tempdiff,trt,sd)

rabbits.mod1 <- lm(tempdiff~trt)

summary(rabbits.mod1)
anova(rabbits.mod1)
predict(rabbits.mod1)

dev.off()