move (version 1.0)

brownian.bridge.dyn: Creates a DBBMM object

Description

The brownian.bridge.dyn function uses a Move object (see Move-class) to calculate the utilization distribution, UD, of the given track. It uses the dynamic Brownian Bridge Movement Model (dBBMM) to do so. The dBBMM has the advantage over the other Brownian Bridge Movement Model that changes in behavior are accounted for. It does so by using the behavioral change point analysis in a sliding window. For details see references.

Usage

brownian.bridge.dyn(object, raster, dimSize, location.error,
                 margin=11, window.size=31, ext=.3, bbox=NA,...)

Arguments

object

an object of the Move-class

raster

a RasterLayer object or numeric value. A numeric value for raster is interpreted as the edge size of the square raster cells (in map units); the according raster will be calculated internally. If a RasterLayer is provided the brownian.b

dimSize

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. Default is 10

location.error

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.

margin

The margin used for the behavioral change point analysis.

window.size

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.

ext

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). Defaul

bbox

vector with 4 numbers defining a bounding box for the raster

...

for additional arguments to program against

Details

There are four ways to launch the brownian.bridge.dyn function which are as follows: 1. Use a raster A RasterLayer object is set for the raster argument which is then used to calculate the UD. 2. Set the cell size To set the cell size, set a numeric value for the raster argument without providing dimSize. The numeric raster argument is used as the cell sizes of the raster. 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. 4. Using default raster When there are no values set, the default raster value is used to calculate and create a RasterLayer object, which is returned to the same function. Note: depending on the size of the area of interest, the default cell size value can result in a large number of cells which may take a very long time to calculate!

The function prints an estimate of the size of the computational task ahead. This can give an indication of how long the computation is going to take. It should scale roughly linearly with the duration of the computations. In our experience 10e9 takes about a minute with an average laptop.

There is one further argument that can be given: time.step. It correspond to the size of the timer intervals taken for every integration step (in minutes). If left NULL 15 steps are taken in the shortest time interval.

References

Kranstauber, B., Kays, R., LaPoint, S. D., Wikelski, M. and Safi, K. (2012), A dynamic Brownian bridge movement model to estimate utilization distributions for heterogeneous animal movement. Journal of Animal Ecology. doi: 10.1111/j.1365-2656.2012.01955.x

Examples

## create a move object
data <- move(system.file("extdata","leroy.csv.gz",package="move"))

## change projection method to aeqd and center the coordinate system to the track
data2 <- spTransform(data, CRSobj="+proj=aeqd", center=TRUE)

## create a DBBMM object
dbbmm <- brownian.bridge.dyn(object=data2[1:100,], location.error=12, dimSize=45, ext=.3, time.step=600)
load(system.file("extdata", "move.RData", package="move"), .GlobalEnv)
  leroy <- leroy[1:50]
  ricky <- ricky[1:50]
  stack <- moveStack(list(leroy,ricky))
  t <- brownian.bridge.dyn(spTransform(stack), location.error=23.5, dimSize=25, time.step=6000, ext=1.2)