R/study-coverage_mcmc.R

In wlandau/fbseqStudies: fbseqStudies

#' @include study-coverage_init.R
NULL

#' @title Function \code{coverage_mcmc}
#' @description MCMC of the coverage study (credible intervals)
#' @export
#' @return path to simulated objects
#' @param path to directory to save simulations and results
#' @param zeronormfactors TRUE/FALSE value: option to set normalization constants h to 0 
coverage_mcmc = function(path = newdir(), zeronormfactors = T){
  path = newdir(path)
  coverage_init(path)
  fit(path, benchmarks = "edgeR", fbseq_methods = "fullybayes", zeronormfactors = zeronormfactors, 
    priors = c("normal", special_beta_priors()))
  path
}