Description of the scripts used for simulations used for ABC model selection and to perform Posterior Predictive Check (PPC) analyses. Simulate data for ABC: runms_Paur.py - A python script to draw values for the parameters tau, theta, ancient migration, current migration, bottleneck and population expansion. This file contains the 7 models tested in this study and outputs the ms' simulation results for both cpDNA and SSR markers in the species P. aurisetus. runms_Pmac.py - A python script to draw values for the parameters tau, theta, ancient migration, current migration, bottleneck and population expansion. This file contains the 7 models tested in this study and outputs the ms' simulation results for both cpDNA and SSR markers in the species P. machrisii. microsat2arp.py - A python script that converts the microsatellite data from ms' accessory program microsat to the arlequin format. simulate data for PPC: runPPC_Paur.py - Retrieve parameter values from the parameter estimates of the selected in the P. aurisetus dataset and use them to simulate new data to perform PPC analysis. runPPC_Pmac.py - Retrieve parameter values from the parameter estimates of the selected in the P. machrisii dataset and use them to simulate new data to perform PPC analysis.