twostrain: Gradient-function for the two-strain SIR model
In epimdr2: Functions and Data for "Epidemics: Models and Data in R (2nd Edition)"

twostrain

R Documentation

Gradient-function for the two-strain SIR model

Description

Gradient-function for the two-strain SIR model

Usage

twostrain(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 two-strain SIR system

Value

A list of gradients

Examples

require(deSolve)
times  = seq(0, 30, by=1/200)
paras  = c(mu = 0.02, N = 1, beta1=500, beta2=750, gamma = 365/5, Theta=0.15, Xi=0.15, Pi=0.8)
start = c(S = 0.999, I1 = 0.001, I2 = 0.00, R1=0, R2=0, J1=0, J2=0, R = 0)
out = as.data.frame(ode(start, times, twostrain, paras))

epimdr2 documentation built on Dec. 28, 2022, 2:23 a.m.