set.seed(98765)
{require(mclust)}
iriss <- iris[sample(150,20),-5]
irisBIC <- mclustBIC(iriss)
siris <- summary.mclustBIC(irisBIC,iriss)
round(mixdens(siris$modelName,iriss,siris$parameters),digits=2)
Run the code above in your browser using DataLab