R/collect_files_n_posteriors.R

In richelbilderbeek/wiritttea: WIRITTTE analysis

#' Collects the number of posteriors of all phylogenies belonging to a
#' multiple parameter file in the melted/uncast/long form
#' @param filenames names of the parameter file
#' @return A dataframe with all number of sampled posteriors of all files
#' @examples
#'  filenames <- c(
#'    find_path("toy_example_1.RDa"),
#'    find_path("toy_example_3.RDa")
#'  )
#'  df <- collect_files_n_posteriors(filenames)
#'  testit::assert(names(df) == c("filename", "n_posteriors"))
#'  testit::assert(nrow(df) == length(filenames))
#'  testit::assert(df$n_posteriors == c(2, 8))
#' @export
collect_files_n_posteriors <- function(filenames) {

  if (length(filenames) < 1) {
    stop(
      "collect_files_n_posteriors: ",
      "there must be at least one filename supplied"
    )
  }

  # posteriors
  n_posteriors <- NULL
  for (filename in filenames) {
    this_n_posteriors <- NULL
    tryCatch(
      this_n_posteriors <- collect_file_n_posteriors(
        filename = filename
      ),
      error = function(msg) {} # nolint
    )
    if (is.null(this_n_posteriors)) {
      this_n_posteriors <- data.frame(
        n_posteriors = NA
      )
    }
    if (!is.null(n_posteriors)) {
      n_posteriors <- rbind(n_posteriors, this_n_posteriors)
    } else {
      n_posteriors <- this_n_posteriors
    }
  }
  df <- data.frame(
    filename = basename(filenames),
    n_posteriors = n_posteriors$n_posteriors
  )
  testit::assert(nrow(df) == length(filenames))
  testit::assert(names(df) == c("filename", "n_posteriors"))
  return(df)
}