qsurvival.nllik {fracprolif} | R Documentation |
Negative Log-Likelihood of a Survival Function
Description
Using a survival function on a specified distribution with a given set of parameters, computes the negative log-likelihood
Usage
qsurvival.nllik(dist, complete.lifespans, censored.lifespans, Q, ...)
Arguments
dist |
A text description of the distribution to use, e.g. emg or norm |
complete.lifespans |
A vector of observations of complete lifespans |
censored.lifespans |
A vector of incomplete observations of lifespans |
Q |
The percent of population in quiescence or non-dividing |
... |
Specific parameters to the |
Value
A numeric value of the negative log-likelihood for the given parameters.
Author(s)
Darren Tyson, Shawn Garbett
See Also
mle-class
q.mle.norm.estimate
q.mle.emg.estimate
q.rates
Examples
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
mle(function(mean, sd, Q){
qsurvival.nllik("norm", mitotic.lifespans, censored.lifespans, Q, mean, sd)
},
method='L-BFGS-B',
lower=list(mean=8, sd=0.1, Q=0.01),
upper=list(mean=30, sd=20, Q=0.9),
start=list(mean=mean(mitotic.lifespans),
sd=sd(mitotic.lifespans),
Q = 0.5))
[Package fracprolif version 1.0.7 Index]