R/genLibSizes.R

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

#' Generate library sizes from log-normal
#'
#' @param fit.ln the library size model
#' @param exprmt.design the design
#'
#' @return The generated libray sizes per batch and group
#' @importFrom stats rlnorm
genLibSizes = function(fit.ln, exprmt.design){
  n.batch = rowSums(exprmt.design$exprmt.config)
  lapply(unique(exprmt.design$sub.batchs), function(b){
    rlnorm(n.batch[b], fit.ln[[b]][["meanlog"]], fit.ln[[b]][["sdlog"]])
  })
}