在 R 中打印 - PCA 旋转组件

Question

我在 R 中做了一个 PCA，我正在尝试打印旋转分量。我几乎是在试图理解我在网上找到的一个片段，如果有人能帮助我，我将不胜感激。请看下面我在网上找到的片段：

    require(stats) 
    prcomp(top2, scale=TRUE) 
    summary(prcomp(top2, scale=TRUE)) 

    for (i in 1:15) { 
    top4[[i]] <- sort(survey.prcomp$rotation[,i], decreasing=TRUE)[1:4]} 
    top4

我正在尝试打印前 15 个主要成分，我得到了 "top4 object not found error"。我是 R 的新手，如果有人能解释一下，我将不胜感激。

可以在 https://www.casact.org/pubs/forum/10spforum/Francis_Flynn.pdf

找到代码片段

非常感谢！

Answer 1

您找到的代码段无效，因为没有声明的 "survey.prcomp" 对象。 "top4" 也不见了。我假设作者错过了这一行：

survey.prcomp <- prcomp(top2, scale=TRUE)

还有这个：

top4 <- list()

然后，如果您的目标是获得前 15 个旋转向量，您可以使用 survey.prcomp$rotation[,1:15]

您粘贴的代码段有所不同。它 returns，对于 15 个主要主成分中的每一个，对载荷（旋转）影响最大的前 4 个变量。

Answer 2

在您粘贴的代码段中，访问了一系列变量但从未分配，因此出现错误。

top2、survey.prcomp 和 top4 从未分配，在您所附的文档中，作者似乎省略了这些行。

在 R 中打印 - PCA 旋转组件

Printing in R - PCA rotation components

r

pca