Calculate equilibrium state variables, RCD model (no referral) — get_equilibrium_states_vivax_rcd_no

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

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

Arguments

I: proportion of infectious individuals at equilibrium (U0+Ul+T0+Tl)
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
sigma: delay to treatment (1/duration of the delay)
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 (U0, Ul, T0, Tl, S0 and Sl)