SEIQRjump: Simulates one jump from a simple stochastic SEIR epidemic...

In jessseliu/epidmod: epidmod

Description

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

Usage

SEIQRjump(x, pars)

Arguments

x	A numeric vector: (time, S, E, I, Q, R, N) at the previous jump.
pars	A numeric vector: (a.rate, e.rate, i.rate, q.rate, r.rate, pE, pSR, pIm).

Value

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

Examples

A <- c(0,1000,0,0,0,0,1000)
B <- c(1/12, 1, 1/4, 0.15, 1/5, 0.3, 0.2, 0.7)
SEIQRjump(A,B)

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