sirwmod: Gradient-function for the SIRWS model

Description Usage Arguments Value Examples

View source: R/rfuncs.R

Description

Gradient-function for the SIRWS model

Usage

1
sirwmod(t, logy, parms)

Arguments

t

Implicit argument for time

logy

A vector with values for the log(states)

parms

A vector with parameter values for the SIRWS system

Value

A list of gradients (in log-coordinates)

Examples

1
2
3
4
5
require(deSolve)
times  = seq(0, 26, by=1/10)
paras  = c(mu = 1/70, p=0.2, N = 1, beta = 200, omega = 1/10, gamma = 17, kappa=30)
start = log(c(S=0.06, I=0.01, R=0.92, W = 0.01))
out = as.data.frame(ode(start, times, sirwmod, paras))

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

Related to sirwmod in epimdr...