Given an mle estimate of the distribution parameters transform into rates of growth and quiescence. The transformation into a growth rate is the exact method of Powell. The quiescence rate uses Poisson estimation. Standard errors are computed using the numerical delta method.
This distribution to use for estimation, e.g. "emg" or "norm"
An object of class
An object of
class q_rate containing (d) the growth rate and (q) the quiescence rate as a numeric vector, as well as additional attributes of stderr, n and df used.
Darren Tyson, Shawn Garbett
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.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
data(ca1d.erlotinib) mitotic.lifespans <- subset(ca1d.erlotinib, !End.of.Expt & !Death & !is.na(Lifespan))$Lifespan censored.lifespans <- subset(ca1d.erlotinib, End.of.Expt & !Death & !is.na(Lifespan))$Lifespan r <- q.rates("emg", q.mle.emg.estimate(mitotic.lifespans, censored.lifespans)) summary(r) r['d'] # Rate of division r['q'] # Rate of quiescense
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.