Single-subject data (two runs) in native space

The example shows the analysis of an SPM dataset, with two conditions: viewing a face image or a scrambled face image.

This example takes a lot of time because the input are lists of 3D images sampled in different positions (encoded by different affine functions).

See also

For more information see the dataset description.

Fetch and inspect the data

Fetch the SPM multimodal_faces data.

from nilearn.datasets import fetch_spm_multimodal_fmri

subject_data = fetch_spm_multimodal_fmri()

[wrapper] Dataset created in /home/runner/nilearn_data/spm_multimodal_fmri
[wrapper] Missing 390 functional scans for session 1.
[wrapper] Data absent, downloading...
[wrapper] Downloading data from
https://www.fil.ion.ucl.ac.uk/spm/download/data/mmfaces/multimodal_fmri.zip ...
[wrapper] Downloaded 3022848 of 134263085 bytes (2.3%%,   46.3s remaining)
[wrapper] Downloaded 7061504 of 134263085 bytes (5.3%%,   37.2s remaining)
[wrapper] Downloaded 11124736 of 134263085 bytes (8.3%%,   34.3s remaining)
[wrapper] Downloaded 15130624 of 134263085 bytes (11.3%%,   32.3s remaining)
[wrapper] Downloaded 19382272 of 134263085 bytes (14.4%%,   30.3s remaining)
[wrapper] Downloaded 23797760 of 134263085 bytes (17.7%%,   28.5s remaining)
[wrapper] Downloaded 27975680 of 134263085 bytes (20.8%%,   27.1s remaining)
[wrapper] Downloaded 32129024 of 134263085 bytes (23.9%%,   25.9s remaining)
[wrapper] Downloaded 35495936 of 134263085 bytes (26.4%%,   25.5s remaining)
[wrapper] Downloaded 39124992 of 134263085 bytes (29.1%%,   24.8s remaining)
[wrapper] Downloaded 43204608 of 134263085 bytes (32.2%%,   23.6s remaining)
[wrapper] Downloaded 47431680 of 134263085 bytes (35.3%%,   22.4s remaining)
[wrapper] Downloaded 51585024 of 134263085 bytes (38.4%%,   21.2s remaining)
[wrapper] Downloaded 55951360 of 134263085 bytes (41.7%%,   19.9s remaining)
[wrapper] Downloaded 60112896 of 134263085 bytes (44.8%%,   18.8s remaining)
[wrapper] Downloaded 63332352 of 134263085 bytes (47.2%%,   18.2s remaining)
[wrapper] Downloaded 66961408 of 134263085 bytes (49.9%%,   17.3s remaining)
[wrapper] Downloaded 71024640 of 134263085 bytes (52.9%%,   16.3s remaining)
[wrapper] Downloaded 75268096 of 134263085 bytes (56.1%%,   15.1s remaining)
[wrapper] Downloaded 79437824 of 134263085 bytes (59.2%%,   14.0s remaining)
[wrapper] Downloaded 83779584 of 134263085 bytes (62.4%%,   12.8s remaining)
[wrapper] Downloaded 87654400 of 134263085 bytes (65.3%%,   11.9s remaining)
[wrapper] Downloaded 91013120 of 134263085 bytes (67.8%%,   11.1s remaining)
[wrapper] Downloaded 94863360 of 134263085 bytes (70.7%%,   10.1s remaining)
[wrapper] Downloaded 99024896 of 134263085 bytes (73.8%%,    9.0s remaining)
[wrapper] Downloaded 102973440 of 134263085 bytes (76.7%%,    8.0s remaining)
[wrapper] Downloaded 106520576 of 134263085 bytes (79.3%%,    7.1s remaining)
[wrapper] Downloaded 109846528 of 134263085 bytes (81.8%%,    6.3s remaining)
[wrapper] Downloaded 113369088 of 134263085 bytes (84.4%%,    5.4s remaining)
[wrapper] Downloaded 116465664 of 134263085 bytes (86.7%%,    4.7s remaining)
[wrapper] Downloaded 119218176 of 134263085 bytes (88.8%%,    4.0s remaining)
[wrapper] Downloaded 122200064 of 134263085 bytes (91.0%%,    3.2s remaining)
[wrapper] Downloaded 125263872 of 134263085 bytes (93.3%%,    2.4s remaining)
[wrapper] Downloaded 128294912 of 134263085 bytes (95.6%%,    1.6s remaining)
[wrapper] Downloaded 130752512 of 134263085 bytes (97.4%%,    1.0s remaining)
[wrapper] Downloaded 133398528 of 134263085 bytes (99.4%%,    0.2s remaining)
[wrapper]  ...done. (38 seconds, 0 min)

[wrapper] Extracting data from
/home/runner/nilearn_data/spm_multimodal_fmri/sub001/multimodal_fmri.zip...
[wrapper] .. done.

[wrapper] Downloading data from
https://www.fil.ion.ucl.ac.uk/spm/download/data/mmfaces/multimodal_smri.zip ...
[wrapper] Downloaded 3547136 of 6852766 bytes (51.8%%,    1.0s remaining)
[wrapper]  ...done. (3 seconds, 0 min)

[wrapper] Extracting data from
/home/runner/nilearn_data/spm_multimodal_fmri/sub001/multimodal_smri.zip...
[wrapper] .. done.

Let’s inspect one of the event files before using them.

import pandas as pd

events = [subject_data.events1, subject_data.events2]

events_dataframe = pd.read_csv(events[0], sep="\t")
events_dataframe["trial_type"].value_counts()

trial_type
scrambled    86
faces        64
Name: count, dtype: int64

We can confirm there are only 2 conditions in the dataset.

from nilearn.plotting import plot_event, show

plot_event(events)

show()

Resample the images: this is achieved by the concat_imgs function of Nilearn.

import warnings

from nilearn.image import concat_imgs, mean_img, resample_img

# Avoid getting too many warnings due to resampling
with warnings.catch_warnings():
    warnings.simplefilter("ignore")
    fmri_img = [
        concat_imgs(subject_data.func1, auto_resample=True),
        concat_imgs(subject_data.func2, auto_resample=True),
    ]
affine, shape = fmri_img[0].affine, fmri_img[0].shape
print("Resampling the second image (this takes time)...")
fmri_img[1] = resample_img(
    fmri_img[1], affine, shape[:3], copy_header=True, force_resample=True
)

Resampling the second image (this takes time)...

Let’s create mean image for display purposes.

mean_image = mean_img(fmri_img, copy_header=True)

Fit the model

Fit the GLM for the 2 runs by specifying a FirstLevelModel and then fitting it.

# Sample at the beginning of each acquisition.
slice_time_ref = 0.0
# We use a discrete cosine transform to model signal drifts.
drift_model = "cosine"
# The cutoff for the drift model is 0.01 Hz.
high_pass = 0.01
# The hemodynamic response function
hrf_model = "spm + derivative"

from nilearn.glm.first_level import FirstLevelModel

print("Fitting a GLM")
fmri_glm = FirstLevelModel(
    smoothing_fwhm=None,
    t_r=subject_data.t_r,
    hrf_model=hrf_model,
    drift_model=drift_model,
    high_pass=high_pass,
)


fmri_glm = fmri_glm.fit(fmri_img, events=events)

Fitting a GLM

View the results

Now we can compute contrast-related statistical maps (in z-scale), and plot them.

from nilearn.plotting import plot_stat_map

print("Computing contrasts")

Computing contrasts

We actually want more interesting contrasts. The simplest contrast just makes the difference between the two main conditions. We define the two opposite versions to run one-tailed t-tests.

contrasts = ["faces - scrambled", "scrambled - faces"]

Let’s store common parameters for all plots.

We plot the contrasts values overlaid on the mean fMRI image and we will use the z-score values as transparency, with any voxel with | Z-score | > 3 being fully opaque and any voxel with 0 < | Z-score | < 1.96 being partly transparent.

plot_param = {
    "vmin": 0,
    "display_mode": "z",
    "cut_coords": 3,
    "black_bg": True,
    "bg_img": mean_image,
    "cmap": "inferno",
    "transparency_range": [0, 3],
}

# Iterate on contrasts to compute and plot them.
for contrast_id in contrasts:
    print(f"\tcontrast id: {contrast_id}")

    results = fmri_glm.compute_contrast(contrast_id, output_type="all")

    plot_stat_map(
        results["stat"],
        title=contrast_id,
        transparency=results["z_score"],
        **plot_param,
    )

• 
• 

        contrast id: faces - scrambled
/home/runner/work/nilearn/nilearn/examples/04_glm_first_level/plot_spm_multimodal_faces.py:135: UserWarning: One contrast given, assuming it for all 2 runs
  results = fmri_glm.compute_contrast(contrast_id, output_type="all")
        contrast id: scrambled - faces
/home/runner/work/nilearn/nilearn/examples/04_glm_first_level/plot_spm_multimodal_faces.py:135: UserWarning: One contrast given, assuming it for all 2 runs
  results = fmri_glm.compute_contrast(contrast_id, output_type="all")

We also define the effects of interest contrast, a 2-dimensional contrasts spanning the two conditions.

import numpy as np

contrasts = np.eye(2)

results = fmri_glm.compute_contrast(contrasts, output_type="all")

plot_stat_map(
    results["stat"],
    title="effects of interest",
    transparency=results["z_score"],
    **plot_param,
)

show()

/home/runner/work/nilearn/nilearn/examples/04_glm_first_level/plot_spm_multimodal_faces.py:152: UserWarning: One contrast given, assuming it for all 2 runs
  results = fmri_glm.compute_contrast(contrasts, output_type="all")
/home/runner/work/nilearn/nilearn/examples/04_glm_first_level/plot_spm_multimodal_faces.py:152: UserWarning: F contrasts should have 20 colmuns, but it has only 2. The rest of the contrast was padded with zeros.
  results = fmri_glm.compute_contrast(contrasts, output_type="all")
/home/runner/work/nilearn/nilearn/examples/04_glm_first_level/plot_spm_multimodal_faces.py:152: UserWarning: Running approximate fixed effects on F statistics.
  results = fmri_glm.compute_contrast(contrasts, output_type="all")

Based on the resulting maps we observe that the analysis results in wide activity for the ‘effects of interest’ contrast, showing the implications of large portions of the visual cortex in the conditions. By contrast, the differential effect between “faces” and “scrambled” involves sparser, more anterior and lateral regions. It also displays some responses in the frontal lobe.

Estimated memory usage: 976 MB