分数不是从princomp生成的，无法生成biplot

Question

当我想使用非princomp本身生成的协方差或相关矩阵时，我遇到biplot的问题，特别是princomp。为简单起见，我将使用比我正在处理的数据集小得多的数据集。

cr <- cov.wt(USArrests)
biplot(princomp(data = USArrests, covmat = cr))

给我错误

Error in biplot.princomp(princomp(data = USArrests, covmat = cr)) : 
  object 'princomp(data = USArrests, covmat = cr)' has no scores

看起来像这里简单的事情，但谷歌到目前为止一直没有用。

Answer 1

data中的princomp参数只能由＆＃34; S3方法用于类＆＃39;＆＃39;＆＃34;。因此，您需要通过以下两种方式之一指定princomp来电：

cr <- cov.wt(USArrests)
pr1 <- princomp(x = USArrests, covmat = cr)
pr2 <- princomp(formula = ~ ., data = USArrests, covmat = cr)

biplot(pr1)
biplot(pr2)