EBSeq/R/RankNorm.R

   1
   2 RankNorm=function(Data){
   3         RankData=apply(Data, 2, rank)
   4         SortData=apply(Data, 2, sort)
   5         SortMean=rowMeans(SortData)
   6         SortMean[SortMean==0]=1
   7         NormMatrix=sapply(1:ncol(Data), function(i)Data[,i]/(SortMean[RankData[,i]]))
   8         NormMatrix[NormMatrix==0]=1
   9         NormMatrix
  10         }
  11