This helper functions provide the same functionality as the prolif_tot set of functions, but with the additional assumption that rates change at some breakpoint during the dataset.
split_tot(t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br=20)
split_d( t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br=20)
split_q( t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br=20)
split_fd( t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br=20)
split_fq( t, x0, q0, d1, d2, q1, q2, ad1, aq1, ad2, aq2, br=20)A vector time points to return expected population size.
The size of the initial propagating (or dividing) population at t=0.
The size of the quiescent or non-propagating population at t=0.
The rate of growth of the propagating population in the first section.
The rate of growth of the propagating population in the second section.
The rate at which propagating population members join the quiescent population in the first section.
The rate at which propagating population members join the quiescent population in the second section.
The rate of death from the dividing or propagating population in the first section.
The rate of death from the quiescent population in the first section.
The rate of death from the dividing or propagating population in the second section.
The rate of death from the quiescent population in the second section.
The time breakpoint for the change in rate
A numerical vector of expected (sub)population totals (or fraction)
The method of determining growth rate is taken from Powell E.0. (1956). Growth Rate and Generation Time of Bacteria, with Special Reference to Continuous Culture. J.Gen.Microbial V15,492-511. This makes a robust estimator in the presence of skewed distributions.