calcZeroInflation: Calculates the zero-inflation part of a hurdle exponential...

In MGDrivE: Mosquito Gene Drive Explorer

Description

Given the probability of an adult mosquito to stay in the same patch throughout its whole lifespan, and its mortality, it calculates the height of the pulse-density part of the hurdle kernel.

Usage

calcZeroInflation(stayThroughLifespanProbability, adultMortality)

Arguments

stayThroughLifespanProbability	Probability of a mosquito to spend its whole lifespan in the same node
adultMortality	Adult mortality rate

Examples

# setup distance matrix
# two-column matrix with latitude/longitude, in degrees
latLong = cbind(runif(n = 5, min = 0, max = 90),
                runif(n = 5, min = 0, max = 180))

# Vincenty Ellipsoid  distance formula
distMat = calcVinEll(latLongs = latLong)

# get hurdle height
# Lets assume 80% stay probs and adult mortality of 0.1
hHeight <- calcZeroInflation(stayThroughLifespanProbability = 0.80,
                             adultMortality = 0.1)

# calculate hurdle exponential distribution over distances
kernMat = calcHurdleExpKernel(distMat = distMat, rate = 10, p0 = hHeight)

MGDrivE documentation built on Oct. 23, 2020, 7:28 p.m.