「統計学関連なんでもあり」の過去ログ--- 048

No.22567　信頼性係数　KR-20　　【やまと】　2018/06/10(Sun) 18:02

よろしくお願いします。

RでKuder RichardsonのKR20を計算するためのスクリプトは公開されているのでしょうか。教えていただけますでしょうか。

No.22568　Re: 信頼性係数　KR-20　　【青木繁伸】　2018/06/10(Sun) 21:13

検索しましたか？

https://rdrr.io/github/cddesja/validateR/man/kr20.html
などもあります。

ここに書かれていることが難しければ，ともかく以下のようにすればよいでしょう。

devtools がまだインストールされていないなら，以下の 1 行を一度だけ

install.packages("devtools")

インストルされていた（インストールした）なら，以下の 2 行を一度だけlibrary(devtools)
install_github("cddesja/validateR")
そして，kr20 を使うときには R を立ち上げるたびに，以下の 1 行を 一度だけlibrary(validateR)
さて，準備が整ったら，データを以下のように準備する
例では 7 人の被検者に 6 項目の変数を 0/1 で準備x = matrix(c(
  1,1,0,1,1,0,
  1,0,0,1,1,1,
  0,1,1,0,1,1,
  0,0,0,1,1,0,
  0,1,0,1,0,1,
  0,0,0,1,1,0,
  0,0,1,0,0,0), byrow=TRUE, ncol=6)
その後，kr20(x)
とすれば，結果は>  kr20(x)
[1] 0.1382488
となる。

＝＝＝＝

なお，やっていることは簡単なので，自分で以下のような関数を書けば，面倒なことは一切なしで結果が得られる。
計算方法につては Wikipedia にさえ書かれているが，例えば，
https://www.radford.edu/.../KR21%20&%20KR20%20examples.xls
には，Excel での計算例が書かれているので，計算方法を理解すれば以下のような関数を書くのは容易である（R がそこそこできればではあるが）KR20 = function(x) {
  cm = colMeans(x)
  k = ncol(x)
  k/(k-1)*(1-sum(cm*(1-cm))/var(rowSums(x)))
}

> KR20(x)
[1] 0.1382488

No.22569　Re: 信頼性係数　KR-20　　【やまと】　2018/06/11(Mon) 07:38

青木先生

大変詳細な説明をいただきありがとうございます。

無事解決できました。