Demo: Running PHLASH on Simulated Data

Running mini_phlash on simulated genomic data.

../../_images/fig_demo_phlash.png — PHLASH demo on msprime-simulated data. Panel A: Observed vs expected SFS under different effective population sizes. Panel B: SFS log-likelihood surface – the maximum identifies the best-fit Ne. Panel C: Random time discretization grid used by PHLASH’s stochastic approach. Panel D: SFS residuals showing the quality of fit.

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