R/binPurity.R

#' @importFrom stats sd

binPurity <- function(mz,intensity,dp = 2){
  intensity <- ceiling(intensity)
  if (sum(intensity) > 1) {
    mzs <- unlist(map2(mz,intensity,~{rep(.x,.y)}))
    purity <- 1 - (sd(mzs)/(1*10^-dp))
  } else {
    purity <- 1
  }
  return(purity)
}

Try the binneR package in your browser

Any scripts or data that you put into this service are public.

binneR documentation built on May 1, 2019, 6:34 p.m.