U f»{j³ ã@s(ddlZdd„Zd dd„Zd dd „ZdS) éNcCs t|ƒ}t d|d¡t|ƒS)z.No frills empirical cdf used in fdrcorrection.é)ÚlenÚnpÚarangeÚfloat)ÚxZnobs©rúS/var/www/html/bwave/python/venv/lib/python3.8/site-packages/mne/stats/multi_comp.pyÚ_ecdf sr çš™™™™™©?Úindepc Cst |¡}|j}| ¡}t |¡}||}| ¡}|dkrDt|ƒ}n<|dkrxt dt dt|ƒd¡¡}t|ƒ|}nt dƒ‚|||k} |   ¡r¨t t  | ¡dƒ} nd} d| d| …<||} tj  | ddd …¡ddd …} d| | dk<| | |¡} | | |¡} | | fS) aTP-value correction with False Discovery Rate (FDR). Correction for multiple comparison using FDR :footcite:`GenoveseEtAl2002`. This covers Benjamini/Hochberg for independent or positively correlated and Benjamini/Yekutieli for general or negatively correlated tests. Parameters ---------- pvals : array_like Set of p-values of the individual tests. alpha : float Error rate. method : 'indep' | 'negcorr' If 'indep' it implements Benjamini/Hochberg for independent or if 'negcorr' it corresponds to Benjamini/Yekutieli. Returns ------- reject : array, bool True if a hypothesis is rejected, False if not. pval_corrected : array P-values adjusted for multiple hypothesis testing to limit FDR. References ---------- .. footbibliography:: )Úir ÚpZposcorr)ÚnZnegcorrçð?rz&Method should be 'indep' and 'negcorr'rTNéÿÿÿÿ)rÚasarrayÚshapeZravelZargsortr ÚsumrrÚ ValueErrorÚanyÚmaxZnonzeroÚminimumÚ accumulateZreshape) ZpvalsÚalphaÚmethodZ shape_initZ pvals_sortindZ pvals_sortedZ sortrevindZ ecdffactorÚcmÚrejectZ rejectmaxZpvals_corrected_rawZpvals_correctedrrr Úfdr_corrections.       rcCs4t |¡}|t|jƒ}|jdd}||k}||fS)aˆP-value correction with Bonferroni method. Parameters ---------- pval : array_like Set of p-values of the individual tests. alpha : float Error rate. Returns ------- reject : array, bool True if a hypothesis is rejected, False if not. pval_corrected : array P-values adjusted for multiple hypothesis testing to limit FDR. r)r)rrrÚsizeZclip)ZpvalrZpval_correctedrrrr Úbonferroni_correctionMs   r )r r )r )Únumpyrr rr rrrr Ús <