mixproj {mclust}R Documentation

Displays one standard deviation of an MVN mixture classification.

Description

mixproj displays one standard deviation of an MVN mxiture classification along with data in for selected pairs of coordinates. Function mvn2plot is used to plot the confidence ellipses.

Usage

mixproj(data, ms, partition, scale = F, k = 15, ...)
mvn2plot(mu, sigma, k)

Arguments

data a matrix of onservations.
ms The result of an mstep calculation (a list consisting of mu and sigma).
partition A integer vector giving an initial classification for each observation.
dimens A vector of length two giving the two variables of the data to be plotted.
scale A logical variable telling whether or not the same scale should be used for both variables so as to preserve geometry. The default does not use the same scale.
k Number of subdivisions for plotting segments of ellipsoids. Default: 8.
... Further arguments to plot.

See Also

mstep, clpairs

Examples

data(iris)
cl <- mhclass(mhtree(iris[,1:4], modelid = "VVV"),3)
z <- me( iris[,1:4], modelid = "VVV", ctoz(cl))
pars <- mstep(iris[,1:4], modelid="VVV", z)
mixproj(iris[,1:4], ms=pars, partition=ztoc(z), dimens=c(1,2))