chainSIR: Gradient-function for the chain-SIR model

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

Description

Gradient-function for the chain-SIR model

Usage

chainSIR(t, logx, params)

Arguments

t	Implicit argument for time
logx	A vector with values for the log-states
params	A vector with parameter values for the chain-SIR system

Value

A list of gradients

Examples

require(deSolve)
times  = seq(0, 10, by=1/52)
paras2  = c(mu = 1/75, N = 1, beta =  625, gamma = 365/14, u=5)
xstart2 = log(c(S=.06, I=c(0.001, rep(0.0001, paras2["u"]-1)), R = 0.0001))
out = as.data.frame(ode(xstart2, times, chainSIR, paras2))

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