Compute the distribution function for \(B_\bday(a)\)

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

Usage

d_parasite_density(
  meshX,
  a,
  FoI_a,
  bday = 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

FoI_a

a cohort trace function

bday

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