Calculate the log likelihood of the negative binomial function (and its derivatives)
llikNbinomMu(x, size, mu, full = FALSE)
x
: Number of successessize
: Size of trialmu
: mu parameter for negative binomialfull
: Add the data frame showing x, mean, sd as well as the fx and derivativesdata frame with fx
for the pdf value of with dProb
that has the derivatives with respect to the parameters at the observation time-point
In an rxode2()
model, you can use llikNbinomMu()
but you have to use all arguments. You can also get the derivative of mu
with llikNbinomMuDmu()
llikNbinomMu(46:54, 100, 40) llikNbinomMu(46:54, 100, 40, TRUE) et <- et(46:54) et$size <- 100 et$mu <- 40 model <- function() { model({ fx <- llikNbinomMu(time, size, mu) dProb <- llikNbinomMuDmu(time, size, mu) }) } rxSolve(model, et)
Matthew L. Fidler
Useful links