AUTHOR=Bucciarelli Vittoria , Vogel Dorian , Wårdell Karin , Coste Jérôme , Blomstedt Patric , Lemaire Jean-Jacques , Guzman Raphael , Hemm Simone , Nordin Teresa 
  
TITLE=Towards robust probabilistic maps in Deep Brain Stimulation: exploring the impact of patient number, stimulation counts, and statistical approaches
  
JOURNAL=Frontiers in Computational Neuroscience
  
VOLUME=Volume 19 - 2025
  
YEAR=2026
  
URL=https://www.frontiersin.org/journals/computational-neuroscience/articles/10.3389/fncom.2025.1699192
  
DOI=10.3389/fncom.2025.1699192
  
ISSN=1662-5188
  
ABSTRACT=IntroductionProbabilistic Stimulation Maps (PSMs) are increasingly employed to identify brain regions associated with optimal therapeutic outcomes in Deep Brain Stimulation (DBS). However, their reliability and generalizability are challenged by the limited size of most patient cohorts and the inherent variability introduced by different statistical methods and input data configurations. This study aimed to investigate the geometrical variability of Probabilistic Sweet Spots (PSS) as a function of both the number of patients (nPat) and the number of stimulations per patient (nStim), and to model a stability boundary defining the minimum data requirements for obtaining geometrically stable PSS.MethodsThree statistical approaches–Bayesian t-test, Wilcoxon test with False Discovery Rate (FDR) correction, and Wilcoxon test with nonparametric permutation correction–were applied to two patient cohorts: a primary cohort of 36 patients undergoing DBS for Parkinson’s Disease (PD), and a secondary cohort of 61 patients treated for Essential Tremor (ET), used to assess generalizability. Stimulation test data was collected intra-operatively for the first cohort and post-operatively for the second one. Geometric stability was evaluated based on variability in PSS volume extent and centroid location.ResultsThe analysis revealed a non-linear trade-off between nPat and nStim to yield stable PSS. A stability boundary was defined, representing the minimum combinations of nPat–nStim required for anatomically robust PSS. Among the tested methods, the Bayesian t-test achieved stability with smaller sample sizes (∼15 patients) and demonstrated a consistent performance across both cohorts. In contrast, the Wilcoxon-based methods showed variable behavior between cohorts, which differed in symptom type and testing phase (intra-operative testing vs. post-operative screening).DiscussionThe proposed PSS stability boundary provides a practical reference for designing DBS studies and stimulation screening protocols aimed at probabilistic mapping. The Bayesian t-test emerged as a reliable method across both cohorts, supporting its potential in studies with limited sample sizes and scenarios where the method needs to be readily generalized to varying symptoms. These findings underscore the importance of considering both cohort size and stimulation count in probabilistic DBS mapping and call for further investigation into method-specific sensitivities to clinical and procedural factors.