在 R 中使用随机投影的奇异协方差和 NA 的 z 值

Question

我在以下代码中得到 NA 的 z 值。我知道这是由于奇异协方差。我已经运行对 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)

谢谢，

基因

Answer 1

看看这个

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")

在 R 中使用随机投影的奇异协方差和 NA 的 z 值

Singular Covariance and NA's for z value using Random Projection in R

r

machine-learning