I have a huge SNP dataset (hundreds of individuals and about 50K markers) including several populations. I need to calculate the value of the Fis (inbreeding coefficient) index for each population, including the level of statistical significance of the index value.
The software and scripts I found so far only return the index value but no statistical significance.
Does anyone know a package, R script or algorithm able to do this?
Thank you in advance.
Licia Colli, PhD
Istituto di Zootecnica – Facolt di Agraria
Universit Cattolica del S. Cuore
via Emilia Parmense, 84
29122 Piacenza (PC)
Tel: +39 – 0523599205
Fax: +39 – 0523599276