str(Bond)
options(contrasts = c(unordered = "contr.SAS", ordered = "contr.poly"))
if (require("lme4", quietly = TRUE, character = TRUE)) {
## compare with output 1.1 on p. 6
print(fm1Bond <- lmer(pressure ~ Metal + (1|Ingot), Bond))
print(anova(fm1Bond))
}
Run the code above in your browser using DataLab