powered by
Negative log-likelihood function for the chain-binomial model
llik.cb(S0, beta, I)
a scalar with value for S0
a scalar with value for beta
a vector incidence aggregated at serial interval
the negative log-likelhood for the model
# NOT RUN { twoweek=rep(1:15, each=2) niamey_cases1=sapply(split(niamey$cases_1[1:30], twoweek), sum) llik.cb(S0=6500, beta=23, I=niamey_cases1) # }
Run the code above in your browser using DataLab