読者です 読者をやめる 読者になる 読者になる

草薙の研究ログ

英語教育関係。でも最近は統計(特にR)ネタが中心。

Cronbachのα係数の信頼区間

Cronbachのα係数の信頼区間って,あんま考えないのだけど,Rではいろんな関数でだしてくれる。で,気づいたんだけど,パッケージごとで値が違う。そりゃそうだ,これ,っていう式がないだろうな。ブートストラップもできるし。

 

いろんな論文があるらしい。

 

Feldt, L. S., Woodruff, D. J., & Salih, F. A. (1987). Statistical inferences for coefficient alpha. Applied Psychological Measurement, 11, 93-103.

 

これが一番よく知られるものってことでいいのかな。

他にも,こういうのもあるけど。

 

Feldt, L.S. (1965). The approximate sampling distribution of Kuder-Richardson reliability coefficient. Psychometrika, 30, 357-370.

Feldt, L.S. (1969). A test of the hypothesis that Cronbach's alpha or Kuder-Richardson coefficient is the same for two tests. Psychometrika, 34, 363-373.

 

psychometricパッケージのalpha.CI関数はこれで求めているようだ。

Feldt, L. S., Woodruff, D. J., & Salih, F. A. (1987). Statistical inferences for coefficient alpha. Applied Psychological Measurement, 11, 93-103.

 

多分一番使うであろう,psychのalpha関数の基本はこれ。

Duhachek, A. and Iacobucci, D. (2004). Alpha's standard error (ase): An accurate and precise confidence interval estimate. Journal of Applied Psychology, 89(5):792-808.

psychは最近bootstrapでも求めるようになったらしい。n.iterでBを指定する。方法はわかんない。パーセンタイル?

 

ltmパッケージのcronbach.alpha関数はブートストラップ。

cronbach.alpha(dat, CI = T, B = 1000)

とやるとB = 1000でブートストラップ。デフォルトは95%。方法はパーセンタイルだな。