当我执行 mclustModel 时,给我的输出与通常的输出不同。
通常输出如下:
$z[,1] [,2] [,3] [,4] [,5]
1 1 0 0 0 0
2 0 1 0 0 0
3 0 0 1 0 0
4 1 0 0 0 0
5 0 0 1 0 0
等等。
然而,这次输出是这样的:
$z [,1] [,2] [,3]
1 1.000000e+00 0.0000000 0
2 1.000000e+00 0.0000000 0
3 1.000000e+00 0.0000000 0
4 1.000000e+00 0.0000000 0
5 1.000000e+00 0.0000000 0
6 5.058816e-65 1.0000000 0
7 2.042567e-72 0.0000000 1
下面的每个数字%z代表一个人,而其他列的集群,在第二个输出中看起来像同一个人可能在两个不同的集群中。这是对的还是有些错误?
我使用了以下代码:
library(foreign)
library(mclust)
idd<- read.dta(data)
idd_ss <-subset(idd,select=c(variable1, variable2, variabel3, variable4, variable5))
iddBIC<- mclustBIC(idd_ss)
cl <- mclustModel(idd_ss,iddBIC)
cl