How to calculate LD in 2 direction?
0
0
Entering edit mode
9.5 years ago

Hi,

I have list of snps in vcf format and I need to calculate LD in such a way that window should be +_90 either side.

I tried both plink and vcf-tools but both works in forward direction(i.e 3rd column value is greater than 2nd column value).

vcftools --vcf testset --hap-r2 --ld-window-bp 90 --out OUTFILE

Sample Output

CHR POS1 POS2 N_CHR R^2 D Dprime
11 60250005 60250039 1008 1.19644e-05 -1.18103e-05 -1
11 60250005 60250058 1008 6.84498e-05 -6.39723e-05 -1
11 60250039 60250058 1008 0.000830469 -0.000767668 -1
11 60250039 60250096 1008 -nan 0 -nan
11 60250039 60250125 1008 -nan 0 -nan
11 60250058 60250096 1008 -nan 0 -nan
11 60250058 60250125 1008 -nan 0 -nan
11 60250058 60250140 1008 -nan 0 -nan
11 60250096 60250125 1008 -nan 0 -nan
11 60250096 60250140 1008 -nan 0 -nan
11 60250096 60250157 1008 -nan 0 -nan
11 60250125 60250140 1008 -nan 0 -nan
11 60250125 60250157 1008 -nan 0 -nan
11 60250140 60250157 1008 -nan 0 -nan
11 60250140 60250230 1008 -nan 0 -nan
11 60250157 60250230 1008 -nan 0 -nan
11 60250230 60250267 1008 9.84931e-05 -9.64506e-05 -1
11 60250230 60250275 1008 0.0004662 -0.000434028 -1
11 60250267 60250275 1008 0.000938967 -0.000868056 -1

How can I get desired output?

SNP • 1.9k views
ADD COMMENT

Login before adding your answer.

Traffic: 1664 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6