tsirSpat {epimdr2}

Function to simulate the spatially-extended seasonally-forced TSIR

Description

Function to simulate the spatially-extended seasonally-forced TSIR for a patchily distributed host population. Coupling is assumed to be global and according to a commuter model so with p patches and a coupling of c, local transmission is reduced by a fraction (1-c*p).

Usage

tsirSpat(beta, alpha, B, N, p, c, inits, type = "det")

Arguments

beta

the seasonal transmission coefficients

alpha

the exponent on I

B

a vector of Births (the length of which determines the length of the simulation)

N

the population size

p

the number of patches

c

the the spatial coupling

inits

a list containing a vector of initial S and a vector of initial I for each patch

type

an argument "det" or "stoc" that determines whether a deterministic or stochastic simulation is done

Value

A list with time series of simulated infected and susceptible hosts

Examples

"see chapter 15 in book"