calcOmega {MGDrivE}

Solve for Omega (additional genotype-specific mortality)

Description

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

Usage

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)