#' 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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.