Protocol Online logo
Top : New Forum Archives (2009-): : Bioinformatics and Biostatistics

How to calculate D' and R2 for linkage disequilibrium - (Jun/11/2013 )

Supposed i have 2 polymorphisms, SNPA and SNPB.
Individuals who carry the wildtype genotype for SNPA also carry the wildtype genotype for SNPB.
Individuals who carry the heterozygous genotype for SNPA also carry the heterozygous genotype for SNPB.

Individuals who carry the mutant genotype for SNPA also carry the mutant genotype for SNPB.

It is pretty obvious that the two polymorphisms are in complete linkage disequilibrium, right?
So the D' should equal to 1.

But, how do I calculate that?

I found a web-based software which seems to be easy to use, but I don't know what number should I key in onto each cell --> http://www.oege.org/software/cubex/

Now supposed I have 1000 wildtype, 800 heterozygous and 600 mutants for both polymorphisms, what number shall I put in the cells?

I tried the following (based on the frequencies of my genotypes):
2000---1800---1600
1800---1600---1400
1600---1400---1200

But I got a D' value of -0.011 instead of 1.

Can anyone help?

You may recommend other softwares or simple methods for calculating the D' and r2.

Thanks in advance.

-jamestoon1-

Nobody answered, but never mind, I already know the answer.
Supposed to fill as follows:

1000---0---0
0---800---0
0---0---600

-jamestoon1-