neurosynth/analysis/stats.py
def fdr(p, q=.05):
    """ Determine FDR threshold given a p value array and desired false
    discovery rate q. """
    s = np.sort(p)
    nvox = p.shape[0]
    null = np.array(range(1, nvox + 1), dtype='float') * q / nvox
    below = np.where(s <= null)[0]
    return s[max(below)] if len(below) else -1
def fdr(p, q=.05):
    s = np.sort(p)
    nvox = p.shape[0]
    null = np.array(range(1, nvox + 1), dtype='float') * q / nvox
    below = np.where(s <= null)[0]
    return s[max(below)] if len(below) else -1
def fdr ( p , q = .05 ) : s = np . sort ( p ) nvox = p . shape [ 0 ] null = np . array ( range ( 1 , nvox + 1 ) , dtype = 'float' ) * q / nvox below = np . where ( s <= null ) [ 0 ] return s [ max ( below ) ] if len ( below ) else - 1
Determine FDR threshold given a p value array and desired false
    discovery rate q.
Determine FDR threshold given a p value array and desired false discovery rate q .