Probabilistic TFCE: A generalized combination of cluster size and voxel intensity to increase statistical power

&NA; The threshold‐free cluster enhancement (TFCE) approach integrates cluster information into voxel‐wise statistical inference to enhance detectability of neuroimaging signal. Despite the significantly increased sensitivity, the application of TFCE is limited by several factors: (i) generalisation to data structures, like brain network connectivity data is not trivial, (ii) TFCE values are in an arbitrary unit, therefore, P‐values can only be obtained by a computationally demanding permutation‐test. Here, we introduce a probabilistic approach for TFCE (pTFCE), that gives a simple general framework for topology‐based belief boosting. The core of pTFCE is a conditional probability, calculated based on Bayes' rule, from the probability of voxel intensity and the threshold‐wise likelihood function of the measured cluster size. In this paper, we provide an estimation of these distributions based on Gaussian Random Field theory. The conditional probabilities are then aggregated across cluster‐forming thresholds by a novel incremental aggregation method. pTFCE is validated on simulated and real fMRI data. The results suggest that pTFCE is more robust to various ground truth shapes and provides a stricter control over cluster “leaking” than TFCE and, in many realistic cases, further improves its sensitivity. Correction for multiple comparisons can be trivially performed on the enhanced P‐values, without the need for permutation testing, thus pTFCE is well‐suitable for the improvement of statistical inference in any neuroimaging workflow. Implementation of pTFCE is available at https://spisakt.github.io/pTFCE. Graphical abstract Figure. No caption available. HighlightsPTFCE is a generalized probabilistic approach for threshold‐free cluster enhancement (TFCE).PTFCE overcomes the limitations of cluster inference.The method is robust to the cluster topology and provides a strict control over false positives.Enhanced p‐values can trivially be corrected for multiple comparisons, without permutation test.easy to integrate into any neuroimaging analysis workflow.

[1]  William A. Cunningham,et al.  Type I and Type II error concerns in fMRI research: re-balancing the scale. , 2009, Social cognitive and affective neuroscience.

[2]  Thomas E. Nichols,et al.  Adjusting the effect of nonstationarity in cluster-based and TFCE inference , 2011, NeuroImage.

[3]  Karl J. Friston,et al.  Detecting Activations in PET and fMRI: Levels of Inference and Power , 1996, NeuroImage.

[4]  K. Worsley,et al.  Unified univariate and multivariate random field theory , 2004, NeuroImage.

[5]  Daniel Rueckert,et al.  Tract-based spatial statistics: Voxelwise analysis of multi-subject diffusion data , 2006, NeuroImage.

[6]  Karl J. Friston,et al.  Assessing the significance of focal activations using their spatial extent , 1994, Human brain mapping.

[7]  Satrajit S. Ghosh,et al.  Nipype: A Flexible, Lightweight and Extensible Neuroimaging Data Processing Framework in Python , 2011, Front. Neuroinform..

[8]  Karl J. Friston,et al.  A unified statistical approach for determining significant signals in images of cerebral activation , 1996, Human brain mapping.

[9]  Christian Genest,et al.  Combining Probability Distributions: A Critique and an Annotated Bibliography , 1986 .

[10]  Mark W. Woolrich,et al.  Advances in functional and structural MR image analysis and implementation as FSL , 2004, NeuroImage.

[11]  K. Wiech,et al.  The Effect of Treatment Expectation on Drug Efficacy: Imaging the Analgesic Benefit of the Opioid Remifentanil , 2011, Science Translational Medicine.

[12]  Thomas E. Nichols,et al.  Controlling the familywise error rate in functional neuroimaging: a comparative review , 2003, Statistical methods in medical research.

[13]  M. Stone The Opinion Pool , 1961 .

[14]  Stephen M. Smith,et al.  Threshold-free cluster enhancement: Addressing problems of smoothing, threshold dependence and localisation in cluster inference , 2009, NeuroImage.

[15]  Karl J. Friston,et al.  Statistical parametric maps in functional imaging: A general linear approach , 1994 .

[16]  Anjali Krishnan,et al.  Cluster-extent based thresholding in fMRI analyses: Pitfalls and recommendations , 2014, NeuroImage.

[17]  Jonathan D. Cohen,et al.  Improved Assessment of Significant Activation in Functional Magnetic Resonance Imaging (fMRI): Use of a Cluster‐Size Threshold , 1995, Magnetic resonance in medicine.

[18]  CM Bennett,et al.  Neural correlates of interspecies perspective taking in the post-mortem Atlantic Salmon: an argument for multiple comparisons correction , 2009, NeuroImage.

[19]  Krzysztof J. Gorgolewski,et al.  Preprocessed Consortium for Neuropsychiatric Phenomics dataset , 2017, F1000Research.

[20]  Krzysztof J. Gorgolewski,et al.  A phenome-wide examination of neural and cognitive function , 2016, Scientific Data.

[21]  Volker Schmid,et al.  Working with the DICOM and NIfTI Data Standards in R , 2011 .

[22]  H. Pashler,et al.  Puzzlingly High Correlations in fMRI Studies of Emotion, Personality, and Social Cognition 1 , 2009, Perspectives on psychological science : a journal of the Association for Psychological Science.

[23]  D. Chakraborty,et al.  Free-response methodology: alternate analysis and a new observer-performance experiment. , 1990, Radiology.

[24]  Thomas E. Nichols Multiple testing corrections, nonparametric methods, and random field theory , 2012, NeuroImage.

[25]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[26]  Klaus Scheffler,et al.  Inflated false negative rates undermine reproducibility in task-based fMRI , 2017, bioRxiv.

[27]  A. Szalay,et al.  The statistics of peaks of Gaussian random fields , 1986 .