This function solves the ODE in Eq. (2.7)
1 | SIR.ODE.full(Tmax, dt = 0.01, beta, gamma, ic)
|
Tmax |
cutoff time of epidemic |
dt |
time increment of ODE |
beta |
present values of beta at the present step at MCMC |
gamma |
present values of gamma at the present step at MCMC |
ic |
initial values of ODE set to c(1.0, rho, 0.0), rho set the present values of beta and gamma at the present step at MCMC |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.