Kernel Home Range Estimation for ArcGIS, using VBA - Fish and ...
Kernel Home Range Estimation for ArcGIS, using VBA - Fish and ...
Kernel Home Range Estimation for ArcGIS, using VBA - Fish and ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
⎛ − x'<br />
x ⎞<br />
⎜ ⎟<br />
⎠<br />
⎝ 2<br />
e<br />
K ( x)<br />
=<br />
(10)<br />
2π<br />
Where x’x is the distance between the evaluation point <strong>and</strong> another point in the distribution,<br />
divided by the smoothing parameter (x’x = distanceP P/hP P). Given equation 9, the convolution of the<br />
multivariate normal density function, KP<br />
(2)<br />
P, becomes:<br />
2<br />
2<br />
⎛ − x'<br />
x ⎞<br />
⎜ ⎟<br />
⎝ 4 ⎠<br />
( 2 ) e<br />
K ( x)<br />
=<br />
(11)<br />
4π<br />
The least-squares cross-validation score now becomes (Worton, 1995; Rodgers <strong>and</strong> Carr, 1998):<br />
M<br />
1<br />
∑∑<br />
i<br />
j<br />
⎛ ⎛ − x'<br />
x ⎞<br />
⎛ − x'<br />
x ⎞<br />
⎛ ⎜ ⎟ ⎞ ⎛ ⎜ ⎟ ⎞⎞<br />
⎜⎜<br />
⎝ 4 ⎠<br />
2<br />
e ⎟ ⎜ ⎝ ⎠<br />
e ⎟⎟<br />
⎜⎜<br />
⎟ − 2⎜<br />
⎟⎟<br />
4π<br />
2π<br />
⎜⎜<br />
⎟ ⎜ ⎟<br />
⎟<br />
⎝⎝<br />
⎠ ⎝ ⎠⎠<br />
n h<br />
1<br />
πnh<br />
( h) =<br />
+<br />
2 2<br />
2<br />
(12)<br />
ABODE uses the minimization process, “Routine GOLDEN” (Sprott, 1991) to minimize MB1B(h), <strong>and</strong><br />
hence find the optimum value <strong>for</strong> h <strong>for</strong> a bivariate normal density. The seed value used in the<br />
minimization process is the hBrefB value (Equation 3), which is subsequently multiplied by 0.01 to get<br />
the lower bound of the search values. For unit variance st<strong>and</strong>ardization, hBrefB is calculated from<br />
the st<strong>and</strong>ardized dataset (i.e. σ = 1). When it is to be used as the smoothing parameter, the hBrefB<br />
value is adjusted <strong>for</strong> the biweight kernel. A(K) is the constant used to make the adjustment. For<br />
conversion to a biweight kernel, A(K) = 2.04 (Silverman, 1986). As stated earlier, this<br />
discrepancy is unexplained.<br />
h opt<br />
A(<br />
K)<br />
σ<br />
= (13)<br />
6<br />
n<br />
Once an h has been estimated from LSCV, this value is adjusted <strong>for</strong> the biweight kernel by<br />
multiplying by the constant A(K) = 2.04 (<strong>for</strong> use as a biweight kernel in the actual density<br />
estimation) (Seaman <strong>and</strong> Powell, 1996; Silverman, 1986:87). It should be noted that ABODE<br />
follows Silverman (1986) by <strong>using</strong> 2.04 as the constant, <strong>and</strong> not the 2.78 reported by Seaman<br />
<strong>and</strong> Powell (1996).<br />
20