All the scripts used are available on Dryad. https://datadryad.org/stash/dataset/doi:10.5061/dryad.t76fk80
I used R to do the blockbootstrap strategy (there is a source code to estimate them):
confidence Intervals
The blockbootstrapped SFS to make demographic inferences are computeed using the functions matsplitter and booti from SFS_bootstrap.R
If you want to estimate confidence intervals of the SFS you wiil have to use the functions matsplitter as well as booti.2pop or booti.allpop
And the demographic inferences files are these ones :
Supplementary file - settings files for demographic inferences
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this
article to respond.