myscratchpad

# 2件法の01データの場合はKuder-Richardsonの公式というのを使うらしい。全く知らなかった。知らないことが多いものである。日本語でぐぐると99件しかヒットしない。皆知らないのだろう。
 # 最近4件法や5件法は順序尺度として扱うべしとの意見があるが、信頼性を出すときはどうするんだろう。クロンバックでいいのだろうか
 # このサイトに少し解説。以下適当な訳と覚書
  # 1. 多くの研究でそうしているように、クロンバックのアルファを使う。その仮定が満たされている限り
  # 2. ポリコリック相関行列からクロンバックのアルファ的なものを計算する (B D Zumbo, A M Gadermann, and C Zeisser. Ordinal versions of coefficients alpha and theta for likert rating scales. Journal of Modern Applied Statistical Methods, 6: 21-29, 2007.)
 # 3. 01データならKuder-Richardsonの公式を使う
 # 4. 他にもいろんな信頼性指標がある (R E Zinbarg, W Revelle, I Yovel, and W Li. Cronbach’s α, Revelle’s β, and McDonald’s ωh: Their relations with each other and two alternative conceptualizations of reliability. Psychometrika, 70(1): 123-133, 2005.)
 # A good review is Bruce Thompson. Score Reliability. Contemporary Thinking on Reliability issues. Sage Publications, 2003. だそうだ。まあとりあえずはクロンバックを使おう。気が向いたらRevelle先生のやつを読もう
# 関連記事：Rで級内相関

# Kuder-Richardsonの公式
 # 以下参考URL。感謝
 # http://www.eflclub.com/elvin/publications/2003/itemanalysis.html
 # http://www.gifted.uconn.edu/siegle/research/Instrument%20Reliability%20and%20Validity/Reliability.htm
 # http://www.jstage.jst.go.jp/article/rika/17/3/17_181/_article/-char/ja

 dat <- data.frame(
 q001 = c(1,1,0,0,1,0,0,0,0),
 q002 = c(1,1,0,1,1,1,0,0,0),
 q003 = c(1,1,1,0,1,0,0,0,0),
 q004 = c(1,1,1,1,1,0,0,0,0),
 q005 = c(1,1,1,1,1,1,0,1,0),
 q006 = c(1,0,1,1,1,0,0,1,0),
 q007 = c(1,1,1,0,0,0,1,0,0),
 q008 = c(1,1,1,1,0,1,1,0,0),
 q009 = c(1,1,1,1,0,0,0,0,0),
 q010 = c(1,1,0,0,0,0,0,0,0)
 )
 dat

 (ni <- ncol(dat)) # 項目数
 (ns <- nrow(dat)) # 参加者数
 (mt <- mean(rowSums(dat))) # 合計得点の平均点
 (sdt <- sd(rowSums(dat))) # 合計得点の標準偏差
 (sds <- sqrt(var(rowSums(dat))*((ns-1)/ns))) # 合計得点の標本標準偏差

 pc <- colSums(dat)/ns # 各項目ごとの正答率 (1の個数/全参加者数)
 pw <- 1-pc # 各項目ごとの誤答率 (0の個数/全参加者数)
 pcpw <- pc*pw
 spcpw <- sum(pcpw)


# Kuder-Richardson 20
 (ni/(ni-1)) * (((sds*sds)-spcpw)/(sdt*sdt)) # 不偏標準偏差を使ったとき
 (ni/(ni-1)) * (((sds*sds)-spcpw)/(sds*sds)) # 標本標準偏差を使ったとき

# Kuder-Richardson 21
# おそらく、ローデータがないときに使う簡易版
ni/(ni-1) * (1-(mt-mt*mt/ni)/(sdt*sdt)) # 不偏標準偏差を使ったとき
ni/(ni-1) * (1-(mt-mt*mt/ni)/(sds*sds)) # 標本標準偏差を使ったとき

# よくわからんけど標本標準偏差を使うのかな



# おまけ。カッパ係数
 # 参考:http://www.mizumot.com/stats/kappa.htm
a <- c(2,1,3,2,2,3,2,1,3,3,3,3,2,1,2,3,1,3,1,2)
b <- c(2,1,3,2,3,3,3,1,3,3,3,3,2,1,2,3,2,3,1,2)
x <- table(a,b) # データフレームではなくtableで集計表にする
x

library(psych)
print.digits(cohen.kappa(x), digits=5)
 # 普通はunweighted kappaのestimateを見る

# psychパッケージ以外。まあpsychパッケージが一番いいけど
library(vcd)
Kappa(x)