Calculating adjusted p-values in Python

python statistics p-value q-value
erikfas picture erikfas · Aug 7, 2014 · Viewed 27k times · Source

So, I've been spending some time looking for a way to get adjusted p-values (aka corrected p-values, q-values, FDR) in Python, but I haven't really found anything. There's the R function p.adjust, but I would like to stick to Python coding, if possible. Is there anything similar for Python?

If this is somehow a bad question, sorry in advance! I did search for answers first, but found none (except a Matlab version)... Any help is appreciated!

Answer

Josef picture Josef · Aug 8, 2014

It is available in statsmodels.

http://statsmodels.sourceforge.net/devel/stats.html#multiple-tests-and-multiple-comparison-procedures

http://statsmodels.sourceforge.net/devel/generated/statsmodels.sandbox.stats.multicomp.multipletests.html

and some explanations, examples and Monte Carlo http://jpktd.blogspot.com/2013/04/multiple-testing-p-value-corrections-in.html

Related questions

Finding Two-Tailed P Value from t-distribution and Degrees of Freedom in Python

How do I determine the P Value of a t-distrobution with n degrees of freedom. Research on this subject points me to this stack exchange answer: https://stackoverflow.com/a/17604216 I assume np.abs(tt) is the T-value, but how …

Read More
Calculating arithmetic mean (one type of average) in Python

Is there a built-in or standard library method in Python to calculate the arithmetic mean (one type of average) of a list of numbers?

Read More
How to normalize an array in NumPy?

I would like to have the norm of one NumPy array. More specifically, I am looking for an equivalent version of this function def normalize(v): norm = np.linalg.norm(v) if norm == 0: return v return v / norm Is there …

Read More