R/create_razzo_gen_experiment.R

In richelbilderbeek/razzo: The Error if Nature is MBD

#' Create a generative experiment that follows the razzo
#' naming conventions and article
#' @author Richèl J.C. Bilderbeek
#' @inheritParams default_params_doc
#' @param rng_seed RNG seed used in inference
#' @export
create_razzo_gen_experiment <- function(
  folder_name = peregrine::get_pff_tempfile(),
  rng_seed = 1
) {
  model_type <- "generative"
  pirouette::create_gen_experiment(
    inference_conditions = pirouette::create_inference_conditions(
      model_type = model_type,
      run_if = "always",
      do_measure_evidence = TRUE
    ),
    inference_model = beautier::create_inference_model(
      tree_prior = beautier::create_bd_tree_prior(),
      mcmc = razzo::get_razzo_mcmc(
        model_type = model_type,
        folder_name = folder_name
      ),
      mrca_prior = create_razzo_mrca_prior()
    ),
    beast2_options = create_razzo_beast2_options(
      model_type = model_type,
      folder_name = folder_name,
      rng_seed = rng_seed
    ),
    est_evidence_mcmc = create_razzo_ns_mcmc(
      folder_name = folder_name,
      model_type = model_type
    ),
    errors_filename = get_errors_filename(
      folder_name = folder_name,
      model_type = model_type
    )
  )
}