MCC <- function(t) {
  P1 <- colSums(t)[1]
  N1 <- colSums(t)[2]
  P2 <- rowSums(t)[1]
  N2 <- rowSums(t)[2]
  f  <- (t[1,1]*t[2,2]-t[1,2]*t[2,1])/sqrt(P1*N1*P2*N2)
return (f)
}