Code from: Source–sink dynamics explains the coexistence of the invasive pest Dryocosmus kuriphilus and its biological control agent Torymus sinensis across French Eastern Pyrenees
Data files
Dec 15, 2025 version files 14.48 KB
-
README.md
1.31 KB
-
Two_sites.R
8.26 KB
-
Within_site.R
4.91 KB
Abstract
This dataset includes the codes we developed to model the spread of a gall-forming hymenopteran parasite of chestnut trees, Dryocosmus kuriphilus, and its control agent, Torymus sinensis, across 23 natural forest sites located in the French Eastern Pyrenees. The model of D. kuriphilus - T. sinensis interactions was derived from the founding host-parasitoid model proposed by Nicholson-Bailey and adapted to account for the effects of the bottom-up and top-down factors regulating the invasion of D. kuriphilus. The first code then corresponds to the modelling of D. kuriphilus and T. sinensis population dynamic within one site, while the second code is an extension of the within-site modelling into a two-site model, allowing for the study of the effects of the pest and control agent dispersal between two sites on their population dynamics.
Dataset DOI: 10.5061/dryad.p8cz8wb3r
Description of the data and file structure
We developed an extension of the founding host-parasitoid mathematical model proposed by Nicholson-Bailey to account for the effects of the bottom-up and top-down factors regulating the invasion of D. kuriphilus, and its control by T. sinensis. To identify these effects and quantify their impact on the D. kuriphilus -T. sinensis. population dynamics , we performed a 2 years field study and combined ecological, genomic, statistical and epidemiological approaches. We then extended our model to investigate the effect of the pest and the control agent dispersal between two sites on their dynamics interaction.
Files and variables
File: Two_sites.R
Description: Code for the Two-sites dynamical model
File: Within_site.R
Description: Code for the within site dynamical model
Code/software
Codes of the within site and two-sites model can be read using R version 4.5.1 with the "tidyverse" package. Both codes are annotated to make them easier to understand.