# R: Gradient-function for the SIRWS model

sirwmod {epimdr} | R Documentation

## Gradient-function for the SIRWS model

### Description

Gradient-function for the SIRWS model

### Usage

    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

    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))

[Package epimdr version 0.6-4 Index]