R/D2K.R

In aplantin/MiRKATS: Microbiome Regression-based Kernal Association Test for Survival (MiRKAT-S)

Defines functions D2K

Documented in D2K

D2K <-
function(D){
  n <- nrow(D)
  centerM <- diag(n) - 1/n
  K <- -0.5*centerM %*% (D*D) %*% centerM
  eK <- eigen(K, symmetric=TRUE)
  K <- eK$vector %*% diag(pmax(0,eK$values)) %*% t(eK$vector)
  return(K)
}