Any Simple Script To Calculate Fst Value Between Populations From Genomewide Snp Data (Available In Plink Format)
3
6
Entering edit mode
10.7 years ago
User 0422 ▴ 150

Hi All,

I am looking for a simple script which calculate the fst between different populations, taking input as PLINK files (either ped or bed).

Thanks

snp ped bed fst plink • 15k views
ADD COMMENT
0
Entering edit mode

I'm surprised this hasn't been answered yet... want it too!

ADD REPLY
2
Entering edit mode
10.2 years ago
Zach Stednick ▴ 660

Why not just use an existing program such as Genepop and write your own script to generate input files? Genepop provides plenty of example input files and once you have input your data it runs very fast.

ADD COMMENT
1
Entering edit mode
10.2 years ago

Because Fst calculations are so much less complicated than iHS, we take a very programmatic approach to Fst calculations for our SNP panels. We have written a program to do this in a high-throughput manner. My colleague wrote a script in R or Perl, but this is so simple that it could even be done in Excel.

ADD COMMENT
0
Entering edit mode

I'm looking to do some quick Fst calculations for a single SNP. I will use excel or a simple R script. I am looking for a good reference to an equation to calculate pair-wise Fst...

ADD REPLY
1
Entering edit mode
9.9 years ago
Davy ▴ 210

The R package hierfstat can do this. As well as computing Fst values, it can group your populations and compute Fst values for the groups. Very flexible. For in the input file I think plink ped files may need some reworking as it's been sometime since I've used it but the vignette can help you with that. Also has some cool functions for generating sample SNP data to vary mutation and migration rates and things like that.

Documentation here: http://cran.r-project.org/web/packages/hierfstat/hierfstat.pdf

ADD COMMENT

Login before adding your answer.

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