# NOT RUN {
data(piepho.cocksfoot)
dat <- piepho.cocksfoot
dat$env <- factor(dat$year)
require(lattice)
densityplot(~date|env, data=dat, main="piepho.cocksfoot - heading date") # Gaussian, not gamma distn
# }
# NOT RUN {
require(lme4)
# A simple mixed model gives results similar to Piepho, Table 2.
m1 <- lmer(date ~ -1 + gen + (1|env), data=dat)
summary(m1)
## Fixed effects:
## Estimate Std. Error t value
## genG01 48.951 1.744 28.07
## genG02 42.270 1.755 24.09
## genG03 45.293 1.744 25.98
## genG04 45.467 1.744 26.08
## genG05 44.215 1.755 25.20
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab