Final From journal on website
- 1. BioOne sees sustainable scholarly publishing as an inherently collaborative enterprise connecting authors, nonprofit publishers, academic institutions, research libraries, and research funders in the common goal of maximizing access to critical research. The utility of auxiliary data in statistical population reconstruction Author(s): Michael V. Clawson, John R. Skalski & Joshua J. Millspaugh Source: Wildlife Biology, 19(2):147-155. Published By: Nordic Board for Wildlife Research DOI: http://dx.doi.org/10.2981/12-076 URL: http://www.bioone.org/doi/full/10.2981/12-076 BioOne (www.bioone.org) is a nonprofit, online aggregation of core research in the biological, ecological, and environmental sciences. BioOne provides a sustainable online platform for over 170 journals and books published by nonprofit societies, associations, museums, institutions, and presses. Your use of this PDF, the BioOne Web site, and all posted and associated content indicates your acceptance of BioOne’s Terms of Use, available at www.bioone.org/page/terms_of_use. Usage of BioOne content is strictly limited to personal, educational, and non-commercial use. Commercial inquiries or rights and permissions requests should be directed to the individual publisher as copyright holder.
- 2. Original articleWildl. Biol. 19: 147-155 (2013) DOI: 10.2981/12-076 Ó Wildlife Biology, NKV www.wildlifebiology.com The utility of auxiliary data in statistical population reconstruction Michael V. Clawson, John R. Skalski & Joshua J. Millspaugh Although statistical population reconstruction (SPR) provides a flexible framework for estimating demographics of harvested populations using age-at-harvest data, that information alone is insufficient. Auxiliary data are needed to ensure all model parameters are estimable and to improve the precision and accuracy of the estimates. We examined the influence of two types of auxiliary information, independent estimates of annual abundance and annual harvest mortality from radio-telemetry studies, on the stability and precision of abundance estimates from SPR. Further, we evaluated whether the timing of auxiliary studies in the reconstruction affected the precision of abundance estimates. Monte Carlo studies simulated auxiliary data with precision levels defined by the coefficients of variation (CV) of 0.05, 0.125, 0.25 and 0.50 corresponding to the three levels of precision suggested by Robson & Regier (1964) for accurate research, accurate management and rough management and a minimal information scenario. For comparable levels of precision, radio-telemetry studies used to estimate harvest mortality stabilized the reconstructed population trends better than independent abundance surveys. However, independent abundance surveys were superior at improving the precision of reconstructed abundance estimates. We found that the timing of auxiliary studies did not influence the stability of SPR estimates, which has important implications for managers designing studies to collect auxiliary data. Our research highlights that different types and quality of auxiliary studies affects the precision and stability of SPR models differently. Key words: abundance estimation, age-at-harvest, harvest mortality, population reconstruction, precision, statistical population reconstruction Michael V. Clawson, School of Environmental and Forest Sciences, University of Washington, 1325 Fourth Avenue, Suite 1820, Seattle, Washington 98101-2509, USA - e-mail: michael.v.clawson@gmail.com John R.Skalski,School of AquaticandFishery Sciences, University of Washington,1325Fourth Avenue, Suite 1820, Seattle, Washington 98101-2509, USA - e-mail: skalski@uw.edu Joshua J. Millspaugh, Department of Fisheries and Wildlife Sciences, University of Missouri, 302 Natural Resources Building, Columbia, Missouri 65211, USA - e-mail: millspaughj@missouri.edu Corrresponding author: John R. Skalski Received 25 July 2012, accepted 14 January 2013 Associate Editor: Erlend B. Nielsen In recent years, statistical population reconstruction (SPR) has emerged as a feasible method for estimat- ing the demographics of harvested wildlife over large geographic areas using age-at-harvest data which are commonly collected by wildlife agencies (Gove et al. 2002, Skalski et al. 2007, 2011, Broms et al. 2010). Its origin can be found in the vast history of quantitative stock assessment in fisheries (Quinn & Deriso 1999). In wildlife science, where methods such as the Downing (1980) method and sex-age-kill (Mill- spaugh et al. 2009) are still the norm, there are advantages in using SPR which offers flexibility and robustness not available in traditional techniques. Further, SPR allows for the simultaneous estimation of survival, abundance, recruitment and harvest mortality whereas traditional reconstruction tech- niques estimate only total abundance. Age-at-harvest data provide the primary source of information in SPR models. However, age-at-har- vest data alone are insufficient to reconstruct popu- lation demographics using SPR. In addition to age- at-harvest data, one or more sets of auxiliary data are needed to estimate one or more of the parameters from the age-at-harvest likelihood, either survival Ó WILDLIFE BIOLOGY 19:2 (2013) 147
- 3. rates, harvest rates, recruitment or abundance. The joint likelihood structure of SPR models is flexible enough to incorporate almost any form of auxiliary data. In the past, it has been common to include catch-effort data as auxiliary information for SPR models. Skalski et al. (2007) calibrated a black-tailed deer Odocoileus hemionus SPR model using catch- effort data. That was possible because hunter effort was deliberately manipulated to produce a strong contrast in harvest rates with alternative levels of effort. Similarly, Skalski et al. (2011) used a five-fold change in trapping effort over time to construct a catch-effort relationship and reconstruct the abun- dance of American martens Martes americana in the upper peninsula of Michigan, USA. In other popu- lations where hunter or trap effort may be relatively constant over time, catch-effort data will likely be an insufficient form of auxiliary information. Thus, although catch-effort data has been a staple in SPR analyses other types of auxiliary data may be more useful. Several other types of auxiliary data have been used in previous SPR analyses including radio- telemetry information on harvest rates and indepen- dentestimatesofabundance.Bromsetal.(2010)used radio-telemetry to help reconstruct greater sage- grouse Centrocercus urophasianus abundance in Oregon, USA. Radio-telemetry data were used to estimate vulnerability coefficients associated with harvest mortality of greater sage-grouse in Oregon. Fieberg et al. (2010) also used radio-telemetry to estimateharvest rates and reconstruct the abundance of a Minnesota black bear Ursus americanus popu- lation. Other researchers have used independent estimates of total abundance as auxiliary data. For example, Gast et al. (submitted) used independent mark-recapture estimates of total abundance to help calibrate an SPR model of elk Cervus elaphus in the upper peninsula of Michigan. Alternatively, Fieberg et al. (2010) chose not to use independent DNA mark-recapture estimates of abundance when recon- structing the Michigan black bear population. In- stead, they chose to use that information as an independent source of confirmatory information. In contrast, Skalski et al. (2007) found that an indepen- dent browse damage index of deer abundance had little or no benefit in reconstructing a black-tailed deer herd. Such index data can help characterize the trend of a populationbut notits absoluteabundance. These examples illustrate the flexibility of SPR, but also raise questions about the potential utility of various types and quality of auxiliary data. Uncer- tainty remains as to the effect of auxiliary data on model stability and how much auxiliary data are needed. Also, it is unclear how the precision of the auxiliary data (i.e. coefficients of variation (CV)) changes the precision of the reconstructed abun- dance estimates. Additionally, it is unclear if the timingofauxiliarystudieschangestheireffectiveness. The purpose of our paper is to provide game managers with guidance on how best to incorporate auxiliary studies in SPR. We evaluate whether precision of the auxiliary study is the sole consider- ation, or whether the types of parameters being estimated are also important. To this end, we compare the performance of SPR where independent estimates of abundance or harvest mortality are available.Wealsoconsiderthetimingoftheauxiliary studies in relation to the duration of an SPR and the relative benefit of more than one auxiliary study in the precision of SPR. Methods Overview of SPR SPR is based on age-at-harvest data collected over time by game management agencies. The observed counts hij (i ¼ l,...,Y; j ¼ 0,..., A) are modeled as a function of the initial abundance of a cohort and the subsequent natural survival and harvest over time and perhaps probabilities of reporting and age determination. Skalski et al. (2007) modeled the diagonals of this age-at-harvest matrix as indepen- dent multinomial distributions (Lij) where the joint likelihood can be written as LAge-at-harvest ¼ YA j¼0 L1j Á YY i¼1 Li0: This likelihood is often accompanied by a likeli- hood model describing the annual probabilities of harvested animals being report and/or aged, i.e. LReporting ¼ YY i¼1 Li; where Li are binomial sampling models for the fractions of animals harvested in year i being reported and/or aged. This likelihood can be omitted if there is 100% reporting and aging of all harvested animals. Together, these two likelihoods are incapable of estimating the demographic parameters of interest, 148 Ó WILDLIFE BIOLOGY 19:2 (2013)
- 4. annual abundance, recruitment, natural survival and harvest mortality. Gove et al. (2002) proved that at least one demographic parameter must be estimated independent of the age-at-harvest data for SPR to be possible.Sometimestheremaybe,say,kindependent auxiliary studies contributing to the reconstruction; hence, the joint likelihood model may be written as L ¼ LAge-at-harvest Á LReporting Á Yk i¼1 LAuxillaryi ð1Þ: Skalski et al. (2007) suggested using a catch-effort likelihood as an auxiliary where the annual harvest numbers are modeled as binomial, random variables as a function of the unknown total abundance in year i (i.e. Ni) and hunter effort (fi) where LCatch-effort ¼ YY i¼1 Ni XA j¼0 hij 0 B @ 1 C ApðfiÞ XA j¼0 hij 3 ð1 - pð fiÞÞ Ni - XA j¼0 hij ð2Þ; and where p(f) is the probability of harvest modeled as a function of effort. A common parameterization for the probability of harvest is pi ¼ 1 - e -cfi ; where c is the vulnerability coefficient (Seber 1982:296, Quinn & Deriso 1999:40). Unless hunter effort has varied dramatically over time, this catch- effort auxiliary may not be adequate to support SPR. Auxuliary likelihoods We considered two alternative forms of auxiliary likelihoods in our evaluation. In one case, an unbiased annual abundance estimate ˆNi was as- sumed to be available with estimated standard error ˆri The auxiliary likelihood then assumed the estimate was asymptotically normally distributed with the likelihood LAux ¼ 1 ffiffiffiffiffiffiffiffiffiffi 2pri p e -1 2 ðNi- ˆNiÞ2 r2 i ð3Þ: We took this approach to allow the auxiliary likelihood to be independent of the form of the abundance surveyand simply reflect surveyprecision (i.e. CV¼ri/Ni). The second auxiliary approach was based on a hypothetical radio-telemetry study to estimateharvestmortalityduringthehuntingseason. Here we used a binomial of the form LAux ¼ Ti di ð1 - e -cfi Þ di ðe -cfi Þ Ti -di ð4Þ; where Ti ¼ number of animals tagged and at risk of harvest, di ¼ number of tagged animals harvested, and where harvest mortality was reparameterized in terms of a vulnerability coefficient (c) and annual year-specific hunter effort (equation 1). Precision in the case of the radio-telemetry study was expressedin terms of CV ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi pið1 - piÞ Ti s =pi: Monte Carlo simulations A Monte Carlo simulation study was used to determine the precision of population reconstruction estimates based on the amount, type and timing of auxiliary studies. A stochastic Leslie matrix model was used to generate age-at-harvest data for popu- lations with different levels of natural survival rates and harvest rates. Recruitment levels were adjusted to produce populations with stationary abundance of approximately 6,000 animals in expectation. Recruitment was generated using a Poisson process, and natural survival and harvest were modeled as binomial processes. Ineachsimulation, 20yearsofdataweregenerated to establish demographic trends with years 21-44 used in the population reconstruction analysis. The full age-class data were generated and used in standard population reconstruction models. The same data were also reanalyzed after pooling the adult age-at-harvest data (i.e. 2.5þ year olds) using the pooled adult reconstruction of Skalski et al. (2012). Demographic scenarios were performed to repre- sent a range of scenarios expected for harvested large mammal populations. Natural survival probabilities were simulated at 0.75 or 0.90 and harvest rates at 0.10 or 0.25. To minimize the number of scenarios investigated, survivalandharvest rateswereassumed constant across all age classes. Auxiliary data were simulated to estimate either the annual abundance ( ˆNi) or a harvest probability (ˆPi) with CV equal to 0.05, 0.125, 0.25 or 0.50. The CVs of 0.05, 0.125 and 0.25 correspond to precision levels described by Robson Regier (1964) as appropriate for accurate research, accurate management and rough manage- Ó WILDLIFE BIOLOGY 19:2 (2013) 149
- 5. ment. A fourthCVof0.50was simulated torepresent a minimum information scenario. At this level of precision, a parameter is estimated within6 100% of the true value 95% of the time. The effect of timing of the auxiliary data was tested for each parameter combination by staging the auxiliary study at either the beginning, middle, near end (i.e. year 23 of the reconstruction) or end (i.e. year 24) of the recon- struction. Average measurement error of the reconstructed abundance estimates was estimated from the vari- ance component expression Eðs 2 ˆNj - s 2 Nj Þ ¼ ðr 2 Nij þ Varð ˆNijjNijÞÞ - r 2 Nij ¼ Varð ˆNijjNijÞ; where s2 ˆNj ¼ P10;00 i¼1 ð ˆNij -N^ ijÞ2 /(10,000 - 1) is the empirical variance among the abundance estimates in the jth year of reconstruction and s2 Nj ¼P10;00 i¼1 ðNij - NijÞ2 /(10,000 - 1) is the empirical vari- ance among the true abundance values in the jth year of the reconstruction. A total of 10,000 simulations per scenario were used to obtain precise estimates of s2 ˆNj and s2 Nj . This approach provides a model-inde- pendent estimate of measurement error. Across the 24 years of reconstruction, reported precision was calculated in terms of median CV of measurement error. Black-tailed deer sensitivity analysis The previous Monte Carlo simulation studies looked at the relationship between the precision of abun- dance estimates from population reconstruction and the use of auxiliary data. This section examines the effect of auxiliary data on the stability of recon- structedpopulationtrends for oneparticular realized data set. The black-tailed deer reconstruction of Skalski et al. (2007) was selected for illustration because no auxiliary likelihood was incorporated in the original population reconstruction. Over a 24- year period, the abundance ranged between 1,500 and 3,500 does. Only catch-effort data were used to calibrate the model. The example is therefore conve- nient for illustrating the relative merits of population reconstruction without and with auxiliary data of varying degrees. Skalski et al. (2012) recommended using point- deletion techniques to determine the stability of population reconstruction to varying amounts of historical information. For a reconstruction to be reliable, the estimated abundance trends should be relatively insensitive to the amount of historical data used in the demographic analysis. They recommend- ed determining how stable the reconstruction abun- dance estimates were when 0, 1, 2, . . . years of the historical data were sequentially eliminated for the analysis. Following the adviceofSkalski etal. (2012), simulated survey data to estimate abundance ( ˆNi) and harvest probability (ˆPi) were added to the original population reconstruction with Cvs ¼ 0.05, 0.125, 0.25 or 0.50. One such survey was assumed to have occurred either at the middle (i.e. 1991) or at the end (2002) of the 24-year population reconstruction (1979-2002). Stability was measured by the relative absolute deviation (RAD) in abundance defined as RAD ¼ 1 y Xy i¼1 jNik - Nij Ni 3 100%; whereNi¼abundanceestimateinyearifromoriginal population reconstruction using all years of data, Nik ¼ abundance estimate in year i from a recon- struction with k historical years of age-at-harvest data deleted and y¼number of years in the truncated reconstruction. The number of years deleted ranged from k¼0, 2, 4, 6, 8, 10 and 12 of the original 24 years of reconstruction. Results Monte Carlo simulations As the precision of the auxiliary studies increased, precision of the reconstructed abundance estimates increased roughly proportional. With minimal in- formation from anauxiliary studywitha CV of 50%, the median CVs of the abundance estimates were intolerably large, usually . 100% (Table 1). As the CV of an auxiliary study used in estimating abun- dance went from 0.25 to 0.05, the median CV of the reconstructed abundance estimates was reduced by morethanhalfandrangedfrom0.354-0.798to0.082- 0.354, respectively. The use of auxiliary data in conjunction with pooled age-class data (i.e. age classes 0.5, 1.5 and 2.5þ) had the same pattern of improvement in precision as occurred for full age- classreconstruction(seeTable1).Theonlydifference was a slight additional reduction in the anticipated CVs. Auxiliary abundance studies had a greater influ- ence on the precision of reconstructed abundance estimates than auxiliary harvest probability studies for equal precision (see Table 1). For example, when 150 Ó WILDLIFE BIOLOGY 19:2 (2013)
- 6. the auxiliary harvest study had a CV of 12.5% (i.e. with 6 25% of the true value 95% of the time), the resulting median CVs for the reconstructed abun- dance estimates ranged from 19.8 to 58.3% (see Table 1). For a similar level of precision in an auxiliary harvest mortality study, the population reconstruction estimates had CVs in the range 57.4- 92.8%. In general, reducing the CV of an auxiliary study produced a commensurate reduction in the CVs of the reconstructed abundance estimates. Timing of the auxiliary studies was generally not important.Thesameimprovementintheprecisionof the population reconstruction estimates occurred regardless of whether the auxiliary study was con- ducted at the beginning, middle or near the end of the investigation. The only exception was when the auxiliarystudywasperformedinthelast(i.e. current) year of reconstruction. Under this circumstance, precision of the population reconstruction willnot be measurably improved until a year thereafter. Incorporating multiple auxiliary studies, 1/3 and 2/3 of the way through the reconstruction period, increased precision of the resulting abundance esti- mates (see Table 1). The effect is most substantial with auxiliary studies of minimal precision (see Table 1). However, a single auxiliary study with a CV of 0.125 results in better precision than two auxiliary studies with a CV of 0.250 each (see Table 1). Black-tailed deer example Augmenting the original black-tailed deer data with an auxiliary abundance survey or a telemetry study, to estimate the vulnerability coefficient, greatly improved the stability of the population reconstruc- tion when the amount of historical data was reduced. Without any auxiliary data, the black-tailed deer reconstruction was very sensitive when four or more years of harvest data were omitted from the analysis (Table 2 and Fig. 1). With six of 24 years of historical data removed, the RAD exceeded 150%. With even more years of data deleted, the reconstruction virtually disintegrated (see Table 2). However, the presence of rough auxiliary studies with a CV¼25% (i.e. 6 50% of the true value 95% of the time) resulted in reasonable stability of the reconstructed population trends. When the auxiliary data estimat- ed the harvest probability (ˆPi) with a CV 0.25, the RAD 7.48% with as many as 12 years of historical data eliminated (see Table 2, Fig. 2). When the auxiliary data provided an abundance estimate within a CV ¼ 0.25, the RAD 22.25% with as many as 12 years of data deleted (see Table 2, Fig. 3). Table 1. Median coefficient of variation (CV) of measurement error in simulated population reconstruction models when including either abundanceorharvestprobabilityauxiliarystudiesatvaryinglevels ofprecisioninthe middleofthereconstruction.Further,weconsideredone (single abundance auxiliary) and two abundance auxiliary studies (double abundance auxiliary). Populations were simulated at high and low levels of harvest and natural survival probabilities. Survival probability Harvest probability Auxiliary CV CV( ˆN) Full model Pooled model Single abundance auxiliary Double abundance auxiliary Single harvest probability auxiliary Single abundance auxiliary 0.75 0.10 0.050 0.082 0.044 0.250 0.119 0.75 0.10 0.125 0.198 0.124 0.574 0.188 0.75 0.10 0.250 0.369 0.280 0.880 0.317 0.75 0.10 0.500 0.953 0.474 2.038 0.895 0.75 0.25 0.050 0.344 0.231 0.573 0.248 0.75 0.25 0.125 0.583 0.467 0.928 0.372 0.75 0.25 0.250 0.798 0.672 1.333 0.530 0.75 0.25 0.500 1.877 0.876 2.441 1.758 0.90 0.10 0.050 0.090 0.068 0.229 0.120 0.90 0.10 0.125 0.221 0.151 0.536 0.191 0.90 0.10 0.250 0.354 0.239 0.795 0.318 0.90 0.10 0.500 0.790 0.404 1.741 0.954 0.90 0.25 0.050 0.325 0.232 0.563 0.218 0.90 0.25 0.125 0.550 0.467 0.883 0.415 0.90 0.25 0.250 0.719 0.642 1.467 0.578 0.90 0.25 0.500 0.894 0.813 2.708 0.861 Ó WILDLIFE BIOLOGY 19:2 (2013) 151
- 7. By the time the auxiliary studies had a precise level suitable for accurate management purposes (CV ¼ 0.125; Robson Regier 1964), the RADs 6.32% with 12 years of data deleted (i.e. 50%) when auxiliary telemetry data were available. Smaller RADs were obtained using radio-telem- etry auxiliary data to estimate harvest probabilities rather than auxiliary abundance surveys for equal levels of precision (see Table 2, Table 3). Stability of the population reconstruction was not affected by whether the auxiliary study was conducted in the middle or near the end of the time series (see Table 3). Discussion Our simulation and sensitivity analyses illustrate several important trade-offsinthe qualityand type of auxiliary studies used in SPR modeling. First, the value of even rather imprecise auxiliary data (e.g. CVs ¼ 25%) on the precision and stability of population reconstruction were evident. This finding demonstrates the general utility of auxiliary infor- mation in SPR and can be used by managers as a guide on the required quality of future auxiliary studies. Second, different types of auxiliary data have different benefits to SPR estimates and the choice of auxiliarydata components ultimately depends on the goal of the resource manager. If the primary concern is precision of the abundance estimates, auxiliary abundance studies are more beneficial than auxiliary radio-telemetry studies for comparable levels of precision. However, managers might be more inter- ested in the inter-annual stability of abundance estimates when designing harvest regulations. If so, auxiliary radio-telemetry studies are more beneficial than abundance auxiliary studies for comparable levels of precision. Another important question in SPR modeling relates to the timing of auxiliary data collection, and our findings suggested that improvements in precision are expected when auxiliary data are collected at any point during the reconstruction Figure 1. Alternative population reconstruc- tions of the black-tailed deer population with 0, 2, 4 or 6 years of historical data removed in the absence of any auxiliary data. Table 2. Relative absolute deviation ðRADÞ in annual abundance estimates from point-deletion sensitivity analyses performed on a statistical populationreconstruction of female black-tailed deer(Skalski etal.2007). Models hadeithernoauxiliarydataor auxiliarydata that estimated abundance or the vulnerability coefficient. Auxiliary studies had a coefficient of variation (CV) of 0.05, 0.125, 0.25 or 0.50 and were simulated in the final year of study. Years removed No auxiliary data Abundance auxiliary data Vulnerability coefficient auxiliary data CV ¼ 0.05 CV ¼ 0.125 CV ¼ 0.25 CV ¼ 0.50 CV ¼ 0.05 CV ¼ 0.125 CV ¼ 0.25 CV ¼ 0.50 2 7.50% 3.18% 3.40% 4.39% 6.12% 0.93% 0.96% 1.49% 3.65% 4 24.94% 3.35% 1.76% 3.07% 9.71% 1.80% 1.55% 1.74% 5.82% 6 152.35% 0.34% 2.77% 10.81% 27.10% 1.25% 0.95% 2.71% 14.21% 8 23943.66% 1.52% 4.33% 22.25% 69.35% 2.89% 1.75% 7.48% 35.81% 10 14386.01% 2.77% 5.52% 15.83% 40.43% 1.79% 2.33% 1.79% 17.91% 12 66998.46% 2.12% 6.32% 21.97% 69.13% 2.16% 1.88% 4.44% 31.29% Mean 13513.50% 2.17% 3.75% 11.70% 32.57% 3.05% 1.48% 1.71% 16.04% 152 Ó WILDLIFE BIOLOGY 19:2 (2013)
- 8. except in the last year. This finding is important for wildlife managers who might have historical age- at-harvest data and are considering conducting a contemporary auxiliary study. These results sug- gest the continued value of collecting auxiliary data after the collection of age-at-harvest data has begun. In other words, age-at-harvest and auxilia- ry data do not need to be collected simultaneously from the start of the study. Further, because auxiliary data can be collected at nearly any point during the reconstruction, a manager can be less concerned if data collection during an auxiliary study proves unsuccessful. For example, if radio- collars malfunction, a manager can attempt a telemetry study at a later date or alter plans and collect another type of auxiliary data without losing the ability to use SPR. Our analyses also indicate that multiple auxiliary studies can further enhance the precision of popula- tion reconstruction estimates. In an adaptive man- Figure 2. Annual abundance trends from a point-deletion sensitivity analysis, with his- toric data removed, on a statistical popula- tion reconstruction of female black-tailed deer, with a simulated auxiliary study to estimate harvest probability in 2002 with a CV of A) 0.05, B) 0.125 and C) 0.250. Ó WILDLIFE BIOLOGY 19:2 (2013) 153
- 9. agement framework where decisions are updated as more and more information becomes available through time, this finding is relevant. SPR is not onlyflexibleenoughtohandlemultipleauxiliarydata types, but we can expect precision of the resulting estimates to improve. However, if a manager is considering whethertocompleteoneor twoauxiliary studies that estimate abundance, it is important to consider that one study with high precision will improve precision of SPR estimates compared with Figure 3. Annual abundance trends from a point-deletion sensitivity analysis, with his- toric data removed, on a statistical popula- tion reconstruction of female black-tailed deer, with a simulated auxiliary study to estimate abundance in 2002 with a CV of A) 0.05, B) 0.125 and C) 0.250.) Table 3. Relative absolute deviation ðRADÞ of abundance estimates from a point-deletion sensitivity analysis of female black-tailed deer comparing auxiliary studies simulated at the end (i.e. 2002) or the center of the reconstruction (i.e. 1991). CV Abundance auxiliary data Vulnerability coefficient auxiliary data 2002 1991 2002 1991 0.050 2.17% 1.76% 1.71% 1.75% 0.125 3.75% 3.76% 1.48% 1.52% 0.250 11.70% 11.77% 3.05% 3.36% 154 Ó WILDLIFE BIOLOGY 19:2 (2013)
- 10. two auxiliary studies with low precision. Although precision of demographic estimates is generally not considered in traditional models of population reconstruction (Millspaugh et al. 2009), we encour- age managers to carefully consider the quality of auxiliary data and howitultimatelyaffectstheresults of SPR. Our assessment, however, is void of any cost- precisioncomparison forradio-telemetryvsauxiliary abundance surveys. The feasibility for each type of survey will vary by species, geographic factors and labor cost. Given that precision of population reconstruction estimates improve roughly propor- tional to the improvement in precision (i.e. reduction in CV) of auxiliary studies, it should be fairly straightforward to perform a cost-benefit analysis. Field investigators should consider their end goals and perform a cost-benefit analysis when planning auxiliary studies to augment SPR models. Such an approach is likely to produce the most useful population estimates and gain support from admin- istrators and stakeholders holding managers ac- countableforcostsandreliabilityofmodelingresults. References Broms, K.M., Skalski, J.R., Millspaugh, J.J., Hagen, C.A. Schulz, J.H. 2010: Using statistical population reconstruc- tion to estimate demographic trends in small game populations. - Journal of Wildlife Management 74: 310- 317. Downing, R.L. 1980: Vital statistics of animal populations. - Wildlife management techniques manual: 247-267. Fieberg, J.R., Shertzer, K.W., Conn, P.B., Noyce, K.V. Garshelis, D.L. 2010: Integrated population modeling of black bears in Minnesota: Implications for monitoring and management. - PLoS ONE 5 (8): e12114. Gast, C.M., Skalski, J.R. Beyer, D.E., Jr. Submitted: Evaluation of fixed- and random-effects models and multistate estimation procedures in statistical population reconstruction. - Journal of Wildlife Management. Gove, N.E., Skalski, J.R., Zager, P. Townsend, R.L. 2002: Statisticalmodelsforpopulationreconstructionusingage- at-harvestdata.-JournalofWildlifeManagement66:310- 320. Millspaugh,J.J., Skalski,J.R.,Townsend, R.L., Diefenbach, D.R., Boyce, M.S., Hansen, L.P. Kammermeyer, K.M. 2009: An evaluation of the sex-age-kill (SAK) model performance. - Journal of Wildlife Management 73: 442- 451. Quinn, T.J., II Deriso, R.B. 1999: Quantitative fish dynamics. - Oxford University Press, Oxford, UK, 542 pp. Robson, D.S. Regier, H.A. 1964: Sample size in Petersen mark-recapture experiments. - Transactions of the Amer- ican Fisheries Society 93: 215-226. Seber, G.A.F. 1982: The estimation of animal abundance and related parameters. - Macmillan, New York, New York, USA, 654 pp. Skalski, J.R., Clawson,M.V. Millspaugh,J.J.2012: Model evaluation in population reconstruction. - Wildlife Biol- ogy 18(3): 225-234. Skalski, J.R., Millspaugh, J.J., Clawson, M.V., Belant, J.L., Etter, D.R., Frawley, B.J. Friedrich, P.D. 2011: Abundance trends of American martens in Michigan basedonstatisticalpopulationreconstruction.- Journalof Wildlife Management 75: 1767-1773. Skalski, J.R., Townsend, R.L. Gilbert, B.A. 2007: Calibrating population reconstruction models using catch-effort and index data. - Journal of Wildlife Man- agement 71: 1309-1316. Ó WILDLIFE BIOLOGY 19:2 (2013) 155