data("Nclus", package="flexmix")
plot(Nclus)
cl1 = stepFlexclust(Nclus, K=2:7, FUN=cclust)
cl1
plot(cl1)
getModel(cl1, 4)
opar=par("mfrow")
par(mfrow=c(2,2))
for(k in 3:6){
image(getModel(cl1, as.character(k)), data=Nclus)
title(main=paste(k, "clusters"))
}
par(opar)Run the code above in your browser using DataLab