r - 绘制超过 2 个因子

Question

假设我进行了因子分析并得到了 5 个相关因子。现在，我想以图形方式表示这些因素对变量的加载。谁能告诉我该怎么做。我可以使用 2 个因素。但是当因子数超过2时就无法做到。

“现代应用统计与 S”中给出了 2 因子绘图，图 11.13。我想创建类似的图表，但有两个以上的因素。请找到上面提到的图的快照:

X 和 y 轴是 2 个因素。

问候，阿里

Answer 1

注意:不是您要找的答案，也可能不正确，这是我的主观想法。

我认为您遇到了在二维屏幕/纸上绘制多个维度的问题。我会说绘制更多因子或 PC 的载荷是没有意义的，但如果你真的坚持:显示前两个(基于特征值)或只创建 2 个因子。或者您也可以通过其他方法(例如 MDS)来减小尺寸。

在 3 维图表中显示 3 个因素的负载将很难清楚，更不用说考虑更多因素了。

---

更新:我有一个梦想是尝试更贴近主题:)

正如@joran 指出的那样，您可以轻松地显示每个 pairs 因素的预测(我在这里不处理轮换):

f <- factanal(mtcars, factors=3)
pairs(f$loadings)

这样你可以显示更多的因素，也可以调整情节，例如:

f <- factanal(mtcars, factors=5)
pairs(f$loadings, col=1:ncol(mtcars), upper.panel=NULL, main="Factor loadings")
par(xpd=TRUE) 
legend('topright', bty='n', pch='o', col=1:ncol(mtcars), attr(f$loadings, 'dimnames')[[1]], title="Variables")

当然，您也可以通过自定义下部三角形，或在上部三角形中显示并在右侧/下方附加图例等方式添加旋转矢量。

或者，如果您的因子不超过 3 个，则只需将变量指向 3D 散点图:

library(scatterplot3d)
f <- factanal(mtcars, factors=3)
scatterplot3d(as.data.frame(unclass(f$loadings)), main="3D factor loadings", color=1:ncol(mtcars), pch=20)

注意:变量名称不应该作为标签放在图上，但在我的拙见中可能会成为一个独特的传说，特别是对于 3D 图。