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))