# calcLognormalKernel: Calculate Lognormal Stochastic Matrix In MGDrivE: Mosquito Gene Drive Explorer

## Description

Given a distance matrix from calcVinEll, calculate a stochastic matrix where one step movement probabilities follow a lognormal density.

## Usage

 1 calcLognormalKernel(distMat, meanlog, sdlog) 

## Arguments

 distMat Distance matrix from calcVinEll meanlog Log mean of Lognormal distribution sdlog Log standard deviation of Lognormal distribution

## Details

The distribution and density functions for the lognormal kernel are given below:

F(x)=\frac{1}{2} + \frac{1}{2} \mathrm{erf}[\frac{\mathrm{ln}x-μ}{√{2}σ}]

f(x)=\frac{1}{xσ√{2π}}\mathrm{exp}≤ft( -\frac{(\mathrm{ln}x-μ)^{2}}{2σ^{2}} \right)

where μ is the mean on the log scale, and σ is the standard deviation on the log scale.

## Examples

  1 2 3 4 5 6 7 8 9 10 11 # 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) # calculate lognormal distribution over distances # mean and standard deviation are just for example kernMat = calcLognormalKernel(distMat = distMat, meanlog = 100, sdlog = 10) 

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