感谢大家阅读我的问题。我试图理解使用R的主成分分析。
我使用这个例子，它出现在factanal（）文档中
首先，我创建我的变量：
v1 <- c(1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,4,5,6)
v2 <- c(1,2,1,1,1,1,2,1,2,1,3,4,3,3,3,4,6,5)
v3 <- c(3,3,3,3,3,1,1,1,1,1,1,1,1,1,1,5,4,6)
v4 <- c(3,3,4,3,3,1,1,2,1,1,1,1,2,1,1,5,6,4)
v5 <- c(1,1,1,1,1,3,3,3,3,3,1,1,1,1,1,6,4,5)
v6 <- c(1,1,1,2,1,3,3,3,4,3,1,1,1,2,1,6,5,4)
然后我使用数据框
df <- data.frame(v1,v2,v3,v4,v5,v6)
我使用factalal commmand
factanal(df,3)
我经常使用这个命令，总是根据我使用的书给出结果。但我试图在sjPlot库中使用sjt.pca命令，因为cronbach的alpha自动计算和html输出：
library(sjPlot)
sjt.pca(df)
这是我不理解唱歌（+或 - ）的地方，例如：
唱歌＆＃34; - ＆＃34;在组件3中，变量v3和v4使用sjt.pca()。这在factanal()中不存在，其中所有加载都是正数且不存在于：
library(psych)
psych::pca(r = df,nfactors = 3)
使用此库，所有加载均为正数。
第二个问题是：为什么他们有不同的订单？ factanal()中的因子1是sjt.pca()中的第2个成分。
据我所知，两者都默认使用varimax旋转。也许这些是基本问题，我不是专家。请帮忙。