  Lymphocyte populations, stimulated in vitro or in vivo, grow as cells divide.
Stochastic models are appropriate because some cells undergo multiple rounds of
division, some die, and others of the same type in the same conditions do not
divide at all. If individual cells behave independently, each can be imagined
as sampling from a probability density of times to division. The most
convenient choice of density in mathematical and computational work, the
exponential density, overestimates the probability of short division times. We
consider a multi-stage model that produces an Erlang distribution of times to
division, and an exponential distribution of times to die. The resulting
dynamics of competing fates is a type of cyton model. Using Approximate
Bayesian Computation, we compare our model to published cell counts, obtained
after CFSE-labelled OT-I and F5 T cells were transferred to lymphopenic mice.
The death rate is assumed to scale linearly with the generation (number of
divisions) and the number of stages of undivided cells (generation $0$) is
allowed to differ from that of cells that have divided at least once
(generation greater than zero). Multiple stages are preferred in posterior
distributions, and the mean time to first division is longer than the mean time
to subsequent divisions.
