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, dimens, scale = F, k = 15, xlim, ylim,
xlab, ylab, col=partition, pch=partition, ...)
mvn2plot(mu, sigma, k, add=T, col=NULL, xlim=NULL, ylim=NULL, xlab=NULL,
ylab=NULL, ...)
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. If not provided, a random selection of two variables is made.
|
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.
|
xlim, ylim, xlab, ylab, col, pch |
Standard arguments to plot with default values.
|
add |
Add ellipses from mvn2plot to the existing plot. |
... |
Further arguments to plot .
|
mu |
Cluster centers.
|
sigma |
Cluster covariance matrices. |
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))
### Plot only one cluster
mvn2plot(pars$mu[1:2,1], pars$sigma[1:2, 1:2, 1], 25, F, col="red")