Computes F statistics according to Weir and Cockerham's estimators. Missing values are allowed but as of today, the NA code is only treated as an extra allele which might bias the result. This function should not be used on haploid data.

1 | ```
Fstat.output(coordinates,genotypes,ploidy,burnin,path.mcmc)
``` |

`coordinates` |
Matrix with one line per individual and two columns |

`genotypes` |
Genotypes of individuals. A matrix with one line per individual and 2 columns per locus |

`burnin` |
Integer: number of saved iterations to discard. |

`ploidy` |
Integer: 1 or 2 (default is 2). Do not use for haploid data. |

`path.mcmc` |
Path to output files directory |

A list with components

`Pairwise.Fis ` |
A matrix of real numbers estimating the pairwise Fis |

`Pairwise.Fst ` |
A matrix of real numbers estimating the pairwise Fst |

Arnaud Estoup for original code in Turbo Pascal.

Translation in Fortran and interface with R by Gilles Guillot

Weir, B.S. and C.C. Cockerham, Estimating F-statistics for the analysis of population structure, Evolution, 1984, vol. 38, 1358-1370.

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.