# Correction to: A cross-validation-based approach for delimiting reliable home range estimates

- Eric R. Dougherty
^{1}Email author, - Colin J. Carlson
^{1}, - Jason K. Blackburn
^{2, 3}and - Wayne M. Getz
^{1, 4}

**Received: **7 November 2017

**Accepted: **7 November 2017

**Published: **5 December 2017

The original article was published in Movement Ecology 2017 5:19

## Correction to: Movement Ecology (2017) 5:19 DOI:10.1186/s40462-017-0110-4

### Original text

*k*and

*s*values on the grid. This analysis entailed the creation of local convex hulls with

*k*nearest neighbors and a scaling factor of

*s*. In all subsequent analyses, we assume that the scaling of time follows a linear formulation; however, when movement patterns more closely exemplify diffusion dynamics, an alternative equation for the TSD may be more accurate [1]. The test points were then laid upon the resulting hulls, and the probability of each was calculated as the proportion of the total number of hulls (equivalent to the total number of points in the training dataset) that contained the test point (Figure 1). Test points that were not contained within any hulls were assigned a probability equal to the inverse of the total number of points in the full movement path divided by 100, effectively penalizing any hull sets that did not include each of the test points. Though an arbitrary selection, the choice of a consistent penalty term across individuals will serve to standardize the procedure. A larger penalty will likely result in a higher optimal

*k*value and bear a closer resemblance to the MCP. The natural log of the probability was calculated and information criterion values analogous to Akaike’s Information Criterion (AIC) were derived using the equation:

The choice of 2*k* as the penalty term was made to maintain a structure analogous to the AIC equation. Given the expansive literature concerning the performance and behavior of AIC under various scenarios, maintaining this structure may offer insight into similar strengths and weaknesses of the proposed approach. Ultimately, without such a penalty, all movement paths would tend towards a *k* equal to the number of points in the training set, such that each individual point was assigned a probability of one. It should be noted that this penalty term is specific to the *k* (nearest neighbors) method, but the underlying cross-validation procedure could very easily be extended for the optimization of the *a* (adaptive parameter) method if an appropriate penalty term is selected. An ideal penalty term would likely result in a increase of the information criterion value by a similar magnitude as in the *k*-based formulation above (i.e., ranging from approximately 10^{0} to 10^{2}).

### Revised text

*k*and

*s*values on the grid. This analysis entailed the creation of local convex hulls with

*k*nearest neighbors and a scaling factor of

*s*. In all subsequent analyses, we assume that the scaling of time follows a linear formulation; however, when movement patterns more closely exemplify diffusion dynamics, an alternative equation for the TSD may be more accurate [1]. The test points were then laid upon the resulting hulls, and the probability of each was calculated as the proportion of the total number of hulls (equivalent to the total number of points in the training dataset) that contained the test point (Figure 1). Test points that were not contained within any hulls were assigned a probability equal to the inverse of the total number of points in the full movement path divided by 100, effectively penalizing any hull sets that did not include each of the test points. Though an arbitrary selection, the choice of a consistent penalty term across individuals will serve to standardize the procedure. A larger penalty will likely result in a higher optimal

*k*value and bear a closer resemblance to the MCP. The natural log of the probability was calculated and information criterion values analogous to the

**Bayesian**Information Criterion

**(BIC)**were derived using the equation:

The choice of k
**∗ ln(**
P
**) as the overall penalty term** was made to maintain a structure analogous to the **BIC** equation. Given the expansive literature concerning the performance and behavior of **BIC** under various scenarios, maintaining this structure may offer insight into similar strengths and weaknesses of the proposed approach. Ultimately, without such a penalty, all movement paths would tend towards a *k* equal to the number of points in the training set, such that each individual point was assigned a probability of one. **An alternative method akin to Akaike’s Information Criterion can also be applied, but the penalty term (2∗**
k
**) does not scale with the total number of test points (in turn, a function of the total length of the movement path) and will likely result in higher optimal**
k
** values than the BIC analogue.** It should also be noted that this penalty term is specific to the *k* (nearest neighbors) method, but the underlying cross-validation procedure could very easily be extended for the optimization of the *a* (adaptive parameter) method if an appropriate penalty term is selected. An ideal penalty term would likely result in a increase of the information criterion value by a similar magnitude as in the *k*-based formulation above (i.e., ranging from approximately 10^{0} to **10**
^{
3
}).

### Explanation of correction

The only difference between the equation here and the one in the original article is the penalty term. In the equation above, increases in the *k* value are penalized more heavily than the simpler 2∗*k* term. The additional benefit of this equation, and the primary reason for its use in the analysis in [2], is that the penalty term scales (in a non-linear fashion) with the total number of test points, offering more flexibility when considering trajectories of varying lengths.

Despite this issue, the fundamental principles underlying the cross-validation method remain sound, and both the original IC equation and the one presented here can be used with confidence. The logic for utilizing a BIC analogue is the same as that for formulating an AIC analogue; the correction outlined here simply enables the replication of the results in the article. The equation and the text in bold above have been altered from the original version of the paper.

## Notes

## Declarations

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Lyons AJ, Turner WC, Getz WM. Home range plus: a space-time characterization of movement over real landscapes. Mov Ecol. 2013; 1(1):2.View ArticlePubMedPubMed CentralGoogle Scholar
- Dougherty ER, Carlson CJ, Blackburn JK, Getz WM. A cross-validation-based approach for delimiting reliable home range estimates. Mov Ecol. 2017; 5(1):19.View ArticlePubMedPubMed CentralGoogle Scholar