If neither n_per_chain
and n
are provided, all iterations are
used.
tidy_steady_states(nm, mcmc, n_per_chain = NULL, n = NULL)
A tidy table containing the mcmc iterations (chain, iteration, parameters), the grouping variables from the network model and the steady state sizes.
A networkModel
object.
The corresponding output from run_mcmc
.
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).
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).
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.