SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Direction of assembled genome/output? bio_informatics Bioinformatics 10 03-25-2015 07:10 AM
Direction! I could use some! NYGen Introductions 2 08-22-2014 05:07 AM
Canonical direction and assembly jlewi De novo discovery 0 09-20-2011 06:23 AM
Sequence direction in alignment newbietonextgen Bioinformatics 1 01-18-2011 11:41 AM
Single direction reads for variants Moggs General 11 06-22-2010 06:55 AM

Reply
 
Thread Tools
Old 04-08-2015, 03:01 AM   #1
man4ish
Junior Member
 
Location: OK

Join Date: Jan 2010
Posts: 7
Default How to calculate LD in 2 direction?

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 10000 --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?

Last edited by man4ish; 04-08-2015 at 05:05 AM.
man4ish is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 12:46 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO