# REML analysis as in O'Neill, p. 130.
dat <- steel.soybeanmet
require(lme4)
m1 <- lmer(yield ~ loc + (1|gen) + (1|gen:loc) + (1|loc:block), data=dat)
# Are gen:loc effects significant? No.
m2 <- lmer(yield ~ loc + (1|gen) + (1|loc:block), data=dat)
anova(m1,m2)
# h^2 deviations
ranef(m2)$genRun the code above in your browser using DataLab