Mclust (5.4.7) 和随后的 predict.Mclust 给出的分类和不确定性有所不同。我不明白为什么？
library(mclust)
library(MASS)
set.seed(123)
c1<-mvrnorm(100,mu=c(0,0),Sigma=matrix(c(2,0,0,2),ncol=2))
c2<-mvrnorm(200,mu=c(3,3),Sigma=matrix(c(3,0,0,3),ncol=2))
d<-rbind(c1,c2)
m<-Mclust(d)
u1<-m$uncertainty
u2<-apply(predict.Mclust(m,d)$z,1,function(x){1-max(x)})
plot(u1,u2)
