This PR implements support for estimation of pi0 and lfdr/PEP by Storey's methods. The implementation follows the Bioconductor/qvalue R code, although with kernel density estimators and spline fitting functions provided by SciPy introducing slight differences.
This PR implements support for estimation of pi0 and lfdr/PEP by Storey's methods. The implementation follows the Bioconductor/qvalue R code, although with kernel density estimators and spline fitting functions provided by SciPy introducing slight differences.