Epidemic_fsMCMC MCMC algorithm to make inference on panel data of an epidemic through forward simulation
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | Epidemic_fsMCMC(
N,
a,
x,
beta0,
gamma0,
kernel = NULL,
no_draws,
s,
T_obs,
k,
lambda,
V,
no_its,
burn_in = 0,
lag_max = NA,
thinning_factor = 1
)
|
N |
Total size of closed population |
x |
panel data observed. Follows a random sample of n individuals and observes them at k timepoints |
beta0 |
beta starting value (Infection Rate Parameter) |
gamma0 |
gamma starting value (Removal rate parameter) |
no_draws |
How many rexp(1) and runif(1) draws to make for use in the Gillespie algorithm. |
T_obs |
the period for which the epidemic is observed. |
k |
how many equally spaced observations take place |
lambda |
RWM proposal parameter |
V |
RWM proposal covariance matrix |
no_its |
The number of MCMC iterations |
burn_in |
How many of the MCMC iterations are thrown away as burn in (Convergence to Stationary Distn) |
lag_max |
When plotting the estimated ACF of samples, what will be the maximum lag estimated/plotted |
thinning_factor |
Create Storage Matrix Proposal Acceptance Counter Propose new beta and gamma using Multiplicative RW propsal Draw New Random Variables Store State Calculating Summary Statistics for samples |
initial_infective |
Number of individuals who are initially infected in the population |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.