lme(formula, data, random, correlation, weights, subset,
method, na.action, control, model, x)
~
operator and the terms, separated by +
operators, on
the right.formula
, random
, correlation
, weights
, and
subset
. By default the variables are taken from the
environment from which ~x1+...+xn | g1/.../gm
, with x1+...+xn
specifying the model for the random effects and g1/.../gm
the
grouping structure (m
mcorStruct
object describing the
within-group correlation structure. See the documentation of
corClasses
for a description of the available corStruct
classes. Defaults to NULL
,
corvarFunc
object or one-sided formula
describing the within-group heteroscedasticity structure. If given as
a formula, it is used as the argument to varFixed
,
corresponding to fixed variance weights. See the dodata
that should be used in the fit. This can be a logical
vector, or a numeric vector indicating which observation numbers are
to be included, or a character vector of th"REML"
the model is fit by
maximizing the restricted log-likelihood. If "ML"
the
log-likelihood is maximized. Defaults to "REML"
.NA
s. The default action (na.fail
) causes
lme
to print an error message and terminate if there are any
incomplete observations.lmeControl
.
Defaults to an empty list.TRUE
the corresponding
components of the fit (the model frame, the model matrices)
are returned.data(bdf, package = "nlme")
fm <- lme(langPOST ~ IQ.ver.cen + avg.IQ.ver.cen, data = bdf,
random = ~ IQ.ver.cen | schoolNR)
summary(fm)
Run the code above in your browser using DataLab