Characterizing cross‐subject spatial interaction patterns in functional magnetic resonance imaging studies: A two‐stage point‐process model

Abstract

We develop a two-stage spatial point process model that introduces new characterizations of activation patterns in multi-subject functional Magnetic Resonance Imaging (fMRI) studies. Conventionally multi-subject fMRI methods rely on combining information across subjects one voxel at a time in order to identify locations of peak activation in the brain. The two-stage model that we develop here addresses shortcomings of standard methods by explicitly modeling the spatial stucture of functional signals and recognizing that corresponding cross-subject functional signals can be spatially misaligned. In our first stage analysis, we introduce a marked spatial point process model that captures the spatial features of the functional response and identifies a configuration of activation units for each subject. The locations of these activation units are used as input for the second stage model. The point process model of the second stage analysis is developed to characterize multi-subject activation patterns by estimating the strength of cross-subject interactions at different spatial ranges. The model uses spatial neighborhoods to account for the cross-subject spatial misalignment in corresponding functional units. We applied our methods to an fMRI study of 21 individuals who performed an attention test. We identified four brain regions that are involved in the test and found that our model results agree well with our understanding of how these regions engage with the tasks performed during the attention test. Our results highlighted that cross-subject interactions are stronger in brain areas that have a more specific function in performing the experimental tasks than in other areas.

Discussion of “Marked Spatial Point Processes: Current State and Extensions to Point Processes on Linear Networks”

“…In the area-interaction process, each point is replaced by a disc and the strength of interaction is measured by the area covered by the discs. Multitype versions of the saturation process have been suggested in Eckel et al (2009), Lee et al (2017) and Rajala et al (2018), and of the area-interaction model e.g. in Picard et al (2009).…”