R/create_mbd_paramses.R

In richelbilderbeek/razzo: The Error if Nature is MBD

#' Create a list of MBD parameters
#'
#' Each element has a unique MBD parameter combination
#' @inheritParams default_params_doc
#' @return a list of \code{mbd_params}
#' @seealso use \link{create_mbd_params} to create one
#'   MBD parameter combination
#' @author Richèl J.C. Bilderbeek
#' @examples
#' create_mbd_paramses()
#' @export
create_mbd_paramses <- function(
  n_replicates = get_razzo_n_replicates()
) {
  df <- mbd.SimTrees::create_params_table(
    n_replicates = n_replicates
  )

  n_paramses <- nrow(df)
  testit::assert(n_paramses > 0)
  mbd_paramses <- list()
  for (i in seq(1, n_paramses)) {
    mbd_params <- razzo::create_mbd_params(
      lambda = df$lambda[i],
      mu = df$mu[i],
      nu = df$nu[i],
      q = df$q[i],
      cond = df$cond[i],
      crown_age = df$crown_age[i],
      seed = df$seed[i]
    )
    mbd_paramses[[i]] <- mbd_params
  }
  mbd_paramses
}