limix.stats.multipletests

limix.stats.multipletests(pvals, alpha=0.05, method='hs', is_sorted=False)[source]

Test results and p-value correction for multiple tests.

Parameters

  • pvals (array_like) – Uncorrected p-values.

  • alpha (float) – FWER, family-wise error rate, e.g. 0.1.

  • method (string) –

    Method used for testing and adjustment of pvalues. Can be either the full name or initial letters. Available methods are

    `bonferroni` : one-step correction
`sidak` : one-step correction
`holm-sidak` : step down method using Sidak adjustments
`holm` : step-down method using Bonferroni adjustments
`simes-hochberg` : step-up method  (independent)
`hommel` : closed method based on Simes tests (non-negative)
`fdr_bh` : Benjamini/Hochberg  (non-negative)
`fdr_by` : Benjamini/Yekutieli (negative)
`fdr_tsbh` : two stage fdr correction (non-negative)
`fdr_tsbky` : two stage fdr correction (non-negative)

  • is_sorted (bool) – If False (default), the p_values will be sorted, but the corrected pvalues are in the original order. If True, then it assumed that the pvalues are already sorted in ascending order.

Returns

  • reject (ndarray, boolean) – True for hypothesis that can be rejected for given alpha.

  • pvals_corrected (ndarray) – P-values corrected for multiple tests.

  • alphacSidak (float) – Corrected alpha for Sidak method.

  • alphacBonf (float) – Corrected alpha for Bonferroni method.

Notes

This is a wrapper around a function from the statsmodels package.