llik.pc | R Documentation |
This function uses binomial likelihoods to estimate the picewise-constant FoI model from age-incidence data
llik.pc(par, age, num, denom, up)
par |
a vector with initial guesses |
age |
a vector with the ages |
num |
a vector with number infected by age |
denom |
a vector with number tested by age |
up |
a vector with upper age-bracket cut-offs |
The negative log-likelihood for a candidate piecewise constant catalytic model
x=c(1,4,8,12,18,24) para=rep(.1,length(x)) ## Not run: optim(par=log(para),fn=loglikpc, age=rabbit$a, num=rabbit$inf, denom=rabbit$n, up=x)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.