## multiple groups, common covariance matrix
if (require ("mvtnorm")){
pcov <- pooled.cov (chondro, chondro$clusters)
rnd <- rmmvnorm (rep (10, 3), mean = pcov$mean, sigma = pcov$COV)
plot (rnd, col = rnd$.group)
}
Run the code above in your browser using DataLab