我正在使用该psych包进行因子分析。我想在fa()对象中或使用fa.diagram().
例如,玩具数据:
require(psych)
n <- 100
choices <- 1:5
df <- data.frame(a=sample(choices, replace=TRUE, size=n),
b=sample(choices, replace=TRUE, size=n),
c=sample(choices, replace=TRUE, size=n),
d=sample(choices, replace=TRUE, size=n))
model <- fa(df, nfactors=2, fm="pa", rotate="promax")
model
Factor Analysis using method = pa
Call: fa(r = df, nfactors = 2, rotate = "promax", fm = "pa")
Standardized loadings (pattern matrix) based upon correlation matrix
PA1 PA2 h2 u2 com
a 0.45 -0.49 0.47 0.53 2.0
b 0.22 0.36 0.17 0.83 1.6
c -0.02 0.20 0.04 0.96 1.0
d 0.66 0.07 0.43 0.57 1.0
Run Code Online (Sandbox Code Playgroud)
我想改变PA1和PA2对FactorA和FactorB,要么通过改变model物体本身,或者在输出调节标签fa.diagram():
的文档fa.diagram有一个labels论点,但没有例子,到目前为止我所做的实验并没有取得成果。非常感谢任何帮助!
随着str(model)我发现$loadings属性,fa.diagram()用来渲染图。修改colnames()的model$loadings成功了。
colnames(model$loadings) <- c("FactorA", "FactorB")
fa.diagram(model)
Run Code Online (Sandbox Code Playgroud)