 Research
 Open Access
 Published:
Modified home range kernel density estimators that take environmental interactions into account
Movement Ecology volume 7, Article number: 16 (2019)
Abstract
Background
Kernel density estimation (KDE) is a major tool in the movement ecologist toolbox that is used to delineate where geotracked animals spend their time. Because KDE bandwidth optimizers are sensitive to temporal autocorrelation, statisticallyrobust alternatives have been advocated, first, datathinning procedures, and more recently, autocorrelated kernel density estimation (AKDE). These yield asymptotically consistent, but very smoothed distributions, which may feature biologically unrealistic aspects such as spilling beyond impassable borders.
Method
I introduce a semiparametric variant of AKDE designed to extrapolate more realistic home range shapes by incorporating movement mechanisms into the bandwidth optimizer and into the base kernels. I implement a first approximative version based on the step selection framework. This method allows accommodating land cover selection, permeability of linear features, and attraction for select landscape features when delineating home ranges.
Results
In a plains zebra (Equus quagga), the reluctance to cross a railway, the avoidance of dense woodland, and the preference for grassland when foraging created significant differences between the estimated home range contours by the new and by previous methods.
Conclusion
There is a tradeoff to find between fully parametric density estimators, which can be very realistic but need to be provided with a good model and adequate environmental data, and nonparametric density estimators, which are more widely applicable and asymptotically consistent, but whose details are bandwidthlimited. The proposed semiparametric approach attempts to strike this balance, but I outline a few areas of future improvement. I expect the approach to find its use in studies that compare extrapolated resource availability and interpolated resource use, in order to discover the movement mechanisms that we need to improve the extrapolations.
Introduction
Many researchers use kernel density estimators (KDE) to extrapolate where a geotracked animal spends its time, often using the 95% extrapolated isopleth as the home range contour [1,2,3]. KDE work by approximating the stationary utilization distribution p(r) of the animal, i.e., its time budget with respect to location r, with a sum of “kernels”, i.e., unimodal distributions κ centered around each recorded location {r_{i}} [4, 5].
The parameter σ_{B}, termed the bandwidth, controls the spread of each kernel around each recorded location, and therefore ultimately the degree of smoothing of the resulting distribution [4, 6, 7]. A small bandwidth yields a distribution with numerous peaks around each cluster of recorded locations; a large bandwidth smooths out these peaks and yields a more spreadout distribution [2, 5]. The choice of an appropriate bandwidth is therefore critical, and indeed usually trumps the influence of the actual shape of the kernels, i.e., the analytical form of function κ [4]. At least three categories of optimal bandwidth selection routines have been developed to inform and automatize this decision [8], but, because of temporal autocorrelation in movement data [9,10,11], only one, the reference function approximation, is recommended for animal tracking applications [12]. In the reference function approximation approach, one optimizes the bandwidth by minimizing an approximated mean integrated square error criterion (MISE).
p_{REF}(r) is the reference function, usually chosen for its mathematical properties. In a nonapproximate MISE, there should be a p(r) term instead of p_{REF}(r), but, since p(r) is the unknown we want to estimate, we need to replace it with the reference function. Ω represents the domain of possible locations, the ∫ ∙ dr notation denotes integration over one realization of the stochastic movement process in space, and the 〈∙〉 notation denotes integration over all realizations of the movement process.
If Eq. 2 did not account for the temporal autocorrelation structure of the movement process, it would introduce a bias in the bandwidth estimate, that can severely impedes comparative inference [7, 9, 11,12,13,14,15]. For example, the optimizer would converge towards a zero bandwidth when the sampling resolution would increase or the amount of temporal autocorrelation in the animal movement would increase, yielding increasingly smaller home ranges [12]. The amount of bias depends on both the sampling resolution and the movement rates of the animal. There are two known ways to deal with temporal autocorrelation in Eq. 2. First, one may subsample the data so that successive records are approximatively independent [9, 16]. The recommended best practice is to keep one record every 3τ where τ is the autocorrelation time (the rate at which the proximity between two records declines with the time lag between them). I hereafter refer to this as “robust KDE” (KDEr). The other way to deal with temporal autocorrelation is to keep all the data in, but specify a temporal autocorrelation model in the kernels that make up \( \widehat{p} \) and in the reference function, yielding a different MISE that is minimized by a different value of the bandwidth [12]. This approach is termed “autocorrelated kernel density estimation” (AKDE). In both cases, the approximation of p(r) by p_{REF}(r) still leads to a “reference function approximation bias” [17]. This bias is usually positive. It can be corrected a posteriori [17], in which case I use the recommended notation “c”, e.g., AKDEc.
Both these robustizing protocols increase the bandwidth compared to standard or naive KDE. The KDEr option also requires the user to discard potentially massive amounts of data. As a consequence, the resulting distribution, although statistically robust and asymptotically consistent, may look oversmoothed and biologically irrelevant [16, 18]. For example, the 95% isopleth might intersect impassable barriers such as coastlines [19]. The issue of boundaries is indeed recurrent in KDE applications, including in other fields than animal tracking data analysis (review in [20, 21]). For home range estimation, the most common response is to clip the distributions at known barriers [19]. This intuitive practice is equivalent to introducing a dose of mechanism in the otherwise fully nonparametric KDE methodology. By specifying where to clip and how to redistribute the weight, we in effect inform the model that the barrier was unpassable. However, instead of providing that information in an ad hoc way at the very end of the process, we could uptake it from the start. The MISE would then become less sensitive to an increase in bandwidth that would otherwise have caused the kernel distributions to spill beyond the barrier, yielding a different optimal bandwidth. We would also incorporate the barrier into the reference function, making it more realistic, and thereby suppressing the reference function approximation bias. Lastly, we would apply the effect of the barrier to each kernel, yielding an estimated utilization distribution that is truncated at the barrier by construct. More generally, in addition to barriers, we can incorporate in this way any movement mechanism that can be formalized using a step selection function [22, 23], e.g., land cover selection [22] and permeability of linear features like roads [24].
Material and methods
Step 1: fitting a mechanistic movement model
At time t, the position r_{t} of the animal was assumed to be drawn from a step selection kernel g_{u} made of the product of an availability function g_{a}, conditioned on the movement path prior to t, R_{t − 1}, and especially the last known position, r_{t − 1}, and of a weight function W, both defined over the movement domain Ω [22, 23, 25].
K_{t} is a scaling constant so that g_{u} sums to one.
The availability function g_{a} was modelled using the OrnsteinUhlenbeck process (OU), a continuoustime stochastic movement model that represents home range behavior as a tendency to revert back to a mean location following random deviations from that mean [26]. The weight function W described environmental interactions: selection of some land covers over others, attraction or repulsion towards fixed landscape features such as human settlements, and barrier permeability, i.e., the rate at which animals avoid crossing linear features such as roads and rivers [25]. Following established practice [22, 27, 28], the analytical form of the weight function W was:
x(r_{t}) describes the environment at location r_{t}. The k^{th} elements of x(r_{t}) (k = 1, … K_{1}) contains a 0 or 1 coding for the presence of land cover type k at location r_{t}. The next K_{2} elements contain the distances from r_{t} to fixed landscape features (or the value of continuous environmental covariates such as climate or vegetation density). δ(r_{t} r_{t − 1}) is a vector containing a 1 in the l^{h} cell if a barrier of type l is crossed when going straight from r_{t − 1} to r_{t}. Note that this straightline permeability model is valid only for small steps that make longer detours extremely unlikely.
The model parameters are included in vectors α, the selection coefficients, and λ, the permeability coefficients. The kth element of α (1 < k ≤ K_{1}) quantifies how much land cover type k is preferred over land cover type 1. The lth element of λ quantifies the reluctance to cross linear features of type l, zero meaning that the linear feature has no effect on movement. Both α and λ were considered constant through time and space in the present application.
For the sake of simplicity, and because my focus here was on the postfitting treatment of α and λ estimates, rather than on the estimation itself, I used a relatively fast but approximate procedure to estimate the parameters in W. Following Johnson, Hooten & Kuhn [29], I reformulated the movement model as a spacetime point process [30, 31]. This meant that the availability function g_{a} was approximated by a Brownian availability window at each time step. However, in the postfitting treatment, I used the OrnsteinUhlenbeck model for g_{a} as announced above. This means that different contradictory models were in practice used to estimate the weight function and the availability function. The detail of the spacetime point process implementation is provided in Additional file 2.
Step 2: bandwidth optimization
To incorporate environmental interactions, I replaced the kernels of Eq. 1 by weighed multivariate Gaussian distributions [20, 21].
φ(r, r_{i}, σ_{B}) denotes the multivariate Gaussian distribution of mean r_{i} and variancecovariance matrix σ_{B}, and \( {K}_i^B \) is a scaling constant so that each kernel sums to one (more details in Additional file 1). Following Fleming et al. [12], I simplified Eq. 5 using σ_{B} = σ_{B} ∙ σ_{0}, where σ_{0} is the variancecovariance matrix of the availability function. This means that the direction of the smoothing was driven by the anisotropy of the movement process.
Importantly, these kernels feature a simplistic permeation model that imply a straightline move from r to r_{i} (Eq. 4). In the zebra case study below, the violation of that straightline assumption was almost without consequence, because the σ_{B} value was moderate and the linear feature exhibited no convolution at all, meaning that even if the path from r to r_{i} was not straight, the chance to cross the linear feature was similar to that of a straight move. However in other applications, researchers may need to consider alternative jackknifing methods, e.g., diffusive permeation kernels [32]. In addition, with Eq. 5 we assume a regular redistribution of the discounted weight across the whole availability domain. In applications where the animals remain near the linear features when they encounter them [19], alternative redistribution rules may be warranted, e.g., reflected kernels [20].
Next, for the reference function, I also used weighted multivariate Gaussian distributions, summed over recorded locations to represent the expected equilibrium distribution of the movement process.
\( {K}_i^0 \) is a scaling constant so that each element of the reference function sums to one.
Like Eq. 5, Eq. 6 features a simplistic permeation model that assumes a straight line move from μ_{0} to r. However, now that simplistic permeation model is applied to the whole home range. Any movement bottleneck, e.g., a constrained corridor between two sections of the home range, would be overly discounted in the resulting reference function. In the zebra case study, we did not have to deal with any such feature. However, in other applications, users may have to consider alternative formulations. Individualbased simulations (IBS; see “Zebra case study” below) could in this case prove particularly useful to generate a “pilot estimate” to use instead of Eq. 6. In addition to offering a straightforward way to incorporate abovementioned movement bottlenecks into the reference function, IBS can be set up so that the simulation step length is short enough that the straightline permeability model remains realistic at all stages of the simulations. However, this option is not implemented yet, and probably warrants further investigation pertaining to sensitivity to simulation parameters (simulated duration, step lengths, etc.).
When using Eq. 5, 6, there is an analytical form for the MISE, derived in Additional file 1. These details expand on the proof of Fleming et al. [12] demonstrating how temporal autocorrelation is incorporated into the KDE framework. The resulting MISE optimization algorithm is however prohibitively timeconsuming. Thus, pending algorithmic improvement, the new approach remains mostly theoretical and exploratory. I however developed a fasterrunning simplified version in the following section.
Simplified version
Because the MISE optimization was so prohibitively timeconsuming, I simplified the protocol by bypassing Step 2 entirely, or more precisely, replacing Step 2 with the corresponding step of the AKDE analytical protocol [33]. In other words, the bandwidth is optimized while taking temporal autocorrelation into account, but without taking environmental interactions into account. The weights from Step 1 are then only applied when eventually computing the distribution (Step 3). I propose the notation EAKDE for the full version where environmental interactions are incorporated in Step 2, and SEAKDE for the simplified version where environmental interactions are not incorporated at Step 2. Importantly, because of the assumption from the straightline permeability model, EAKDE is not compulsorily always more reliable than SEAKDE.
Step 3: computing the kernel density and correcting for the remaining reference function approximation bias
From the estimated σ_{B}, I computed isopleths of \( \widehat{p}\left(\boldsymbol{r}\right) \) using Eq. 5. I then applied the reference function bias correction routine of Fleming and Calabrese [17] to the isopleths, but only when implementing SEAKDE, hence the SEAKDEc notation hereafter. When implementing EAKDE, I considered that by changing the reference function (Eq. 6), I got rid of the reference function approximation bias. This is certainly a strong assumption that was however supported by the empirical results.
Zebra case study and comparison with alternatives to KDE
I analyzed data collected from a plains zebra (Equus quagga) in and near Hwange national park, Zimbabwe (26.861E, − 18.624 N). The study individual (individual local identifier: Ganda) was monitored between Jan 2011 and Sept 2012 with a collarmounted GPS that recorded one location every hour [34]. I rescaled the recent Hwange vegetation map [35] at a 150 m resolution and pooled vegetation classes into 4 categories to reduce computing time for this illustration case. Other landscape features known to influence zebra space use included water holes, a railway with adjacent road that marks the northern border of the park, and a town. These were all included in the step selection model (Eq. 3).
I compared four variants of the KDE methodology: a naïve reference functionbased bandwidth optimizer that did not account for temporal autocorrelation (KDE), the robustized approach where the data were subsampled before analysis (KDEr), the AKDE approach where the autocorrelation structure was incorporated in the reference function (AKDEc), and finally the new methods (EAKDE and SEAKDEc). For KDE and KDEr I used the kde2d function in Rpackage MASS. The position autocorrelation time τ required to subsample the data for KDEr was estimated from the ctmm.fit routine in Rpackage ctmm [33]. For AKDEc, I used the akde function in ctmm within the recommended analytical protocol [33]. For EAKDE, the algorithm and mathematical justification are described (with words) in Additional file 1. For SEAKDEc, I used the AKDE bandwidth but then applied the environmental interaction weights when computing the distribution.
I also implemented three nonKDE methodologies. First, I used the asymptotic distribution of the fitted OrnsteinUhlenbeck model to draw the ellipses that most closely approximated the home range and core area. Second, I implemented the movementbased kernel density estimator (MKDE), in which the kernels are replaced by step selection functions, i.e., the g_{u}(r_{i + 1} R_{i}) estimated at Step 1 [19, 36]. MKDE, like the Brownian bridge, computes the probability that a location was used between records [36]. The method focuses on one realization of the movement path and on the process uncertainty around the interpolated path between records [37]. By contrast, KDE extrapolators average the utilization distribution across realizations of the movement path. Lastly, I implemented an individualbased simulation procedure (IBS) [38, 39] to generate 1000 1monthlong tracks with one record per hour, each track starting from a randomly selected recorded location, and moving stochastically according to the model described in Eq. 1. This yielded a cloud of 720,000 simulated records, from which I computed the density of records per pixel, thereby obtaining a rasterized cumulative utilization distribution that quantified the time budget under the fitted mechanistic movement model.
For each estimator, I computed the home range area (95% isopleth) and the core area (50% isopleth), as well as the home range scale computed as the root mean square distance of the distribution to its centroid. I also computed the “amplitude” of the core area and of the home range as the longest distance between two points on the isopleth. These represent different ways to measure the home range. In particular, the home range area of the asymptotic OU distribution is proportional to the movement variance or home range scale.
Results and discussion
Comparison between KDE variants in the zebra case study
The standard KDE yielded the smallest home range area and amplitude among all KDE variants (Fig. 1 and Table 1). As reviewed in the introduction, this small estimated home range size is partly caused by an unwanted property of the standard KDE in the presence of temporal autocorrelation [7, 9, 11,12,13,14,15]. The KDEr version indeed yielded a much larger estimated home range. AKDEc yielded a smaller home range estimate than KDEr, with a notably smaller core area leading to a large estimated home range scale. Neither AKDEc nor KDEr provided any information about landcover selection or reluctance to cross the railway (Fig. 1), as expected by construction. The most visually compelling effect of using SEAKDEc and EAKDE was that the space to the east of the railway was weighted down compared to AKDEc and KDEr. In addition, the core area was markedly irregular in shape, reflecting the avoidance of the densest woodland cover type and the preference for pure grassland. The composition of the home range remained similar across all methodologies. In particular, while c.30% of the raw data was recorded in grassland, only 1–6% of the estimated home range was estimated to be grassland. This result stems from the scarcity and patchiness of grassland in the area, meaning that any extrapolation was bound to incorporate more nongrassland than grassland. In addition, bushland and woodland are sometimes actively selected by zebras, e.g., after a predation event [34], meaning that the fitted landcover selection model did not strictly discount these landcover types and that it might be interesting to fit a timevarying landcover selection model in this context. Lastly, the estimated EAKDE bandwidth (0.20) was smaller than that of AKDE (0.31) yielding a smaller estimated home range from EAKDE than SEAKDEc. The decrease in bandwidth suggests that the methodology successfully took up the information that some of the movement variance was caused by resource selection, not stochastic diffusion. These empirical results overall suggest that the new reference function suppressed the reference function approximation bias.
Extrapolation/interpolation, parametric/nonparametric
As outlined by several authors previously, MKDE does not measure the same thing as KDE [33, 36]. MKDE implements an interpolation. The potential for confusion has led some authors to recommend against using the terminology of “utilization distribution” for MKDE and other interpolative methods, and to restrict the use of that phrase to extrapolated distributions [33]. Accordingly, in the zebra case study, the “home range” area was much smaller when estimated with MKDE than with other methods (Table 1, Fig. 1). The isopleth of the MKDE distribution quantifies the process uncertainty around the interpolated path, and the interpolated time budget during the study. By contrast, KDEr and AKDEc are designed to delineate asymptotically consistent, statistically robust, conservative buffers around activity centers.
Like KDEr and AKDEc, SEAKDEc, and EAKDE extrapolate the utilization probability. But, contrary to KDEr and AKDEc, they include a dose of mechanism (Eq. 1) into the base kernels of the extrapolation (Eq. 5). This yields what can be termed a semiparametric extrapolation. Compared to AKDEc, the dose of mechanism modifies what the method considers a plausible realization of the movement process. The objective is to combine the asymptotic consistence of AKDE with the biological realism of fully parametric methodologies. One of the main criticism of fully parametric extrapolations [40] is indeed their sensitivity to the goodness of fit of the underlying mechanistic movement model, which the semiparametric approach partly relaxes. In terms of biological inference, SEAKDEc, and EAKDE extrapolate the potentially accessible resources under a known set of movement rules and under the constraint that the movement paths must all pass through the recorded locations. Comparing extrapolations from different models can help infer new movement mechanisms or assess model parsimony. Finally, compared to the IBS approach, SEAKDEc and EAKDE provide three key advantages: no tuning parameters, full conditioning on the recorded locations, and asymptotic consistency. As highlighted above, a way to better articulate the IBS approach with EAKDE would be to use IBS to generate the pilot estimate upon which to base the reference function.
Conclusion
The key message is that it should soon be possible to make the statistically robust, asymptotically consistent alternatives to KDE less bandwidthlimited than they currently are, and make them yield more realistic, less ovoid home range shapes. I introduced new semiparametric methodologies, based on the stepselection framework [23]. I outlined several avenues for future improvement. I expect EAKDE to function as part of an iterative process by which semimechanistic extrapolations are compared to realized resource use until no significant improvement can be made by adding new movement rules in the extrapolation process. The point would be to give less importance to the time spent at a given location, and more importance to the ratio between availability and use.
Abbreviations
 AKDE:

autocorrelated kernel density estimator
 EAKDE:

autocorrelated kernel density estimator with environmental interactions.
 KDE:

kernel density estimator
 MKDE:

movementbased kernel density estimator
References
 1.
Fieberg J, Börger L. Could you please phrase “home range” as a question? J Mammal. 2012;93:890–902.
 2.
Kie JG, Matthiopoulos J, Fieberg J, Powell RA, Cagnacci F, Mitchell MS, et al. The homerange concept: are traditional estimators still relevant with modern telemetry technology? Philos Trans R Soc B Biol Sci. 2010;365:2221–31.
 3.
Powell RA, Mitchell MS. What is a home range? J Mammal. 2012;93:948–58.
 4.
Silverman BW. Density estimation for statistics and data analysis. London: Chapman and Hall; 1986.
 5.
Worton B. Kernel methods for estimating the utilization distribution in homerange studies. Ecology. 1989;70:164–8.
 6.
Kie JG. A rulebased ad hoc method for selecting a bandwidth in kernel homerange analyses. Anim Biotelemetry. 2013;1:13.
 7.
Börger L, Franconi N, De Michele G, Gantz A, Meschi F, Manica A, et al. Effects of sampling regime on the mean and variance of home range size estimates. J Anim Ecol. 2006;75:1393–405.
 8.
Silverman BW. Spline smoothing: the equivalent variable kernel method. Ann Stat Institute of Mathematical Statistics. 1984;12:898–916.
 9.
Hansteen TL, Andreassen HP, Ims RA. Effects of spatiotemporal scale on autocorrelation and home range estimators. J Wildl Manag. 1997;61:280–90.
 10.
Dray S, RoyerCarenzi M, Calenge C. The exploratory analysis of autocorrelation in animalmovement studies. Ecol Res. 2010;25:673–81.
 11.
Swihart RK, Slade NA. Testing for Independence of observations in animal movements. Ecology. 1985;66:1176–84.
 12.
Fleming CH, Fagan WF, Mueller T, Olson KA, Leimgruber P, Calabrese JM. Rigorous homerange estimation with movement data: a new autocorrelated kerneldensity estimator. Ecology. 2015;96:1182–8.
 13.
Girard I, Ouellet JP, Courtois R, Dussault C, Breton L. Effects of sampling effort based on GPS telemetry on homerange size estimations. J Wildl Manag. 2002;66:1290–300.
 14.
Laver PN, Kelly MJ. A critical review of home range studies. J Wildl Manag. 2008;72:290–8.
 15.
Seaman DE, Millspaugh JJ, Kernohan BJ, Brundige GC, Raedeke KJ, Gitzen RA. Effects of sample size on kernel home range estimates. J Wildl Manag. 1999;63:739.
 16.
De Solla SR, Bonduriansky R, Brooks RJ. Eliminating autocorrelation reduces biological relevance of home range estimates. J Anim Ecol. 1999;68:221–34.
 17.
Fleming CH, Calabrese JM. A new kernel density estimator for accurate homerange and speciesrange area estimation. Methods Ecol Evol. 2017;8:571–9.
 18.
Slaght JC, Horne JS, Surmach SG, Gutiérrez RJ. Home range and resource selection by animals constrained by linear habitat features: an example of Blakiston’s fish owl. J Appl Ecol. 2013;50:1350–7.
 19.
Benhamou S, Cornélis D. Incorporating movement behavior and barriers to improve kernel home range space use estimates. J Wildl Manag. 2010;74:1353–60.
 20.
Jones MC. Simple boundary correction for kernel density estimation. Stat Comput. 1993;3:135–46.
 21.
Tenreiro C. Boundary kernels for distribution function estimation. Revstat Stat J. 2013;11:169–90.
 22.
Forester JD, Im HK, Rathouz PJ. Accounting for animal movement in estimation of resource selection functions: sampling and data analysis. Ecology. 2009;90:3554–65.
 23.
Lele SR, Keim JL. Weighted distributions and estimation of resource selection probability functions. Ecology. 2006;87:3021–8.
 24.
Bischof R, Steyaert SMJG, Kindberg J. Caught in the mesh: roads and their networkscale impediment to animal movement. Ecography. 2017;40:1369–80.
 25.
Beyer HL, Gurarie E, Börger L, Panzacchi M, Basille M, Herfindal I, et al. “You shall not pass!”: quantifying barrier permeability and proximity avoidance by animals. J Anim Ecol. 2016;85:43–53.
 26.
Dunn JE, Gipson PS. Analysis of radio telemetry data in studies of home range. Biometrics. 1977;33:85–101.
 27.
Johnson DS, Thomas DL, Ver Hoef JM, Christ A. A general framework for the analysis of animal resource selection from telemetry data. Biometrics. 2008;64:968–76.
 28.
Fieberg J, Matthiopoulos J, Hebblewhite M, Boyce MS, Frair JL. Correlation and studies of habitat selection: problem, red herring or opportunity? Philos Trans R Soc B Biol Sci. 2010;365:2233–44.
 29.
Johnson DS, Hooten MB, Kuhn CE. Estimating animal resource selection from telemetry data using point process models. J Anim Ecol. 2013;82:1155–64.
 30.
Cressie NAC. Statistics for spatial data. New York: Wiley; 1993.
 31.
Berman M, Turner TR. Approximating point process likelihoods with GLIM. J R Stat Soc C Appl Stat. 1992;41:31–8.
 32.
Robb WL. Thin silicone membranes. Their permeation properties and some applications. Ann N Y Acad Sci. 1968;146:119–37.
 33.
Calabrese JM, Fleming CH, Gurarie E. Ctmm: an r package for analyzing animal relocation data as a continuoustime stochastic process. Methods Ecol Evol. 2016;7:1124–32.
 34.
Courbin N, Loveridge AJ, Macdonald DW, Fritz H, Valeix M, Makuwe ET, et al. Reactive responses of zebras to lion encounters shape their predatorprey space game at large scale. Oikos. 2016;125:829–38.
 35.
Arraut EM, Loveridge AJ, ChamailléJammes S, Fox HV, Macdonald DW. The 20132014 vegetation structure map of Hwange National Park, Zimbabwe, produced using free images and software. KOEDOE  African Prot Area Conserv Sci. 2018;60:a1497.
 36.
Benhamou S. Dynamic approach to space and habitat use based on biased random bridges. PLoS One. 2011;6:e14592.
 37.
Fleming CH, Fagan WF, Mueller T, Olson KA, Leimgruber P, Calabrese JM. Estimating where and how animals travel: an optimal framework for path reconstruction from autocorrelated tracking data. Ecology. 2016;97:576–82.
 38.
Signer J, Fieberg J, Avgar T. Estimating utilization distributions from fitted stepselection functions. Ecosphere. 2017;8:e01771.
 39.
Wang M, Grimm V. Home range dynamics and population regulation: an individualbased model of the common shrew Sorex araneus. Ecol Model. 2007;205:397–409.
 40.
Moorcroft PR, Lewis MA. Mechanistic home range analysis. Princeton: Princeton University Press; 2004.
Acknowledgements
S. ChamailléJammes provided the zebra data and constructive comments on a previous draft. I thank C. Fleming for a detailed and helpful review, J. Calabrese for insightful discussions of previous drafts, and an anonymous reviewer. V. Miele advised on the incorporation of C++ code.
Funding
The zebra project was funded by ANR08BLAN0022 and ANR16CE02–000101.
Availability of data and materials
The data have been uploaded to MoveBank (www.movebank.org) under the study name “Plains zebra ChamailléJammes Hwange NP”.
Author information
Affiliations
Contributions
The author read and approved the final manuscript.
Corresponding author
Correspondence to Guillaume Péron.
Ethics declarations
Ethics approval and consent to participate
The zebra project operates under Zimbabwe Parks and Wildlife Management Authority permit # 23(1)(c)(ii)03/2009.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional files
Additional file 1:
Appendix A: Description of the EAKDE bandwidth optimizer. (PDF 271 kb)
Additional file 2:
Appendix B: Additional method elements: The spacetime point process likelihood and Approximate routine to determine whether a linear feature is crossed. (PDF 241 kb)
Rights and permissions
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.
About this article
Received
Accepted
Published
DOI
Keywords
 Movement ecology
 AKDE
 Temporal autocorrelation
 Step selection function
 Resource selection
 Point process pattern
 Semiparametric