###################################################
### chunk number 1:
###################################################
 options(width=60)


###################################################
### chunk number 2:
###################################################
 library(MASS)     #needed for whiteside data set

 # see a description of data

 ?whiteside       #documentation
 str(whiteside)   #a description
 head(whiteside)  #first few observations


###################################################
### chunk number 3:
###################################################
 by(whiteside, whiteside$Insul, FUN=summary)


###################################################
### chunk number 4:
###################################################
plot(Gas ~ Temp, data=whiteside, col=Insul)


###################################################
### chunk number 5:
###################################################
 fit <- lm(Gas ~ Temp, data=whiteside)
 fit


###################################################
### chunk number 6:
###################################################
 plot(Gas ~ Temp, data=whiteside, col=Insul)
 abline(fit$coef)


###################################################
### chunk number 7:
###################################################
 plot(fit$fitted, fit$resid, col=whiteside$Insul)
 abline(h=0)

###################################################
### verbatim chunk for lattice
###################################################

 library(lattice)
 xyplot(Gas ~ Temp | Insul, data=whiteside)

###################################################
### verbatim chunk for lattice
###################################################
panel <- function(x, y) {
    panel.xyplot(x, y)
    panel.lmline(x, y)
    }

xyplot(Gas ~ Temp | Insul, data=whiteside, panel=panel)

###################################################
### chunk number 8:
###################################################
 fitBefore <- lm(Gas ~ Temp, data=whiteside, subset=(Insul=="Before"))
 fitBefore
 summary(fitBefore)
 anova(fitBefore)


###################################################
### chunk number 9:
###################################################
 plot(Gas ~ Temp, data=whiteside, col=Insul)
 abline(fitBefore$coef)


###################################################
### chunk number 10:
###################################################
 fitBefore <- lm(Gas ~ Temp, data=whiteside, subset=(Insul=="Before"))

 fitAfter <- update(fitBefore, subset=(Insul=="After"))

 fitBefore
 fitAfter


###################################################
### chunk number 11:
###################################################
 plot(Gas ~ Temp, data=whiteside, col=Insul)
 abline(fitBefore$coef)
 abline(fitAfter$coef, col=2)
 legend("topright", inset=.01, c("Before","After"), col=1:2, lty=1)

###################################################
### more verbatim chunks
###################################################

 fit2 <- lm(Gas ~ Insul * Temp, data=whiteside)
 fit2
 summary(fit2)$coef

 fit3 <- lm(Gas ~ Insul+Temp, data=whiteside)
 anova(fit3, fit2)