R Script For Fst (Reynolds Et Al. 1983)

Question

[user]

Anybody has a simple R script for Fst (Reynolds et al. 1983). thanks

[Neilfws]

In future a little more detail would be appreciated, such as title of paper and a link to it. In fact, a simple Google search using the title returns an answer, see below.

Answer 1

[Neilfws]

Looks like the function:

dist.genpop(... , method = 3)

from the package adegenet will calculate this for you.

[Neilfws]

By the way, that was the 3rd Google hit when searching using the title of the Reynolds paper.

[user]

Thanks man, Indeed this package outputs dist values over a bunch of markers between populations.

I have a input file with SNP, AFpop1, AFpop2 (with different population size) and what I need is an estimate of Fst for each SNP position. A link for the formula is here http://bioinformatics.org/~tryphon/populations/formules/reynold_drw.xhtml

[Jorge Amigo]

I had to develop a program that had to deal with several different situations for computing Fst values. I found it easier to program it myself using the stripped formulas on this worked example:

http://www.uwyo.edu/dbmcd/popecol/maylects/fst.html