conc <- rep(1:3,each=2,times=6) day <- rep(1:3,each=12) run <- rep(1:6,each=6) glucose <- c(41.2,42.6,135.7,136.8,163.2,163.3, 41.2,41.4,143.0,143.3,181.4,180.3, 39.8,40.3,132.4,130.3,173.6,173.9, 41.5,43.0,134.4,130.0,174.9,175.6, 41.9,42.7,137.4,135.2,166.6,165.5, 45.5,44.7,141.1,139.1,175.0,172.0) spectro <- data.frame(conc,day,run,glucose) attach(spectro) conc <- factor(conc) day <- factor(day) run <- factor(run) glucose.mod1 <- aov(glucose ~ conc + day + conc:day + day/run + conc:day/run) summary(glucose.mod1) glucose.mod2 <- aov(glucose ~ conc + day + conc:day + day/run + Error(conc:day/run)) summary(glucose.mod2) library(nlme) glucose.mod3 <- lme(fixed=glucose~conc,random=~1|conc/day/run) summary(glucose.mod3) anova(glucose.mod3) library(lme4) glucose.mod4 <- lmer(glucose~conc+(1|conc:day)+(1|run:day)+(1|day/conc:run)) summary(glucose.mod4) anova(glucose.mod4)