Data and code from: Temporal replication and sampling efficiency in high-frequency autonomous and manual eDNA sampling
Data files
Jul 06, 2026 version files 1.94 GB
-
Final_Files.zip
1.94 GB
-
README.md
5.18 KB
Abstract
Environmental DNA (eDNA) sampling is a valuable tool for ecological research, but its reliability as a quantitative approach can be improved by accounting for unexplained spatiotemporal variability and autocorrelation in eDNA concentrations. Sampling designs with temporal replication can separate sampling variation from temporal variation, but obtaining true temporal replicates is difficult when samples are collected sequentially. Autosamplers may facilitate high-frequency eDNA sampling, but their sequential sampling structure raises questions about when closely spaced samples can reasonably be treated as temporal replicates. This 19.5-day study evaluated temporal variation in invasive Round Goby (Neogobius melanostomus) eDNA concentration at a single large-river site with short-term hydrologic variability. We used two synchronized autosamplers and a backpack sampler to collect samples as frequently as every 15 min using polyethersulfone (PES) or mixed cellulose ester (MCE) filters. We fit hierarchical temporal models that decomposed sampling and temporal variation, evaluated evidence for residual temporal correlation under alternative priors for the effective range, and used simulation to assess bias when closely spaced sequential samples were treated as simultaneous temporal replicates. We found evidence of residual temporal correlation at short lag scales, but the effective range of temporal correlation was weakly identified and sensitive to prior assumptions. Hydrologic covariates had limited explanatory support once temporal structure was accounted for. Simulations showed that treating sequential samples collected over 15–60 min windows as simultaneous temporal replicates can bias estimates of sampling and temporal variation, especially when temporal correlation is relatively weak. MCE filters produced higher expected recovered copy numbers than PES filters, while evidence for systematic differences among sampler types was weaker. These results emphasize that eDNA temporal replication should be collected as close together in time as possible, and that uncertainty in temporal correlation should be explicitly considered when designing and analyzing high-frequency eDNA surveys.
Dataset DOI: 10.5061/dryad.wstqjq2zz
Description of the data and file structure
This repository contains the data and analysis files that accompany the paper with the same name as this repository where the methods and results are fully described. Briefly, this data consists of eDNA copy number observations of round goby collected at a single site in NY on the Hudson River. eDNA samples were collected with 2 autosamplers and a backpack sampler.
Files and variables
File: Final_Files.zip
Empirical folder
Empirical folder, Data Processing Files subfolder:
Data Processing Files: Files to process the raw field data for analysis
1. 24-008_attributes.csv: Contains sample-level information used for analysis: dates, times, sampler types, filter types, sample volumes
2. 24-008_techrep new.csv: Contains replicate-level qPCR copy number data used for analysis
3. Metadata.csv: Contains metadata for 24_008_attributes.csv and 24-008_techrep new.csv
4. Process Data.R: This file processes the two csv files above to build the dataset for analysis. We also obtain the streamgage covariates here.
5. gage.covs.RData: Contains streamgage measurements for discharge, turbidity, specific conductance, and temperature that are obtained in the Process Data.R script
6. conductivity.data.RData: Conductivity data from the streamgage that is used in the simulation study. This is built in Process Data.R
7. Robot_data.RData: The processed data we use to fit models
Empirical folder, Model Fitting Files subfolder:
1. Robot_data.RData: The processed data we use to fit models, same as 6 in the Data Processing Folder
There are four models: noExp (no temporal correlation), spExp_d0 (GP model, uniform prior on d0), spExp_logRho (GP model, uniform prior on log(rho)), spExp_phi (GP model, unform prior on phi).
2. Fit X RJ multi.R: file to fit model X, with X being one of the four models above
3. NimModel Temporal X RJ.R: the model file for model X
4. output_X_RJ.RData: The posteriors for model X
Empirical folder, Posterior Processing Files subfolder:
1. Robot_data.RData: The processed data we use to fit models, same as 6 in the Data Processing Folder
2. output_X_RJ.RData: The posteriors for model X, same as in the model fitting files folder
3. Process Posterior Functions noExp.R: Function to get pivotal discrepancy measures and their p-values, get point and interval estimates for latent variables
4. Process Posterior Functions spExp.R: Same as 3 but for the spExp models
5. Process Posteriors X RJ.R: For model X, call the appropriate Process Posterior Functions file above (noExp or spExp), summarize output, get parameter estimates (excluding latent variables), make plots and tables
6. table1_X.RData: RJMCMC parameter estimates from model X
7. table2_X.RData: parameter estimates from model X excluding RJMCMC parameters
8. auto_probs_X.RData: autocorrelation pivotal discrepancy measure p-values for model X
9. moment_pvals_X.RData: central moment pivotal discrepancy measure p-values for model X
10. dispIndex.X.RData: dispersion index pivotal discrepancy measure p-value for model X
11. diff_table_X.RData: Efficiency comparison table from model X
12. PPCtable_X.RData: Table of posterior predictive check summaries for model X
13. PPcalc.R: Function to run posterior predictive check
14. PPcalc2.R: Function to run replicate-level mixed predictive check
15. PPcalc3.R: Function to run sample-level mixed predictive check
16. PPcalc4.R: Function to run time-level mixed predictive check
17. PPchecks_X.R: Runs posterior and mixed predictive checks
18. Combine Tables.R: File to combine the tables across all four models
19. prior post plots: Plots the priors and posteriors for the temporal models
Simulation folder
1. simSpExp.R is the data simulator
2. conductivity.data.R is used to simulate data in the data simulator
3. Simulate Datasets X.R where X = S0a, S0b, S1a, S1b, S2a, S2b, S3a, S3b are the data simulation scenarios with different parameter values. The spacing of within burst spacing between samples for S0, S1, S2, and S3 are 0, 0.25, 0.5, and 1 hour, respectively. Then, "a" scenarios have an effective range of 3.56 hours and "b" scenarios have an effective range of 48.38 hours.
4. Fit_Scenarios.R runs the files to fit models for each scenario
5. Fit_X_single where X = S0a, S0b, S1a, S1b, S2a, S2b, S3a, S3b fits models for each scenario
6. NimModel Temporal spExp tempCov.R is the nimble model file used for all scenarios
7. Process Fits.R processes the output from fitted models
8. Folders X_datasets contain the simulated data set for scenario X
9. Folders X_fits contain the posteriors for each simulated data set for scenario X
Code/software
Data can be viewed with any program that can load a .csv. The analysis files are for Program R.
Access information
Other publicly accessible locations of the data:
- NA
Data was derived from the following sources:
- NA
