tidy_steady_states {isotracer}

Build a tidy table with the calculated steady states for each iteration

Description

If neither n_per_chain and n are provided, all iterations are used.

Usage

tidy_steady_states(nm, mcmc, n_per_chain = NULL, n = NULL)

Arguments

nm

A networkModel object.

mcmc

The corresponding output from run_mcmc.

n_per_chain

Integer, number of iterations randomly drawn per chain. Note that iterations are in sync across chains (in practice, random iterations are chosen, and then parameter values extracted for those same iterations from all chains).

n

Integer, number of iterations randomly drawn from mcmc. Note that iterations are *not* drawn in sync across chains in this case (use n_per_chain if you need to have the same iterations taken across all chains).

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 tidy table containing the mcmc iterations (chain, iteration, parameters), the grouping variables from the network model and the steady state sizes.