calcOmega function

Solve for Omega (additional genotype-specific mortality)

Solves for root of equation of geometrically-distributed lifespan for value of omega.

calcOmega(mu, lifespanReduction)

Arguments

• mu: Daily mortality probability (discrete-time hazard, called muAd in code)
• lifespanReduction: Target reduced lifespan, between 0 and 1 (target average lifespan will be $\frac{1}{\mu_{Ad}} \times lifespanReduction$)

Examples

# reduce lifespan by 10%
#  Example mu is an average for Aedes
newOmega <- calcOmega(mu = 0.11, lifespanReduction = 0.90)