How to generate a smoothing plot for Fst along the genome from Lositan output?

Hi All,

I want to make a genome wide Fst plot using the Lositan output along with P-Values calculated by lositan.

Any inputs will be appreciated.

Thanks

snp lositan fst

• 3,907 views • link

updated 11.5 years ago by Zev.Kronenberg

• written 11.7 years ago by Anurag •

1 answer

I don't know how to do it from Lositan, but I have all the details how to do it using my code:

https://github.com/jewmanchue/vcflib/wiki

Check out the library, run wcFst:

https://github.com/jewmanchue/vcflib/wiki/Association-testing-with-GPAT

smooth it:

https://github.com/jewmanchue/vcflib/wiki/Smoothing-with-GPAT

• 0 views • link

written 11.5 years ago by Zev.Kronenberg

More posts like this

Compute P-Values for SNP Fst from Pairwise Fst Calculations (--weir-fst-pop)
written 1.6 years ago by Armin •

Hi Everyone, I've calculated pairwise Fst values for two populations using PLINK's --weir-fst-pop option, and I now have Fst values for each SNP. How can …
manhattan plot
written 5.0 years ago by anithanagaraj93 •

I calculated fst values by using following command, "vcftools --vcf a.vcf --weir-fst-pop pop1.txt --weir-fst-pop pop2.txt --out ex" and got negative and -nan Fst values. if …
Fst Calculation and Visualization
written 6.2 years ago by safiq713 •

Dear all, I am working with ddRAD data, and I have allele and genotype frequencies of different SNPs in different populations. I have calculated Fst …
Fst with p-values
written 8.8 years ago by Simo •

Hi all, I've been working with genome-wide data for about 200 samples. I have 150.000 SNPs and I already calculated the Fst using Eigenstrat. I …
Lositan using 2 populations only? Loci falling in straight line where Fst=He
written 10.8 years ago by lauren.c.white89 •

Hello, I'm using Lositan with 2 populations and around 200 SNP loci. Most of my loci fall in a straight line on the plot where …
Calculating Fst from simulations
written 10.8 years ago by Rubal

Hi All, I calculated Weir and Cockerham's Fst using VCF tools and want to use simulations to have a feeling whether my high Fst regions …
Calculating statistically significant outlier for Pairwise Fst obtained from VCFTools
written 11.8 years ago by Anurag •

Hi, I calculated pairwise Fst using VCFTools: vcftools --vcf input_data.vcf --weir-fst-pop population_1.txt --weir-fst-pop population_2.txt --out pop1_vs_pop2 What method I should use for statistical significance to …
Plotting Fst
written 12.6 years ago by thjnant

Hello, I have a huge SNP dataset of 25 millions SNP genome wide. I have divided it into different chromosomes and this FST plot …
Mapping Fst Values To A Manhattan Plot By Rs#
written 13.1 years ago by confusedious

Hello everyone, I have genome-wide pairwise Fst values for ~1 million SNPs and would like to put them onto a Manhattan plot. My data files …
R: Manhattan Plot Of Fst Values Instead Of -Log(P)
written 14.4 years ago by Abdel

Hi all, I have a tab delimited file with the following columns: <pre><code>CHR SNP bp Fst </code></pre> where Fst is the <a href='http://en.wikipedia.org/wiki/Fixation_index'>fixation index</a>. Does …