在以下代码中，我的z值出现了NA值。我知道这是由于奇异协方差引起的。我在iris数据集上运行了类似的测试，并没有遇到同样的错误。是我做错了什么，还是这是数据的特性所致？

R语言代码如下：

install.packages('mclust')library('mclust')mydata <- read.table("http://archive.ics.uci.edu/ml/machine-learning-databases/wine-quality/winequality-red.csv", sep=";", header=TRUE);est <- meVVV(mydata[,-12], unmap(mydata[,12]))randProj(mydata[,-12], seeds=200, parameters = variance, z = est$z, truth = iris[,5],         what = "errors", identify = TRUE, scale = TRUE)

谢谢，

@[隐藏人名]

回答：

请看这里

table(mydata[, 12]) # 3   4   5   6   7   8 # 10  53 681 638 199  18

这意味着你在使用一个包含10个观测值的聚类，而你的模型中有11个特征（变量），因此会出现警告信息和NA值。

如果你想指定模型为VVV，请使用以下代码：

 fm <- Mclust(mydata[, -12],  modelNames="VVV")