Question: Using Ldhat To Estimate Effective Population Size
2
gravatar for GabrielMontenegro
5.1 years ago by
United Kingdom
GabrielMontenegro510 wrote:

Hi everyone:

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.

snps recombination • 2.7k views
ADD COMMENTlink modified 5.1 years ago by Giovanni M Dall'Olio26k • written 5.1 years ago by GabrielMontenegro510
1
gravatar for Giovanni M Dall'Olio
5.1 years ago by
London, UK
Giovanni M Dall'Olio26k wrote:

You may want to read how they did it in this paper: Prado-Martinez et al, Great ape genetic diversity and population history. In particular, see Figure 3 of the paper, and page 199 of the Supplementary Materials.

You can estimate the effective population size, and how it varied over time, using the psmc software.

ADD COMMENTlink written 5.1 years ago by Giovanni M Dall'Olio26k

Dear Giovanni. Thanks for the reply I will take a look at those references. Lets see how it goes.. Thanks again.

ADD REPLYlink written 5.1 years ago by GabrielMontenegro510

For what I see they estimate Ne based on the population mutation rate, and not the recombination rate. And is it posible to use the psmc software on genotype phased data? I don't have sequences.

ADD REPLYlink written 5.1 years ago by GabrielMontenegro510
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1072 users visited in the last hour