Susceptible-Infected-Recovered-Susceptible Model with Vaccination

SIRS_vaccination_ode(t, x, params)

Arguments

t

The timestep over which to calculate derivatives

x

A numeric vector of compartment populations.

params

A named vector of parameter values.

Value

A vector of derivatives

Examples

##Model Input
S_u_0 <- 989
I_u_0 <- 1
R_u_0 <- 0
S_v_0 <- 0
I_v_0 <- 0
R_v_0 <- 0
beta <- 3
chi <- 0.5
tau <- 2
lambda <- 0.7
dt <- 1

parameters <- c(beta = beta, tau = tau,
                chi = chi, lambda = 0.7)
inits <- c(S_u = S_u_0, I_u = I_u_0, R_u_0 = R_u_0,
           S_v = S_v_0, I_v = I_v_0, R_v_0 = R_v_0)

SIRS_vaccination_ode(1, inits, parameters)
#> [[1]]
#>        S_u        S_u        I_u        S_v        S_v        I_v 
#> -2.9969697  0.9969697  2.0000000  0.0000000  0.0000000  0.0000000 
#>