245 R -- クロンバックのα信頼性係数 青木繁伸 2002/01/13 (日) 16:13
x はデータ行列。合計値は含まない。 alpha <- function(x) { k <- ncol(x) k/(k-1)*(1-sum(apply(x,2,var))/var(apply(x,1,sum))) } # example x <- matrix(c(49,44,37,54, 36,36,36,42, 42,51,45,35, 47,67,54,40, 54,59,68,54, 51,55,59,67, 45,36,48,46, 72,49,50,58), byrow=TRUE, ncol=4) alpha(x) # returned value # [1] 0.7403366
alpha <- function(x) { k <- ncol(x) k/(k-1)*(1-sum(apply(x,2,var))/var(apply(x,1,sum))) } # example x <- matrix(c(49,44,37,54, 36,36,36,42, 42,51,45,35, 47,67,54,40, 54,59,68,54, 51,55,59,67, 45,36,48,46, 72,49,50,58), byrow=TRUE, ncol=4) alpha(x) # returned value # [1] 0.7403366
● 「統計学関連なんでもあり」の過去ログ--- 017 の目次へジャンプ● 「統計学関連なんでもあり」の目次へジャンプ ● 直前のページへ戻る