I am trying to use the LDhat package to estimate the population recombination rate of my samples (phased data based on genowide SNPs) to then estimate the effective population size (Ne). I have seen that you can do that (based on some publications) but I couldn't figure out the way to do that. I think that there was only some minimal algebra involved , but again I am not very sure...
I think it goes something like this:
rho = 4 Ne x r
Ne : Effective population size r : recombination rate
If I can get the recombination rate of my region, say for chromosome 1, then I could use that estimate to solve for Ne. Is there some way to get that rate? I've been exploring the deCODE website, where there is information on genetic map, but there is only the recombination rate for pair of SNPs and not for the region.
Any help would be appreciated.
Thanks in advance.