fb4package (version 2.0.0)

plot_sensitivity.fb4_result: Plot temperature sensitivity analysis for a Bioenergetic object

Description

Runs analyze_growth_temperature_sensitivity and plots the result. Sensitivity analysis requires a Bioenergetic object (not an fb4_result), because it re-runs the model across a grid of temperatures and p_values. Use plot(bio_obj, type = "sensitivity") as the primary interface; this function is the underlying implementation.

Usage

plot_sensitivity.fb4_result(
  bio_obj,
  temperatures = seq(4, 20, by = 2),
  p_values = seq(0.3, 1, by = 0.1),
  simulation_days = 365,
  color_scheme = "grayscale",
  add_annotations = TRUE,
  verbose = FALSE,
  ...
)

Value

Called for its plotting side-effect. Invisibly returns NULL.

Arguments

bio_obj

Bioenergetic object with species parameters, temperature profile, diet, and simulation settings.

temperatures

Numeric vector of absolute temperatures (°C) to test. Default seq(4, 20, by = 2).

p_values

Numeric vector of p_values (proportion of Cmax) to evaluate. Must be in (0, 1]. Default seq(0.3, 1.0, by = 0.1).

simulation_days

Number of simulation days. Default 365.

color_scheme

Color scheme for the plot. Default "grayscale".

add_annotations

Add optimal temperature annotations. Default TRUE.

verbose

Show analysis progress. Default FALSE.

...

Additional arguments passed to plot_growth_temperature_sensitivity().

Examples

# \donttest{
data(fish4_parameters)
sp   <- fish4_parameters[["Oncorhynchus tshawytscha"]]$life_stages$adult
info <- fish4_parameters[["Oncorhynchus tshawytscha"]]$species_info
bio  <- Bioenergetic(
  species_params     = sp,
  species_info       = info,
  environmental_data = list(
    temperature = data.frame(Day = 1:30, Temperature = rep(12, 30))
  ),
  diet_data = list(
    proportions = data.frame(Day = 1:30, Prey1 = 1.0),
    energies    = data.frame(Day = 1:30, Prey1 = 5000),
    prey_names  = "Prey1"
  ),
  simulation_settings = list(initial_weight = 100, duration = 30)
)
bio$species_params$predator$ED_ini <- 5000
bio$species_params$predator$ED_end <- 5500
plot_sensitivity.fb4_result(
  bio_obj         = bio,
  temperatures    = c(10, 14),
  p_values        = c(0.4, 0.7),
  simulation_days = 30,
  verbose         = FALSE
)
# }