R/geometric_mean.R

In SCIL-leuven/qPCR: qPCR analysis from viia7

#' Geometric mean
#'
#' Normalize to housekeeping genes
#' @param x numeric
#' @return numeric
#' @export
gm_mean = function(x, na.rm=TRUE){
  exp(sum(log(x[x > 0]), na.rm=na.rm) / length(x))
}