This README.txt was generated on [2022-09-25] by [Daniel R. Eacker] GENERAL INFORMATION 1. Title of Dataset: Data and code for: Summer elk calf survival in a partially migratory population. 2. Author Information Co-investigator 1 Name: Jodi E. Berg (Corresponding Author) ORCID: 0000-0003-0678-8137 Institution: University of Alberta, Edmonton Current affiliation: Alta Science & Engineering, Moscow, ID 83843 Address: 116 St & 85 Ave, Edmonton, AB T6G 2R3, Canada Email: jberg@ualberta.ca Co-investigator 2 Name: Daniel R. Eacker ORCID: 0000-0003-0194-0726 Institution: Taurus Wildlife Consulting Address: 1017 Edwin Pl, Juneau, AK 99801 USA Email: daniel.eacker@tauruswildlifeconsulting.com Co-investigator 3 Name: Mark Hebblewhite ORCID: 0000-0001-5382-1361 Institution: University of Montana, Missoula Address: 32 Campus Drive, Missoula, MT 59812 USA Email: mark.hebblewhite@umontana.edu Co-investigator 4 Name: Evelyn H. Merrill ORCID: 0000-0001-7737-958X Institution: University of Alberta, Edmonton Address: 116 St & 85 Ave, Edmonton, AB T6G 2R3, Canada Email: emerrill@ualberta.ca 3. Date of data collection: 2004-05-14 to 2016-08-31 4. Geographic location of data collection: 51°38’ N, -115°30’ W 5. Funding sources that supported the collection of the data: Long-term Research in Environmental Biology grants (LTREB) to M. Hebblewhite and E. H. Merrill (1556248, 2038704) 6. Recommended citation for this dataset: Berg et al. (2022), Data and R code from: Summer elk calf survival in a partially migratory population, Dryad, Dataset and R code. DATA & FILE OVERVIEW 1. Description of dataset These data and computer code (written in R, https://www.r-project.org) were created to statistically evaluate a suite of intrinsic and extrinsic risk factors related to calf elk and their mothers body condition and age. Specifically, known-fate data were collected from 94 elk calves monitored from 2013-2016 in a partially migratory elk (Cervus canadensis) population in Alberta, Canada. Along with adult female data on pregnancy status, age, and body condition, we created a time-to-event dataset that allowed us to analysis calf mortality risk in a time-to-event approach. We also estimated pooled survivorship and cause-specific mortality, as well as stratifying these metrics by migration tactic (resident vs. eastern migrant). Cox proportional hazards models were used to evaluate calf mortality risk in terms of forage biomass (kg/ha), bear predation risk (from an RSF), and other factors that varied between migration tactics. We tested for differences in a number of maternal reproductive parameters (e.g., pregnancy status) and for calf explanatory variables between migrant and resident elk segments. We also use cumulative incidence functions to estimates cause-specific mortality in this multiple carnivore system. Ultimately, we hope that this work helps wildlife managers anticipate how elk calf survival and partial migration dynamics are affected by grizzly bear predation, and our study builds on a long-term partial migration study at the Ya Ha Tinda Ranch in Alberta, Canada. 2. File List: File 1 Name: CumIncidence.R File 1 Description: Function to estimate cumulative incidence functions File 2 Name: SurvivalData.csv File 2 Description: Elk calf survival records and associated intrinsic variables File 3 Name: RiskCovariates_26day.csv File 3 Description: Calf risk covariate data for 26-day calf-rearing area File 4 Name: AdultFemaleData.csv File 4 Description: Biological data collected from adult female elk during capture File 5 Name: tte_day2.csv File 5 Description: Processed time-to-event data with daily survival timescale File 6 Name: MotherCalfBaseAnalysis.pdf File 6 Description: First section of R analysis code File 7 Name: CalfK-MSurvivorship.pdf File 7 Description: Second section of R analysis code File 8 Name: CalfCauseSpecificMortality.pdf File 8 Description: Third section of R analysis code File 9 Name: CalfRiskRegression.pdf File 9 Description: Fourth section of R analysis code METHODOLOGICAL INFORMATION We used a continuous time-to-event approach on a daily mortality timescale to assess factors influencing mortality risk for elk calves over a 90-day period with birth as the origin. We right-censored calves with failed tags (n = 1) and calves that survived past the 90-day period in their year of birth. We used a generalized Kaplan-Meier (KM) estimator to provide estimates of overall calf survival during the 90-day post-calving period. We stratified estimates by sex, year, and migration tactic, then used log-rank tests to test for significant differences in survival among these groups. We report 95% confidence intervals (CI) for KM survival estimates on the complementary log-log scale, which provide better estimates of uncertainty near the boundaries of 0 or 1. We used cumulative incidence functions to estimate cause-specific pooled mortality rates as well as stratified by migration tactic. We assessed factors influencing mortality risk using the Andersen-Gill formulation of the Cox proportional hazards (PH) model because it accommodated time-dependent risk covariates. We evaluated a suite of nested candidate models that included both intrinsic and extrinsic variables. Further, because of the importance of understanding differences in migratory tactics, we conducted Welch’s two-sample t-tests to evaluate differences in terms of adult female elk age and body condition, and calf sex and birth mass; we used a nonparametric Kolmogorov-Smirnov test to investigate phenological differences in calf birth dates. We also tested for differences between migratory tactics across extrinsic variables using a Welch’s two sample t-test for human disturbance metrics and linear mixed-effects (LME) models with random intercepts to account for correlated observations within calves for forage biomass and predation risk. In summary, grizzly bear predation risk was considered a time-dependent covariate in our mortality risk analysis similar to forage biomass that changed with variation in plant phenology, while other covariates were static through time. We standardized (by subtracting the mean and dividing by 2 SDs) all continuous predictor variables before analysis to allow relative coefficient comparisons among categorical and continuous predictor variables . We also screened variables for collinearity avoided including variables in the same model that had a correlation coefficient of >|0.5|. We further screened out clearly unimportant variables (i.e., P > 0.05) before conducting model selection. As a conservative measure, we used Bayesian Information Criterion (BIC) for model selection. We also used model averaging, given the complexity of the full model relative to our sample size, and reported unconditional standard errors and naturally averaged model coefficients (i.e., non-shrinkage estimates). We present only models in 95% of the cumulative weight. We used Monte Carlo sampling to predict survivorship and the 95% uncertainty interval from model-averaged coefficient estimates and the Breslow estimate of the cumulative baseline hazard function . For brevity, we provide the time-to-event dataset and analysis code rather than include all of the code, GIS, etc. used to estimate calf-rearing areas and extract risk covariates for each individual. Rather, we provide the data and all code to reproduce all results presented in the manuscript. The analysis code is provided as a series of 4 pdf documents that were created in Rstudio using Rmarkdown. The analysis code is broken into the following sections the mirror the results in the manuscript: I) summary of adult female reproductive parameters and other summary statistics and tests, II) survivorship estimates from Kaplan-Meier along with log rank tests, III) calf cause-specific morality estimation, and IV) Cox proportional hazards for calf mortality risk analysis. DATA-SPECIFIC INFORMATION FOR: SurvivalData.csv 1. Number of variables: 17 2. Number of cases/rows: 94 3. Variable List: cowid: adult female individual identifier [character] Age: adult female age (years); 3-18 [integer] collar: type of collar; GPS/VHF [character] calfid: calf individual identifier [character] Capt_Date: calf capture date [character] birthdate: calf birth date [character] censordate: calf right-censoring date [character] Mort_End: mortality date [character] yr: year of study; 2013-2016 [integer] Mort_Status: mortality status; Censor/Mort [character] mortcause: cause of death [character] migstatus: migration tactic; Resident/Eastern [character] sex: calf sex; Male/Female [character] birthwt: calf weight at birth (kg); 12.2-23.0 [numeric] ifbf: percent ingesta-free body fat of adult females; 3.9-17.3 [numeric] julian: Julian birth date of calf adjusted for calendar year; 134-177 [numeric] fdate: Timing of calf birth; After/Before/Peak [character] 4. Missing data codes: None 5. Abbreviations used: N/A; not applicable 6. Other relevant information: Note that we do some processing of this dataset in the code to covert some character classes to factors. DATA-SPECIFIC INFORMATION FOR: RiskCovariates_26day.csv 1. Number of variables: 20 2. Number of cases/rows: 4189 3. Variable List: calfid: calf individual identifier [character] date: calendar date [character] birthdate: calf birth date [character] censordate: calf right-censoring date [character] yr: year of study; 2013-2016 [integer] migstatus: migration tactic; Resident/Eastern [character] forreg: percent cut (logged) forest; 0-38.7 [numeric] elev: elevation (m); 1,438-2,055 [numeric] aspect: aspect (degrees); 90.5-238.7 [numeric] slope: slope (degrees); 2.3-21.8 [numeric] hill: hillshade; 145.8-196.8[numeric] wolf: wolf predation risk; 0-0.4 [numeric] bear: unstandardized bear predation risk; 0.4-4.9 [numeric] lion: cougar predation risk; 0.5-0.8 [numeric] prcntopn: percent open forest; 4.1-95.2 [numeric] prcntcls: percent closed forest; 0.2-93.4[numeric] totfor: forage biomass (g/m2); 25.9-402.58 [numeric] burn: percent burned; 0-89.7 [numeric] line: density of linear features (m/km2); 366.1-5,446.2 [numeric] ndvi: normalized difference vegetation index; 0.5-0.9 [numeric] 4. Missing data codes: None 5. Abbreviations used: N/A; not applicable 6. Other relevant information: Note that this covariate data is combined in the analysis code with SurvivalData.csv to create the daily time-to-event dataset. DATA-SPECIFIC INFORMATION FOR: AdultFemaleData.csv 1. Number of variables: 6 2. Number of cases/rows: 154 3. Variable List: cowid: adult female individual identifier [character] yr: year of study; 2013-2016 [integer] pregnant: indicator of pregnancy status; TRUE/FALSE [logical] ifbf: % ingesta-free body fat of adult females; 1.5-17.3 [numeric] migstatus: migration tactic; Resident/Eastern [character] Age (years): adult female age; 2-21 [integer] 4. Missing data codes: None 5. Abbreviations used: N/A; not applicable 6. Other relevant information: This data includes adult female elk that are not included in the SurvivalData.csv file because some of the females did not produce offspring during the study. DATA-SPECIFIC INFORMATION FOR: tte_day2.csv 1. Number of variables: 38 2. Number of cases/rows: 4,189 3. Variable List: calfid: calf individual identifier [character] day: number of days after calf was born; 1-91 [integer] cowid: adult female individual identifier [character] Age: adult female age (years); 3-18 [integer] collar: type of collar; GPS/VHF [character] Capt_Date: calf capture date [character] birthdate: calf birth date [character] censordate: calf right-censoring date [character] Mort_End: mortality date [character] yr: year of study; 2013-2016 [integer] Mort_Status: mortality status; Censor/Mort [character] mortcause: cause of death [character] migstatus: migration tactic; Resident/Eastern [character] sex: calf sex; Male/Female [character] birthwt: calf weight at birth (kg); 12.2-23.0 [numeric] ifbf: percent ingesta-free body fat of adult females; 3.9-17.3 [numeric] julian: Julian birth date of calf adjusted for calendar year; 134-177 [numeric] fdate: Timing of calf birth; After/Before/Peak [character] enter: individual entry time; 0-90 [integer] exit: individual exit time; 1-91 [integer] event: fate status; 0/1 [integer] date: calendar date [character] forreg: percent cut (logged) forest; 0-38.7 [numeric] elev: elevation (m); 1,438-2,055 [numeric] aspect: aspect (degrees); 90.5-238.7 [numeric] slope: slope (degrees); 2.3-21.8 [numeric] hill: hillshade; 145.8-196.8[numeric] wolf: wolf predation risk; 0-0.4 [numeric] bear: unstandardized bear predation risk; 0.4-4.9 [numeric] lion: cougar predation risk; 0.5-0.8 [numeric] prcntopn: percent open forest; 4.1-95.2 [numeric] prcntcls: percent closed forest; 0.2-93.4[numeric] totfor: forage biomass (g/m2); 25.9-402.58 [numeric] burn: percent burned; 0-89.7 [numeric] line: density of linear features (m/km2); 366.1-5,446.2 [numeric] ndvi: normalized difference vegetation index; 0.5-0.9 [numeric] bear_rsc: standardized bear predation risk; 0-0.47 [numeric] totfor_kgha: forage biomass (kg/ha); 258.8-4,025.8 [numeric] 4. Missing data codes: None 5. Abbreviations used: N/A; not applicable 6. Other relevant information: Note that these variables are not centered and scaled by two standard deviations as this is done in the analysis code.