Calculate equilibrium state variables, in model with delay to treatment — get_equilibrium_states_vivax

Calculates equilibrium state variables U0, Ul, T0, Tl, S0 and Sl based on I and model parameters

get_equilibrium_states_vivax_delay(
  I,
  lambda,
  r,
  gamma,
  f,
  alpha,
  beta,
  sigma,
  rho,
  delta,
  omega
)

Arguments

I: proportion of infectious individuals at equilibrium (U0+Ul)
lambda: transmission rate
r: blood clearance rate
gamma: liver clearance rate
f: relapse frequency
alpha: effective treatment probability
beta: probability of liver clearance
sigma: delay to treatment (1/duration of the delay)
rho: observation rate
delta: importation rate
omega: vector control intensity. If vector control is not the purpose of the analysis, omega can be fixed to 1.

A list with the equilibrium states (U0, Ul, S0 and Sl)