Analyses of RADseq data are possible given these files:
i. ipyrad params: this is the parameter file associated with our processing of RADseq data in ipyrad. 

ii. bellmerge.vcf: After processing reads in ipyrad, the resulting data are stored in the variant call format (VCF) in bellmerge.vcf. These data can be analyzed with 3 R scripts to generate results presented in the paper.

iii. diversity.R: extracts estimates of the population mutation parameter and Tajima's D

iv. inbreeding coefficient.R: calculates observed and expected heterozygosity from genotypic data and the inbreeding coefficient of populations.

v. TDoA.R: calculates directionality statistics and identifies the geographic location of a recent range expansion. To properly run this file, one must set a filepath containing the files in the TDoA dependencies folder.