<<目次へ


クラメールの連関係数



 クラメールの連関係数を求めるための関数です。

cramer.coe <- function(dat){
        row.sum <- apply(dat, 1, sum)
        col.sum <- apply(dat, 2, sum)
        row.length <- length(dat[,1])
        col.length <- length(dat[1,])
        my.sum = 0
        for(i in 1:row.length){
                for(j in  1:col.length){
                        num <- dat[i,j]^2 / (row.sum[i] * col.sum[j])
                        my.sum <- my.sum + num
                }
        }
        coe <- sqrt((my.sum - 1) / (row.length - 1))
        coe
}


〜使用例〜

> dat <- matrix(c(18,0,4,3,7,16,8,12,8),ncol=3,byrow=TRUE)
> dat
     [,1] [,2] [,3]
[1,]   18    0    4
[2,]    3    7   16
[3,]    8   12    8
> cramer.coe(dat)
[1] 0.4624463