Runs the variability assessment for a Manly mixture model.

1 |

`X ` |
dataset matrix (n x p) |

`model ` |
Manly mixture model |

Returns the estimated variance-covariance matrix for model parameter estimates.

Manly.EM

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ```
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)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.