steady {mizer} | R Documentation |
Set initial values to a steady state for the model
Description
The steady state is found by running the dynamics while keeping reproduction,
resource and other components constant until the size spectra no longer
change much (or until time t_max
is reached, if earlier).
Usage
steady(
params,
t_max = 100,
t_per = 1.5,
dt = 0.1,
tol = 0.1 * dt,
return_sim = FALSE,
preserve = c("reproduction_level", "erepro", "R_max"),
progress_bar = TRUE
)
Arguments
params |
A MizerParams object |
t_max |
The maximum number of years to run the simulation. Default is 100. |
t_per |
The simulation is broken up into shorter runs of |
dt |
The time step to use in |
tol |
The simulation stops when the relative change in the egg
production RDI over |
return_sim |
If TRUE, the function returns the MizerSim object holding
the result of the simulation run, saved at intervals of |
preserve |
Specifies whether the |
progress_bar |
A shiny progress object to implement a progress bar in a shiny app. Default FALSE. |
Details
If the model use Beverton-Holt reproduction then the reproduction parameters
are set to values that give the level of reproduction observed in that
steady state. The preserve
argument can be used to specify which of the
reproduction parameters should be preserved.
Value
A MizerParams or a MizerSim object
Examples
params <- newTraitParams()
species_params(params)$gamma[5] <- 3000
params <- steady(params)
plotSpectra(params)