coyne2 {epimdr2}

Edition 2 Gradient-function for Coyne et al's rabies model

Description

Edition 2 Gradient-function for Coyne et al's rabies model

Usage

coyne2(t, logx, parms)

Arguments

t

Implicit argument for time

logx

A vector with values for the log-states

parms

A vector with parameter values for the dynamical system

Value

A list of gradients for the log system

Examples

require(deSolve)
times  = seq(0, 50, by=1/520)
paras  = c(gamma = 0.0397, b = 0.836, a = 1.34, sigma = 7.5, 
alpha = 66.36, beta = 33.25, c = 0, rho = 0.8)
start = log(c(S=12.69/2, E1=0.1, E2=0.1, I = 0.1, R = 0.1))
out = as.data.frame(ode(start, times, coyne, paras))