we used the resulting theta (also inferred by moments) and a mutation rate106 of 1.66 × 10−8 to generate corresponding effective population sizes (Ne). To convert generations to years, we assumed a generation time of 25 years. The 95% confidence intervals were generated by resampling the site frequency spectrum 1,000 times and using the Godambe information matrix to generate parameter uncertainties107. A more detailed description is available in Supplementary Information 1.10.
