day <- rep(1:3,each=6)
run <- rep(1:6,each=3)
glucose <- c(42.5,43.3,42.9,42.2,41.4,41.8,
             48.0,44.6,43.7,42.0,42.8,42.8,
             41.7,43.4,42.5,40.6,41.8,41.8)

clin <- data.frame(day,run,glucose)

attach(clin)

day <- factor(day)
run <- factor(run)

glucose.mod1 <- aov(glucose ~ day + day/run)

summary(glucose.mod1)

glucose.mod2 <- aov(glucose ~ day + Error(run))

summary(glucose.mod2)

library(nlme)
glucose.mod3 <- lme(fixed=glucose~1,random=~1|day/run)

summary(glucose.mod3)
intervals(glucose.mod3)


library(lme4)
glucose.mod4 <- lmer(glucose~1+(1|day/run))
summary(glucose.mod4)