Re-Evaluating Neonatal-Age Models for Ungulates: Does Model Choice Affect Survival Estimates?

New-hoof growth is regarded as the most reliable metric for predicting age of newborn ungulates, but variation in estimated age among hoof-growth equations that have been developed may affect estimates of survival in staggered-entry models. We used known-age newborns to evaluate variation in age estimates among existing hoof-growth equations and to determine the consequences of that variation on survival estimates. During 2001–2009, we captured and radiocollared 174 newborn (≤24-hrs old) ungulates: 76 white-tailed deer (Odocoileus virginianus) in Minnesota and South Dakota, 61 mule deer (O. hemionus) in California, and 37 pronghorn (Antilocapra americana) in South Dakota. Estimated age of known-age newborns differed among hoof-growth models and varied by >15 days for white-tailed deer, >20 days for mule deer, and >10 days for pronghorn. Accuracy (i.e., the proportion of neonates assigned to the correct age) in aging newborns using published equations ranged from 0.0% to 39.4% in white-tailed deer, 0.0% to 3.3% in mule deer, and was 0.0% for pronghorns. Results of survival modeling indicated that variability in estimates of age-at-capture affected short-term estimates of survival (i.e., 30 days) for white-tailed deer and mule deer, and survival estimates over a longer time frame (i.e., 120 days) for mule deer. Conversely, survival estimates for pronghorn were not affected by estimates of age. Our analyses indicate that modeling survival in daily intervals is too fine a temporal scale when age-at-capture is unknown given the potential inaccuracies among equations used to estimate age of neonates. Instead, weekly survival intervals are more appropriate because most models accurately predicted ages within 1 week of the known age. Variation among results of neonatal-age models on short- and long-term estimates of survival for known-age young emphasizes the importance of selecting an appropriate hoof-growth equation and appropriately defining intervals (i.e., weekly versus daily) for estimating survival.


Introduction
Survival of young ungulates often drives annual fluctuations in population growth because adult survival is relatively constant in comparison to young [1][2][3].Therefore, determining factors that influence the ecology and mortality of young ungulates is important for population management and understanding how pre-hunt survival affects sustainable harvest rates [4,5].Yet, collecting data on neonates can be challenging because of their cryptic coloration and inactivity during the first month of life; thus, making capture difficult and survival information costly to collect [4].
Regional and seasonal variation in survival rates and causespecific mortality of young ungulates with respect to sex, age, animal density, and environmental conditions further complicates the interpretation and limits the application of results of such studies [6][7][8][9][10].Additionally, survival and cause-specific mortality of young vary as a function of age, with marked changes occurring within the first few weeks of life [5,[11][12][13][14].This strongly suggests accurate information on date of birth is critical to understanding age-dependent patterns of survival and risk to specific sources of mortality.Inaccuracies in age estimates of captured neonates could affect estimates of survival because age-at-capture determines the interval that an individual enters and exits staggered-entry models [15][16][17].
To estimate age-at-capture, researchers have developed regression models to predict age (y-axis) of newborn white-tailed deer (Odocoileus virginianus; [18][19][20]), mule deer (O.hemionus; [21]), and pronghorn (Antilocapra americana; [22]) from measurements of new-hoof growth (x-axis) of captive animals.The results of such models, however, may not be consistent with estimates of age from hoof-growth measurements collected from free-ranging animals [17].Although slopes among published equations describing the relationship between age (days) as a function of the independent variable new-hoof growth (mm) often are similar, intercept terms differ, thereby producing different age estimates for the same neonate.For example, intercepts of published growth equations ranged from 28.29 to 0.66 for white-tailed neonates [17][18][19][20] and from 26.30 to 5.29 for mule deer neonates [17,21].
Our first objective was to evaluate existing models for estimating age of neonates from measurements of hoof growth of known-age (#24-hours old), wild newborns of 3 species (white-tailed deer, mule deer, and pronghorn) from study sites in California, South Dakota, and Minnesota, USA.We expected age estimates of newborns to vary among hoof-growth equations because of the variation in slopes and intercepts among models.Our second objective was to quantify the consequences of estimating age-atcapture of neonate ungulates on the resulting estimates of survival at both 30 and 120 days.We predicted that survival estimates would differ among models.Specifically, differences in estimation of age (days) at the time of capture would influence the time interval at which a neonate entered and exited the population atrisk in survival analyses.

Study Area
Our study was conducted in Minnesota (white-tailed deer), South Dakota (white-tailed deer and pronghorn), and California (mule deer), USA (Fig. 1).We captured pronghorn in western South Dakota (Fall River, Harding, and Custer counties).Wind Cave National Park (WCNP), our study area in Custer County, was located in the southern Black Hills and was enclosed by a 2.5m woven-wire fence with cattle guards to prevent passage by ungulates [23].Our study areas in Fall River and Harding counties were characterized by a mosaic of mixed-grass prairie interspersed with shrubs (e.g., sagebrush [Artemisia spp.]) and patches of forest.Topography was rolling, with occasional buttes and intermittent streams [24,25].
We captured white-tailed deer fawns in study areas in northcentral and eastern South Dakota, and in southwestern Minnesota, USA.North-central South Dakota (Edmunds and Faulk counties) was characterized by previously glaciated, rolling prairie interspersed with pothole wetland complexes, cultivated agricultural land, intermittent streams, and river floodplains [26].The eastern South Dakota study area (Brookings County) lies within the Prairie Coteau Region formed by the Wisconsin Glaciation [27].The Prairie Coteau historically contained numerous wetlands and in eastern South Dakota, approximately 35% were drained through anthropogenic modifications (e.g., agriculture; [25,28]).The southwest Minnesota study area was characterized by flat-torolling topography [29].Lincoln and Pipestone counties, Minnesota, are within the Prairie Coteau physiographic region, whereas Redwood and Renville counties occur within the Minnesota River Valley physiographic region.The Minnesota River Valley was a linear corridor that was heavily forested with small interspersed grassland remnants and adjacent lands comprised primarily of cultivated crops [13,30].Habitat in the region was fragmented and dominated by intense row-crop agriculture [31].
In California, our study area was located in the Sierra Nevada (a high-elevation mountain range) and focused on a migratory mule deer population that wintered on the western edge of the Great Basin that included Fresno, Inyo, Madera, and Mono counties.Summer range for mule deer in the Sierra Nevada occurred on both sides of the Sierra crest at elevations ranging from 2,200 to .3,000m [9,32].The western slope of the summer range was substantially more mesic and dominated by the upper montane and mixed conifer vegetation zones, whereas the xeric, eastern slope of the Sierra Nevada (,2,100 m) was dominated by the sagebrush vegetation zone [9,32].

Adult Capture
We captured adult female white-tailed deer and mule deer on winter range using a hand-held net gun fired from a helicopter [33].In Minnesota, female white-tailed deer were immobilized with ketamine hydrochloride (5 mg/kg IM; Ketaset; Fort Dodge Laboratories, Fort Dodge, IA, USA) and xylazine hydrochloride (1 mg/kg IM; Xylaject, Phoenix Pharmaceutical Inc., St. Joseph, MO, USA).We used yohimbine hydrochloride (0.125 mg/kg IV; Yobine, Ben Venue Laboratories, Inc., Bedford, OH, USA) as an antagonist [34].We did not use chemical immobilization during adult mule deer capture in California [32].White-tailed deer and mule deer females were blindfolded and hobbled prior to transport via helicopter to a central processing station.
We assumed all captured adult white-tailed deer in Minnesota were pregnant, given the high pregnancy rates for adult deer inhabiting these agricultural regions [35,36].We assumed that not all captured mule deer at the California site were pregnant, and determined pregnancy status and litter size of mule deer with trans-abdominal ultrasonography [37].We captured 14 adult female white-tailed deer in Minnesota and 263 adult female mule deer in California.All white-tailed deer and 178 mule deer were fitted with vaginal implant transmitters (VIT; Model M3930, Advanced Telemetry Systems, Isanti, MN, USA) using standard procedures [30,38,39].

Neonate Capture
To spatially focus our search efforts for newborns, we used VITs, pre-and post-partum behavior of radiocollared dams and post-partum behavior of unmarked dams to locate birth sites.We captured neonatal white-tailed deer during late May and early June 2001-2004 in eastern South Dakota and southwestern Minnesota [13,30,31], and 2007-2009 in north-central South Dakota [40][41][42].We captured pronghorn neonates during late May and early June 2002-2005 in western South Dakota [43], and neonatal mule deer in California during mid-June through mid-July 2006-2008 [14,44].
We monitored VITs 1-3 times daily using a very high frequency (VHF) receiver and a truck-mounted, null-peak antenna system with an electronic digital compass (C100 Compass Engine; KVH Industries, Inc., Middletown, RI, USA; [45]), Cessna 182 and 185 aircraft (Cessna Aircraft Co., Wichita, KS, USA) fitted with 2, 2element Yagi antennas or a hand-held Yagi antenna (Advanced Telemetry Systems).When a VIT indicated a birth had occurred, we immediately (,4 hours in Minnesota, ,12 hours in California) located the unit and radiocollared female.We used the location of the VIT and radiocollared female to focus a search area, and implemented intensive ground searches to locate neonates.Because we checked all VIT signals a minimum of once per day during the fawning period, time of birth was known for newborns captured #12 hours after receiving a signal the VIT had been expelled.Because VITs may be expelled pre-partum [39,46], we used other indicators of age for newborns captured .12hours after the VIT was shed and for newborns captured from dams without VITs.We classified neonates as #24-hrs old if they met all of the following criteria: 1) tips of hooves were soft and supple, often white, and a soft semi-gelatinous sulphur-yellow pad was still present on the bottoms of hooves and tips of dewclaws [18], 2) umbilicus was still wet and fresh in appearance [18], and 3) alarm bradycardia (includes a motionless behavioral response and decreased heart rate) and a high level of passiveness were displayed [20,47].We removed all neonates .24-hrsold from our analyses because we were interested in evaluating the performance of new-hoof growth equations based on known-age newborns.
We did not use VITs for white-tailed deer in eastern and northcentral South Dakota, or for pronghorn in western South Dakota.Instead, we relied on post-partum behavior of females [13,35,42,[48][49][50] or direct observations of neonates [51] to facilitate capture.Nevertheless, we used the same 3 criteria used for females with VITs when ascertaining if neonates were #24hrs-old.
We manually restrained neonates of all 3 species and carefully measured new-hoof growth (the distance from the hairline to the growth-ring line on the outer edge of the hoof [20]) to the nearest 0.1 mm using a dial caliper, and used average new-hoof growth of the 2 front hooves in subsequent analyses.The same researcher measured or trained technicians to measure new-hoof growth at each site throughout the study period.We fitted each neonate with an expandable radiocollar (model M4210, Advanced Telemetry Solutions, Isanti, MN, USA; model CB-6, Telonics Inc., Mesa, AZ, USA; model TS-37, Telemetry Solutions, Concord, CA, USA) to monitor survival.We relocated neonates daily for the first 30 days of life and then $1 time per week through at least 120 days old using a truck mounted null-peak antenna, Cessna 182 and 185 aircraft, or a hand-held Yagi antenna.When we detected a mortality signal, ground personnel located the collar within 12 hours (range 0.5-12 hours).We examined physical evidence at the collar-recovery site to ascertain whether the device was shed prematurely or the neonate died.We used evidence of struggle, blood (on ground or collar), or remains as indications of cause of neonate mortality.Collars attached to fencing, trees, or shrubs, and collars with all folds expanded with no evidence of mortality in the immediate area (50-m radius) were considered prematurely shed [44].We estimated date of collar loss as the mid-point between date last known alive and date of mortality signal.
Animal handling methods used in this project followed guidelines recommended by the American Society of Mammalogists [52] and were approved by the Institutional Animal Care and Use Committee at South Dakota State University (Approval nos.04-A009, 00-A038, 02-A037, 02-A043, 02-A001, 02-A002) and an Independent Institutional Animal Care and Use Committee at Idaho State University (Protocol number 650-0410).Data collection was authorized by California Department of Fish and Wildlife, South Dakota Department of Game, Fish and Parks, and Minnesota Department of Natural Resources.All data were collected on private land with permission of individual landowners as well as public land with permission of California Department of Fish and Wildlife, South Dakota Department of Game, Fish and Parks, Minnesota Department of Natural Resources, and the U.S. Fish and Wildlife Service.Field studies did not involve endangered or protected species.

Statistical analysis
For white-tailed deer, we evaluated hoof-growth equations of Haugen and Speake (hereafter Haugen and Speake; [18]), Sams [17]).For pronghorn, we evaluated the Tucker and Garner (hereafter Tucker and Garner; [22]) equation (Table 1).We evaluated hoof-growth equations by: 1) comparing and contrasting age of known-age (#24-hours old) newborns with age estimates using each hoof-growth equation based on new-hoof growth measurements, and 2) comparing age-specific estimates of 30-day and 120day survival from known ages against age estimated from each hoof-growth equation.We analyzed each species separately and used analysis of variance (ANOVA) and Tukey-Kramer pairwise comparisons to determine differences in estimates of new-hoof growth among study sites and years, and controlled for study site by blocking to account for possible measurement bias by investigators.We conducted statistical tests using SAS version 9.3 [53], with an experiment-wide error rate of 0.05.
We evaluated the effect of age estimates from the various hoofgrowth equations proposed for each species on survival estimates using known-fate models in Program MARK [54] with the logitlink function to model survival of neonates to 30 and 120 days-ofage.The focus of known-fate models is the probability of surviving an interval between sampling occasions.The sampling or detection probabilities are assumed to be 1 because the animal is radio-tagged [55].The model is a product of simple binomial likelihoods and the precision of known-fate models is typically high [55].We selected 30-day survival because research on neonates ,1 month of age provides critical information regarding reproduction, sex ratios, mortality, movements, and behavior [48].We selected 120-day survival to approximate summer survival of young ungulates.For the known-aged model, all neonates entered the analysis during the first interval because only neonates known to be ,24 hrs old were included in our analysis.In contrast, when estimating age using hoof-growth equations, newborns that were estimated to be 0-24 hrs old at capture entered the survival analysis during the first interval, and newborns estimated to be 24-48 hrs old at capture entered the analysis in the second interval, and so on.When using hoof-growth equations with negative intercepts [19][20][21], we truncated estimated age of neonates to zero when estimates were ,0-days old.We right-censored all animals that prematurely shed collars because censoring likely was independent of the fate of the animal [44].During daily monitoring, we right-censored animals on the day the collar was retrieved.During less intense monitoring (i.e., $1 time per week), we right-censored animals at the mid-point between last relocation and day the collar was retrieved.
We evaluated survival models that assumed daily survival was constant throughout the 30-day or 120-day period, and models that allowed daily survival to vary randomly among days (daily survival) or among weeks (weekly survival; [54,55]).For each newborn, we created an encounter history using the neonate's actual age as well as an encounter history for the neonate's age estimated using each hoof-growth equation.We evaluated potential differences in survival as a function of estimated age by comparing models with different group assignments.For example, for mule deer, we contrasted models that considered survival varied among groups (known-age [KA], Robinette [R], and Haskell [H] models; KA, R, H), with a model that considered survival of known-age neonates different from those estimated by hoof-growth equations (assumed survival differed between the 2 groups with survival equal between Robinette and Haskell equations; KA, R = H), with a model that considered survival was constant (assumes survival was similar between the groups; KA = R = H), and so forth.For mule deer and pronghorn, we modeled all combinations of group assignments for known age newborns with ages estimated from hoof-growth equations.Because there were 4 hoof-growth equations for white-tailed deer, we selected model combinations based on similarities in mean age estimates calculated from the models.For instance, model KA,B,S,HS,H tested the hypothesis that survival differed among the 5 groups (known-age, Brinkman [B], Sams [S], Haugen and Speake [HS], and Haskell [H]), whereas KA = B = S, HS = H assumed survival to be the same for known-age neonates and those aged using the Brinkman and Sams equations but differed from those aged with the Haugen and Speake and Haskell equations (which were similar).We used Akaike's Information Criterion corrected for small sample size (AIC c ) to select models that best described the data.We used model weights (w i ) as an indication of relative support for each model, and considered models with DAIC c values ,2 from that of the top-ranked models to be viable alternatives [56].We used Program CONTRAST to test for differences in survival among cohorts in our top-ranked models [57].

Results
During May-July 2001-2009, we captured 407 neonates.Of these, 76 of 204 white-tailed deer, 61 of 119 mule deer, and 37 of 142 pronghorn were estimated at #24-hrs old and thus, were used in subsequent analyses.We captured 18 white-tailed deer and all 61 mule deer with the aid of VITs.All others were captured via postpartum behavior of females or direct observation of the birth of neonates.New-hoof growth was similar (t 34 = 21.44,P = 0.26) between white-tailed deer neonates located using VITs and those located using more traditional methods.Estimated age of newborns generated from existing hoof-growth equations varied by .15days for white-tailed deer, .20 days for mule deer, and .10 days for pronghorn (Table 2).Based on known-age of newborns, accuracy (i.e., the proportion of neonates assigned to the correct age) was greatest for the Brinkman and Sams equations for aging neonatal white-tailed deer (Table 3).For mule deer, the Robinette model was more accurate than the Haskell model (Table 3).Accuracy of the Tucker and Garner model was 0.0% for known-age newborn pronghorns (   Models are those of Brinkman [20], Haskell [17], Sams [19], Haugen and Speake [18], Robinette [21], or Tucker and Garner [22].doi: The top-ranked 30-day survival model (KA = B = S, H = HS_weekly survival; w i = 0.96) for white-tailed deer indicated there were differences in survival estimates among hoof-growth equations (Table 4).The model indicated survival for known-age neonates and those with age estimated using Brinkman or Sams equations were similar (0.80, SE = 0.03, 95% CI = 0.74-0.85),but differed marginally (x 2 1 = 2.56, P = 0.11) from those with age estimated using Haskell and Haugen and Speake (0.72, SE = 0.04, 95% CI = 0.65-0.79)equations (Fig. 2).All remaining 30-day survival models for white-tailed deer were .2DAIC c from the topranked model.Conversely, the top-ranked survival model for 120 days post-birth (KA = B = S = H = HS_weekly survival; w i = 0.99) indicated similar survival among hoof-growth equations and known ages of fawns (Table 5).All remaining white-tailed deer survival models up to 120-days old were .2DAIC c from the selected model.
The top-ranked 30-day (KA = TG; w i = 0.69) and 120-day (KA = TG_weekly survival; w i = 1.00) models for pronghorn a Models were obtained from Brinkman [20], Haskell [17], Sams [19], Haugen and Speake [18], Robinette [21], or Tucker and Garner [22].b Accuracy was defined as proportion of neonates assigned to the correct age and classified within 1, 2, and 3 days of their true age.doi:10.1371/journal.pone.0108797.t003 indicated no difference in survival between known-age newborns and those with ages derived from the Tucker and Garner equation (Table 8, 9).All other pronghorn survival models were .2DAIC c from the top-ranked 30-and 120-day models.

Discussion
Variability in estimates of birth date and neonatal age affected estimates of survival generated from staggered-entry models [16,17].This was especially true when high rates of perinatal mortality occurred because estimates of age-at-capture determine the intervals at which neonates entered and exited survival analyses.Our results indicated that variability in estimates of ageat-capture based on neonatal-age models that were applied to known-age newborns (#24-hours old) affected estimates of 30-day survival for both white-tailed deer and mule deer, and estimates of 120-day survival for the latter species.Model results indicated no variability for short-or long-term survival estimates for pronghorn.The top-ranked model indicated similar survival among known-age (KA) neonates and those aged using the Brinkman (B) and Sams (S) hoof-growth equations, which differed from neonates aged using the Haskell (H) and Haugen and Speake (HS) hoof-growth equations; survival was similar between Haskell and Haugen and Speake equations.Equations were obtained from Brinkman [20], Sams [19], Haskell [17] and Haugen and Speake [18].doi:10.1371/journal.pone.0108797.g002[20], H = Haskell [17], S = Sams [19], or HS = Haugen and Speake [18] neonatal-age model.b Temporal scale represents constant, daily, or weekly survival among intervals.However, the relatively low DAIC c value of the second-ranked 30day survival model suggests a similar influence on this species.Age estimates for newborns varied by .10days for white-tailed deer and .15days for mule deer compared with .6 days for pronghorn.The marked differences in age estimates for neonates among neonatal-age models when combined with high rates of early mortality likely led to the significant differences in 30-day survival estimates for both deer species and 120-day survival estimates for mule deer.Over longer time periods, survival estimates of white-tailed deer and pronghorn were not influenced greatly by methods of age estimation.Some hoof-growth equations used to estimate age of neonate ungulates were imprecise, which suggests the need for caution in modeling survival at too fine a temporal scale, especially if true birth dates are unknown.Most hoof-growth equations predicted ages within 1 week of the known age of newborns.Hence, potential inaccuracies in age estimates of newborns of up to 20 days may contribute to negative bias in survival estimates and decreased fit of survival models based on a daily encounter history.Although improvements in technology and performance of VITs facilitate the capture of known-age  neonates and provide the opportunity to evaluate temporal patterns in ecological phenomena, our results indicate that lack of precision in neonatal-age models should be taken into consideration when birthdate is unknown.
Variability in estimates of survival among neonatal-age models for white-tailed deer and mule deer was related to intercepts of those models (Table 1).Models with negative intercepts (i.e., Brinkman, Sams, and Robinette) and expected positive values of new-hoof growth for newborns at 0 days-of-age yielded estimates of survival up to 30 and 120 days-of-age that were similar to the known-age neonate group.Conversely, survival was lower for neonates aged using models (i.e., Haskell, Haugen and Speake) with positive intercepts.The tendency for neonatal-age models with negative intercepts to predict survival estimates similar to  120-day survival.Consequently, intercepts of equations to estimate age based on new-hoof growth play a key role in accurate estimation of age for neonatal ungulates.Additionally, differences in intercepts and growth rates (slopes) among neonatal-age models and variation in new-hoof growth among study sites and years in our study support the hypothesis that relationships between newhoof growth and age may be population-and time-specific [17].
Variation in intercepts among growth models may be a function of sampling variance among biologists when measuring new-hoof growth.For example, based on data from the Brinkman equation, a 0.75-mm change (error) in measurement of new-hoof growth at birth would result in a 40.4% change in estimate of the regression intercept.This would result in an intercept similar to the Sams equation.Measurement error could be reduced within studies by having a single person conduct all new-hoof growth measurements, or at a minimum have a single person train all personnel.However, measurements across studies likely result in variance caused by observer bias that contribute to differences in intercepts.Also, variation in intercepts may possibly be explained by biological variance such as gestation length and physical condition of the mother [58][59][60].Increasing birth weight is related to longer gestation period [58] while increasing body mass is correlated with greater hoof growth [19].Gestation period can range dramatically for ungulate species [58][59][60] and may be shortened or lengthened depending upon nutritional status of the female and environmental conditions [14,[61][62][63].Additionally, new-hoof growth of whitetailed deer differed between females fed a high-and low-protein diet; neonates born to females on a low-protein diet had shorter hoof-growth measurements than those born to females on a highprotein diet [19].Other ancillary variables such as maternal nutritional condition, birth mass, and litter size also may contribute to differences in new-hoof growth and thus, intercepts in equations for estimating age.
Probability of survival can fluctuate markedly with age during the first weeks of life for neonatal ungulates and generally is thought to occur because of size, agility, activity, and vulnerability of neonates [5,[12][13][14]64].Our results indicate that discrepancies in estimates of age from the true age of wild-captured neonates can alter results of temporal survival patterns and thus, interpretation of factors influencing survival.For example, initial observations suggested the greatest period of vulnerability for white-tailed deer in Minnesota and South Dakota aged using the Brinkman equation occurred during the first 2 weeks-of-life [13].These results were consistent with Rohm et al. [5] who aged fawns using the Haugen and Speake equation and attributed the greatest period of mortality to changes in habitat availability and coyote (Canis latrans) behavior.Conversely, the same neonates aged using the Haskell equation would not support these conclusions but rather supported Nelson and Woolf [11], who observed that neonate mortality was highest during 2-8 weeks-of-life.They [11] hypothesized that neonates were safe from predation because of their sedentary behavior when 0-2 weeks old, were most vulnerable when neonates became active during 2-8 weeks old, and could evade predators when .8 weeks old.Understanding the behavior of wild, young ungulates, their vulnerability to mortality, and assessing the influence of management actions on probability of survival could be undermined by using inaccurately aged neonates in survival analyses.
The advent of powerful modeling techniques available in Program MARK permits the use of individual covariates such as birth weight to examine their influence on survival [5,13,14,42].Weight at birth is a key factor affecting probability of survival for young ungulates because it is associated with strength and viability of neonates [10,19,[64][65][66].Estimates of birth weight calculated using weight at capture and estimated age to back-calculate that metric [5,41,67], would vary depending on the hoof-growth equation used.However, this should have little effect on model results using estimated birth weight as a covariate, unless overestimation of ages is large enough to yield estimates of birth weight that are truncated at 0 kg (i.e., Haskell for white-tailed and mule deer).Without this truncation, the relationship between survival and estimated birth weight will remain stable because larger neonates will have a greater probability of survival than smaller neonates even if absolute values for estimated birth weight are biased upward or downward for all individuals.Consequently, caution should be used if management objectives include identifying a threshold in birth weight below which fawn survival is compromised, because thresholds will be biased low when ages at capture are overestimated.
Timing of parturition coincides with the flush of nutrients during the spring to support the costs of lactation [10,68], allow sufficient growth and accumulation of body reserves of young before winter [69], potentially avoid high predation pressure [12,70], and enhance survival of young.Although other methods for determining peak parturition are available including evaluation of movement with GPS data or use of VITs [38,39,71], most  studies rely on back-calculated birth dates from estimated age-atcapture to determine parturition dates [10].As with estimated birth weights, relationships between parturition dates and other characteristics will remain relative, but identifying actual dates for peak parturition or thresholds in those mathematical relationships could be biased by inaccurate neonatal-age models.For example, Lomas and Bender [10] observed an 18-29 day shift in mean birth dates of mule deer in north-central New Mexico between the 1980s and the early 2000s in response to a marked decline in habitat quality in the region.Similar apparent shifts could be noted between studies that use different models to estimate ages of neonates even if little change had actually occurred.

Conclusions
Survival of young drives annual population trajectories for ungulates and influences management strategies for harvest, habitat treatments, and predator control [67].Survival of neonates is routinely estimated through capture and collaring during the first few weeks of life, with subsequent monitoring for survival.Neonatal-age models based on new-hoof growth have been regarded as the most accurate method to back-calculate date of birth, but as we demonstrated, choice of model can have a profound effect on age-dependent patterns of mortality and shortterm estimates of survival.Our results indicated that estimates of summer survival were more robust to variation in estimates of ageat-capture and support the reliability of most previously reported estimates of survival that used neonatal-age models.We encourage researchers to use caution, however, when interpreting estimates of survival, birth weights, and parturition dates when age is estimated based on hoof-growth equations because some models perform better than others.In most studies, a portion of wild-captured neonates may be confidently identified as newborn either through observation of birth or robust criteria such as that used in our study.Therefore, we suggest testing for differences in birthdates and estimated birth weights between known-age neonates and those whose ages are estimated to assess the potential for bias associated with those hoof-growth equations [14].Alternatively, researchers could estimate their own growth models [17] if sufficient data were available.Finally, our analyses indicate that modeling survival in daily intervals is too fine a temporal scale when birth date is unknown because of the potential inaccuracies among models available to estimate age of neonates.We suggest that weekly survival intervals are more appropriate because most hoof-growth models accurately predicted neonatal age within one week.

Figure 2 .
Figure 2. Candidate models of 30-day survival of white-tailed deer (Odocoileus virginianus) neonates in South Dakota and Minnesota, USA.Survival estimates up to 30 days from top-ranked model (KA = B = S, H = HS) of white-tailed deer (Odocoileus virginianus) neonates (n = 71), South Dakota and Minnesota, USA, 2001-2009.The top-ranked model indicated similar survival among known-age (KA) neonates and those aged using the Brinkman (B) and Sams (S) hoof-growth equations, which differed from neonates aged using the Haskell (H) and Haugen and Speake (HS) hoof-growth equations; survival was similar between Haskell and Haugen and Speake equations.Equations were obtained from Brinkman[20], Sams[19], Haskell[17] and Haugen and Speake[18].doi:10.1371/journal.pone.0108797.g002

c
Akaike's Information Criterion corrected for small sample size[55].d Difference in the AIC c value of the top-ranked model and that of the model under consideration.e Akaike weight [55].f Number of parameters.doi:10.1371/journal.pone.0108797.t004

Table 1 .
Intercept and slope for published neonatal-age models used to estimate age (days) of neonatal white-tailed deer, mule deer, and pronghorn based on measurements of new hoof growth (mm).
Age estimates were generated from published neonatal-age models.

Table 4 .
Survival models using Program MARK for white-tailed deer from birth to 30 days using known-age newborn fawns and ages estimated from published neonatal-age models based on new-hoof growth, South Dakota and Minnesota,USA, 2001USA,  -2009.   .
a KA = known-age newborn fawns, B = Brinkman

Table 5 .
Survival models using Program MARK for white-tailed deer from birth to 120 days using known-age newborn fawns and ages estimated from published neonatal-age models based on new-hoof growth, South Dakota and Minnesota,USA, 2001USA,  -2009.   .

Table 6 .
Survival models using Program MARK for fawn mule deer from birth to 30 days using known-age newborn fawns and ages estimated from published neonatal-age models based on new-hoof growth, California, USA, 2005-2007.

Table 9 .
Survival models using Program MARK for pronghorn from birth to 120 days using known-age newborn fawns and ages estimated from published neonatal-age models based on new-hoof growth, western South Dakota,USA, 2002USA,  -2005.   .