seirmod {epimdr2}

Gradient-function for the SEIR model

Description

Gradient-function for the SEIR model

Usage

seirmod(t, y, parameters)

Arguments

t

Implicit argument for time

y

A vector with initial values for the states

parameters

A vector with parameter values for the SEIR system

Value

A list of gradients

Examples

require(deSolve)
times  = seq(0, 10, by=1/120)
paras  = c(mu = 1/50, N = 1, beta =  1000, sigma = 365/8, gamma = 365/5)
start = c(S=0.06, E=0, I=0.001, R = 0.939)
out=ode(start, times, seirmod, paras)