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错误”。我是R的新手，如果有人可以解释一下，我会很感激。

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

找到该代码段

非常感谢！

Answer 1

您找到的代码段不起作用，因为没有声明＆＃34; survey.prcomp＆＃34;宾语。＆＃34; TOP4＆＃34;也不见了。我认为作者错过了这一行：

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

还有这一个：

top4 <- list()

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

进行此操作

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

Answer 2

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

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

在R - PCA旋转组件中打印

2 个答案: