SEIQHRFjump: Simulates one jump from a simple stochastic SEIQHRF epidemic...

Description Usage Arguments Value Examples

Description

SEIQHRFjump returns states of next jump given the states right now.

Usage

1
SEIQHRFjump(x, pars)

Arguments

x

A numeric vector: (time, S, E, I, Q, H, R, F, N) at the previous jump.

pars

A numeric vector: (a.rate, e.rate, i.rate, h1.rate, q.rate, h2.rate, r1.rate, r2.rate, r3.rate, f.rate, pE, pSR, pIim, pQim, pHim).

Value

A numeric vector: (time, S, E, I, Q, H, R, F, N) immediately after the next jump.

Examples

1
2
3
A <- c(0,1000,0,0,0,0,0,0,1000)
B <- c(4,2,1,2,1,3,1,2,1,2,0.9, 0.3, 0.4,0.1,0.1)
SEIQHRFjump(A,B)

jessseliu/epidmod documentation built on Sept. 15, 2020, 8:53 p.m.