R/estLibSizeDistr.R

In SPsimSeq: Semi-parametric simulation tool for bulk and single-cell RNA sequencing data

#' Estimate log-normal distribution for the library sizes 
#'
#' @param LS observed library sizes
#' @param batch batches
#' 
#' @return Estimated log-normal parameter library sizes
#' @importFrom fitdistrplus fitdist
estLibSizeDistr = function(LS, batch){
  tapply(LS, batch, function(L.b){
      fitdist(L.b, distr = "lnorm")$estimate 
  })
}