densityplot produces plots of the densities. The
function automatically separates the observed and imputed
data. The functions extend the usual features of
## S3 method for class 'mids':
densityplot(x, data, na.groups = NULL, groups = NULL,
as.table = TRUE, plot.points = FALSE, theme = mice.theme(),
mayreplicate = TRUE, thicker = 2.5, allow.multiple = TRUE,
outer = TRUE,
drop.unused.levels = lattice.getOption("drop.unused.levels"),
panel = lattice.getOption("panel.densityplot"),
default.prepanel = lattice.getOption("prepanel.default.densityplot"), ...,
subscripts = TRUE, subset = TRUE)mids object, typically created by
mice() or mice.mids().is.na(x$data).
groups arguments
in na.groups
because it evaluates in the completed data
data.frame(complete(x, "long", inc=TRUE)) (as
usual), whereas na.groupsdensityplot
that signals whether the points should be plotted.mice.theme
produces a short list of default colors, line width, and
so on. The extensive list may be obtained from
trellis.par.get(). Global graphicaldensityplot. Multiplication
factor of the line width of the observed density.
thicker=1 uses the same thickness for the observed
and imputed data.xyplot.xyplot.xyplot.xyplot.xyplot.xyplot.xyplot.xyplot."trellis". The
update method can be
used to subsequently update components of the object, and
the print method
(usually called by default) will plot it on an appropriate
plotting device.na.groups may be used to specify
(combinations of) missingness in any of the variables. The
argument groups can be used to specify groups based
on the variable values themselves. Only one of both may be
active at the same time. When both are specified,
na.groups takes precedence over groups.Use the subset and na.groups together to
plots parts of the data. For example, select the first
imputed data set by by subset=.imp==1.
Graphical paramaters like col, pch and
cex can be specified in the arguments list to alter
the plotting symbols. If length(col)==2, the color
specification to define the observed and missing groups.
col[1] is the color of the 'observed' data,
col[2] is the color of the missing or imputed data.
A convenient color choice is col=mdc(1:2), a
transparent blue color for the observed data, and a
transparent red color for the imputed data. A good choice
is col=mdc(1:2), pch=20, cex=1.5. These choices can
be set for the duration of the session by running
mice.theme().
van Buuren S and Groothuis-Oudshoorn K (2011). mice:
Multivariate Imputation by Chained Equations in R.
Journal of Statistical Software, 45(3), 1-67.
mice, xyplot,
stripplot, bwplot,
Lattice for an overview of the package, as
well as densityplot,
panel.densityplot,
print.trellis,
trellis.par.setrequire(lattice)
imp <- mice(boys, maxit=1)
### density plot of head circumference per imputation
### blue is observed, red is imputed
densityplot(imp, ~hc|.imp)
### All combined in one panel.
densityplot(imp, ~hc)Run the code above in your browser using DataLab