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)