Demo: Running discoal on Simulated Data

Running mini_discoal on simulated genomic data.

discoal demo of selective sweep simulation. Panel A: Deterministic sweep trajectories for strong and weak selection. Panel B: Diversity reduction profile around a hard sweep – simulation (green) vs the Maynard Smith & Haigh approximation r/(r+s) (dashed red), showing a clear valley of reduced diversity centred on the selected site. Panel C: TMRCA distribution under a sweep (linked site) vs neutral coalescence – sweeps dramatically reduce coalescence times. Panel D: Sweep duration vs selection coefficient, matching the theoretical prediction.

This figure was generated by running the mini_discoal implementation on data simulated with msprime. The script is at figures/fig_demo_discoal.py.