Question: Weir-Fst Vcftools
1
gravatar for thjnant
5.4 years ago by
thjnant90
Germany
thjnant90 wrote:

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 • 2.9k views
ADD COMMENTlink modified 5.4 years ago by Adam990 • written 5.4 years ago by thjnant90
0
gravatar for Adam
5.4 years ago by
Adam990
United States
Adam990 wrote:

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

Thanks,

Adam

ADD COMMENTlink written 5.4 years ago by Adam990

good comment, I check it now

ADD REPLYlink written 5.4 years ago by thjnant90

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

ADD REPLYlink written 5.4 years ago by thjnant90

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

ADD REPLYlink written 5.4 years ago by Adam990

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 REPLYlink written 5.4 years ago by thjnant90

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 REPLYlink written 5.4 years ago by Adam990
Please log in to add an answer.

Help
Access

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