Converts p-values in q-values, see (Storey and Tibshirani, 2003)

The only function that needs to be called is 'estimate()'. It will accept a numpy array of pvalues and will return a numpy array of qvalues. For instance, given some uniform p-values:

import numpy as np 
import qvalue

pv = np.random.uniform(0.0, 1.0, size = (1000,))

it's possible to convert them in q-values by calling

 qv = qvalue.estimate(pv)

estimate() also includes a memory-efficient (but not CPU-efficient) procedure for the case when it's not possible to store both the p-values and the q-values in memory at the same time. In addition it's possible to change the number of hypotheses tested or some parameters of the procedure (like pi0).

We refer to the documentation (in ipython, 'qvalue.estimate?') for more informations.