inhomoModPoissonTLforSIWR: Modified inhomogeneous Poisson tau-leaping algorithm for SIWR...

In vnminin/bayessir: Particle marginal Metropolis-Hastings (PMMH) algorithm for time-varying Susceptible-Infectious-Recovered-Susceptible (SIRS) model

Description

Uses a tau-leaping algorithm to simulate S, I, W, and R forward in time when value of environmental force of infection changes

Usage

inhomoModPoissonTLforSIWR(startState, popSize, Wmax, startTime, intervalLength,
  parameters, alphas, allbreaks, deltatint, ncrit)

Arguments

startState	Integer vector of size three containing the initial values for S, I, and W
popSize	Population size
startTime	Start time of simulation
intervalLength	Length of simulation interval
parameters	Numeric vector of size six containing the current values of the infectious contact rate, betaW, the recovery rate, the rate at which immunity is lost, kappa, and eta
alphas	Numeric vector containing the values of alpha
allbreaks	Numeric vector containing the times at which the value of alpha changes
deltatint	Initial value to use for tau in tau-leaping algorithm
ncrit	Critical size for modified tau-leaping algorithm; if the population of a compartment is lower than this number a single step algorithm is used until the population gets above the critical size.

Value

Vector containing values of S, I, and W at the end of the simulation interval

vnminin/bayessir documentation built on May 3, 2019, 6:17 p.m.