Examining Temporal Sample Scale and Model Choice with Spatial Capture-Recapture Models in the Common Leopard Panthera pardus

Joshua F. Goldberg; Tshering Tempa; Nawang Norbu; Mark Hebblewhite; L. Scott Mills; Tshewang R. Wangchuk; Paul Lukacs

doi:10.1371/journal.pone.0140757

Abstract

Many large carnivores occupy a wide geographic distribution, and face threats from habitat loss and fragmentation, poaching, prey depletion, and human wildlife-conflicts. Conservation requires robust techniques for estimating population densities and trends, but the elusive nature and low densities of many large carnivores make them difficult to detect. Spatial capture-recapture (SCR) models provide a means for handling imperfect detectability, while linking population estimates to individual movement patterns to provide more accurate estimates than standard approaches. Within this framework, we investigate the effect of different sample interval lengths on density estimates, using simulations and a common leopard (Panthera pardus) model system. We apply Bayesian SCR methods to 89 simulated datasets and camera-trapping data from 22 leopards captured 82 times during winter 2010–2011 in Royal Manas National Park, Bhutan. We show that sample interval length from daily, weekly, monthly or quarterly periods did not appreciably affect median abundance or density, but did influence precision. We observed the largest gains in precision when moving from quarterly to shorter intervals. We therefore recommend daily sampling intervals for monitoring rare or elusive species where practicable, but note that monthly or quarterly sample periods can have similar informative value. We further develop a novel application of Bayes factors to select models where multiple ecological factors are integrated into density estimation. Our simulations demonstrate that these methods can help identify the “true” explanatory mechanisms underlying the data. Using this method, we found strong evidence for sex-specific movement distributions in leopards, suggesting that sexual patterns of space-use influence density. This model estimated a density of 10.0 leopards/100 km² (95% credibility interval: 6.25–15.93), comparable to contemporary estimates in Asia. These SCR methods provide a guide to monitor and observe the effect of management interventions on leopards and other species of conservation interest.

Citation: Goldberg JF, Tempa T, Norbu N, Hebblewhite M, Mills LS, Wangchuk TR, et al. (2015) Examining Temporal Sample Scale and Model Choice with Spatial Capture-Recapture Models in the Common Leopard Panthera pardus. PLoS ONE 10(11): e0140757. https://doi.org/10.1371/journal.pone.0140757

Editor: Christof Markus Aegerter, University of Zurich, SWITZERLAND

Received: April 6, 2015; Accepted: September 30, 2015; Published: November 4, 2015

Copyright: © 2015 Goldberg et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All data are available from the Dryad Database (accession number 10.5061/dryad.mr1pt)

Funding: National Aeronautics and Space Administration (grant # NNX11OA47G) to MH; Ugyen Wangchuck Institute for Conservation and Environment Jigme Singye Wangchuck Research and Training Fund to TT; Bhutan Foundation research grant to TT. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Across the globe, large carnivores face threats from habitat loss and fragmentation, prey depletion, and poaching [1]. These risk factors present an important conservation challenge because large carnivores can play critical ecological roles as keystone and umbrella conservation species [1–5]. The combination of the ecological impact of carnivores and risks to their persistence merit increased vigilance of their population status and trends to support effective management [6, 7]. However, many of these species challenge conventional monitoring methods, as large carnivores often occur at low densities in dense cover habitats, and have a shy, solitary, elusive, cryptic and/or nocturnal nature [8–12]. In short, large carnivores are frequently difficult to detect.

To address the imperfect detectability of large carnivores, researchers have adopted non-invasive, remote sampling techniques, such as automatically triggered camera traps, hair snares or scat surveys that likely reduce animal avoidance of target sites. For many species, these data can identify unique individuals based on genetics or unique markings, facilitating a capture-mark-recapture (CMR) analysis that explicitly incorporates detection probability into the estimate of abundance and/or density. Spatial capture-recapture (SCR) models further refine these techniques, by incorporating the spatial distribution of an individual’s movement inferred from observed captures [13, 14], which may be particularly appropriate given the territorial habits of many large carnivores. Furthermore, SCR models address a critical weakness within conventional CMR methods by linking estimates of abundance to a well-defined area, eliminating the need to estimate density with an ad-hoc estimate of area [14–17]. Simulations confirm that SCR analysis reduces bias in estimates of abundance or density under a range of conditions [12, 14, 18–20]. Additionally, SCR models can integrate ecological and behavioral factors into density estimates using individual, trap or spatial covariates, which can further improve the accuracy of estimates [12, 20–22].

Despite these advances, SCR methods raise new questions about sampling design and its relationship to population estimates and precision, particularly as these techniques are applied to new species and systems, including large carnivores. Russell et al. [23], Sollman et al. [18] and Sun et al. [19] examine model sensitivity to the spatial configuration of SCR study design. Wilton et al. [20] consider the arrangement of traps with reference to black bears, while Tobler and Powell [12] address issues of camera spacing and sampling area using parameter values specific to jaguars. Potential differences related to temporal variation in detections within SCR models remain untested, and guidelines developed for conventional CMR analysis may create confusion for practitioners of SCR methods. Conventional guidelines advise collapsing data temporally to improve detection probability [24, 25], despite the risk of lost detections and bias, and the potential for continuous-time analyses in an SCR framework [12, 26]. While these continuous-time models may provide the most technically sound approach, their recent development and complexity may limit their application in many conservation settings. Thus, understanding the relationship between sampling interval duration and the precision of density estimates in discrete-time models still hold practical value in maximizing the effectiveness and applicability of these approaches in a management setting.

Beyond sampling structure, the ability to incorporate covariates into density estimates raises important questions about model selection, particularly when estimates disagree. While maximum likelihood-based SCR estimation techniques have adopted an Akaike information criterion (AIC) model selection framework, no robust methods have been widely applied in Bayesian approaches to SCR analysis. These limitations in assessing the relative merits of different models restrict the usefulness of these methods, especially for data-poor species like large carnivores.

Here, we develop two methodological advances that may improve sampling design and analysis of capture-recapture studies for a variety of taxa: (1) we examine the impact of the temporal resolution of the data on SCR model estimates and precision and (2) we implement a Bayes factor approach to model selection with multiple covariates. We address these issues through a simulation study and using a remote camera-trapping data set for the common leopard (Panthera pardus) from the lower, subtropical foothills of Royal Manas National Park (RMNP), Bhutan. The common leopard occurs from sub-Saharan Africa to the Russian Far East, as well as on the islands of Sri Lanka and Java [27, 28] and faces threats similar to many other large carnivores. In Bhutan, human-wildlife conflict over livestock depredation has resulted in large numbers of leopards being killed [29–31]. Despite this threat to the common leopard population, knowledge of leopard populations is limited in spatial extent and relies on conventional CMR techniques [32]. Our application within this system and our simulations show that different sample interval lengths have a substantial impact on the precision of SCR density estimates and that Bayes factors can help discriminate between models including sex-specific capture probabilities and sex-specific distributions of movement. Given these results, we make recommendations for other capture-recapture studies, particularly those of large carnivores.

Materials and Methods

Study Area

RMNP is located in the southern foothills of Bhutan (90^° 57’37.61” E, 26^° 47’31.27” N) and borders India’s Manas Tiger Reserve, forming a trans-boundary conservation landscape. RMNP covers 1057 km², with elevations from 90 m in the southern foothills to 2900 m in the north (Fig 1). RMNP experiences hot, humid summers followed by cool, dry winters with annual maximum temperatures ranging from 20°C to 34°C. RMNP has diverse vegetation communities including; cool, moist and warm broadleaf forests, wetlands, subtropical dry forests, and subtropical scrub and grasslands, as well as agricultural fringe areas. Common prey species for wild felids include: sambar (Cervus unicolor), barking deer (Muntiacus muntjak), wild pigs (Sus scrofa), serow (Capricornis thar), goral (Naemorhedous goral), Himalayan crestless porcupine (Hystrix brachyuran), langurs (Trachypithecus spp), macaques (Macaca spp.) and several bird species.

Download:

Fig 1. Spatial Density Estimate of Common Leopards (Panthera pardus) from the Best-Supported Model with Inset Study Area Map.

The posterior spatial density estimate of common leopards/100 km² from the best-supported spatial capture-recapture model, σ_sex, in the lower foothills of Royal Manas National Park (RMNP), Bhutan for sampling carried out during 2010–2011. The 162 km² sampling area is displayed with the solid black line, the RMNP boundary with broken black line, camera-trapping stations with leopard detections with crosses and camera-trapping stations without leopard detections black circles. Each station represents a pair of cameras. Inset: RMNP (light gray) in Bhutan with the location of the 162 km² gridded study area (black) for common leopards in 2010–2011.

https://doi.org/10.1371/journal.pone.0140757.g001

Ethics Statement

The Ministry of Agriculture and Forests, Royal Government of Bhutan, approved this study as a part of long-term predator monitoring in Royal Manas National Park. The Ugyen Wangchuck Institute for Conservation and Environment (UWICE) conducted this study in collaboration with the University of Montana. UWICE is a governmental research institute with the mandate to conduct studies in the protected areas of Bhutan. Since we did not handle any animals and used non-invasive, remote camera traps, no animal handling, care nor ethics permits were necessary.

Field Methods

We conducted a camera trapping study in the lower foothills of the Manas administrative range of RMNP. This area was chosen due to high expected wild felid diversity and to minimize logistical and security constraints [33]. We placed pairs of camera traps (Reconyx Inc., Holmen, WI, U.S.A.) in 29 locations within this area from mid-November to mid-February for a total of 91 days. To distribute camera traps, we first created a grid of 2.5 km x 2.5 km cells that covered 162 km², although international and topographic boundaries made some irregular (Fig 1). We then distributed available camera pairs among the 56 total cells, but logistical constraints sometimes prevented camera placement in accordance with the grid, such that three grid cells had two pairs of cameras and one camera pair was placed outside of the grid.

We placed cameras along walking paths and animal trails that showed signs of wildlife use to maximize the probability of photographing target felid species. Cameras were mounted on trees or cut poles and placed two to four meters from the focal movement pathway. We cleared vegetation between the camera and the trail. To minimize elephant damage to cameras, we placed fresh elephant dung on and around cameras. We did not move camera traps during the 91-day sampling period and visited cameras as needed (every two to four weeks) to change batteries and memory cards. We identified individuals based on unique pelage patterns and determined the sex of these individuals as possible from sex-specific cues, such as visible genitalia or the presence of young. We performed this identification visually without the aid of any pattern-recognition software. A single author (TT) verified all leopard identifications to minimize any observer effects.

Analytical Methods

We estimated the abundance and density of leopards within our study area using a spatial capture-recapture model [14, 21, 23, 34]. We followed the hierarchical model formulation described by Royle et al. [35] for a camera-trapping grid. The model relates the observations, y_ijk, of individual i in trap j during sampling interval k to the latent distribution of activity centers. We let the observation, y_ijk, take the value of one for a capture and zero if not captured to produce a capture history for all individuals in all traps over all sampling intervals. We treated multiple detections of an individual in a particular trap during the same sampling interval as a single capture. Individuals could be captured in multiple traps during a sampling interval. We followed the formulation of the observation process used by Gardner et al. [36] and Russell et al. [23]. Within this model structure, we allowed the shape of the movement distribution to vary between an exponential and half-normal model (i.e., 1 ≤ θ ≤ 2, where θ represents the exponent of the kernel from the exponential family of distributions).

We fit a suite of 16 models to the data to assess the impact of varying the duration of sampling intervals and different covariates. We used 4 different sampling interval durations, corresponding to daily (91 sampling occasions), weekly (13 sampling occasions), monthly (three sampling occasions) and quarterly (one sampling occasion) sample periods. For each of these durations, we fit four different covariate combinations: (1) a model where detectability depended only on distance between a trap and the individual’s activity center (Distance), (2) a model with sex as a covariate on baseline detectability (Sex), (3) a model with sex as a covariate on the scale of the movement distribution (σ_sex), and (4) a model with sex as a covariate on baseline detectability and the scale of the movement distribution (Sex + σ_sex). We fit these models using an area that buffered the trapping grid by 10 km, yielding a state space of 1,551 km². This buffer incorporates individuals with activity centers outside of the trapping grid, but whose movement range extends into the trapping grid, allowing for a robust estimate of density.

Bayesian Model Analysis

We used a Bayesian approach to model analysis using data augmentation [14, 37, 38], which has been successfully implemented in a number of recent SCR models [21, 23, 36]. This technique adds a sufficiently large number of all-zero (un-encountered) capture histories to create a dataset of size M individuals. We determined the augmentation to be large enough, when the number of augmented individuals did not limit posterior estimates of population size. We chose a uniform prior distribution from [0, M] on population size.

We fit our models using Markov chain Monte Carlo (MCMC) methods in R [39], using the SCRbayes package (available at: https://sites.google.com/site/spatialcapturerecapture/scrbayes-r-package; S1 File). We ran models for 20,000 iterations, discarded the first 5,000 iterations as burn-in and further thinned the chain by skipping every other iteration to reduce autocorrelation, leaving 7,500 iterations in our posterior sample. We assessed the convergence of the MCMC samples by examining trace plots and histograms for each parameter. From these converged samples, we computed the mean, median and 95% credibility intervals for the model parameters.

Model Selection

The flexibility of SCR models to incorporate covariates at different levels allows estimates to include added complexity, but also presents ecologists with potential model selection problems when using Bayesian methods. Estimates and precision may vary among different models creating uncertainty regarding the population status and preventing insight into the mechanisms generating patterns of abundance [12, 23, 40]. These model selection problems can be especially troublesome in an applied context, where managers must make decisions regarding where and how to distribute limited resources.

To assess the relative merit of the different model parameterizations and bridge this methodological gap, we computed the posterior model probabilities for models with daily sampling intervals. Let each model parameterization be represented by M_l for l = 0,…,3, for the data, y. From Bayes’s theorem, the relative posterior probabilities among the models in the candidate set are given by Where Pr(y|M_l) is the likelihood of the data given the model and Pr(M_l) is the prior probability of a model being correct. For our evaluation, we chose an uninformative prior and set Pr(M_l) = 0.25 for all models. Under this assumption, Bayes’s theorem then reduces to the proportion of the total likelihood attributed to a particular model. Gelfand and Dey [41] proposed an approximation to this quantity using the posterior distribution of the model parameters, such that Where f(·) is a probability density function of the same dimension as ω_l, the parameters in the model, and ω_l⁽ⁱ⁾ is a sample from the posterior distribution for a particular iteration of the MCMC, i = 1,…, c. For our models, c took the value of 7,500, the number of posterior samples. The denominator is the likelihood multiplied by the density at the value of the parameters for any given iteration of the chain. We chose f(·)~t₇₅₀₀(, ∑), a multivariate t-distribution with 7,500 degrees of freedom, centered at the mean of the posterior samples and with a scale matrix, ∑, the correlation matrix of the samples. The Gelfand and Dey [41] approximation compares favorably to alternative approximations of the marginal likelihood [42, 43]. Substituting these marginal likelihoods into Bayes’s theorem above gives the relative probability of each model given the data.

We assessed the support of the additional sex-specific parameters for detectability and scale of the movement distribution using Bayes factors [44, 45]. To compute the Bayes factor between the distance-only model, M₀, and each of the sex-specific models, M_l (l = 1, 2, 3), we used which describes the ratio of the posterior odds to prior odds for M_l and M₀ [42]. The Bayes factor provides a Bayesian analog to the likelihood ratio test, although Bayes factors differ in two important respects. Bayes factors use the integrated marginal likelihood instead of the maximized likelihood, such that model uncertainty is incorporated into the test statistic, and Bayes factors do not require that hypotheses have nested forms.

Simulations

We conducted a simulation study to corroborate our findings with respect to the effects of sampling interval on estimated density and precision, and test our model selection techniques. We simulated 89 SCR datasets for a population of size N = 250 (150 female, 100 male) distributed over an 11 x 11 unit continuous statespace for a density of 1.98 individuals/unit². Our simulations exposed these individuals to an 8 x 8 trap grid centered within the statespace for 90 (daily) sampling occasions. We generated capture histories using a Binomial observation model with a half-normal hazard rate detection function and sex-specific baseline detection probabilities, i.e. the Sex model. To parameterize the data generation, we used estimates from models fit to our field data. We set σ = 1 unit, λ_{0, female} = 0.05 and β_sex = -1.61. For each simulation, we then subsampled the data to have 1 (quarterly), 3 (monthly) and 13 (weekly) sampling occasions. Finally, we fit the Distance, Sex, σ_sex, and Sex + σ_sex models to all temporal resolutions of each simulated dataset. To fit these models, we discretized the statespace into an 11.5 x 11.5 grid of points and fixed θ = 1 (a half-normal hazard rate detection function). We ran all models with 5,000 MCMC iterations, discarding the first 1,000 iterations as burn-in. We use these results to assess the accuracy and precision of model estimates and the fidelity of our Bayesian model selection methods.

Results

Camera Trapping

We captured 22 individual common leopards a total of 82 times between 17 November 2010 and 15 February 2011 in 2,639 camera-trapping nights. The median number of captures per individual was 2 (range: 1–16 captures) and individuals were captured in between one and eight different traps in daily sampling opportunities. We identified six of these individuals as males and the remainder (16 leopards) as female based on photographic evidence of genetalia and/or offspring. We considered camera stations to have operated continuously for the duration of the study, since both cameras at a trap location never malfunctioned simultaneously.

Model Results

We fit 16 models to the camera captures of common leopard we obtained from sampling. The density estimates varied from a minimum of 4.1 leopards/100 km² in the Distance-only model with monthly sampling periods to 15.8 leopards/100 km² in the sex + σ_sex model with a single, quarterly sampling interval (Table 1). For all sampling interval lengths, the Distance model estimated a lower density of animals than any of the models with sex-specific parameters, such that the Distance model median did not fall within the 95% credibility intervals (CI) of the estimates from the other models. The models with sex-specific parameters provided comparable density estimates and none differed at the 95% credibility level (Table 1; Fig 2). Sampling interval length affected the density estimates of models including a sex-specific baseline detection probability, as the density estimates of the Sex and Sex + σ_sex model declined with shorter sampling intervals (Table 1; Fig 2). These trends, however, were not significant at the 95% credibility level. The precision of the density estimates also varied with sampling interval length. Shorter sampling intervals typically yielded narrower 95% CIs for all models, especially the Sex and Sex + σ_sex models, although we note that the Distance model with a weekly sampling interval did not follow this pattern (Table 1; Fig 2). Sampling interval length also produced a distinct pattern in the baseline detection probability (λ₀) and sex-specific covariate on this parameter (β_sex). Shorter sampling intervals had lower baseline detection probabilities (for both sexes, in sex-specific models) (S2 Table). This result matches the intuitive expectation that the detection probability declines as the interval length becomes shorter.

Download:

Table 1. Common Leopard Population Density Estimates, 95% Credibility Intervals and Standardized 95% Credibility Interval Width for All Sampling Intervals and Model Parameterizations.

https://doi.org/10.1371/journal.pone.0140757.t001

Download:

Fig 2. Variation in Median Leopard Density Estimates and 95% Credibility Interval Width for Different Sampling Intervals.

Effects of changing sample period duration on median estimated densities (in leopards/100 km²; A) and 95% credibility interval width (in leopards/100 km²; B) of common leopards in Royal Manas National Park, Bhutan during winter 2010–2011 for 16 spatial capture-recapture models. The 16 models include all possible combinations of the four sampling periods—Quarterly, Monthly, Weekly and Daily—and four model specifications—Distance (dark gray), Sex (medium gray), σ_sex (light gray) and Sex + σ_sex (white). Error bars represent 95% credibility intervals.

https://doi.org/10.1371/journal.pone.0140757.g002

In models with a daily sampling interval, the σ_sex model received 0.971 of the posterior model probability among the four candidate model parameterizations. The Bayes factor comparison to the distance model showed strong support for this model, as it was over 33 times as likely as the distance model and had positive values of log₁₀ B_l0 and 2 ln B_l0 (Table 2). This model estimated a median density of 10.0 leopards/100 km² (95% CI of 6.25–15.93 leopards/100 km²) (Fig 1), which corresponds to an abundance of 155 leopards (95% CI: 97–247 leopards) within the statespace.

Download:

Table 2. Bayesian Model Selection Results for Daily Sampling Interval Models for Common Leopards.

https://doi.org/10.1371/journal.pone.0140757.t002

Simulations

The simulated datasets represented capture histories for a median of 223 unique individuals (95% CI: 212.0–232.6) recorded a median of 2479 times (95% CI: 2266.0–2752.2) over 90 sampling occasions, 2325 times (95% CI: 2111.8–2569.0) over 13 sampling occasions, 1864 times (95% CI: 1703.2–2053.0) over 3 sampling occasions and 1255 times (95% CI: 1142.8–1360.8) over a single sampling interval. On average, capture events recorded 144 female (95% CI: 138.2–148.0) and 79 male individuals (95% CI: 70.4–86.0).

The median of median density estimates across simulated datasets for each model and sample interval length varied from a minimum of 1.84 individuals/unit² with the Distance model, using either daily or weekly sampling periods, to a maximum of 2.35 individuals/unit² with the Sex + σ_sex model with a quarterly sampling period (Table 3). The Distance model produced significantly lower density estimates than models that incorporated sex-specific parameters, but density estimates from models that included sex-specific parameters did not differ at the 95% credibility level (Table 3). This pattern reflects a negative bias of estimates from the Distance model and a positive bias of estimates from models with sex-specific detection parameters.

Download:

Table 3. Simulated Population Density Estimates, 95% Credibility Intervals and Standardized 95% Credibility Interval Width for All Sampling Intervals and Model Parameterizations.

https://doi.org/10.1371/journal.pone.0140757.t003

Precision of the density estimates depended upon sampling interval length. The Distance, Sex and Sex + σ_sex models showed improved precision when moving from quarterly sampling interval to finer temporal resolutions (monthly, weekly or daily sampling intervals) (Table 3). Precision did not differ substantially or predictably among monthly, weekly or daily sampling intervals (Table 3).

The Sex model received the greatest support across all simulations for all sampling intervals. The Bayes factors model selection process identified the Sex model as having the highest posterior probability at a minimum of 60.7% of simulations with a weekly sampling interval and as much as 71.9% of simulations with quarterly sampling intervals (Table 4). As the most frequent top model choice, the Sex model had the highest median posterior probability estimate in all sampling intervals, receiving a minimum of 94% median posterior model probability with weekly sampling intervals (Table 4). These posterior probabilities varied significantly across simulations for models with sex-specific baseline detection probabilities, as the 95% credibility intervals of the Sex and Sex + σ_sex models ranged from 0 to 1. However, the Distance and σ_sex model received relatively little support across all simulations and sampling intervals with the σ_sex model identified as the top model in 2.2% of simulation replicates with quarterly sampling intervals.

Download:

Table 4. Bayesian Model Selection Results for Daily Sampling Interval Models Fit to Simulated Data.

https://doi.org/10.1371/journal.pone.0140757.t004

Discussion

We provide two advances in the implementation and use of SCR models through a simulation study and our analysis of the common leopard population in the lower foothills RMNP, Bhutan. First, our investigation of different sample interval lengths provides some guidance to reduce the variability of SCR population estimates, particularly for difficult to detect species, like large carnivores. Second, we apply Bayes factors in a novel context to perform SCR model selection. These results suggest new insights and avenues for future development of these methods within population ecology.

Sample Interval Length and SCR Estimates

Sample interval length affected the precision of density estimates in our SCR analysis of the leopard data and our simulation study. For the common leopard, we found that analyzing our data with daily sample intervals provided the best precision for all models, especially those with covariates (Table 1; Fig 2). The advantages of daily sampling intervals for model precision did not hold in our simulations, where the principle improvement in precision occurred when moving from a quarterly sampling interval to any finer temporal resolution (Table 3). We believe these divergent results largely depend upon the size of the leopard versus the simulated datasets. Even at our finest temporal resolution (daily sampling intervals), we recorded only 82 leopard captures, a relatively sparse dataset by any measure. In contrast, our simulated data included 1255 capture events (95% CI: 1142.8–1360.8) at the coarsest, quarterly scale. Despite these relatively rich datasets, we still observed improved precision with shorter sampling intervals that increased the number of unique capture events in our simulation. However, we note that these gains in precision for our simulations were of smaller magnitude than those observed in the leopard models. The differences in the size of the leopard and simulated datasets likely occurred due to a higher density of individuals and the larger, more regular trapping grid of our simulations [18–20], even though we used estimates from our leopard analysis to parameterize the simulations. Thus, the impact of sampling interval length may have an asymptotic relationship with model precision that further depends upon the size of the gathered dataset and by proxy the study design.

These results suggest some guidelines for the implementation of SCR models. For camera trapping surveys of elusive species, where samples have accurate time signatures and limited captures take place, we recommend using daily (or finer, if recapture rates are high) subdivisions of the data [26, 46]. While these short sample intervals do come with added computation time in both likelihood-based and Bayesian analysis, this specificity maximizes the potential number of recaptures of individuals and allows for better estimation of the movement distributions central to the SCR framework. Of course, some sampling techniques, such as hair snares or scat surveys, may not allow such temporal specificity of samples [23, 36]. For these alternative methods, our results demonstrate substantive gains in precision when moving from quarterly sample intervals to monthly or weekly intervals, especially with high capture rates as in our simulations (Tables 1 and 3; Fig 2). Taken together, our study clearly cautions against highly collapsed sampling intervals in SCR studies.

Bayesian Model Selection

Our simulation study largely validated our Bayes factor model selection methods. The true data generating process in our simulations incorporated a sex-specific baseline detection probability and models that incorporated this sex-specific baseline detection probability received decisive support in almost all simulations (Table 4). The Sex model (with only sex-specific baseline detection probabilities) received the strongest and most frequent support, while the Sex + σ_sex model (with an additional parameter for sex-specific movement distributions) represented a clear runner-up across the four different sampling intervals. Moreover, our selection methods more frequently identified the Sex + σ_sex model as most probable with weekly, and to a lesser extent monthly and daily sampling intervals (Table 4). While these sampling intervals produced datasets rich enough to demonstrate improved precision (Table 3), they may also have contained enough information to support additional, extraneous model complexity and over-fit the available data. We find support for this view in the overlapping confidence intervals for the σ_male and σ_female parameters in this model (S2 Table). The σ_sex component of the Sex + σ_sex model did not capture true sex-specific differences, but rather accounted for small, happenstance variation in some datasets [47]. Thus, our Bayes factors model selection may show some bias towards more complex models with increasing sample size, but these selection errors appear identifiable from examining model estimates.

We may have further improved model precision and reduced the observed bias of the density estimates by refining our simulation approach. Across all sampling interval lengths in our simulations, we observed a negative bias in the median densities estimated by the Distance model and a positive bias in the sex-specific covariate models. We feel that these inaccuracies in the model estimates, particularly the Sex model, which generated the data, likely arise from using a statespace that did not adequately capture the area that contains the activity centers of all detectable individuals [35, 48]. Thus, the model does not account for individuals at the edge of effective sampling area of the trapping array [48]. We find support for this view in the generally negative bias of the σ parameter estimates, especially for more abundant females in the σ_sex model (S2 Table). This spurious model parameter likely reduced the bias of the density estimate by accounting for the diminished range for individual movements at the edge of the statespace. Furthermore, this pattern may have contributed to the observed support for the Sex + σ_sex model. More generally, this observation raises larger questions about the sensitivity of SCR model results to the definition of the statespace.

Having verified our implementation of the Bayes factors approximation in our simulations, we interpret the model selection results from models of common leopards in RMNP. Among models fit with data from the most precise (daily) sampling interval, median density estimates varied from 4.5 in the Distance model to 11.0 in the Sex + σ_sex model (Table 1). Thus, model selection represents a non-trivial matter. On the basis of our Bayes factors approach, we selected the σ_sex model, which allows the size of male and female movement ranges to vary. This model received an overwhelming proportion of the posterior model probability (0.971) and was over 33 times more likely than the null, Distance-only model (Table 2). We estimated male movement ranges to be nearly twice the size of female movement ranges (Table 5), which is supported by our recapture of males over a wider range than females. The relatively small movement ranges of females led the model to estimate a relatively large proportion of undetected females, as shown by the strong female bias in sex ratio, ψ_sex (Table 5).

Download:

Table 5. Parameter Estimates and 95% Credibility Intervals from Spatial Capture-Recapture Models with Daily Sampling Intervals for Common Leopards.

https://doi.org/10.1371/journal.pone.0140757.t005

These ecological covariates may entail some cost to the precision of the density estimate. In our simulation and application to the leopard system, sex-specific models, including our top models, had wider 95% CIs than those of the Distance model. For the leopard analysis, we may have improved the precision of the model by arranging our traps in clusters to provide more spatial recaptures of females, as has been shown in simulations and other systems [18–20]. To the contrary, the relative precision of male and female movement ranges indicates that our sampling may have provided a sufficient number of spatial recaptures to adequately estimate this parameter for females. The size of our trapping grid was approximately 2×σ_female, which falls within recommendations for trap spacing relative to movement distributions from simulations [18, 19]. Overall, this uncertainty underscores the need to maximize both the number of individuals captured and the locations at which they are recaptured in spatial capture-recapture studies [40, 49].

More broadly, the approximation to Bayes factors that we adopt represent only one approach to Bayesian model selection. Hooten and Hobbs [50] provide a timely review of this emerging topic in ecological research, although the appropriateness of the available techniques for SCR models varies. They outline two principle approaches to the problem of Bayesian model selection in the literature. The first class of methods attempts to minimize out-of-sample predictive error by applying some penalty to a measure of model goodness-of-fit. AIC represents the most frequently used method of this type in maximum likelihood analyses. Bayesian analogs to this approach include the deviance information criterion (DIC) [51] and Watanabe-Akaike information criterion (WAIC) [52], however the assumptions of these techniques seem a poor fit for the nature of SCR analyses. Related methods, such as extensions of the WAIC measure [53, 54] or posterior predictive loss [55, 56], may be better suited to SCR models, but come with their own computational challenges [50, 53]. In contrast to these predictive measures of model performance, Bayesian model weighting methods, such as Bayes factors, quantify uncertainty within a set of models given the available data. Promising alternatives to Bayes factors in this wider class of selection techniques include Bernoulli indicator variables (e.g., [57, 58]) or reversible-jump MCMC (RJMCMC) [59], although the former method shows a strong sensitivity to prior distributions that require proper tuning [38, 58] and the latter approach can entail significant technical challenges ([50], but see [60]). Regardless of the chosen model selection tool, researchers should attempt to explore how and under what circumstances each candidate model fails to better understand the underlying ecological process [47, 61].

Ecological Implications

The support for sex-specific parameters in our leopard example highlights the need to consider relevant behavioral variation in estimating abundance and/or density for large carnivores. Previous leopard radio tracking [62, 63] and camera trapping [40] studies have found variation in home or movement range sizes between males and females. The significant difference between our top model and the distance model shows that failing to consider this factor would have led to an underestimate of leopard density in our study area. Similarly, Gray and Prum [40] found support for sex-specific differences in baseline detection probability, suggesting these differences may occur because males use trails more than females. We found some evidence for sex-specific differences in baseline detectability, since β_sex had 95% CIs that did not overlap zero (Table 2). Although alternative approaches to model selection would have included this term in the top model [23], these models had extremely low posterior probabilities in our Bayes factor analysis (<0.001). On the basis of these low posterior probabilities, we rejected these models in favor of the σ_sex model, although we note that the point estimates between these models did not differ significantly. In our study, differences in detectability appear to have occurred because males had larger movement ranges than females, not because of inherent differences in baseline detection probability (e.g. from trail use patterns). Furthermore, our model selection techniques gave additional insight into leopard behavior in our system and improved our density estimate.

A comparison of our estimate to other study areas from the Indian sub-continent demonstrates the importance of considering animal ecology in SCR models [49, 64, 65]. Our estimate of 10.0 leopards/100 km² (95% CI: 6.25–15.93 leopards/100 km²) was within the reported range of leopard densities, from 1.0–28.9 leopards/100 km² [65], across their geographic distribution. In the adjoining Indian Manas National Park, Borah et al. [49] reported a SCR-based estimate of 3.4 ± 0.82 SE leopards/100 km², similar to our estimate for the Distance model of 4.5 leopards/100 km² (95% CI: 2.90–6.64). However, they did not consider sex-specific covariates, which we found to increase density in our study. The Bayes factor model selection approach demonstrated strong support for sex-specific movement distributions suggesting the SCR-based estimate of Borah et al. [49] may be biased low. More broadly, SCR methods in combination with Bayes factors may help resolve the discrepancies that commonly occur in estimates of large carnivore density [12, 23, 65].

Looking towards the future, we see a great deal of promise in SCR models incorporating auxiliary ecological information. Our estimate may have been further improved by incorporating a metric of habitat quality [66], although we did not have any such indices available for our study area. In addition to habitat, leopard density may further depend on the community of sympatric large carnivores. Harihar et al. [64] found that leopard density declined from a median estimate of 9.09 leopards/100 km² (95% CI: 4.54–17.98) to 1.44 leopards/100 km² (95% CI: 0.41–6.61), corresponding to an increase in tiger density across years. These contrasting changes in density across an entire study area raise the possibility that conspecific interference or avoidance may affect patterns of abundance or movement within a study at finer spatial and/or temporal scales. These interactions could be integrated into SCR analyses by introducing the presence or abundance of other species as spatial covariates or through multi-species models. More broadly, SCR methods have been integrated with concepts of resource selection and ecological distance to provide inference on factors affecting habitat use and animal movements [22, 67]. These sorts of analyses expand the scope of our understanding, beyond simply density or abundance, of large carnivores, such as leopards.

Conclusions

Many large carnivores, including leopards, have a broad geographic distribution, but increasingly face threats from habitat disturbance, poaching, inadequate prey resources and human-wildlife conflicts. Despite the many risks facing large carnivores, the status of populations across species’ ranges remains poorly understood. SCR methods hold a great deal of promise to help fill this critical gap of knowledge for large carnivores by explicitly incorporating spatial heterogeneity in detections into our estimates of abundance. These methods yield the greatest precision when combined with a high temporal resolution of sampling (daily intervals), although more modest weekly or monthly sampling intervals can be informative, particularly with high capture rates. Moreover, we implement and validate a Bayesian model selection method for use with these models. These methodological advancements may contribute to population assessments, help prioritize conservation efforts and improve understanding of large carnivore ecology as it relates to population biology. These steps towards a more holistic view of conservation may pay dividends for large carnivores and threatened species globally.

Supporting Information

S1 File. R Code to Implement Spatial Capture Recapture Analysis with Bayes Factors Model Selection.

R code to implement spatial capture recapture models with a novel approach to model selection using the SCRbayes package. This script reads the data for the camera trap locations and leopard capture histories available through the Dryad Database (http://dx.doi.org/10.5061/dryad.mr1pt). This data is used to fit four spatial capture-recapture models with the SCRbayes package for R (package available at: https://sites.google.com/site/spatialcapturerecapture/scrbayes-r-package). Finally, we conduct model selection with Bayes factors. Users should consult the package help files and additional resources for further information.

https://doi.org/10.1371/journal.pone.0140757.s001

(DOCX)

S2 File. Internal R Code to Conduct Bayes Factors Model Selection.

This function within the SCRbayes package will operate on standard model output from SCRj.fn to calculate posterior marginal likelihoods and Bayes factors from spatial capture-recapture models implemented in the SCRbayes package. The SCRbayes package can be accessed at: https://sites.google.com/site/spatialcapturerecapture/scrbayes-r-package. This code comes as a function (SCR.bf) within the package. Users should consult the package help files and additional resources for further information.

https://doi.org/10.1371/journal.pone.0140757.s002

(DOCX)

S1 Table. Parameter Estimates and 95% Credibility Intervals from All Spatial Capture-Recapture Models for Common Leopards.

Median parameter estimates with 95% credibility intervals in parentheses from spatially explicit capture-recapture models of common leopards in Royal Manas National Park during 2010–2011 for all 16 combinations of covariates and data. λ₀ gives the baseline capture probability at an individual’s activity center per sample interval per camera station. β_sex denotes the effect of sex on detection probability on the log scale. The σ parameters describe the scale of an individual’s movement distribution in km, which varies by sex in some models. ψ_sex estimates the proportion of the population that is male. θ represents the shape parameter of the individual’s movement distribution, where 0.5 is exponential and 1.0 is Gaussian.

https://doi.org/10.1371/journal.pone.0140757.s003

(DOCX)

S2 Table. Median Parameter Estimates and 95% Credibility Intervals from All Spatial Capture-Recapture Models from Simulations.

Median of median model parameter estimates with median 95% credibility intervals in parentheses from spatially explicit capture-recapture models fit to all simulated datasets for all 16 combinations of covariates and data. λ₀ gives the baseline capture probability at an individual’s activity center per sample interval per camera station. β_sex denotes the effect of sex on detection probability on the log scale. The σ parameters describe the scale of an individual’s movement distribution in units, which varies by sex in some models. ψ_sex estimates the proportion of the population that is male. True values used to generate the simulated were: λ₀ = 0.05 (with daily sampling intervals), β_sex = -1.61, σ_male = σ_female = 1.0, ψ_sex = 0.4.

https://doi.org/10.1371/journal.pone.0140757.s004

(DOCX)

Acknowledgments

We are grateful to Royal Manas National Park staff for their armed escort and logistical support and Pema Dendup of UWICE. We also thank J. Andrew Royle and Robin E. Russell for analytical support. We thank His Excellency, Honorable Minister, Honorable Secretary, and Honorable Director General of Forests for their continued support and commitment to conservation in Bhutan. We appreciate the constructive criticism of four reviewers.

Author Contributions

Conceived and designed the experiments: JFG; TT; MH; LSM. Performed the experiments: TT. Analyzed the data: JFG; TT. Contributed reagents/materials/analysis tools: JFG; TT; NN; MH; LSM; TRW; PL. Wrote the paper: JFG; TT; MH; LSM; PL.

References

1. Ripple WJ, Estes JA, Beschta RL, Wilmers CC, Ritchie EG, Hebblewhite M, et al. Status and Ecological Effects of the World’s Largest Carnivores. Science. 2014;343(6167):1241484. pmid:24408439
- View Article
- PubMed/NCBI
- Google Scholar
2. Terborgh J. Maintenance of Diversity in Tropical Forests. Biotropica. 1992;24(2):283–92.
- View Article
- Google Scholar
3. Ross J, Hearn AJ, Johnson PJ, Macdonald DW. Activity patterns and temporal avoidance by prey in response to Sunda clouded leopard predation risk. J Zool. 2013;290(2):96–106.
- View Article
- Google Scholar
4. Paviolo A, De Angelo CD, Di Blanco YE, Di Bitetti MS. Jaguar Panthera onca population decline in the upper Parana Atlantic forrest of Argentina and Brazil. Oryx. 2008;42(4):554–61.
- View Article
- Google Scholar
5. Datta A, Anand MO, Naniwadekar R. Empty forests: Large carnivore and prey abundance in Namdapha National Park, north-east India. Biol Cons. 2008;141(5):1429–35. WOS:000257008400024.
- View Article
- Google Scholar
6. Mills LS. Conservaton of Wildlife Populations: Demography, Genetics and Management. 2 ed. West Sussex, UK: Wiley-Blackwell; 2013. 342 p.
7. Williams BK, Nichols JD, Conroy MJ. Analysis and management of animal populations. San Diego, CA, U.S.A.: Academic Press; 2002.
8. Sanei A, Zakaria M, Yusof E, Roslan M. Estimation of leopard population size in a secondary forest within Malaysia's capital agglomeration using unsupervised classification of pugmarks. Trop Ecol. 2011;52(2):209–17.
- View Article
- Google Scholar
9. O'Connell AF, Nichols JD, Karanth KU. Camera traps in animal ecology: methods and analyses. Tokyo: Springer; 2011.
10. Karanth KU. Estimating tiger Panthera tigris populations from camera-trap data using capture-recapture models. Biol Cons. 1995;71(3):333–8. WOS:A1995QH95100015.
- View Article
- Google Scholar
11. Karanth KU, Nichols JD. Estimation of tiger densities in India using photographic captures and recaptures. Ecology. 1998;79(8):2852–62.
- View Article
- Google Scholar
12. Tobler MW, Powell GV. Estimating jaguar densities with camera traps: Problems with current designs and recommendations for future studies. Biol Cons. 2013;159:109–18.
- View Article
- Google Scholar
13. Efford M. Density estimation in live—trapping studies. Oikos. 2004;106(3):598–610.
- View Article
- Google Scholar
14. Royle JA, Young KV. A hierarchical model for spatial capture-recapture data. Ecology. 2008;89(8):2281–9. pmid:18724738
- View Article
- PubMed/NCBI
- Google Scholar
15. Efford MG, Borchers DL, Byrom AE. Density estimation by spatially explicit capture—recapture: likelihood-based methods. In: Thomson DL, Cooch EG, Conroy MJ, editors. Modeling demographic processes in marked populations. New York, NY, U.S.A.: Springer; 2009. p. 255–69.
16. Borchers DL, Efford M. Spatially explicit maximum likelihood methods for capture—recapture studies. Biometrics. 2008;64(2):377–85. pmid:17970815
- View Article
- PubMed/NCBI
- Google Scholar
17. Royle JA, Nichols JD, Karanth KU, Gopalaswamy AM. A hierarchical model for estimating density in camera—trap studies. J Appl Ecol. 2009;46(1):118–27.
- View Article
- Google Scholar
18. Sollman R, Gardner B, Belant JL. How does spatial study design influence density estimates from spatial capture-recapture models? PLOS ONE. 2012;7(4):e34575. pmid:22539949
- View Article
- PubMed/NCBI
- Google Scholar
19. Sun CC, Fuller AK, Royle JA. Trap configuration and spacing influences parameter estimates in spatial capture-recapture models. PLOS ONE. 2014;9(2):e88025. pmid:24505361
- View Article
- PubMed/NCBI
- Google Scholar
20. Wilton CM, Puckett EE, Beringer J, Gardner B, Eggert LS, Belant JL. Trap Array Configuration Influences Estimates and Precision of Black Bear Density and Abundance. PLOS ONE. 2014;9(10):e111257. pmid:25350557
- View Article
- PubMed/NCBI
- Google Scholar
21. Sollmann R, Furtado MM, Gardner B, Hofer H, Jácomo AT, Tôrres NM, et al. Improving density estimates for elusive carnivores: accounting for sex-specific detection and movements using spatial capture—recapture models for jaguars in central Brazil. Biol Cons. 2011;144(3):1017–24.
- View Article
- Google Scholar
22. Royle JA, Chandler RB, Sun CC, Fuller AK. Integrating resource selection information with spatial capture-recapture. Method Ecol Evol. 2013;4(6):520–30. WOS:000320117200003.
- View Article
- Google Scholar
23. Russell RE, Royle JA, Desimone R, Schwartz MK, Edwards VL, Pilgrim KP, et al. Estimating abundance of mountain lions from unstructured spatial sampling. J Wildl Manage. 2012;76(8):1551–61. WOS:000310252300003.
- View Article
- Google Scholar
24. Foster RJ, Harmsen BJ. A critique of density estimation from camera—trap data. J Wildl Manage. 2012;76(2):224–36.
- View Article
- Google Scholar
25. Maffei L, Noss AJ, Silver SC, Kelly MJ. Abundance/density case study: Jaguars in the Americas. In: O'Connell AF, Nichols JD, Karanth KU, editors. Camera traps in animal ecology: Methods and analyses. New York, NY, U.S.A.: Springer; 2011. p. 119–44.
26. Borchers D, Distiller G, Foster RJ, Harmsen BJ, Milazzo L. Continuous-time spatially explicit capture-recapture models with an application to jaguar camera-trap survey. Method Ecol Evol. 2014;7:656–65.
- View Article
- Google Scholar
27. Bailey T. The African Leopard: Ecology and Behavior of a Solitary Felid. New York, NY, U.S.A.: Columbia University Press; 1993.
28. Uphyrkina O, Johnson WE, Quigley H, Miquelle D, Marker L, Bush M, et al. Phylogenetics, genome diversity and origin of modern leopard, Panthera pardus. Mol Ecol. 2001;10(11):2617–33. pmid:11883877
- View Article
- PubMed/NCBI
- Google Scholar
29. Wang SW, Macdonald D. Livestock predation by carnivores in Jigme Singye Wangchuck National Park, Bhutan. Biol Cons. 2006;129(4):558–65.
- View Article
- Google Scholar
30. Sangay T, Vernes K. Human—wildlife conflict in the Kingdom of Bhutan: patterns of livestock predation by large mammalian carnivores. Biol Cons. 2008;141(5):1272–82.
- View Article
- Google Scholar
31. DoFPS. Departmental report on wildlife crimes in Bhutan. In: Services DoFaP, editor. Thimphu, Bhutan: Ministry of Agriculture and Forests, Royal Government of Bhutan; 2014.
32. Wang SW, Macdonald DW. The use of camera traps for estimating tiger and leopard populations in the high altitude mountains of Bhutan. Biol Cons. 2009;142(3):606–13. WOS:000263742500014.
- View Article
- Google Scholar
33. Tempa T, Hebblewhite M, Mills LS, Wangchuk TR, Norbu N, Wangchuk T, et al. Royal Manas National Park, Bhutan: a hot spot for wild felids. Oryx. 2013;47(02):207–10.
- View Article
- Google Scholar
34. Gardner B, Royle JA, Wegan MT. Hierarchical models for estimating density from DNA mark-recapture studies. Ecology. 2009;90(4):1106–15. WOS:000264753400027. pmid:19449704
- View Article
- PubMed/NCBI
- Google Scholar
35. Royle JA, Karanth KU, Gopalaswamy AM, Kumar NS. Bayesian inference in camera trapping studies for a class of spatial capture-recapture models. Ecology. 2009;90(11):3233–44. WOS:000271457300026. pmid:19967878
- View Article
- PubMed/NCBI
- Google Scholar
36. Gardner B, Royle JA, Wegan MT, Rainbolt RE, Curtis PD. Estimating Black Bear Density Using DNA Data From Hair Snares. J Wildl Manage. 2010;74(2):318–25. WOS:000274243800018.
- View Article
- Google Scholar
37. Royle JA, Dorazio RM, Link WA. Analysis of multinomial models with unknown index using data augmentation. J Comput Graph Stat. 2007;16(1):67–85. WOS:000244878400004.
- View Article
- Google Scholar
38. Royle JA, Dorazio RM. Hierarchical modeling and inference in ecology: the analysis of data from populations, metapopulations and communities. San Diego, CA, U.S.A.: Academic Press; 2008.
39. R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for statistical computing; 2013.
40. Gray TNE, Prum S. Leopard density in post—conflict landscape, Cambodia: Evidence from spatially explicit capture—recapture. J Wildl Manage. 2012;76(1):163–9.
- View Article
- Google Scholar
41. Gelfand AE, Dey DK. Bayesian model choice: asymptotics and exact calculations. J Roy Stat Soc B Met. 1994:501–14.
- View Article
- Google Scholar
42. Kass RE, Raftery AE. Bayes factors. J Am Stat Assoc. 1995;90(430):773–95.
- View Article
- Google Scholar
43. Liu C, Liu Q. Marginal likelihood calculation for the Gelfand—Dey and Chib methods. Econ Lett. 2012;115(2):200–3.
- View Article
- Google Scholar
44. Jeffreys H, editor Some tests of significance, treated by the theory of probability. Proceedings of the Cambridge Philosophical Society; 1935: Cambridge Univ Press.
45. Jeffreys H. The theory of probability. Oxford, UK: Oxford University Press; 1961.
46. Bischof R, Hameed S, Ali H, Kabir M, Younas M, Shah KA, et al. Using time-to-event analysis to complement hierarchical methods when assessing determinants of photographic detectability during camera trapping. Methods Ecol Evol. 2014;5:44–53.
- View Article
- Google Scholar
47. Gelman A, Rubin DB. Avoiding model selection in Bayesian social research. Sociol Methodol. 1995;25:165–74.
- View Article
- Google Scholar
48. Royle JA, Chandler RB, Sollman R, Gardner B. Spatial Capture Recapture. Waltham, MA, USA: Academic Press; 2013.
49. Borah J, Sharma T, Das D, Rabha N, Kakati N, Basumatary A, et al. Abundance and density estimates for common leopard Panthera pardus and clouded leopard Neofelis nebulosa in Manas National Park, Assam, India. Oryx. 2013:1–7.
- View Article
- Google Scholar
50. Hooten M, Hobbs N. A guide to Bayesian model selection for ecologists. Ecol Monogr. 2015;85(1):3–28.
- View Article
- Google Scholar
51. Spiegelhalter DJ, Best NG, Carlin BP, Van Der Linde A. Bayesian measures of model complexity and fit. J Roy Stat Soc B Met. 2002;64(4):583–639.
- View Article
- Google Scholar
52. Watanabe S. Asymptotic equivalence of Bayes cross validation and widely applicable information criterion in singular learning theory. J Mach Learn Res. 2010;11:3571–94.
- View Article
- Google Scholar
53. Gelman A, Hwang J, Vehtari A. Understanding predictive information criteria for Bayesian models. Stat Comput. 2014;24(6):997–1016.
- View Article
- Google Scholar
54. Ando T, Tsay R. Predictive likelihood for Bayesian model selection and averaging. Int J Forecasting. 2010;26(4):744–63. http://dx.doi.org/10.1016/j.ijforecast.2009.08.001.
- View Article
- Google Scholar
55. Gelfand AE, Ghosh SK. Model choice: A minimum posterior predictive loss approach. Biometrika. 1998;85(1):1–11.
- View Article
- Google Scholar
56. Ibrahim JG, Chen M-H, Sinha D. Criterion-based methods for Bayesian model assessment. Stat Sinica. 2001;11(2):419–44.
- View Article
- Google Scholar
57. Kuo L, Mallick B. Variable selection for regression models. Sankhya Ser B. 1998:65–81.
- View Article
- Google Scholar
58. Carlin BP, Chib S. Bayesian model choice via Markov chain Monte Carlo methods. J Roy Stat Soc B Met. 1995:473–84.
- View Article
- Google Scholar
59. Green PJ. Reversible jump Markov chain Monte Carlo computation and Bayesian model determination. Biometrika. 1995;82(4):711–32.
- View Article
- Google Scholar
60. Barker RJ, Link WA. Bayesian Multimodel Inference by RJMCMC: A Gibbs Sampling Approach. Am Stat. 2013;67(3):150–6.
- View Article
- Google Scholar
61. Gelman A, Shalizi CR. Philosophy and the practice of Bayesian statistics. Brit J Math Stat Psy. 2013;66(1):8–38.
- View Article
- Google Scholar
62. Odden M, Wegge P. Spacing and activity patterns of leopards Panthera pardus in the Royal Bardia National Park, Nepal. Wildlife Biol. 2005;11(2):145–52.
- View Article
- Google Scholar
63. Simcharoen S, Barlow ACD, Simcharoen A, Smith JLD. Home range size and daytime habitat selection of leopards in Huai Kha Khaeng Wildlife Sanctuary, Thailand. Biol Cons. 2008;141(9):2242–50. WOS:000260027100008.
- View Article
- Google Scholar
64. Harihar A, Pandav B, Goyal SP. Responses of leopard Panthera pardus to the recovery of a tiger Panthera tigris population. J Appl Ecol. 2011;48(3):806–14. WOS:000290587100034.
- View Article
- Google Scholar
65. Kalle R, Ramesh T, Qureshi Q, Sankar K. Density of tiger and leopard in a tropical deciduous forest of Mudumalai Tiger Reserve, southern India, as estimated using photographic capture—recapture sampling. Acta Theriol. 2011;56(4):335–42.
- View Article
- Google Scholar
66. Proffitt K, Goldberg JF, Hebblewhite M, Russell RE, Jimenez B, Robinson H, et al. Integrating resource selection in spatial capture-recapture models for large carnivores. Ecography. in revision.
67. Royle JA, Chandler RB, Gazenski KD, Graves TA. Spatial capture-recapture for jointly estimating population density and landscape connectivity. Ecology. 2013;94:287–94. pmid:23691647
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Ripple WJ, Estes JA, Beschta RL, Wilmers CC, Ritchie EG, Hebblewhite M, et al. Status and Ecological Effects of the World’s Largest Carnivores. Science. 2014;343(6167):1241484. pmid:24408439
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Terborgh J. Maintenance of Diversity in Tropical Forests. Biotropica. 1992;24(2):283–92.
View Article
Google Scholar

[6] View Article

[7] Google Scholar

[ref3] 3. Ross J, Hearn AJ, Johnson PJ, Macdonald DW. Activity patterns and temporal avoidance by prey in response to Sunda clouded leopard predation risk. J Zool. 2013;290(2):96–106.
View Article
Google Scholar

[9] View Article

[10] Google Scholar

[ref4] 4. Paviolo A, De Angelo CD, Di Blanco YE, Di Bitetti MS. Jaguar Panthera onca population decline in the upper Parana Atlantic forrest of Argentina and Brazil. Oryx. 2008;42(4):554–61.
View Article
Google Scholar

[12] View Article

[13] Google Scholar

[ref5] 5. Datta A, Anand MO, Naniwadekar R. Empty forests: Large carnivore and prey abundance in Namdapha National Park, north-east India. Biol Cons. 2008;141(5):1429–35. WOS:000257008400024.
View Article
Google Scholar

[15] View Article

[16] Google Scholar

[ref6] 6. Mills LS. Conservaton of Wildlife Populations: Demography, Genetics and Management. 2 ed. West Sussex, UK: Wiley-Blackwell; 2013. 342 p.

[ref7] 7. Williams BK, Nichols JD, Conroy MJ. Analysis and management of animal populations. San Diego, CA, U.S.A.: Academic Press; 2002.

[ref8] 8. Sanei A, Zakaria M, Yusof E, Roslan M. Estimation of leopard population size in a secondary forest within Malaysia's capital agglomeration using unsupervised classification of pugmarks. Trop Ecol. 2011;52(2):209–17.
View Article
Google Scholar

[20] View Article

[21] Google Scholar

[ref9] 9. O'Connell AF, Nichols JD, Karanth KU. Camera traps in animal ecology: methods and analyses. Tokyo: Springer; 2011.

[ref10] 10. Karanth KU. Estimating tiger Panthera tigris populations from camera-trap data using capture-recapture models. Biol Cons. 1995;71(3):333–8. WOS:A1995QH95100015.
View Article
Google Scholar

[24] View Article

[25] Google Scholar

[ref11] 11. Karanth KU, Nichols JD. Estimation of tiger densities in India using photographic captures and recaptures. Ecology. 1998;79(8):2852–62.
View Article
Google Scholar

[27] View Article

[28] Google Scholar

[ref12] 12. Tobler MW, Powell GV. Estimating jaguar densities with camera traps: Problems with current designs and recommendations for future studies. Biol Cons. 2013;159:109–18.
View Article
Google Scholar

[30] View Article

[31] Google Scholar

[ref13] 13. Efford M. Density estimation in live—trapping studies. Oikos. 2004;106(3):598–610.
View Article
Google Scholar

[33] View Article

[34] Google Scholar

[ref14] 14. Royle JA, Young KV. A hierarchical model for spatial capture-recapture data. Ecology. 2008;89(8):2281–9. pmid:18724738
View Article
PubMed/NCBI
Google Scholar

[36] View Article

[37] PubMed/NCBI

[38] Google Scholar

[ref15] 15. Efford MG, Borchers DL, Byrom AE. Density estimation by spatially explicit capture—recapture: likelihood-based methods. In: Thomson DL, Cooch EG, Conroy MJ, editors. Modeling demographic processes in marked populations. New York, NY, U.S.A.: Springer; 2009. p. 255–69.

[ref16] 16. Borchers DL, Efford M. Spatially explicit maximum likelihood methods for capture—recapture studies. Biometrics. 2008;64(2):377–85. pmid:17970815
View Article
PubMed/NCBI
Google Scholar

[41] View Article

[42] PubMed/NCBI

[43] Google Scholar

[ref17] 17. Royle JA, Nichols JD, Karanth KU, Gopalaswamy AM. A hierarchical model for estimating density in camera—trap studies. J Appl Ecol. 2009;46(1):118–27.
View Article
Google Scholar

[45] View Article

[46] Google Scholar

[ref18] 18. Sollman R, Gardner B, Belant JL. How does spatial study design influence density estimates from spatial capture-recapture models? PLOS ONE. 2012;7(4):e34575. pmid:22539949
View Article
PubMed/NCBI
Google Scholar

[48] View Article

[49] PubMed/NCBI

[50] Google Scholar

[ref19] 19. Sun CC, Fuller AK, Royle JA. Trap configuration and spacing influences parameter estimates in spatial capture-recapture models. PLOS ONE. 2014;9(2):e88025. pmid:24505361
View Article
PubMed/NCBI
Google Scholar

[52] View Article

[53] PubMed/NCBI

[54] Google Scholar

[ref20] 20. Wilton CM, Puckett EE, Beringer J, Gardner B, Eggert LS, Belant JL. Trap Array Configuration Influences Estimates and Precision of Black Bear Density and Abundance. PLOS ONE. 2014;9(10):e111257. pmid:25350557
View Article
PubMed/NCBI
Google Scholar

[56] View Article

[57] PubMed/NCBI

[58] Google Scholar

[ref21] 21. Sollmann R, Furtado MM, Gardner B, Hofer H, Jácomo AT, Tôrres NM, et al. Improving density estimates for elusive carnivores: accounting for sex-specific detection and movements using spatial capture—recapture models for jaguars in central Brazil. Biol Cons. 2011;144(3):1017–24.
View Article
Google Scholar

[60] View Article

[61] Google Scholar

[ref22] 22. Royle JA, Chandler RB, Sun CC, Fuller AK. Integrating resource selection information with spatial capture-recapture. Method Ecol Evol. 2013;4(6):520–30. WOS:000320117200003.
View Article
Google Scholar

[63] View Article

[64] Google Scholar

[ref23] 23. Russell RE, Royle JA, Desimone R, Schwartz MK, Edwards VL, Pilgrim KP, et al. Estimating abundance of mountain lions from unstructured spatial sampling. J Wildl Manage. 2012;76(8):1551–61. WOS:000310252300003.
View Article
Google Scholar

[66] View Article

[67] Google Scholar

[ref24] 24. Foster RJ, Harmsen BJ. A critique of density estimation from camera—trap data. J Wildl Manage. 2012;76(2):224–36.
View Article
Google Scholar

[69] View Article

[70] Google Scholar

[ref25] 25. Maffei L, Noss AJ, Silver SC, Kelly MJ. Abundance/density case study: Jaguars in the Americas. In: O'Connell AF, Nichols JD, Karanth KU, editors. Camera traps in animal ecology: Methods and analyses. New York, NY, U.S.A.: Springer; 2011. p. 119–44.

[ref26] 26. Borchers D, Distiller G, Foster RJ, Harmsen BJ, Milazzo L. Continuous-time spatially explicit capture-recapture models with an application to jaguar camera-trap survey. Method Ecol Evol. 2014;7:656–65.
View Article
Google Scholar

[73] View Article

[74] Google Scholar

[ref27] 27. Bailey T. The African Leopard: Ecology and Behavior of a Solitary Felid. New York, NY, U.S.A.: Columbia University Press; 1993.

[ref28] 28. Uphyrkina O, Johnson WE, Quigley H, Miquelle D, Marker L, Bush M, et al. Phylogenetics, genome diversity and origin of modern leopard, Panthera pardus. Mol Ecol. 2001;10(11):2617–33. pmid:11883877
View Article
PubMed/NCBI
Google Scholar

[77] View Article

[78] PubMed/NCBI

[79] Google Scholar

[ref29] 29. Wang SW, Macdonald D. Livestock predation by carnivores in Jigme Singye Wangchuck National Park, Bhutan. Biol Cons. 2006;129(4):558–65.
View Article
Google Scholar

[81] View Article

[82] Google Scholar

[ref30] 30. Sangay T, Vernes K. Human—wildlife conflict in the Kingdom of Bhutan: patterns of livestock predation by large mammalian carnivores. Biol Cons. 2008;141(5):1272–82.
View Article
Google Scholar

[84] View Article

[85] Google Scholar

[ref31] 31. DoFPS. Departmental report on wildlife crimes in Bhutan. In: Services DoFaP, editor. Thimphu, Bhutan: Ministry of Agriculture and Forests, Royal Government of Bhutan; 2014.

[ref32] 32. Wang SW, Macdonald DW. The use of camera traps for estimating tiger and leopard populations in the high altitude mountains of Bhutan. Biol Cons. 2009;142(3):606–13. WOS:000263742500014.
View Article
Google Scholar

[88] View Article

[89] Google Scholar

[ref33] 33. Tempa T, Hebblewhite M, Mills LS, Wangchuk TR, Norbu N, Wangchuk T, et al. Royal Manas National Park, Bhutan: a hot spot for wild felids. Oryx. 2013;47(02):207–10.
View Article
Google Scholar

[91] View Article

[92] Google Scholar

[ref34] 34. Gardner B, Royle JA, Wegan MT. Hierarchical models for estimating density from DNA mark-recapture studies. Ecology. 2009;90(4):1106–15. WOS:000264753400027. pmid:19449704
View Article
PubMed/NCBI
Google Scholar

[94] View Article

[95] PubMed/NCBI

[96] Google Scholar

[ref35] 35. Royle JA, Karanth KU, Gopalaswamy AM, Kumar NS. Bayesian inference in camera trapping studies for a class of spatial capture-recapture models. Ecology. 2009;90(11):3233–44. WOS:000271457300026. pmid:19967878
View Article
PubMed/NCBI
Google Scholar

[98] View Article

[99] PubMed/NCBI

[100] Google Scholar

[ref36] 36. Gardner B, Royle JA, Wegan MT, Rainbolt RE, Curtis PD. Estimating Black Bear Density Using DNA Data From Hair Snares. J Wildl Manage. 2010;74(2):318–25. WOS:000274243800018.
View Article
Google Scholar

[102] View Article

[103] Google Scholar

[ref37] 37. Royle JA, Dorazio RM, Link WA. Analysis of multinomial models with unknown index using data augmentation. J Comput Graph Stat. 2007;16(1):67–85. WOS:000244878400004.
View Article
Google Scholar

[105] View Article

[106] Google Scholar

[ref38] 38. Royle JA, Dorazio RM. Hierarchical modeling and inference in ecology: the analysis of data from populations, metapopulations and communities. San Diego, CA, U.S.A.: Academic Press; 2008.

[ref39] 39. R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for statistical computing; 2013.

[ref40] 40. Gray TNE, Prum S. Leopard density in post—conflict landscape, Cambodia: Evidence from spatially explicit capture—recapture. J Wildl Manage. 2012;76(1):163–9.
View Article
Google Scholar

[110] View Article

[111] Google Scholar

[ref41] 41. Gelfand AE, Dey DK. Bayesian model choice: asymptotics and exact calculations. J Roy Stat Soc B Met. 1994:501–14.
View Article
Google Scholar

[113] View Article

[114] Google Scholar

[ref42] 42. Kass RE, Raftery AE. Bayes factors. J Am Stat Assoc. 1995;90(430):773–95.
View Article
Google Scholar

[116] View Article

[117] Google Scholar

[ref43] 43. Liu C, Liu Q. Marginal likelihood calculation for the Gelfand—Dey and Chib methods. Econ Lett. 2012;115(2):200–3.
View Article
Google Scholar

[119] View Article

[120] Google Scholar

[ref44] 44. Jeffreys H, editor Some tests of significance, treated by the theory of probability. Proceedings of the Cambridge Philosophical Society; 1935: Cambridge Univ Press.

[ref45] 45. Jeffreys H. The theory of probability. Oxford, UK: Oxford University Press; 1961.

[ref46] 46. Bischof R, Hameed S, Ali H, Kabir M, Younas M, Shah KA, et al. Using time-to-event analysis to complement hierarchical methods when assessing determinants of photographic detectability during camera trapping. Methods Ecol Evol. 2014;5:44–53.
View Article
Google Scholar

[124] View Article

[125] Google Scholar

[ref47] 47. Gelman A, Rubin DB. Avoiding model selection in Bayesian social research. Sociol Methodol. 1995;25:165–74.
View Article
Google Scholar

[127] View Article

[128] Google Scholar

[ref48] 48. Royle JA, Chandler RB, Sollman R, Gardner B. Spatial Capture Recapture. Waltham, MA, USA: Academic Press; 2013.

[ref49] 49. Borah J, Sharma T, Das D, Rabha N, Kakati N, Basumatary A, et al. Abundance and density estimates for common leopard Panthera pardus and clouded leopard Neofelis nebulosa in Manas National Park, Assam, India. Oryx. 2013:1–7.
View Article
Google Scholar

[131] View Article

[132] Google Scholar

[ref50] 50. Hooten M, Hobbs N. A guide to Bayesian model selection for ecologists. Ecol Monogr. 2015;85(1):3–28.
View Article
Google Scholar

[134] View Article

[135] Google Scholar

[ref51] 51. Spiegelhalter DJ, Best NG, Carlin BP, Van Der Linde A. Bayesian measures of model complexity and fit. J Roy Stat Soc B Met. 2002;64(4):583–639.
View Article
Google Scholar

[137] View Article

[138] Google Scholar

[ref52] 52. Watanabe S. Asymptotic equivalence of Bayes cross validation and widely applicable information criterion in singular learning theory. J Mach Learn Res. 2010;11:3571–94.
View Article
Google Scholar

[140] View Article

[141] Google Scholar

[ref53] 53. Gelman A, Hwang J, Vehtari A. Understanding predictive information criteria for Bayesian models. Stat Comput. 2014;24(6):997–1016.
View Article
Google Scholar

[143] View Article

[144] Google Scholar

[ref54] 54. Ando T, Tsay R. Predictive likelihood for Bayesian model selection and averaging. Int J Forecasting. 2010;26(4):744–63. http://dx.doi.org/10.1016/j.ijforecast.2009.08.001.
View Article
Google Scholar

[146] View Article

[147] Google Scholar

[ref55] 55. Gelfand AE, Ghosh SK. Model choice: A minimum posterior predictive loss approach. Biometrika. 1998;85(1):1–11.
View Article
Google Scholar

[149] View Article

[150] Google Scholar

[ref56] 56. Ibrahim JG, Chen M-H, Sinha D. Criterion-based methods for Bayesian model assessment. Stat Sinica. 2001;11(2):419–44.
View Article
Google Scholar

[152] View Article

[153] Google Scholar

[ref57] 57. Kuo L, Mallick B. Variable selection for regression models. Sankhya Ser B. 1998:65–81.
View Article
Google Scholar

[155] View Article

[156] Google Scholar

[ref58] 58. Carlin BP, Chib S. Bayesian model choice via Markov chain Monte Carlo methods. J Roy Stat Soc B Met. 1995:473–84.
View Article
Google Scholar

[158] View Article

[159] Google Scholar

[ref59] 59. Green PJ. Reversible jump Markov chain Monte Carlo computation and Bayesian model determination. Biometrika. 1995;82(4):711–32.
View Article
Google Scholar

[161] View Article

[162] Google Scholar

[ref60] 60. Barker RJ, Link WA. Bayesian Multimodel Inference by RJMCMC: A Gibbs Sampling Approach. Am Stat. 2013;67(3):150–6.
View Article
Google Scholar

[164] View Article

[165] Google Scholar

[ref61] 61. Gelman A, Shalizi CR. Philosophy and the practice of Bayesian statistics. Brit J Math Stat Psy. 2013;66(1):8–38.
View Article
Google Scholar

[167] View Article

[168] Google Scholar

[ref62] 62. Odden M, Wegge P. Spacing and activity patterns of leopards Panthera pardus in the Royal Bardia National Park, Nepal. Wildlife Biol. 2005;11(2):145–52.
View Article
Google Scholar

[170] View Article

[171] Google Scholar

[ref63] 63. Simcharoen S, Barlow ACD, Simcharoen A, Smith JLD. Home range size and daytime habitat selection of leopards in Huai Kha Khaeng Wildlife Sanctuary, Thailand. Biol Cons. 2008;141(9):2242–50. WOS:000260027100008.
View Article
Google Scholar

[173] View Article

[174] Google Scholar

[ref64] 64. Harihar A, Pandav B, Goyal SP. Responses of leopard Panthera pardus to the recovery of a tiger Panthera tigris population. J Appl Ecol. 2011;48(3):806–14. WOS:000290587100034.
View Article
Google Scholar

[176] View Article

[177] Google Scholar

[ref65] 65. Kalle R, Ramesh T, Qureshi Q, Sankar K. Density of tiger and leopard in a tropical deciduous forest of Mudumalai Tiger Reserve, southern India, as estimated using photographic capture—recapture sampling. Acta Theriol. 2011;56(4):335–42.
View Article
Google Scholar

[179] View Article

[180] Google Scholar

[ref66] 66. Proffitt K, Goldberg JF, Hebblewhite M, Russell RE, Jimenez B, Robinson H, et al. Integrating resource selection in spatial capture-recapture models for large carnivores. Ecography. in revision.

[ref67] 67. Royle JA, Chandler RB, Gazenski KD, Graves TA. Spatial capture-recapture for jointly estimating population density and landscape connectivity. Ecology. 2013;94:287–94. pmid:23691647
View Article
PubMed/NCBI
Google Scholar

[183] View Article

[184] PubMed/NCBI

[185] Google Scholar

Figures

Abstract

Introduction

Materials and Methods

Study Area

Ethics Statement

Field Methods

Analytical Methods

Bayesian Model Analysis

Model Selection

Simulations

Results

Camera Trapping

Model Results

Simulations

Discussion

Sample Interval Length and SCR Estimates

Bayesian Model Selection

Ecological Implications

Conclusions

Supporting Information

S1 File. R Code to Implement Spatial Capture Recapture Analysis with Bayes Factors Model Selection.

S2 File. Internal R Code to Conduct Bayes Factors Model Selection.

S1 Table. Parameter Estimates and 95% Credibility Intervals from All Spatial Capture-Recapture Models for Common Leopards.

S2 Table. Median Parameter Estimates and 95% Credibility Intervals from All Spatial Capture-Recapture Models from Simulations.

Acknowledgments

Author Contributions

References