Differing output: vcftools' --weir-fst-pop and R hierfstats package's varcomp.glob()
1
0
Entering edit mode
12 weeks ago
S ▴ 10

I'm getting different Fst values when I use vcftools --weir-fst-pop vs. the varcomp.glob() function from the R hierfstats package. I'm not sure why because they both calculate Fst based on Weir & Cockerham's 1984 paper. Can someone please explain why?

Vcftools documentation states that this is calculated on a "per-site basis." (link: https://vcftools.github.io/man_latest.html)

Hierfstats' varcomp.glob() documentation states that this function "will produce the estimation of the variance components for each locus and overall, as well as the matrix of hierarchical F - statistics". (link: https://rdrr.io/cran/hierfstat/man/varcomp.glob.html)

Fst, group 1: vcftools = 0.079; hierfstats = 0.

Fst, group 2: vcftools = 0.049; hierfstats = 0.005.

Thanks so much

vcftools hierfstats fst • 324 views
ADD COMMENT
0
Entering edit mode
11 weeks ago
S ▴ 10

Just found this post comparing methods for computing Fst. These differences could be related to differences in how each tool handles missing data, among other things: https://thierrygosselin.github.io/assigner/articles/web_only/fst_comparisons.html

ADD COMMENT

Login before adding your answer.

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