R/kldist.R

In bapred: Batch Effect Removal and Addon Normalization (in Phenotype Prediction using Gene Data)

Defines functions kldist

Documented in kldist

kldist <-
function(xba, batch) {

  xba <- scale(xba)

  allpairs <- combn(levels(batch), 2)

  sum(apply(allpairs, 2, function(y) (sum(batch==y[1]) + sum(batch==y[2]))*kldistTwo(xba[batch==y[1],], xba[batch==y[2],]))/(nrow(xba)*(length(levels(batch))-1)))

}

Try the bapred package in your browser