<<目次へ

四分点相関係数（φ係数）

　四分点相関係数を求める関数です。引数には2行2列のデータ行列（=2×2分割表）を指定してください。

phi.coe <- function(dat){
        x <- c(rep(1,dat[1,1]),rep(1,dat[1,2]),rep(0,dat[2,1]),rep(0,dat[2,2]))
        y <- c(rep(1,dat[1,1]),rep(0,dat[1,2]),rep(1,dat[2,1]),rep(0,dat[2,2]))
        phi <- cor(x,y)
        phi
}

〜使用例〜

> my.table <- matrix(c(4, 1, 2, 3), ncol=2, byrow=TRUE)
> my.table
     [,1] [,2]
[1,]    4    1
[2,]    2    3
> phi.coe(my.table)
[1] 0.4082483