The Mississippi
data frame has 37 rows and 3 columns.
This data frame contains the following columns:
an ordered factor with levels
3
< 5
< 2
< 1
< 4
< 6
a numeric vector
a factor with levels
1
2
3
Littel, R. C., Milliken, G. A., Stroup, W. W., and Wolfinger, R. D. (1996), SAS System for Mixed Models, SAS Institute (Data Set 4.2).
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19  str(Mississippi)
if (require("lattice", quietly = TRUE, character = TRUE)) {
dotplot(drop(influent:Type) ~ y, groups = Type, Mississippi)
}
if (require("lme4", quietly = TRUE, character = TRUE)) {
options(contrasts = c(unordered = "contr.SAS", ordered = "contr.poly"))
## compare with output 4.1, p. 142
print(fm1Miss < lmer(y ~ 1 + (1influent), Mississippi))
## compare with output 4.2, p. 143
print(fm1MLMiss < update(fm1Miss, REML=FALSE))
## BLUP's of random effects on p. 142
ranef(fm1Miss)
## BLUP's of random effects on p. 144
print(ranef(fm1MLMiss))
#intervals(fm1Miss) # interval estimates of variance components
## compare to output 4.8 and 4.9, pp. 150152
print(fm2Miss < lmer(y ~ Type+(1influent), Mississippi, REML=TRUE))
print(anova(fm2Miss))
}

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.