heplot3d(mod, ...)
## S3 method for class 'mlm':
heplot3d(mod, terms, hypotheses, term.labels = TRUE,
hyp.labels = TRUE, variables = 1:3, error.ellipsoid = !add,
factor.means = !add, grand.mean = !add, remove.intercept = TRUE,
type = c("II", "III", "2", "3"), manova, size = c("evidence", "effect.size"),
level = 0.68, alpha = 0.05, segments = 40, col = palette()[-1], lwd=c(1, 4),
bg.col = c("white", "black"), fogtype = c("none", "exp2", "linear", "exp"),
fov = 30, offset = 0.01, xlab, ylab, zlab, add = FALSE, verbose = FALSE,
warn.rank = FALSE, ...)"mlm".TRUE, defaults to all terms;
if FALSE, no terms are plotted.linear.hypothesis function in the car package;
the listTRUE (the default) the names of the terms are used;
if FALSE, term labels are not plotted.TRUE (the default) the names of components of the list of
hypotheses are used; if FALSE, hypothesis labels are not plotted.1:3.TRUE, plot the error ellipsoid; defaults to TRUE,
if the argument add is FALSE (see below).TRUE or FALSE; defaults to TRUE,
if the argument add is FALSE (see below).TRUE, plot the centroid for all of the data;
defaults to TRUE,
if the argument add is FALSE (see below).TRUE (the default), do not plot the
ellipsoid for the intercept even if it is in the MANOVA table."II", "III", "2", or "3", where "II" is
the default (and "2" is a synomym).Anova.mlm object for the model; if absent a
MANOVA is computed. Specifying the argument can therefore save
computation in repeated calls."evidence", the default, the scaling is done so that
a ``significant'' hypothesis elliposid extends outside of the error
ellipsoid; if "effect.size0.68.size="evidence", then the
hypothesis ellipsoid is scaled so that it just touches the error ellipsoid at the
specified alpha level; a larger hypothesis elllipsoid ther40.c(1, 4)."white" or "black", defaulting
to "white"."none". See bg.viewpoint.0.01.TRUE, add to the current plot; the default is FALSE.
If TRUE, the error ellipsoid is not plottedTRUE, print the MANOVA table and details of hypothesis
tests; the default is FALSE.TRUE, do not suppress warnings about the rank of the
hypothesis matrix when the ellipsoid collapses to an ellipse or line;
the default is FALSE.heplot3d returns nothing useful --- it is used for the side effect of
drawing or adding to an rgl graph. The function leaves an object named
.frame in the global environment, containing the rgl object IDs for
the axes, axis labels, and bounding box; these are deleted and the axes, etc.
redrawn if the plot is added to.Anova, linear.hypothesis,
heplot, pairs.mlm, rgl.# Soils data, from car package
soils.mod <- lm(cbind(pH,N,Dens,P,Ca,Mg,K,Na,Conduc) ~ Block + Contour*Depth, data=Soils)
Anova(soils.mod)
heplot(soils.mod, variables=c("Ca", "Mg"))
pairs(soils.mod, terms="Depth", variables=c("pH", "N", "P", "Ca", "Mg"))
heplot3d(soils.mod, variables=c("Mg", "Ca", "Na"))Run the code above in your browser using DataLab