★ R -- κ統計量 ★

 246 R -- κ統計量  青木繁伸  2002/01/13 (日) 16:59


246. R -- κ統計量  青木繁伸  2002/01/13 (日) 16:59
kappa.stat <- function(o)
{
    n <- sum(o)
    e <- crossprod(apply(o,1,sum),apply(o,2,sum))/n
    qo <- 1-(po <- sum(diag(o))/n)
    qe <- 1-(pe <- sum(diag(e))/n)
    kappa <- 1-qo/qe
    sk <- sqrt(po*qo/(n*qe^2))
    sk0 <- sqrt(pe/(n*qe))
    result <- c(kappa, sk, sk0, kappa-qnorm(0.975)*sk, kappa+qnorm(0.975)*sk, kappa/sk0, pnorm(z,lower=F))
    names(result) <- c("kappa", "sigma-kappa", "sigma-kappa-0", "95% lcl", "95% ucl", "Z value", "P value")
    result
}

# example
x <- matrix(c(12,6,1, 3,19,4, 2,5,34), byrow=TRUE, 3)
kappa.stat(x)


● 「統計学関連なんでもあり」の過去ログ--- 017 の目次へジャンプ
● 「統計学関連なんでもあり」の目次へジャンプ
● 直前のページへ戻る