This function creates a utilization distribution according to the Bivariate Gaussian Bridge model.
# S4 method for dBGBvariance,RasterLayer,numeric
dynBGB(move, raster, locErr, timeStep, ...)
# S4 method for .MoveTrackSingle,RasterLayer,numeric
dynBGB(move, raster, locErr, margin, windowSize, ...)
# S4 method for .MoveTrackSingle,numeric,ANY
dynBGB(move, raster, locErr, ext, ...)
# S4 method for .MoveTrackSingle,missing,ANY
dynBGB(move, raster, locErr, dimSize, ext, ...)
It returns an object of the class dynBGB-class.
a move
or dBGBvariance
object. This object must be in a projection different to longitude/latitude (one suitable for euclidean geometry), use spTransform
to transform your coordinates.
a RasterLayer
object or a numeric value. If a RasterLayer
is provided the dynBGB
starts to calculate the UD based on that raster. If a numeric value is provided it is interpreted as the resolution of the square raster cells (in map units); the according raster will be calculated internally.
single numeric value or vector of the length of coordinates that describes the error of the location (sender/receiver) system in map units. Or a character string with the name of the column containing the location error can be provided.
It correspond to the size of the timer intervals taken for every integration step (in minutes). If left NULL
20.1 steps are taken in the shortest time interval. See 'Details'. Optional.
The margin used for the behavioral change point analysis. This number has to be odd.
The size of the moving window along the track. Larger windows provide more stable/accurate estimates of the brownian motion variance but are less well able to capture more frequent changes in behavior. This number has to be odd.
Describes the amount of extension of the bounding box around the animal track. It can be numeric (same extension into all four directions), vector of two (first x, then y directional extension) or vector of four (xmin, xmax, ymin, ymax extension). Only considered in combination with a numeric raster argument or the dimSize argument.
numeric. dimSize is only used if raster
is not set. dimSize
is interpreted as the number of cells along the largest dimension of the track. The according raster will be calculated internally.
Currently no other arguments implemented.
Bart Kranstauber & Anne Scharf
There are three ways to launch the dynBGB
function:
1. Use a raster:
A RasterLayer object is set for the raster argument which is then used to calculate the UD.
(needed arguments: move, raster(=RasterLayer
), locErr, margin, windowSize; optional arguments: timeStep)
2. Set the cell size
To set the cell size, set a numeric value for the raster argument which is used as the cell sizes of the raster.
(needed arguments: move, raster(=numeric
), locErr, margin, windowSize, ext; optional arguments: timeStep)
3. Set the number of cells (col/row)
To set the number of cells along the largest dimension a numeric dimSize argument can be set.
(needed arguments: move, dimSize, locErr, margin, windowSize, ext; optional arguments: timeStep)
timeStep
. The default value is the shortest time interval divided by 20.1. This means, if there is a location recorded e.g. every 40 mins, the function divides each segment into 1.99 mins chunks upon which it does the calculation. If for some reason there is one time interval of 20 secs in the track, each segment of the track will be divided into 1secs chunks, increasing the calculation time immensely. Before calculating the DBBMM, use e.g. min(timeLag(x=
myMoveObject, units="mins"
)) to check which is the duration of the shortest time interval of the track. If the track contains time intervals much shorter than the scheduled on the tag, set the timeStep
e.g. to the scheduled time interval at which the tag was set, divided by 20.1.
Kranstauber, B., Safi, K., Bartumeus, F.. (2014), Bivariate Gaussian bridges: directional factorization of diffusion in Brownian bridge models. Movement Ecology 2:5. doi:10.1186/2051-3933-2-5.
dynBGBvariance, getMotionVariance, getVolumeUD, contour, outerProbability, raster, raster2contour, brownian.bridge.dyn, brownian.motion.variance.dyn
data(leroy)
leroy <- leroy[230:265,]
## change projection method to aeqd and center the coordinate system to the track
dataAeqd <- spTransform(leroy, CRSobj="+proj=aeqd +ellps=WGS84", center=TRUE)
dBGB <- dynBGB(dataAeqd, locErr=9, raster=10, ext=0.5, windowSize=31, margin=15, timeStep=15/20.1)
plot(dBGB, col=hsv(sqrt(1:700/1000)))
lines(dataAeqd)
Run the code above in your browser using DataLab