calculate_steady_state: Calculate steady-state compartment sizes for a network
Description
This is an experimental function. It attempts to calculate steady-state
compartment sizes using the set parameter values and the initial compartment
sizes. Use it with caution!
Usage
calculate_steady_state(nm)
Value
A tibble containing steady-state compartment sizes.
Arguments
nm
A network model, with set parameter values.
Details
Note about how steady state sizes for split compartments are calculated: the
steady size of the active portion is calculated divide it is divided by the
active fraction (portion.act parameter) to get the total size including the
refractory portion. In this case we get a "steady-state" refractory portion,
consistent with steady state size of active fraction and with portion.act
parameter.
m <- aquarium_mod
m <- set_prior(m, constant_p(0), "lambda")
m <- set_params(m, sample_params(m))
proj <- project(m, end = 40)
plot(proj)
z <- calculate_steady_state(m)
z
z$stable_sizes