Calculate hapflk p-values
Documentation utilisateur
07/09/2015
This python script estimates the hapflk chi-squared density and corresponding p-values.
The script requires the numpy, scipy and statsmodels python packages.
The input file must be a .hapflk result file, possibly containing all values from a genome-wide scan. The script creates a new file (suffix .hapflk_sc) which is the original file plus 2 new colums: the rescaled statistic and the corresponding p-values.
Usage :
python scaling_chi2_hapflk.py all.hapflk K N
where K is the number of clusters used and N the number of populations.