[adegenet-forum] Fstat vs pairwise.fst for 2 populations

Manuel A. Morales mmorales at williams.edu
Fri Nov 10 17:20:05 CET 2017


I may have a fundamental misunderstanding of what's happening but it
seems to me that the functions Fst() and pairwise.fst() should give the
same value for the case of two populations, which they do not.

A reproducible example:
data(nancycats)
obj1 <- seppop(nancycats)$P01
obj2 <- seppop(nancycats)$P02
obj3 <- repool(obj1, obj2)
fstat(obj3)
pairwise.fst(obj3)

And output:
> fstat(obj3)
            pop       Ind
Total 0.1307741 0.2804306
pop   0.0000000 0.1721722
> pairwise.fst(obj3)
         1
2 0.080185

Any help would be very much appreciated.

Best,
Manuel



More information about the adegenet-forum mailing list