because they include voxels that do not belong to the average functional cluster defined by the resels. Thus, n = 1 ROIs were used to maximize the number of voxels within the smoothness of the data. Thus, the average BOLD and CM signals in the left and right medial dorsal nuclei of the thalamus (MDTHA; xyz = ±3, −15, 9 mm), cerebellum (CER; culmen; xyz = ±3, −45, −6 mm), and the rostral anterior cingulate cortex (rACC; xyz = ±3, 42, −3 mm) were extracted from the individual BOLD and CM maps in 9×9×9 mm3 ROI volumes containing 27 voxels to using IDL. We selected these ROIs based on our previous fMRI studies [16]–[18], [35]. These ROI masks were created and centered at the precise coordinates listed in Table 2; the coordinates of the ROI masks were kept fixed across subjects and conditions. Note that compared to variable-volume ROIs, fixed-volume (shape and size) ROIs are advantageous in order to minimize statistical confounds. Brain regions were labeled using the Talairach daemon (http://www.talairach.org/) [36] and a query range of 5 mm.
