R/skewdiv.R

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

Defines functions skewdiv

Documented in skewdiv

skewdiv <-
  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]))*skewdivTwo(xba[batch==y[1],], xba[batch==y[2],]))/(nrow(xba)*(length(levels(batch))-1)))
    
  }
  

Try the bapred package in your browser