from statsmodels.stats.multitest import fdrcorrection rejected, q-value = fdrcorrection(my_list_of_p_values)