R/getMedianNormalizationFactors.R

In nicohuttmann/htmnanalysis: Analysis of Omics Data

#' Returns factors from PQN normalization
#'
#' @param data
#'
#' @return
#' @export
#'
#'
#'
getMedianNormalizationFactors <- function(data) {

  factors <- matrix(NA, nrow = nrow(data), ncol = nrow(data))
  rownames(factors) <- colnames(factors) <- rownames(data)

  for(i in 1:nrow(data)) {
    factors[i, ] <- medianFactor(factorizeMatrix(data, i))
  }

  apply(factors, 2, geomean)

}