How to calculate Hudson's Fst at Population and SNP level?
5
1
Entering edit mode
7.6 years ago
bibb77 ▴ 90

Dear BioStars community:

I would like to know if there is a software or script that calculates Hudson's Fst estimator (http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1205159/) at both Population and SNP level, the first to asses genetic distance between population and the second to detect highly divergent SNPs. I have used the EIGENSOFT software but this only produces Fst values at population level (pairwise matrix), but not for single SNPs. I also have used the well known method from Weir and Cockerham implemented in PLINK, but Bhatia et al, 2013 (http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3759727/#B16) recommends to use Hudson's estimator because it behaves better with uneven population sizes.

Thanks in advance.

Fst SNP Population-differentiation • 6.6k views
ADD COMMENT
1
Entering edit mode
6.6 years ago

Using a VCF file as input, SNiPlay implements VCFtools that calculates Fst (Weir and Cockerham’s) for each marker... http://sniplay.southgreen.fr/cgi-bin/analysis_v3.cgi

ADD COMMENT
0
Entering edit mode
7.4 years ago
gideond • 0

Hi there,

Any updates on this? Did you find something useful? Or implement the estimator yourself? I found this: https://bio.cst.temple.edu/~hey/program_files/Sites/SITES_Documentation.htm but it doesn't look well maintained and I am pressed to find actual usage examples in the documentation.

I'm thinking about implementing a simple version myself.

Hope to hear your thoughts.

ADD COMMENT
0
Entering edit mode
7.4 years ago
gideond • 0

The second paper you mentioned actually links to a software package called EIGENSOFT. It is available here: http://www.hsph.harvard.edu/alkes-price/software/. I have to admit that this has not made things any easier for me as I do not understand what a block-jackknife estimator is and there is no mention of Hudson estimators anywhere in the documentation (or source code). I'll continue keeping you posted.

ADD COMMENT
0
Entering edit mode

I was also frustrated by the lack of detail in the documentation--I see no mention of the type of Fst estimator.

ADD REPLY
0
Entering edit mode
6.6 years ago

Try this http://spsmart.cesga.es/

ADD COMMENT

Login before adding your answer.

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