I am wondering how to make a silhouette plot on the results from Mclust packge, where measure of silhouette for each sample in each group is shown?

Using "Fcp" package I can check the average silhouette.

mb = Mclust(iris[,-5], 3)

cs = cluster.stats(dist(iris[1:4]), mb$classification)

#install.packages("fpc")
library("fpc")

cs[c("within.cluster.ss","avg.silwidth")]

Thanks a lot.