R/summarize_simulation_error.R

In Enpeterson/outputsim: Output Simulation Results

summarize_simulation_error <-
  function(data) {
    # tib <-
    #   readRDS(paste(output.dir, '/mcmctib_', simnum, '.RDS', sep = ""))

    tidydat <- data %>%
      select(
        contains("sens"),
        contains("spec"),
        contains("vradj.ct"),
        contains("gamma.truematvr.ct")
      )

    if (ncol(tidydat) < 62) { #for testing purposes
      #only and all selected parameter should be in this set
      stop(
        paste0(
          "At least oneparameter of interest is missing. .\n",
          "Check to make sure mcmctib_ is correct"
        )
      )
      call. = FALSE
    }
    # else {
    #   tidydat <- tidydat %>%
    #     gather(key = "parameter", value = "value") %>%
    #     group_by(parameter) %>%
    #     summarize(
    #       lower = quantile(value, .025),
    #       median = quantile(value, .5),
    #       upper = quantile(value, .975)
    #     ) %>%
    #     print()
    # }
  }