nimare.stats.nlogp_fdr

nlogp_fdr(nlogp, method='bh')

Perform FDR correction on nlogp values.

Added in version 0.21.0.

The log-space counterpart of pymare.stats.fdr(), step for step: the step-up procedure only ever divides a p-value by a positive factor and takes running minima, both of which carry over to logs unchanged. Sorting by nlogp is sorting by p, so the set of tests declared significant at any alpha is identical.

Parameters:

• nlogp (numpy.ndarray) – Natural logarithms of the uncorrected p-values.

• method ({"bh", "by"}, optional) – Either “bh” (Benjamini-Hochberg [1]) or “by” (Benjamini-Yekutieli [2]). Default is “bh”.

Returns:

Natural logarithms of the corrected p-values.

Return type:

numpy.ndarray

References

[1]

Yoav Benjamini and Yosef Hochberg. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal statistical society: series B (Methodological), 57(1):289–300, 1995. URL: https://doi.org/10.1111/j.2517-6161.1995.tb02031.x, doi:10.1111/j.2517-6161.1995.tb02031.x.

[2]

Yoav Benjamini and Daniel Yekutieli. The control of the false discovery rate in multiple testing under dependency. The Annals of Statistics, 29(4):1165 – 1188, 2001. URL: https://doi.org/10.1214/aos/1013699998, doi:10.1214/aos/1013699998.