#pdf("kuehl_07_01.pdf") day <- rep(1:4,each=8) machine <- rep(1:4,each=2,times=4) serum <- c(142.3,144.0,148.6,146.9,142.9,147.4,133.8,133.2, 134.9,146.3,145.2,146.3,125.9,127.6,108.9,107.5, 148.6,156.5,148.6,153.1,135.5,138.9,132.1,149.7, 152.0,151.4,149.7,152.0,142.9,142.3,141.7,141.2) triglyc <- data.frame(day,machine,serum) attach(triglyc) day <- factor(day) machine <- factor(machine) triglyc.mod1 <- aov(serum ~ day*machine) summary(triglyc.mod1) library(lme4) triglyc.mod2 <- lmer(serum~1 + (1|day) + (1|machine) + (1|day:machine)) summary(triglyc.mod2) #dev.off()