Call parasite_density and return the PDF: $$f_B(\xi; a, \tau |h) = \log_{10} \left( \sum_{M_\tau(a|h)>0} 10^{P_\tau(a |h)}\right)$$

d_parasite_density(
  meshX,
  a,
  FoIpar,
  tau = 0,
  hhat = 1,
  r = 1/200,
  RBC_par = par_lRBC_static(),
  Fmu_par = par_Fmu_base(),
  Omega_par = par_Omega_beta(),
  pWda = par_Wda_none()
)

Arguments

meshX

a mesh over parasite densities

a

host cohort age

FoIpar

parameters that define an FoI function

tau

the cohort birthday

hhat

a local scaling parameter for the FoI

r

the clearance rate for a simple infection

RBC_par

parameters to compute log10RBC

Fmu_par

parameters to compute Fmu

Omega_par

parameters to compute d_Omega

pWda

parameters to dispatch Wda

Value

a numeric vector of length meshX