set.seed(98765)
require(mclust)
iriss <- iris[sample(150,20),-5]
irisBIC <- mclustBIC(iriss,G=5,modelNames="VEV")
siris <- summary(irisBIC,iriss)
emp <- extract.mixturepars(siris,2)
emp$pro
round(emp$mean,digits=1)
emp$variance$modelName
round(emp$variance$scale,digits=2)
Run the code above in your browser using DataLab