Calculate equilibrium state variables — get_equilibrium_states_vivax

Calculates equilibrium state variables I0, Il, S0 and Sl based on I and model parameters

get_equilibrium_states_vivax_rcd(
  I,
  lambda,
  r,
  gamma,
  f,
  alpha,
  beta,
  rho,
  delta,
  omega,
  iota_star,
  nu,
  tau,
  eta,
  rho2
)

Arguments

I: proportion of infectious individuals at equilibrium (I0+Il)
lambda: transmission rate
r: blood clearance rate
gamma: liver clearance rate
f: relapse frequency
alpha: effective treatment probability
beta: probability of liver clearance
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.
iota_star: RCD parameter: maximum number of index cases investigated per population per unit of time at equilibrium
nu: RCD parameter: number of secondary cases investigated per index case
tau: RCD parameter: targeting ratio of the RCD
eta: RCD parameter: probability that an investigated cases is detected (symptoms, test sensitivity, etc.)
rho2: RCD parameter: observation rate for reactively-detected cases

Value

A list with the equilibrium states (I0, Il, S0 and Sl)