Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Detecting outlier loci from pairwise Fst values
0
gravatar for Earendil
4.2 years ago by
Earendil • 20
Sweden
Earendil • 20 wrote:

Having calculated pairwise Fst values with vcftools, I now need to find the threshold for outlier loci. I've decided to follow the second suggestion from an answer of this thread: Calculating statistically significant outlier for Pairwise Fst obtained from VCFTools which is:

"2.  Permute your genotypes and re-run Fst many times.  This would be considered an empirical p-value, or probability."

Since I am having no statistical background, would there be a simple explanation on how to implement this?
vcftools fst outlier • 3.0k views
ADD COMMENT • link •
modified 19 months ago by Shangzhe Zhang • 0 • written 4.2 years ago by Earendil • 20
0
gravatar for Shangzhe Zhang
19 months ago by
Shangzhe Zhang • 0
China
Shangzhe Zhang • 0 wrote:

Dear Earendil, I'd like to know that if you solve the problem. I'm new to NGS analysis and I'm stuck in this problem. Hope for your help!
ADD COMMENT • link written 19 months ago by Shangzhe Zhang • 0

Dear Shangzhe,

That was long time ago, I didn't find out how to do that so instead I used the software Bayescan which directly detects Fst outliers, you would might want to take a look at that.

ADD REPLY • link written 19 months ago by Earendil • 20
0
gravatar for Shangzhe Zhang
19 months ago by
Shangzhe Zhang • 0
China
Shangzhe Zhang • 0 wrote:

Dear Earendil,

Thanks for you reply. Coincidently one of my colleage used this method to find the outliers of Fst.

On my opinion this method is useful when your windowed-Fst data normally distributed and its mean is almost 1, which is difficult to indentify the outliers using the p-value against the normal distriobution defined by its mean and standard distribution.

The method is anctually called mutiple testing. For our data, this method aviod the bais due to the windows with few SNP pairs but with the Fst value relatively high.

Maybe this is ambiguous, I copy the link of the article of my colleage and the article of the method below. Hope these can help:

http://www.pnas.org/content/115/2/E236 (see the "method" part) https://www.researchgate.net/publication/40688145_How_does_multiple_testing_correction_work_Nat_Biotechnol

And thank you again. I will see Bayescan.
ADD COMMENT • link modified 19 months ago • written 19 months ago by Shangzhe Zhang • 0
Please log in to add an answer.

Similar posts • Search »

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