phi.coef <-
function(t) {
f <- (t[1,1]*t[2,2]-t[1,2]*t[2,1])/sqrt(colSums(t)[1]*colSums(t)[2]*rowSums(t)[1]*rowSums(t)[2])
return(f)
}