https://github.com/nfusi/qvalue
Raw File
Tip revision: c215bdcfac262d72d1382263c7d507dd90535642 authored by Nicolò Fusi on 14 July 2016, 21:34:32 UTC
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.
back to top