This repository contains all code needed to reproduce results from the forthcoming paper, "Predicting the invasiveness of underdominance gene drives".
- 1D:
- oneD-diffusion-model.slim is for underdominance alleles and homing underdominance gene drives.
- oneD-diffusion-model-tade.slim is for TADE modification and suppression drives.
- 2D:
- twoD-diffusion-model.slim is for underdominance alleles and homing underdominance gene drives.
- twoD-diffusion-model-tade.slim is for TADE modification and suppression drives.
See the SLiM model README for details on running these models.
-
The critical release width (in 1D) or diameter (in 2D) based on deSolve simulations can be found numerically using the
find_critical_rde_release_statfunction in prediction-functions.R. It relies on code from: -
The phat threshold for TADE modification or suppression can also be found numerically using the
find_phat_tadefunction in prediction-functions.R
See deSolve simulation README for a more detailed description on how find_critical_rde_release_stat works and how we simulated 1D and 2D ODEs in deSolve.
For each threshold-dependent system, we varied the release width (in 1D) or diameter (in 2D) in SLiM and recorded outputs in compressed folders of data
summary-statistics.Rmd shows how summary statistics, such as the SLiM critical release size, were calculated from SLiM data. These summary statistics are provided in the following Excel sheets:
- 1D:
- 2D:
- plotting.Rmd shows how each subplot was created
- Subplots were arranged in Adobe Illustrator 2025
For any questions about the code, please email Isabel at ikk23@cornell.edu