CT <- matrix(c(42,8,10,40), nrow=2, byrow=TRUE)

CT

​

#クロス集計表から計算 　p.136参照

Phi <- abs(42*40-8*10)/sqrt((42+10)*(8+40)*(42+8)*(10+40))

Phi

​

​#パッケージpsychを利用するならば

library(psych)

phi(CT)

​