plotMKDE: Make an image plot of an MKDE.

In mkde: 2D and 3D Movement-Based Kernel Density Estimates (MKDEs)

Make an image plot of an MKDE.

Description

Makes an image plot of an MKDE. If the MKDE is 3D, the contours will be based on the entire MKDE, but only one level indexed by z.index will be plotted.

Usage

plotMKDE(mkde.obj, z.index=1, probs=c(0.99, 0.95, 0.90, 0.75, 0.5, 0.0),
cmap=rev(rainbow(length(probs)-1)), add=FALSE, ...)	

Arguments

mkde.obj	A 2D, 2.5D, or 3D MKDE object created with initialize2DMKDE or initialize3DMKDE and density initialized with initializeDensity
z.index	Index for the z-dimension of the density if a 3D MKDE
probs	Probabilities for image contours.
cmap	Color map for image plot.
add	FALSE to make a new plot, TRUE to add to existing plot
...	Additional graphical parameters.

Details

A plot of the density for the specified level in the z-dimension (which should be 1, the defaul value, for a 2D or 2.5D MKDE) is generated.

Value

No value is returned.

Author(s)

Jeff A. Tracey, PhD
USGS Western Ecological Research Center, San Diego Field Station
[e-mail]
James Sheppard, PhD
San Diego Zoo Institute for Conservation Research
[e-mail]

Examples

library(raster)
data(pandadem)
cell.sz <- 0.25*mean(res(pandadem))
ext <- extent(pandadem)
nx <- 4*ncol(pandadem)
ny <- 4*nrow(pandadem)
mkde.obj <- initializeMKDE2D(ext@xmin, cell.sz, nx, ext@ymin, cell.sz,
ny)

# set up movement data
data(panda)
mv.dat <- initializeMovementData(panda$time, panda$x, panda$y, 
sig2obs=25.0, t.max=185.0)

# estimate density
dens.res <- initializeDensity(mkde.obj, mv.dat)
mkde.obj <- dens.res$mkde.obj
mv.dat <- dens.res$move.dat
plotMKDE(mkde.obj)

mkde documentation built on July 9, 2023, 6:41 p.m.