Solve the queuing model \(M/M/\infty\)

A wrapper to solve the queuing model \(M/M/\infty\) (see dMoIda).

The function automatically sets the maximum MoI to be computed, and it sets initial conditions. The equations are solved using deSolve::ode and returned at regular intervals dt from age 0 up to Amax (in days).

Usage

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

Arguments

h

the force of infection

FoIpar

\(h_\tau(a)\), a list formatted to compute FoI

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 list with the orbits by name

See also

dMoIda