calculate_steady_state: Calculate steady-state compartment sizes for a network

In isotracer: Isotopic Tracer Analysis Using MCMC

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)

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.

Value

A tibble containing steady-state compartment sizes.

Examples

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

isotracer documentation built on April 4, 2025, 3:57 a.m.