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%。方法はパーセンタイルだな。