Expected Proliferation Total

Share:

Description

Returns the expected population size over time given quiescent growth model parameters.

Usage

1
prolif_tot(t, x0, q0=0, d, q, ad, aq, tol=1e-6)

Arguments

t

A vector time points to return expected population size.

x0

The size of the initial propagating (or dividing) population at t=0.

q0

The size of the quiescent or non-propagating population at t=0.

d

The rate of growth of the propagating population.

q

The rate at which propagating population members join the quiescent population.

ad

The rate of death from the dividing or propagating population.

aq

The rate of death from the quiescent population.

tol

The tolerance between parameters before switching to limit models. The solution to the quiescent growth model has some boundary cases that would result in a division by zero. These happen when the rate of growth (d) is equal to rate of quiescence (q). When the absolute difference between these parameters is less than the tolerance, the model switches to the limit cases. The default value is good for most common cases. Preferably a very small number, that must be greater than zero.

Value

A numerical vector of expected population totals.

Author(s)

Shawn Garbett

See Also

prolif_q prolif_d prolif_fq prolif_fd q.rates

Examples

1
  prolif_tot(1:100, 0.5, 0.5, 0.04, 0.03, 0.001, 0.001)