calcGammaKernel {MGDrivE} | R Documentation |
Calculate Gamma Stochastic Matrix
Description
Given a distance matrix from calcVinEll
, calculate a
stochastic matrix where one step movement probabilities follow a gamma density.
Usage
calcGammaKernel(distMat, shape, rate)
Arguments
distMat |
Distance matrix from |
shape |
Shape parameter of |
rate |
Rate parameter of |
Details
The distribution and density functions for the gamma kernel are given below:
where is the Gamma function,
is the lower incomplete
gamma function, and
are the shape and rate parameters, respectively.
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)
# calculate gamma distribution over distances
# shape and rate are just for example
kernMat = calcGammaKernel(distMat = distMat, shape = 1, rate = 1)
[Package MGDrivE version 1.6.0 Index]