Computes the expected count of a subpopulation that is propagating inside population undergoing quiescent growth.
Usage
prolif_d(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 dividing subpopulation totals.