Susceptible-High-risk-latent-Low-risk-latent-Infected-Recovered Model with Simple Demographics
SHLIR_demographics_ode(t, x, params)
The timestep over which to calculate derivatives
A numeric vector of compartment populations.
A named vector of parameter values.
A vector of derivatives
##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)#> [] #> S S H H I #> -2.831198 -4.156457 5.000000 1.487654 0.500000 #>