Skip to contents

Compare Multiple Growth Models

Source: R/plot_vital_functions.R
compare_growth_models.Rd

Creates a faceted comparison of multiple growth model fits.

Usage

compare_growth_models(
  ...,
  data = NULL,
  age_col = "age",
  length_col = "length",
  sex_col = "sex",
  sex_labels = c(`1` = "Female", `2` = "Male"),
  age_range = NULL,
  k_based_list = NULL,
  which_model_list = NULL,
  colors = NULL,
  colorblind = FALSE,
  base_size = 12
)

Arguments

...

Named CmdStanMCMC objects from fit_bayesian_growth. Names become facet labels.

data

Optional data.table/data.frame with observed data.

age_col, length_col, sex_col

Column names in data.

sex_labels

Named character vector. Labels for sex in plot legend and facets. Default c("1" = "Female", "2" = "Male").

age_range

Numeric vector of length 2. If NULL, auto-determined.

k_based_list

Logical vector indicating which models are k-based. If NULL, assumes all are maturity-based.

which_model_list

Integer vector of model types. If NULL, assumes VBGM.

colors

Color palette.

colorblind

Logical. Use colorblind-safe palette? Default FALSE.

base_size

Base font size.

Value

A ggplot object with faceted model comparison.

See also

vignette("visualization") for comprehensive plotting examples.

vignette("model_diagnostics") for model comparison via LOO-CV.

Statistical Methods: Model Assessment for LOO-CV theory.

plot_growth_curve, compute_loo, compare_loo

Examples

if (FALSE) { # \dontrun{
compare_growth_models(
  "von Bertalanffy" = vb_fit,
  "Gompertz" = gomp_fit,
  "Logistic" = log_fit,
  data = shark_data
)

# With French labels
compare_growth_models(
  "von Bertalanffy" = vb_fit,
  "Gompertz" = gomp_fit,
  data = shark_data,
  sex_labels = c("1" = "Femelle", "2" = "Mâle")
)
} # }