library(dplyr)
library(ggplot2)
library(mclust)
set.seed(2016)
centers <- data.frame(cluster=factor(1:3), size=c(100, 150, 50),
x1=c(5, 0, -3), x2=c(-1, 1, -2))
points <- centers %>% group_by(cluster) %>%
do(data.frame(x1=rnorm(.$size[1], .$x1[1]),
x2=rnorm(.$size[1], .$x2[1]))) %>%
ungroup()
m = Mclust(points %>% dplyr::select(x1, x2))
tidy(m)
head(augment(m, points))
glance(m)
Run the code above in your browser using DataLab