Plot a ranef.lme Object
## S3 method for class 'ranef.lme' plot(x, form = NULL, omitFixed = TRUE, level = Q, grid = TRUE, control, xlab, ylab, strip, ...)
x |
an object inheriting from class
|
form |
an optional formula specifying the desired type of plot.
See also ‘Details:’. |
omitFixed |
an optional logical value indicating whether
columns with values that are constant across groups should be
omitted. Default is |
level |
an optional integer value giving the level of grouping
to be used for |
grid |
an optional logical value indicating whether a grid should
be added to plot. Only applies to plots associated with two-sided
formulas in |
control |
an optional list with control values for the
plot, when |
xlab, ylab |
axis labels, each with a sensible default. |
strip |
|
... |
optional arguments passed to the Trellis |
If form
is missing, or is given as a one-sided formula, a
Trellis dot-plot (via dotplot()
from pkg
lattice) of the random effects is generated, with a different
panel for each random effect (coefficient). Rows in the dot-plot are
determined by the form
argument (if not missing) or by the row
names of the random effects (coefficients). Single factors (~g
)
or crossed factors (~g1*g2
) are allowed. For a single factor,
its levels determine the dot-plot rows (with possibly multiple dots
per row); otherwise, if form
specifies a crossing of factors, the dot-plot rows are determined by
all combinations of the levels of the individual factors in the
formula.
If form
is a two-sided formula, the left hand side must
be a single random effect (coefficient) and the right hand side is
formed by covariates in x
separated by +
.
An xyplot()
Trellis display is generated,
with a different panel for each variable listed in the right hand side
of form
. Scatter plots are generated for numeric variables and
boxplots are generated for categorical (factor
or
ordered
) variables.
a Trellis plot of the estimated random-effects (coefficients) versus covariates, or groups.
José Pinheiro and Douglas Bates bates@stat.wisc.edu
fm1 <- lme(distance ~ age, Orthodont, random = ~ age | Subject) plot(ranef(fm1)) fm1RE <- ranef(fm1, aug = TRUE) plot(fm1RE, form = ~ Sex) plot(fm1RE, form = age ~ Sex) # "connected" boxplots
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.