Manly.var {ManlyMix} | R Documentation |
Variability assessment of Manly mixture model
Description
Runs the variability assessment for a Manly mixture model.
Usage
Manly.var(X, model = NULL, conf.CI = NULL)
Arguments
X |
dataset matrix (n x p) |
model |
Manly mixture model |
conf.CI |
confidence level, say 95 percent confidence |
Details
Returns the estimated variance-covariance matrix and confidence intervals for model parameter estimates.
Value
V |
variance-covariance matrix. |
CI |
confidence intervals for each parameter. |
See Also
Manly.EM
Examples
set.seed(123)
#Use iris dataset
K <- 3; p <- 4
X <- as.matrix(iris[,-5])
#Use k-means clustering result
#all skewness parameters set to be 0.1 as the initialization of the EM algorithm
id.km <- kmeans(X, K)$cluster
la <- matrix(0.1, K, p)
#Run the EM algorithm with Manly mixture model
M.EM <- Manly.EM(X, id.km, la)
# Run the variability assessment
Manly.var(X, M.EM, conf.CI = 0.95)
[Package ManlyMix version 0.1.15 Index]