The neuroimaging community heavily relies on statistical inference to explain measured brain activity given the experimental paradigm. Undeniably, this method has led to many results, but it is limited by the richness of the generative models that are deployed, typically in a mass-univariate way. Such an approach is suboptimal given the high-dimensional and complex spatiotemporal correlation structure of neuroimaging data.