mkde2Dinteraction: Probability of 2D spatial-temporal interaction.

Description Usage Arguments Details Value Author(s) Examples

View source: R/mkdeFunctions.R

Description

Probability of 2D spatial-temporal interaction.

Usage

1
mkde2Dinteraction(mkde.obj, move.dat0, move.dat1, t.step, d.thresh)	

Arguments

mkde.obj

An MKDE object created with initialize2DMKDE

move.dat0

A move data object for the first individual created with initializeMovementData

move.dat1

A move data object for the second individual created with initializeMovementData

t.step

A time step used for numerical integration over the movement trajectory

d.thresh

The value of the kernel below which its contibrution to the overal density is considered negligible

Details

This function assumes that the two individual animals were observed at the same times. The cell values returned in the mkde.obj can be summed to obtain a global measure of spatio-temporal interaction.

Value

Returns a list with the following elements:

mkde.obj

An updated MKDE object containing the cell-level Bhattacharyya coefficients

move.dat0

A move data object for the first individuals with updated variance parameters

move.dat1

A move data object for the second individuals with updated variance parameters

Author(s)

Jeff A. Tracey, PhD
USGS Western Ecological Research Center, San Diego Field Station
[email protected]
James Sheppard, PhD
San Diego Zoo Institute for Conservation Research
[email protected]
Jun Zhu, PhD
University of Wisconsin-Madison
[email protected]

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
library(raster)
data(panda)
mv.dat0 <- initializeMovementData(panda$time, panda$x, panda$y, 
sig2obs=25.0, t.max=185.0)

n <- nrow(panda)
v <- 20.0 # increase from 0 to increase difference
mv.dat1 <- initializeMovementData(panda$time, panda$x+rnorm(n, 0, v), 
panda$y+rnorm(n, 0, v), sig2obs=25.0, t.max=185.0)

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

res <- mkde2Dinteraction(mkde.obj, mv.dat0, mv.dat1, 10.0, 1e-20)
mkde.obj <- res$mkde.obj
mv.dat0 <- res$move.dat0
mv.dat1 <- res$move.dat1

Example output

Loading required package: sp
2D MKDE Interaction Computation: STARTING
	Processing move step 1 of 146
	Processing move step 2 of 146
	Processing move step 3 of 146
	Processing move step 4 of 146
	Processing move step 5 of 146
	Processing move step 6 of 146
	Processing move step 7 of 146
	Processing move step 8 of 146
	Processing move step 9 of 146
	Processing move step 10 of 146
	Processing move step 11 of 146
	Processing move step 12 of 146
	Processing move step 13 of 146
	Processing move step 14 of 146
	Processing move step 15 of 146
	Processing move step 16 of 146
	Processing move step 17 of 146
	Processing move step 18 of 146
	Processing move step 19 of 146
	Processing move step 20 of 146
	Processing move step 21 of 146
	Processing move step 22 of 146
	Processing move step 23 of 146
	Processing move step 24 of 146
	Processing move step 25 of 146
	Processing move step 26 of 146
	Processing move step 27 of 146
	Processing move step 28 of 146
	Processing move step 29 of 146
	Processing move step 30 of 146
	Processing move step 31 of 146
	Processing move step 32 of 146
	Processing move step 33 of 146
	Processing move step 34 of 146
	Processing move step 35 of 146
	Processing move step 36 of 146
	Processing move step 37 of 146
	Processing move step 38 of 146
	Processing move step 39 of 146
	Processing move step 40 of 146
	Processing move step 41 of 146
	Processing move step 42 of 146
	Processing move step 43 of 146
	Processing move step 44 of 146
	Processing move step 45 of 146
	Processing move step 46 of 146
	Processing move step 47 of 146
	Processing move step 48 of 146
	Processing move step 49 of 146
	Processing move step 50 of 146
	Processing move step 51 of 146
	Processing move step 52 of 146
	Processing move step 53 of 146
	Processing move step 54 of 146
	Processing move step 55 of 146
	Processing move step 56 of 146
	Processing move step 57 of 146
	Processing move step 58 of 146
	Processing move step 59 of 146
	Processing move step 60 of 146
	Processing move step 61 of 146
	Processing move step 62 of 146
	Processing move step 63 of 146
	Processing move step 64 of 146
	Processing move step 65 of 146
	Processing move step 66 of 146
	Processing move step 67 of 146
	Processing move step 68 of 146
	Processing move step 69 of 146
	Processing move step 70 of 146
	Processing move step 71 of 146
	Processing move step 72 of 146
	Processing move step 73 of 146
	Processing move step 74 of 146
	Processing move step 75 of 146
	Processing move step 76 of 146
	Processing move step 77 of 146
	Processing move step 78 of 146
	Processing move step 79 of 146
	Processing move step 80 of 146
	Processing move step 81 of 146
	Processing move step 82 of 146
	Processing move step 83 of 146
	Processing move step 84 of 146
	Processing move step 85 of 146
	Processing move step 86 of 146
	Processing move step 87 of 146
	Processing move step 88 of 146
	Processing move step 89 of 146
	Processing move step 90 of 146
	Processing move step 91 of 146
	Processing move step 92 of 146
	Processing move step 93 of 146
	Processing move step 94 of 146
	Processing move step 95 of 146
	Processing move step 96 of 146
	Processing move step 97 of 146
	Processing move step 98 of 146
	Processing move step 99 of 146
	Processing move step 100 of 146
	Processing move step 101 of 146
	Processing move step 102 of 146
	Processing move step 103 of 146
	Processing move step 104 of 146
	Processing move step 105 of 146
	Processing move step 106 of 146
	Processing move step 107 of 146
	Processing move step 108 of 146
	Processing move step 109 of 146
	Processing move step 110 of 146
	Processing move step 111 of 146
	Processing move step 112 of 146
	Processing move step 113 of 146
	Processing move step 114 of 146
	Processing move step 115 of 146
	Processing move step 116 of 146
	Processing move step 117 of 146
	Processing move step 118 of 146
	Processing move step 119 of 146
	Processing move step 120 of 146
	Processing move step 121 of 146
	Processing move step 122 of 146
	Processing move step 123 of 146
	Processing move step 124 of 146
	Processing move step 125 of 146
	Processing move step 126 of 146
	Processing move step 127 of 146
	Processing move step 128 of 146
	Processing move step 129 of 146
	Processing move step 130 of 146
	Processing move step 131 of 146
	Processing move step 132 of 146
	Processing move step 133 of 146
	Processing move step 134 of 146
	Processing move step 135 of 146
	Processing move step 136 of 146
	Processing move step 137 of 146
	Processing move step 138 of 146
	Processing move step 139 of 146
	Processing move step 140 of 146
	Processing move step 141 of 146
	Processing move step 142 of 146
	Processing move step 143 of 146
	Processing move step 144 of 146
	Processing move step 145 of 146
	Processing move step 146 of 146
	Normalizing Constant = 6.12647e-05(16322.6)
	Maximum cell Bhattacharyya coeff = 0.0150099
	Overall Bhattacharyya coeff = 0.880828
2D MKDE Interaction Computation: DONE

mkde documentation built on May 2, 2019, 6:46 a.m.