Solve a system of differential equations to compute the true PR and the MoI

solve_dpda(h, FoIpar, r = 1/200, tau = 0, Amax = 730, dt = 1)

Arguments

h

a scaling factor on the FoI

FoIpar

parameters that define an FoI function

r

the clearance rate for a simple infection

tau

the cohort birthday

Amax

The maximum runtime (in days)

dt

The output frequency (in days)

Value

a data.frame with the orbits