seirmod2: Gradient-function for the forced SEIR model

In epimdr: Functions and Data for "Epidemics: Models and Data in R"

Description

Gradient-function for the forced SEIR model

Usage

seirmod2(t, y, parms)

Arguments

t	Implicit argument for time
y	A vector with values for the states
parms	A vector with parameter values for the SIR system

Value

A list of gradients

Examples

require(deSolve)
times  = seq(0, 10, by=1/120)
paras  = c(mu = 1/50, N = 1, beta0 = 1000, beta1 = 0.2, sigma = 365/8, gamma = 365/5)
start = c(S=0.06, E=0, I=0.001, R = 0.939)
out=ode(y=start, times=times, func=seirmod2, parms=paras)

epimdr documentation built on March 26, 2020, 7:41 p.m.