Learn R
lme4 (version 1.1-12)

plots.thpr: Mixed-Effects Profile Plots (Regular / Density / Pairs)

Description

Xyplot, Densityplot, and Pairs plot methods for a mixed-effects model profile.

xyplot() draws “zeta diagrams”, also visualizing confidence intervals and their asymmetry.

densityplot() draws the profile densities.

splom() draws profile pairs plots. Contours are for the marginal two-dimensional regions (i.e. using df = 2).

Usage

"xyplot"(x, data = NULL, levels = sqrt(qchisq(pmax.int(0, pmin.int(1, conf)), df = 1)), conf = c(50, 80, 90, 95, 99)/100, absVal = FALSE, scales=NULL, which = 1:nptot, ...)


"densityplot"(x, data, ...)


"splom"(x, data, levels = sqrt(qchisq(pmax.int(0, pmin.int(1, conf)), 2)), conf = c(50, 80, 90, 95, 99)/100,  which = 1:nptot, draw.lower = TRUE, draw.upper = TRUE, ...)

Arguments

x
a mixed-effects profile, i.e., of class "thpr", typically resulting from profile(fm) where fm is a fitted model from lmer (or its generalizations).
data
unused - only for compatibility with generic.
levels
the contour levels to be shown; usually derived from conf.
conf
numeric vector of confidence levels to be shown as contours.
absVal
logical indicating if abs(.)olute values should be plotted, often preferred for confidence interval visualization.
scales
plotting options to be passed to xyplot
which
integer or character vector indicating which parameters to profile: default is all parameters (see profile-methods for details).
draw.lower
(logical) draw lower-triangle (zeta scale) panels?
draw.upper
(logical) draw upper-triangle (standard dev/cor scale) panels?
...
further arguments passed to xyplot, densityplot, or splom from package lattice, respectively.

Value

xyplot:
a density plot, a "trellis" object (lattice package) which when print()ed produces plots on the current graphic device.
densityplot:
a density plot, a "trellis" object, see above.
splom:
a pairs plot, aka scatterplot matrix, a "trellis" object, see above.

See Also

profile, notably for an example.

Examples

Run this code
## see   example("profile.merMod")

Run the code above in your browser using DataCamp Workspace