Legacy effects of four decades of insecticide applications on contemporary benthic macroinvertebrates in eastern Canadian rivers

These files support the manuscript "Legacy effects of four decades of insecticide applications on contemporary benthic macroinvertebrates in eastern Canadian rivers," which has been submitted to Proceedings of the Royal Society B for consideration for publication. This manuscript evaluates the relative importance of historical insecticide applications and modern environmental variables in explaining benthic macroinvertebrate assemblages at 274 sites across the province of New Brunswick, Canada. No new data was generated for the work described in this manuscript. Historical insecticide application data for the province of New Brunswick were taken from Heartz et al. 2023 (https://doi.org/10.1002/ecy.4068), and benthic macroinvertebrate data were downloaded from the publicly available CABIN database (https://open.canada.ca/data, record ID 13564ca4-e330-40a5-9521-bfb1be767147). The authors cannot make their compiled data set publicly available because the subset of CABIN data used for this analysis includes privately maintained data that are only available upon request from respective stakeholders. Please contact the corresponding author for guidance on obtaining these data directly from CABIN. The attached files include analysis_R_scripts: All R codes required to reproduce our analysis. 01_cabin_prep_BMI_data cleans the CABIN benthic macroinvertebrate data for analysis. 02_cabin_prep_env_data cleans the CABIN site-level environmental data for analysis. 03_pesticide_etl_calculation calculates the estimated toxicity load (ETL) in each reach and catchment we studied. 04_ pesticide_apps_per_year calculates the number of years each insecticide was applied in each reach and catchment we studied. 05_site_year_selection_script subsets the CABIN data to include only one observation per site and calculates univariate indices (e.g., %EPT) for each site. 06_data_merge_for_analysis compiles all data (insecticide applications and environmental variables at the reach and catchment scales + CABIN macroinvertebrate and site-level environmental data) 07_univariate_analysis performs multiple regression analysis with the merged data. 08_mulivariate_analysis_familydata performs PCA and RDA analysis with the merged data. 09_misc_figures renders figures used in the manuscript text. 10_site_year_selection_100x loops over script #5 100x to confirm that our results were robust to data subsetting and imputation. 11_data_merge_100x loops over script #6 100x to confirm that our results