Susceptible-High-risk-latent-Low-risk-latent-Infected-Recovered Model with Simple Demographics

Susceptible-High-risk-latent-Low-risk-latent-Infected-Recovered Model with Simple Demographics

SHLIR_demographics_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_0 <- 989
H_0 <- 10
L_0 <- 0
I_0 <- 1
R_0 <- 0
beta = 3 # Rate of transmission
gamma_H = 1/5 # Rate of progression to active symptoms from high risk latent
nu = 1/2 #Rate of progression from high to low risk latent
gamma_L = 1/100 # Rate of progression to active symptoms for low risk latent
tau = 1/2 # Rate of recovery
mu = 1/81 # Rate of natural mortality
parameters <- c(beta = beta, gamma_H = gamma_H, gamma_L = gamma_L, nu = nu, tau = tau, mu = mu)
inits <- c(S = S_0, H = H_0, L = L_0, I = I_0, R_0 = R_0)

SHLIR_demographics_ode(1, inits, parameters)
#> [[1]]
#>         S         S         H         H         I 
#> -2.831198 -4.156457  5.000000  1.487654  0.500000 
#>