We calculated voxel-by-voxel pairwise crosscorrelations between prewhitened resting-state BOLD Functional Magnetic Resonance ImagingFunctional Magnetic Resonance Imaging (fMRI)A functional neuroimaging procedure using MRI technology that measures brain activity by detecting changes associated with blood flow. This technique relies on the fact that cerebral blood flow and neuronal activation are coupled. When an area of the brain is in use, blood flow to that region also increases.[citation needed] The primary form of fMRI uses the blood-oxygen-level dependent (BOLD) contrast, discovered by Seiji Ogawa. This is a type of specialized brain and body scan used to map neural activity in the brain or spinal cord of humans or other animals by imaging the change in blood flow (hemodynamic response) related to energy use by brain cells. Since the early 1990s, fMRI has come to dominate brain mapping research because it does not require people to undergo shots, surgery, or to ingest substances, or be exposed to ionising radiation, etc. time series recorded from 60 cortical areas (30 per hemisphere) in 18 human subjects (nine women and nine men). Altogether, more than a billion-and-a-quarter pairs of BOLD time series were analyzed. For each pair, a crosscorrelogram was computed by calculating 21 crosscorrelations, namely at zero lag ± 10 lags of 2 s duration each. For each crosscorrelogram, in turn, the crosscorrelation with the highest absolute value was found and its sign, value, and lag were retained for further analysis. In addition, the crosscorrelations at zero lag (irrespective of the location of the peak) were also analyzed as a special case. Based on known varying density of anatomical connectivity, we distinguished four general brain groups for which we derived summary statistics of crosscorrelations between voxels within an area (group I), between voxels of paired homotopic areas across the two hemispheres (group II), between voxels of an area and all other voxels in the same (ipsilateral) hemisphere (group III), and voxels of an area and all voxels in the opposite (contralateral) hemisphere (except those in the homotopic area) (group IV). We found the following. (a) Most of the crosscorrelogram peaks occurred at zero lag, followed by ±1 lag; (b) over all groups, positive crosscorrelations were much more frequent than negative ones; (c) average crosscorrelation was highest for group I, and decreased progressively for groups II-IV; (d) the ratio of positive over negative crosscorrelations was highest for group I and progressively smaller for groups II-IV; (e) the highest proportion of positive crosscorrelations (with respect to all positive ones) was observed at zero lag; and (f) the highest proportion of negative crosscorrelations (with respect to all negative ones) was observed at lag = 2. These findings reveal a systematic pattern of crosscorrelations with respect to their sign, magnitude, lag and brain group, as defined above. Given that these groups were defined along a qualitative gradient of known overall anatomical connectivity, our results suggest that functional interactions between two voxels may simply reflect the density of such anatomical connectivity between the areas to which the voxels belong.