The role of the bandwidth matrix in influencing kernel home range estimates for snakes using VHF telemetry data

Context. Despite the diversity of available home range estimators, no single method performs equally well in all circumstances. It is therefore important to understand how different estimators perform for data collected under diverse conditions. Kernel density estimation is a popular approach for home range estimation. While many studies have evaluated different kernel bandwidth selectors, few studies have compared different formulations of the bandwidth matrix using wildlife telemetry data. Additionally, few studies have compared the performance of kernel bandwidth selectors using VHF radio-telemetry data from small-bodied taxa. Aims. In this study, we used eight different combinations of bandwidth selectors and matrices to evaluate their ability to meet several criteria that could be potentially used to select a home range estimator. Methods. We used handheld VHF telemetry data from two species of snake displaying non-migratory and migratory movement patterns. We used subsampling to estimate each estimator's sensitivity to sampling duration and fix rate and compared home range size, the number of disjunct volume contours and the proportion of telemetry fixes not included in those contours among estimators. Key Results. We found marked differences among bandwidth selectors with regards to our criteria but comparatively little difference among bandwidth matrices for a given bandwidth selector. Least-squares cross-validation bandwidths exhibited near-universal convergence failure whereas likelihood cross-validation bandwidths showed high sensitivity to sampling duration and fix rate. The reference, plug-in and smoothed cross-validation bandwidths were more robust to variation in sampling intensity, with the former consistently producing the largest estimates of home range size. Conclusions. Our study illustrates the performance of multiple kernel bandwidth estimators for estimating home ranges with datasets typical of many small-bodied taxa. The reference and plug-in bandwidths with an unconstrained bandwidth matrix generally had the best performance. However, our study concurs with earlier studies indicating that no single home range estimator performs equally well in all circumstances. Implications. Although we did not find strong differences between bandwidth matrices, we encourage the use of unconstrained matrices because of their greater flexibility in smoothing data not parallel to the coordinate axes. We also encourage researchers to select an estimator suited to their study objectives and the life history of their study organism.