因子数の求め方

記念すべき第2回目の記事です．
今回は，因子数の求め方です．

ご参考に〜^ ^

＜サンプルデータ＞
7個の質問に対して7件法で回答したサンプルデータがあります．
Excelに以下のように入力し，ファイル名を「sample2.xlsx」として保存します．

RstudioのConsole画面に以下のように入力します．

install.packages("readxl")
library(readxl)
x<-read_excel("sample2.xlsx")

1行目は，readxlパッケージ(1)のインストールをしています．
                readxlパッケージはexcelをRに読み込む関数です．
2行目は，readxlパッケージを読み込んでいます．
3行目は，xに「sample2.xlsx」のファイルデータを格納しています．

install.packages("psych")
library(psych)

1行目は，psychパッケージ(2)をインストールしています．
2行目は，psychパッケージを読み込んでいます．

VSS.scree(x)

スクリープロットを描写することができます．

この図をみると，因子数が2であることがわかります．

> corx <- cor(x)
> eigen(corx)$values
[1] 3.2879653 1.3244321 0.7527276 0.6224593 0.4322615 0.3624570 0.2176971
> eigen(corx)$values/sum(eigen(corx)$values)
[1] 0.46970933 0.18920459 0.10753251 0.08892276 0.06175165 0.05177957 0.03109959

2行目は固有値，4行目は寄与率を算出しています．
ガットマン基準では固有値が1以上を因子数決定の目安にしています．


そうそう，因子数を決定する際はいくつかの方法がありますからね〜．

文献
(1)Hadley Wickham and Jennifer Bryan (2019). readxl: Read Excel Files. R package version 1.3.1. https://CRAN.R-project.org/package=readxl
(2)Revelle, W. (2018) psych: Procedures for Personality and Psychological Research, Northwestern University, Evanston, Illinois, USA, https://CRAN.R-project.org/package=psych Version = 1.8.12.