sig <- diag (2)
cen <- rep (1, 2)
x <- rbind(MASS::mvrnorm(360, cen * 0, sig),
MASS::mvrnorm(540, cen * 5, sig * 6 - 2),
MASS::mvrnorm(100, cen * 2.5, sig * 50)
)
clus.1 <- tclust(x, k = 2, alpha=0.1, restr.fact=12)
clus.2 <- tclust(x, k = 3, alpha=0.1, restr.fact=1)
dsc.1 <- DiscrFact(clus.1)
summary(dsc.1)
dsc.2 <- DiscrFact(clus.2)
summary(dsc.2)
Run the code above in your browser using DataLab