Spatial Mark-Resight for Categorically Marked Populations with an Application to Genetic Capture-Recapture

The estimation of animal population density is a fundamental goal in wildlife ecology and management, commonly met using mark recapture or spatial mark recapture (SCR) study designs and statistical methods. Mark-recapture methods require the identification of individuals; however, for many species and sampling methods, particularly noninvasive methods, no individuals or only a subset of individuals are individually identifiable. The unmarked SCR model, theoretically, can estimate the density of unmarked populations; however, it produces biased and imprecise density estimates in many sampling scenarios typically encountered. Spatial mark-resight (SMR) models extend the unmarked SCR model in three ways: 1) by introducing a subset of individuals that are marked and individually identifiable, 2) introducing the possibility of individual-linked telemetry data, and 3) introducing the possibility that the capture-recapture data from the survey used to deploy the marks can be used in a joint model, all improving the reliability of density estimates. The categorical spatial partial identity model (SPIM) improves the reliability of density estimates over unmarked SCR along another dimension, by adding categorical identity covariates that improve the probabilistic association of the latent identity samples. Here, we combine these two models into a "categorical SMR" model to exploit the benefits of both models simultaneously. We demonstrate using simulations that SMR alone can produce biased and imprecise density estimates with sparse data and/or when few individuals are marked. Then, using a fisher (Pekania pennanti) genetic capture-recapture data set, we show how categorical identity covariates, marked individuals, telemetry data, and jointly modeling the capture survey used to deploy marks with the resighting survey all combine to improve inference over the unmarked SCR model. As previously seen in an application of the categorical SPIM to a real-world data set, the fisher data set demonstrates that individual heterogeneity in detection function parameters, especially the spatial scale parameter sigma, introduces positive bias into latent identity SCR models (e.g., unmarked SCR, SMR), but the categorical SMR model provides more tools to reduce this positive bias than SMR or the categorical SPIM alone. We introduce the possibility of detection functions that vary by identity category level, which will remove individual heterogeneity in detection function parameters than is explained by categorical covariates, such as individual sex. Finally, we provide efficient SMR algorithms that accommodate all SMR sample types, interspersed marking and sighting periods, and any number of identity covariates using the 2-dimensional individual by trap data in conjunction with precomputed constraint matrices, rather than the 3-dimensional individual by trap by occasion data used in SMR algorithms to date.