Note that this maximization is done using standard numerical methods implemented in GNU Scientific Library. The underlying idea of this approach is to characterize the extreme tail of the null distribution without directly sampling from it, something that would entail a huge computational burden. Finally, we can approximate an adjusted P-value pb from the best nominal P-value pn and from the maximum likelihood (ML) fitted beta distribution with: (5)pb=P(U≤pn)
