01.01.2015 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!