#  first read the data directly from Dr. Oehlerts website

Resin1 <- read.table("http://www.stat.umn.edu/~gary/book/fcdae.data/exmpl3.2",header=TRUE)

Resin1

Resin1$temp <- as.factor(Resin1$temp)      # to make sure that temp is a factor, not numeric

boxplot(y ~ temp, data = Resin1, main = "Resin data")

Resin1.lm <- lm(y ~ temp, data=Resin1)

anova(Resin1.lm)
summary(Resin1.lm)

windows()             # opens a new graphics window under Windows operating system

par(mfrow=c(2,2))  
plot(Resin1.lm)
par(mfrow=c(1,1))