Weir-Fst Vcftools
1
1
Entering edit mode
7.9 years ago
thjnant ▴ 130

Hello,

I am using the weir-fst in vcftools to calculate the pairwise fst. I have also written my own script using the Weir and Cockerham 1984.

The puzzling thing is that I get different estimates, the difference is between 0.004 to 0.01 but this is still puzzling to me, if both calculations are based on the same formula and considering the deterministic nature of weir-fst, I should get exactly the same estimates using exactly the same data.

I was wondering whether it is possible to get access to vcftools pairwise fst source code. I would also appreciate if you could share your ideas about this in case you have similar experience.

Cheers, Lucia

vcftools • 4.3k views
ADD COMMENT
0
Entering edit mode
7.9 years ago
Adam ★ 1.0k

Please make sure you're using the latest version of VCFtools. Older versions had a small bug when calculating Fst.

Thanks,

Adam

ADD COMMENT
0
Entering edit mode

good comment, I check it now

ADD REPLY
0
Entering edit mode

Yes, I am using the version VCFtools-v0.1.10.2, so, I think this should be fine.

ADD REPLY
0
Entering edit mode

Actually, there were a couple of bugs fixed in v0.1.11, so please upgrade. Thank you.

ADD REPLY
0
Entering edit mode

Just installed v0.1.11 and did the analysis again, the results are exactly the same as the previous analysis. To give you an idea about the difference, the number I get for my analysis for example for one pair of populations is 0.3692089 using my own script and 0.36488 using vcf tools weir fst.

ADD REPLY
0
Entering edit mode

The vcftools calculation has been checked against the geneclust R package (at least the version in the SVN), and gets the same results. If you compare to geneclust, do you get the same issue?

ADD REPLY

Login before adding your answer.

Traffic: 1331 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