psych: principal - 载荷分量

Question

我的问题与 psych 包中的 principal() 函数有关。

set.seed(0)
x <- replicate(8, rnorm(10))
pca.x <- principal(x, nf=4, rotate="varimax")

我知道如果我想查看负载 table，我可以使用 loading.x <-loadings(pca.x)，然后我将得到以下结果。

> loading.x
Loadings:
     RC1    RC3    RC4    RC2   
[1,]        -0.892 -0.205  0.123
[2,]  0.154  0.158  0.909       
[3,] -0.660  0.255 -0.249  0.392
[4,] -0.352  0.412  0.614 -0.481
[5,]  0.950 -0.208  0.117       
[6,] -0.302  0.111         0.860
[7,]  0.852        -0.195 -0.358
[8,] -0.109  0.903         0.265

                 RC1   RC3   RC4   RC2
SS loadings    2.323 1.934 1.373 1.342
Proportion Var 0.290 0.242 0.172 0.168
Cumulative Var 0.290 0.532 0.704 0.871

我的第一个困惑是加载对象。从技术上讲，它是一个矩阵，但是看它的维度，它是8 * 4，这意味着不包括下半部分。

基本上，我想要实现的是单独提取这部分：

                 RC1   RC3   RC4   RC2
SS loadings    2.323 1.934 1.373 1.342
Proportion Var 0.290 0.242 0.172 0.168
Cumulative Var 0.290 0.532 0.704 0.871

将它放在 data.frame 或矩阵中，而不是在控制台中查看它。 post Extracting output from principal function in psych package as a data frame 中似乎有 William Revelle 的回答。能够单独提取下半部分，但 print 函数仍然提供了全部内容。

其实我也很好奇开发者是如何构造一个loading对象的（看源码也搞不懂）。另外，我需要的部分在 'pca.x' 列表的其他地方找不到，至少在格式化的 table 中找不到。我在 mac 和 psych 1.5.1.

上使用 Rstudio 版本 0.98.1102、R 3.1.2

提前致谢！

Answer 1

这是部分回答，但由于是我的包裹，我会给出更完整的回答。

PCA 或 FA 因子载荷 tables 的摘要 table 在打印函数中计算。它被返回（通过打印无形地）。但是，它可用作 Vaccounted 对象。

即PCA 或 FA 输出的摘要 table

set.seed(0)
x <- replicate(8, rnorm(10))
pca.x <- principal(x, nf=4, rotate="varimax")
p <- print(pca.x)

round(p$Vaccounted,2)   #shows the summary of the loadings table
                       PC1  PC3  PC4  PC2
SS loadings           2.32 1.93 1.37 1.34
Proportion Var        0.29 0.24 0.17 0.17
Cumulative Var        0.29 0.53 0.70 0.87
Proportion Explained  0.33 0.28 0.20 0.19
Cumulative Proportion 0.33 0.61 0.81 1.00

这也适用于 fa 功能。

psych: principal - 载荷分量

psych: principal - loadings components

r

principal

pca

psych