R/ADJ.R

In mudfold: Multiple UniDimensional unFOLDing

ADJ <-function(X){
  J <- colnames(X)
  K <- ncol(X)
  n <- nrow(X)
  ADJ.mat <- matrix(0,nrow=K,ncol=K)
  ADJ.mat <- t(as.matrix(X)) %*% as.matrix(X) / n
  ADJ.mat[!lower.tri(ADJ.mat)] <- NA
  dimnames(ADJ.mat) <- list(J,J)
  return(ADJ.mat)
}