Question: Calculating Tajima's D from vcf data using bed files
gravatar for spiral01
15 months ago by
spiral0180 wrote:

I need to calculate Tajima's D for vcf data, using a bed file to indicate the windows in which I want to calculate the statistic. I tried to use vcftools like so:

vcftools --gzvcf file.vcf.gz --bed file.bed --TajimaD 500 --stdout > output.txt

However vcftools requires that you indicate the size of the windows for which you wish to calculate the statistic (the 500 above). As such it ends up calculating Tajima's D for each scrolling window, instead of for each of the windows stipulated in my bed file.

Is there anyway to circumvent this, or any tool that is better suited to this task?

snp • 1.2k views
ADD COMMENTlink modified 5 months ago by Gabriel20 • written 15 months ago by spiral0180
gravatar for Gabriel
5 months ago by
Gabriel20 wrote:

Hi @spiral01, I am with the same issue here. Did you find a way to do it? Thanks! Gabriel

ADD COMMENTlink written 5 months ago by Gabriel20
Please log in to add an answer.


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