我正在对我拥有的数据集执行主成分分析。这是我迄今为止使用的代码:
pca1<-prcomp(data1,center=TRUE,scale=TRUE)
summary(pca1)
pca1$rotation
pca1$sdev^2
library("devtools")
library("ggbiplot")
plot1 <- ggbiplot(pca1,obs.scale=1, var.scale=1, groups=data2, ellipse=TRUE)
我得到错误:
Error in `[.data.frame`(x, order(x, na.last = na.last, decreasing = decreasing)) : undefined columns selected
供您参考,“data1”是一个只有数值变量的数据框,“data2”包含一列分类变量,我试图通过这些变量对数值变量进行分组。两个数据集的行数相同。我之前在非常相似的数据集上运行了这个确切的代码,没有任何问题。
关于可能是什么问题的任何建议?