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.

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 Sept. 22, 2023, 1:07 a.m.