R/summarize_runsetting.R

#' Summarize MCMC Samples (across multiple simulations)
#'
#'Using a loop funciton, takes in MCMC samples from each simulation run, combines them into one tibble that list MCMC samples by simulation number.
#'Global parameters and sens.ct, spec.ct, gammatruematvr.ct, and vradj.ct are the parameters included in each tibble. The combined tibble is saved to the output folder.
#' @param output.dir
#' @param simsetting
#' @param nsim
#'
#' @return One tibble consisting of all MCMC samples from all simulations (with parameters listed in description. Note not all parameters have been included)
#' @export
#'
#' @examples
summarize_runsetting <- function(output.dir, simsetting, nsim){
  tiblist <- list()
  for(i in 1:nsim){
    tiblist[[i]] <- readRDS(paste(output.dir, "sumtib_", i, '.RDS', sep = ""))
  }

  simsetting_tib <- do.call("rbind", tiblist)

  saveRDS(simsetting_tib, paste(output.dir, "sumsetting_", simsetting, ".RDS", sep=""))
  return(simsetting_tib)
}
Enpeterson/outputsim documentation built on May 24, 2019, 9:53 a.m.