Probability waves: Adaptive cluster-based correction by convolution of p-value series from mass univariate analysis

• Published in: Journal of neuroscience methods 2021-06, Vol.357, p.109155-109155, Article 109155
• Main Authors: Abramov, Dimitri Marques, Miranda de Sá, Antonio Mauricio F.L.
• Format: Article
• Language: English
• Subjects: Cluster-based statistics; Convolution; False discovery rate; Mass univariate analysis; Multiple comparisons; Type-II error
• ISSN: 0165-0270; 1872-678X
• DOI: 10.1016/j.jneumeth.2021.109155

•The p-values correction from multiple pairwise comparisons of ERP waves is proposed.•The correction was performed by convolving p-values vector with a Gaussian window.•This window length was obtained based on the auto-correlation from the ERP signals.•Simulations indicate our method to be less conservative than False Discovery Rate.•The method seems to correct Type I and II errors in Mass Univariate Analysis. Methods for p-value correction are criticized for either increasing Type II error or improperly reducing Type I error in large exploratory data analysis. This text considers patterns in probability vectors resulting from mass univariate analysis to correct p-values, where clusters of significant p-values may indicate true H0 rejection. We used ERP experimental data from control and ADHD boys to test the method. The Log10 of p-vector was convolved with a Gaussian window whose length was set as the shortest lag above which autocorrelation of each ERP wave may be assumed to have vanished. We realized Monte-Carlo simulations (MC) to (1) evaluate confidence intervals of rejected and non-rejected areas of our data, (2) to evaluate differences between corrected and uncorrected p-vectors or simulated ones in terms of distribution of significant p-values, and (3) to empirically verify the type-I error rate (comparing 10,000 pairs of mixed samples whit control and ADHD subjects). The differences between simulation or raw p-vector and corrected p-vectors were, respectively, minimal and maximal for window length set by autocorrelation in p-vector convolution. Our method was less conservative while FDR methods rejected basically all significant p-values.The MC simulations presented 2.78 ± 4.83% of difference (20 channels) from corrected p-vector, while difference from raw p-vector was 596 ± 5.00% (p = 0.0003). As a cluster-based correction, the present new method seems to be biological and statistically suitable to correct p-values in mass univariate analysis of ERP waves, which adopts adaptive parameters to correction.