https://github.com/nfusi/qvalue
Tip revision: c215bdcfac262d72d1382263c7d507dd90535642 authored by Nicolò Fusi on 14 July 2016, 21:34:32 UTC
Merge pull request #3 from alexisboukouvalas/master
Merge pull request #3 from alexisboukouvalas/master
Tip revision: c215bdc
README.md
qvalue
======
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.