R/sum_across_sim.R

In Enpeterson/outputsim: Output Simulation Results

#' Summarize across runsettings
#'
#' @param runsettings Vector of names of runsettings
#'
#' @return Table of bias and coverage proportions for each parameter in each runsetting
#' @export
#'
#' @examples sum_across_sim(c("alldata","nobreakdown"))
#'

sum_across_sim <- function(runsettings) {
  #takes in a vector of the runsettings to summarize


summaster <- tibble(
  parameter = as.character(),
  lower = as.numeric(),
  median = as.numeric(),
  upper = as.numeric(),
  sim = as.numeric(),
  runsetting = as.character(),
  truth = as.numeric(),
  coverage = as.numeric(),
  bias = as.numeric()
)

  for (i in 1:length(runsettings)) { #reads in each summarized runsetting tibble and joins them
    #output.dir <- paste(getwd(), '/Output/', runsettings[i], sep = "")
    output.dir <- paste(output.toy.dir, runsettings[i], sep = "")
    dat <-
      readRDS(paste(output.dir, "/sumsetting_", runsettings[i], ".RDS", sep = ""))
    summaster <- full_join(summaster, dat)
  }


sumtable <- summaster %>%
  group_by(parameter, runsetting) %>%
  summarize(coverage = mean(coverage),
            bias = (mean(truth) - mean(median)))

saveRDS(sumtable, paste(output.toy.dir, "sumacrossruns.RDS", sep = ""))
return(sumtable)
}