Manly.var: Variability assessment of Manly mixture model

In ManlyMix: Manly Mixture Modeling and Model-Based Clustering

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)

ManlyMix documentation built on April 29, 2023, 9:17 a.m.