#pdf("kuehl_07_02.pdf")

day <- rep(1:4,each=2,times=2)
method <- rep(1:2,each=8)
serum <- c(142.3,144.0,134.9,146.3,148.6,156.5,152.0,151.4,
           142.9,147.4,125.9,127.6,135.5,138.9,142.9,142.3)

chem_exp <- data.frame(day,method,serum)

attach(chem_exp)

day <- factor(day)
method <- factor(method)


chem.mod0 <- aov(serum ~ method + day + method:day)
summary(chem.mod0)

chem.mod1 <- aov(serum ~ method + day + Error(method:day))
summary(chem.mod1)

library(nlme)
chem.mod2 <- lme(fixed=serum~method,random=~1|day/method)

summary(chem.mod2)
intervals(chem.mod2)
anova(chem.mod2)

library(lme4)
chem.mod3 <- lmer(serum~method+(1|day)+(1|method:day))
summary(chem.mod3)
anova(chem.mod3)

#dev.off()