Thread by @PessoaBrain

How do you test 941 voxels without having to correct for multiple comparisons? 🔥🔥
Build 1⃣ Bayesian multilevel model and estimate 1⃣ multidimensional posterior distribution. No correction needed! (see paper below)🧵

Inspired by work with @gangchen6 on Bayesian multilevel models at the ROI level:
link.springer.com/article/10.100…
🧵

Lab effort to apply the framework at the ROI and in a voxelwise manner (just insula for now) to study the effect of controllability on shock-plus-sound stressors:
biorxiv.org/content/10.110…
Looking forward to running whole-brain soon.🧵

If you like the amygdala and the bed nucleus of the stria terminalis (BST), the paper could be of interest. The BST was very robustly modulated by controllability. In contrast the central amygdala did not show much evidence. The periaqueductal gray (PAG) even less so.🔚

Forgot to mention that Bayesian multilevel posterior estimation possible due to #Stan and #rstan. Amazing tools.