Question: Detecting outlier loci from pairwise Fst values
0
gravatar for Earendil
4.2 years ago by
Earendil20
Sweden
Earendil20 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 COMMENTlink modified 19 months ago by Shangzhe Zhang0 • written 4.2 years ago by Earendil20
0
gravatar for Shangzhe Zhang
19 months ago by
China
Shangzhe Zhang0 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 COMMENTlink written 19 months ago by Shangzhe Zhang0

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 REPLYlink written 19 months ago by Earendil20
0
gravatar for Shangzhe Zhang
19 months ago by
China
Shangzhe Zhang0 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 COMMENTlink modified 19 months ago • written 19 months ago by Shangzhe Zhang0
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: 1510 users visited in the last hour