データ
library(psych) library(GPArotation) data(bfi) d1 <- bfi[1:25]
res01 <- fa(d1, nfactors = 5, fm = "ml", rotate = "promax", scores=TRUE) fa.diagram (res01, cut=0.3, simple=FALSE, sort=TRUE, digits=3)
オプション
- cut: abs(loading)>cutの負荷が表示される
- simple: 項目ごとの最大荷重のみを表示
- digits: : 係数の桁数
- sort: 因子負荷量をソートしてから図を表示する
- graphviz: 出力にRgraphvizを使うか
過去のpsychパッケージに関するエントリ
https://ides.hatenablog.com/entry/2019/03/19/093726 https://ides.hatenablog.com/entry/2019/04/08/171145 https://ides.hatenablog.com/entry/2019/07/14/160817 https://ides.hatenablog.com/entry/2019/07/31/031654 https://ides.hatenablog.com/entry/2019/08/02/052958