statistical thresholding of source localization maps

hello, everyone,

If I use MNE or eLORETA methods to obtain source estimates for 10 subjects at ERP peak times, such as data of 10 (subject) * 20484 (source vertex count), how should I determine statistical strategies and thresholds to identify statistically significant brain regions in 10 subjects. In the MNE Python example, conditions are difference. How to perform statistical analysis on one condition to determine the threshold?

Thank you.