提取 psych 包 for.poly 函数中因子的相关矩阵答案

【问题标题】：extract the correlation matrix for the factors in the psych package's fa.poly function提取 psych 包 for.poly 函数中因子的相关矩阵
【发布时间】：2012-07-18 23:33:23
【问题描述】：

我在caracal's great example 工作，对二分数据进行因子分析，现在我正努力从psych 包的fa.poly 函数生成的对象中提取因子。

谁能帮我从fa.poly 对象中提取因子（并查看相关性）？

有关工作示例，请参阅caracal's example。

【问题讨论】：

标签： r data-structures extract

【解决方案1】：

在本例中，您创建一个对象：

faPCdirect <- fa.poly(XdiNum, nfactors=2, rotate="varimax")    # polychoric FA

所以在 faPCdirect 的某个地方有你想要的。我建议使用str() 来检查faPCdirect 的结构

> str(faPCdirect)

List of 5
 $ fa   :List of 34
  ..$ residual    : num [1:6, 1:6] 4.79e-01 7.78e-02 -2.97e-0...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:6] "X1" "X2" "X3" "X4" ...
  .. .. ..$ : chr [1:6] "X1" "X2" "X3" "X4" ...
  ..$ dof         : num 4
  ..$ fit      
      ...skip stuff....
  ..$ BIC         : num 4.11
  ..$ r.scores    : num [1:2, 1:2] 1 0.0508 0.0508 1
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:2] "MR2" "MR1"
  .. .. ..$ : chr [1:2] "MR2" "MR1"
  ..$ R2          : Named num [1:2] 0.709 0.989
  .. ..- attr(*, "names")= chr [1:2] "MR2" "MR1"
  ..$ valid       : num [1:2] 0.819 0.987
  ..$ score.cor   : num [1:2, 1:2] 1 0.212 0.212 1

所以这表示这个对象是一个由五个组成的列表，第一个元素称为fa，其中包含一个称为score.cor 的元素，它是一个 2x2 矩阵。我认为你想要的是非对角线。

> faPCdirect$fa$score.cor
          [,1]      [,2]
[1,] 1.0000000 0.2117457
[2,] 0.2117457 1.0000000

【讨论】：

您知道这种相关性是如何计算的吗？是皮尔逊相关吗？