calcZeroInflation.Rd
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.
calcZeroInflation(stayThroughLifespanProbability, adultMortality)
stayThroughLifespanProbability | Probability of a mosquito to spend its whole lifespan in the same node |
---|---|
adultMortality | Adult mortality rate |
# 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)