The 22 autosomes and X chromosome were split into 833,754 and 18,039 SNP intervals respectively, as defined by the boundaries of all recombination events found in the combined data set. We estimated female and male recombination rates within each of these SNPs intervals. First, we build female and male unrefined maps based on the maximum likelihood estimates of rates (Supplementary Methods). Second, we implemented a MCMC procedure to build the refined maps. We assumed that recombination rate in each interval was independent from every other SNP interval. We modelled the underlying recombination rate within a SNP interval as a random variable drawn from a gamma distribution, with the prior parameters selected via inspection of the unrefined map. We used a Gibbs sampler to iteratively sample the SNP interval location of each recombination event, and then resampled recombination rates consistent with the event assignment. For each SNP interval, we obtained a mean posterior estimate of the recombination rate and the associated CI by running 1.3 million iterations, and discarding the first 300,000 as burn-in. We sampled from the chain every 100 iterations to obtain 10,000 samples. We applied this method separately to female and male recombination events, obtaining separate female and male refined maps. Our method to refine the estimates of recombination rates is available here: github.com/auton1/rMCMC. The recombination maps generated can be downloaded here: https://github.com/cbherer/Bherer_etal_SexualDimorphismRecombination..