Note
Go to the end to download the full example code
Decoding source space data#
Decoding to MEG data in source space on the left cortical surface. Here univariate feature selection is employed for speed purposes to confine the classification to a small number of potentially relevant features. The classifier then is trained to selected features of epochs in source space.
# Author: Denis A. Engemann <denis.engemann@gmail.com>
# Alexandre Gramfort <alexandre.gramfort@inria.fr>
# Jean-Remi King <jeanremi.king@gmail.com>
# Eric Larson <larson.eric.d@gmail.com>
#
# License: BSD-3-Clause
import numpy as np
import matplotlib.pyplot as plt
from sklearn.pipeline import make_pipeline
from sklearn.preprocessing import StandardScaler
from sklearn.feature_selection import SelectKBest, f_classif
from sklearn.linear_model import LogisticRegression
import mne
from mne.minimum_norm import apply_inverse_epochs, read_inverse_operator
from mne.decoding import (cross_val_multiscore, LinearModel, SlidingEstimator,
get_coef)
print(__doc__)
data_path = mne.datasets.sample.data_path()
meg_path = data_path / 'MEG' / 'sample'
fname_fwd = meg_path / 'sample_audvis-meg-oct-6-fwd.fif'
fname_evoked = meg_path / 'sample_audvis-ave.fif'
subjects_dir = data_path / 'subjects'
Set parameters
raw_fname = meg_path / 'sample_audvis_filt-0-40_raw.fif'
event_fname = meg_path / 'sample_audvis_filt-0-40_raw-eve.fif'
fname_cov = meg_path / 'sample_audvis-cov.fif'
fname_inv = meg_path / 'sample_audvis-meg-oct-6-meg-inv.fif'
tmin, tmax = -0.2, 0.8
event_id = dict(aud_r=2, vis_r=4) # load contra-lateral conditions
# Setup for reading the raw data
raw = mne.io.read_raw_fif(raw_fname, preload=True)
raw.filter(None, 10., fir_design='firwin')
events = mne.read_events(event_fname)
# Set up pick list: MEG - bad channels (modify to your needs)
raw.info['bads'] += ['MEG 2443'] # mark bads
picks = mne.pick_types(raw.info, meg=True, eeg=False, stim=True, eog=True,
exclude='bads')
# Read epochs
epochs = mne.Epochs(raw, events, event_id, tmin, tmax, proj=True,
picks=picks, baseline=(None, 0), preload=True,
reject=dict(grad=4000e-13, eog=150e-6),
decim=5) # decimate to save memory and increase speed
Opening raw data file /home/circleci/mne_data/MNE-sample-data/MEG/sample/sample_audvis_filt-0-40_raw.fif...
Read a total of 4 projection items:
PCA-v1 (1 x 102) idle
PCA-v2 (1 x 102) idle
PCA-v3 (1 x 102) idle
Average EEG reference (1 x 60) idle
Range : 6450 ... 48149 = 42.956 ... 320.665 secs
Ready.
Reading 0 ... 41699 = 0.000 ... 277.709 secs...
Filtering raw data in 1 contiguous segment
Setting up low-pass filter at 10 Hz
FIR filter parameters
---------------------
Designing a one-pass, zero-phase, non-causal lowpass filter:
- Windowed time-domain design (firwin) method
- Hamming window with 0.0194 passband ripple and 53 dB stopband attenuation
- Upper passband edge: 10.00 Hz
- Upper transition bandwidth: 2.50 Hz (-6 dB cutoff frequency: 11.25 Hz)
- Filter length: 199 samples (1.325 sec)
[Parallel(n_jobs=1)]: Using backend SequentialBackend with 1 concurrent workers.
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 0.0s remaining: 0.0s
[Parallel(n_jobs=1)]: Done 2 out of 2 | elapsed: 0.0s remaining: 0.0s
[Parallel(n_jobs=1)]: Done 3 out of 3 | elapsed: 0.0s remaining: 0.0s
[Parallel(n_jobs=1)]: Done 4 out of 4 | elapsed: 0.0s remaining: 0.0s
[Parallel(n_jobs=1)]: Done 366 out of 366 | elapsed: 0.9s finished
Not setting metadata
143 matching events found
Setting baseline interval to [-0.19979521315838786, 0.0] sec
Applying baseline correction (mode: mean)
Created an SSP operator (subspace dimension = 3)
4 projection items activated
Using data from preloaded Raw for 143 events and 151 original time points (prior to decimation) ...
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
Rejecting epoch based on EOG : ['EOG 061']
28 bad epochs dropped
Compute inverse solution
snr = 3.0
noise_cov = mne.read_cov(fname_cov)
inverse_operator = read_inverse_operator(fname_inv)
stcs = apply_inverse_epochs(epochs, inverse_operator,
lambda2=1.0 / snr ** 2, verbose=False,
method="dSPM", pick_ori="normal")
366 x 366 full covariance (kind = 1) found.
Read a total of 4 projection items:
PCA-v1 (1 x 102) active
PCA-v2 (1 x 102) active
PCA-v3 (1 x 102) active
Average EEG reference (1 x 60) active
Reading inverse operator decomposition from /home/circleci/mne_data/MNE-sample-data/MEG/sample/sample_audvis-meg-oct-6-meg-inv.fif...
Reading inverse operator info...
[done]
Reading inverse operator decomposition...
[done]
305 x 305 full covariance (kind = 1) found.
Read a total of 4 projection items:
PCA-v1 (1 x 102) active
PCA-v2 (1 x 102) active
PCA-v3 (1 x 102) active
Average EEG reference (1 x 60) active
Noise covariance matrix read.
22494 x 22494 diagonal covariance (kind = 2) found.
Source covariance matrix read.
22494 x 22494 diagonal covariance (kind = 6) found.
Orientation priors read.
22494 x 22494 diagonal covariance (kind = 5) found.
Depth priors read.
Did not find the desired covariance matrix (kind = 3)
Reading a source space...
Computing patch statistics...
Patch information added...
Distance information added...
[done]
Reading a source space...
Computing patch statistics...
Patch information added...
Distance information added...
[done]
2 source spaces read
Read a total of 4 projection items:
PCA-v1 (1 x 102) active
PCA-v2 (1 x 102) active
PCA-v3 (1 x 102) active
Average EEG reference (1 x 60) active
Source spaces transformed to the inverse solution coordinate frame
Decoding in sensor space using a logistic regression
# Retrieve source space data into an array
X = np.array([stc.lh_data for stc in stcs]) # only keep left hemisphere
y = epochs.events[:, 2]
# prepare a series of classifier applied at each time sample
clf = make_pipeline(StandardScaler(), # z-score normalization
SelectKBest(f_classif, k=500), # select features for speed
LinearModel(LogisticRegression(C=1, solver='liblinear')))
time_decod = SlidingEstimator(clf, scoring='roc_auc')
# Run cross-validated decoding analyses:
scores = cross_val_multiscore(time_decod, X, y, cv=5, n_jobs=None)
# Plot average decoding scores of 5 splits
fig, ax = plt.subplots(1)
ax.plot(epochs.times, scores.mean(0), label='score')
ax.axhline(.5, color='k', linestyle='--', label='chance')
ax.axvline(0, color='k')
plt.legend()
[Parallel(n_jobs=1)]: Using backend SequentialBackend with 1 concurrent workers.
0%| | Fitting SlidingEstimator : 0/31 [00:00<?, ?it/s]
3%|3 | Fitting SlidingEstimator : 1/31 [00:00<00:02, 14.71it/s]
6%|6 | Fitting SlidingEstimator : 2/31 [00:00<00:01, 14.77it/s]
10%|9 | Fitting SlidingEstimator : 3/31 [00:00<00:01, 17.88it/s]
13%|#2 | Fitting SlidingEstimator : 4/31 [00:00<00:01, 16.94it/s]
16%|#6 | Fitting SlidingEstimator : 5/31 [00:00<00:01, 18.70it/s]
19%|#9 | Fitting SlidingEstimator : 6/31 [00:00<00:01, 17.83it/s]
23%|##2 | Fitting SlidingEstimator : 7/31 [00:00<00:01, 19.09it/s]
26%|##5 | Fitting SlidingEstimator : 8/31 [00:00<00:01, 18.31it/s]
29%|##9 | Fitting SlidingEstimator : 9/31 [00:00<00:01, 19.30it/s]
32%|###2 | Fitting SlidingEstimator : 10/31 [00:00<00:01, 20.17it/s]
35%|###5 | Fitting SlidingEstimator : 11/31 [00:00<00:00, 20.94it/s]
39%|###8 | Fitting SlidingEstimator : 12/31 [00:00<00:00, 20.04it/s]
42%|####1 | Fitting SlidingEstimator : 13/31 [00:00<00:00, 20.73it/s]
45%|####5 | Fitting SlidingEstimator : 14/31 [00:00<00:00, 21.35it/s]
48%|####8 | Fitting SlidingEstimator : 15/31 [00:00<00:00, 21.92it/s]
52%|#####1 | Fitting SlidingEstimator : 16/31 [00:00<00:00, 21.02it/s]
55%|#####4 | Fitting SlidingEstimator : 17/31 [00:00<00:00, 21.55it/s]
58%|#####8 | Fitting SlidingEstimator : 18/31 [00:00<00:00, 22.05it/s]
61%|######1 | Fitting SlidingEstimator : 19/31 [00:00<00:00, 21.22it/s]
65%|######4 | Fitting SlidingEstimator : 20/31 [00:00<00:00, 21.70it/s]
68%|######7 | Fitting SlidingEstimator : 21/31 [00:01<00:00, 20.97it/s]
71%|####### | Fitting SlidingEstimator : 22/31 [00:01<00:00, 21.43it/s]
74%|#######4 | Fitting SlidingEstimator : 23/31 [00:01<00:00, 21.86it/s]
77%|#######7 | Fitting SlidingEstimator : 24/31 [00:01<00:00, 21.16it/s]
81%|######## | Fitting SlidingEstimator : 25/31 [00:01<00:00, 21.58it/s]
84%|########3 | Fitting SlidingEstimator : 26/31 [00:01<00:00, 20.94it/s]
87%|########7 | Fitting SlidingEstimator : 27/31 [00:01<00:00, 21.35it/s]
90%|######### | Fitting SlidingEstimator : 28/31 [00:01<00:00, 20.75it/s]
94%|#########3| Fitting SlidingEstimator : 29/31 [00:01<00:00, 21.16it/s]
97%|#########6| Fitting SlidingEstimator : 30/31 [00:01<00:00, 21.55it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 21.64it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 21.19it/s]
[Parallel(n_jobs=1)]: Using backend SequentialBackend with 1 concurrent workers.
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 0.1s remaining: 0.0s
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 0.1s finished
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 1.7s remaining: 0.0s
0%| | Fitting SlidingEstimator : 0/31 [00:00<?, ?it/s]
3%|3 | Fitting SlidingEstimator : 1/31 [00:00<00:02, 14.83it/s]
6%|6 | Fitting SlidingEstimator : 2/31 [00:00<00:01, 19.91it/s]
10%|9 | Fitting SlidingEstimator : 3/31 [00:00<00:01, 17.77it/s]
13%|#2 | Fitting SlidingEstimator : 4/31 [00:00<00:01, 19.92it/s]
16%|#6 | Fitting SlidingEstimator : 5/31 [00:00<00:01, 18.52it/s]
19%|#9 | Fitting SlidingEstimator : 6/31 [00:00<00:01, 19.92it/s]
23%|##2 | Fitting SlidingEstimator : 7/31 [00:00<00:01, 18.85it/s]
26%|##5 | Fitting SlidingEstimator : 8/31 [00:00<00:01, 19.92it/s]
29%|##9 | Fitting SlidingEstimator : 9/31 [00:00<00:01, 19.04it/s]
32%|###2 | Fitting SlidingEstimator : 10/31 [00:00<00:01, 19.92it/s]
35%|###5 | Fitting SlidingEstimator : 11/31 [00:00<00:01, 19.16it/s]
39%|###8 | Fitting SlidingEstimator : 12/31 [00:00<00:00, 19.92it/s]
42%|####1 | Fitting SlidingEstimator : 13/31 [00:00<00:00, 20.61it/s]
45%|####5 | Fitting SlidingEstimator : 14/31 [00:00<00:00, 21.24it/s]
48%|####8 | Fitting SlidingEstimator : 15/31 [00:00<00:00, 20.42it/s]
52%|#####1 | Fitting SlidingEstimator : 16/31 [00:00<00:00, 21.00it/s]
55%|#####4 | Fitting SlidingEstimator : 17/31 [00:00<00:00, 21.53it/s]
58%|#####8 | Fitting SlidingEstimator : 18/31 [00:00<00:00, 20.76it/s]
61%|######1 | Fitting SlidingEstimator : 19/31 [00:00<00:00, 21.26it/s]
65%|######4 | Fitting SlidingEstimator : 20/31 [00:00<00:00, 21.73it/s]
68%|######7 | Fitting SlidingEstimator : 21/31 [00:01<00:00, 21.00it/s]
71%|####### | Fitting SlidingEstimator : 22/31 [00:01<00:00, 21.46it/s]
74%|#######4 | Fitting SlidingEstimator : 23/31 [00:01<00:00, 21.89it/s]
77%|#######7 | Fitting SlidingEstimator : 24/31 [00:01<00:00, 22.30it/s]
81%|######## | Fitting SlidingEstimator : 25/31 [00:01<00:00, 21.55it/s]
84%|########3 | Fitting SlidingEstimator : 26/31 [00:01<00:00, 21.96it/s]
87%|########7 | Fitting SlidingEstimator : 27/31 [00:01<00:00, 22.34it/s]
90%|######### | Fitting SlidingEstimator : 28/31 [00:01<00:00, 21.62it/s]
94%|#########3| Fitting SlidingEstimator : 29/31 [00:01<00:00, 22.00it/s]
97%|#########6| Fitting SlidingEstimator : 30/31 [00:01<00:00, 22.37it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 22.50it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 21.75it/s]
[Parallel(n_jobs=1)]: Using backend SequentialBackend with 1 concurrent workers.
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 0.1s remaining: 0.0s
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 0.1s finished
[Parallel(n_jobs=1)]: Done 2 out of 2 | elapsed: 3.3s remaining: 0.0s
0%| | Fitting SlidingEstimator : 0/31 [00:00<?, ?it/s]
3%|3 | Fitting SlidingEstimator : 1/31 [00:00<00:02, 14.83it/s]
6%|6 | Fitting SlidingEstimator : 2/31 [00:00<00:01, 14.85it/s]
10%|9 | Fitting SlidingEstimator : 3/31 [00:00<00:01, 17.98it/s]
13%|#2 | Fitting SlidingEstimator : 4/31 [00:00<00:01, 20.09it/s]
16%|#6 | Fitting SlidingEstimator : 5/31 [00:00<00:01, 18.64it/s]
19%|#9 | Fitting SlidingEstimator : 6/31 [00:00<00:01, 20.03it/s]
23%|##2 | Fitting SlidingEstimator : 7/31 [00:00<00:01, 18.93it/s]
26%|##5 | Fitting SlidingEstimator : 8/31 [00:00<00:01, 20.00it/s]
29%|##9 | Fitting SlidingEstimator : 9/31 [00:00<00:01, 20.91it/s]
32%|###2 | Fitting SlidingEstimator : 10/31 [00:00<00:00, 21.70it/s]
35%|###5 | Fitting SlidingEstimator : 11/31 [00:00<00:00, 22.39it/s]
39%|###8 | Fitting SlidingEstimator : 12/31 [00:00<00:00, 21.22it/s]
42%|####1 | Fitting SlidingEstimator : 13/31 [00:00<00:00, 21.85it/s]
45%|####5 | Fitting SlidingEstimator : 14/31 [00:00<00:00, 22.42it/s]
48%|####8 | Fitting SlidingEstimator : 15/31 [00:00<00:00, 22.93it/s]
52%|#####1 | Fitting SlidingEstimator : 16/31 [00:00<00:00, 21.86it/s]
55%|#####4 | Fitting SlidingEstimator : 17/31 [00:00<00:00, 22.35it/s]
58%|#####8 | Fitting SlidingEstimator : 18/31 [00:00<00:00, 22.81it/s]
61%|######1 | Fitting SlidingEstimator : 19/31 [00:00<00:00, 21.87it/s]
65%|######4 | Fitting SlidingEstimator : 20/31 [00:00<00:00, 22.32it/s]
68%|######7 | Fitting SlidingEstimator : 21/31 [00:00<00:00, 21.50it/s]
71%|####### | Fitting SlidingEstimator : 22/31 [00:01<00:00, 21.94it/s]
74%|#######4 | Fitting SlidingEstimator : 23/31 [00:01<00:00, 21.14it/s]
77%|#######7 | Fitting SlidingEstimator : 24/31 [00:01<00:00, 21.57it/s]
81%|######## | Fitting SlidingEstimator : 25/31 [00:01<00:00, 20.92it/s]
84%|########3 | Fitting SlidingEstimator : 26/31 [00:01<00:00, 21.34it/s]
87%|########7 | Fitting SlidingEstimator : 27/31 [00:01<00:00, 20.73it/s]
90%|######### | Fitting SlidingEstimator : 28/31 [00:01<00:00, 21.15it/s]
94%|#########3| Fitting SlidingEstimator : 29/31 [00:01<00:00, 20.59it/s]
97%|#########6| Fitting SlidingEstimator : 30/31 [00:01<00:00, 20.99it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 21.04it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 21.14it/s]
[Parallel(n_jobs=1)]: Using backend SequentialBackend with 1 concurrent workers.
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 0.1s remaining: 0.0s
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 0.1s finished
[Parallel(n_jobs=1)]: Done 3 out of 3 | elapsed: 5.0s remaining: 0.0s
0%| | Fitting SlidingEstimator : 0/31 [00:00<?, ?it/s]
3%|3 | Fitting SlidingEstimator : 1/31 [00:00<00:02, 14.83it/s]
6%|6 | Fitting SlidingEstimator : 2/31 [00:00<00:01, 19.92it/s]
10%|9 | Fitting SlidingEstimator : 3/31 [00:00<00:01, 17.79it/s]
13%|#2 | Fitting SlidingEstimator : 4/31 [00:00<00:01, 19.92it/s]
16%|#6 | Fitting SlidingEstimator : 5/31 [00:00<00:01, 21.46it/s]
19%|#9 | Fitting SlidingEstimator : 6/31 [00:00<00:01, 19.80it/s]
23%|##2 | Fitting SlidingEstimator : 7/31 [00:00<00:01, 20.95it/s]
26%|##5 | Fitting SlidingEstimator : 8/31 [00:00<00:01, 21.90it/s]
29%|##9 | Fitting SlidingEstimator : 9/31 [00:00<00:00, 22.69it/s]
32%|###2 | Fitting SlidingEstimator : 10/31 [00:00<00:00, 21.30it/s]
35%|###5 | Fitting SlidingEstimator : 11/31 [00:00<00:00, 22.01it/s]
39%|###8 | Fitting SlidingEstimator : 12/31 [00:00<00:00, 22.65it/s]
42%|####1 | Fitting SlidingEstimator : 13/31 [00:00<00:00, 23.21it/s]
45%|####5 | Fitting SlidingEstimator : 14/31 [00:00<00:00, 23.70it/s]
48%|####8 | Fitting SlidingEstimator : 15/31 [00:00<00:00, 22.45it/s]
52%|#####1 | Fitting SlidingEstimator : 16/31 [00:00<00:00, 22.94it/s]
55%|#####4 | Fitting SlidingEstimator : 17/31 [00:00<00:00, 23.39it/s]
58%|#####8 | Fitting SlidingEstimator : 18/31 [00:00<00:00, 22.33it/s]
61%|######1 | Fitting SlidingEstimator : 19/31 [00:00<00:00, 22.78it/s]
65%|######4 | Fitting SlidingEstimator : 20/31 [00:00<00:00, 21.87it/s]
68%|######7 | Fitting SlidingEstimator : 21/31 [00:00<00:00, 22.30it/s]
71%|####### | Fitting SlidingEstimator : 22/31 [00:00<00:00, 22.72it/s]
74%|#######4 | Fitting SlidingEstimator : 23/31 [00:01<00:00, 23.09it/s]
77%|#######7 | Fitting SlidingEstimator : 24/31 [00:01<00:00, 22.22it/s]
81%|######## | Fitting SlidingEstimator : 25/31 [00:01<00:00, 22.61it/s]
84%|########3 | Fitting SlidingEstimator : 26/31 [00:01<00:00, 22.98it/s]
87%|########7 | Fitting SlidingEstimator : 27/31 [00:01<00:00, 22.17it/s]
90%|######### | Fitting SlidingEstimator : 28/31 [00:01<00:00, 21.47it/s]
94%|#########3| Fitting SlidingEstimator : 29/31 [00:01<00:00, 21.85it/s]
97%|#########6| Fitting SlidingEstimator : 30/31 [00:01<00:00, 21.21it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 22.05it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 22.09it/s]
[Parallel(n_jobs=1)]: Using backend SequentialBackend with 1 concurrent workers.
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 0.1s remaining: 0.0s
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 0.1s finished
[Parallel(n_jobs=1)]: Done 4 out of 4 | elapsed: 6.6s remaining: 0.0s
0%| | Fitting SlidingEstimator : 0/31 [00:00<?, ?it/s]
3%|3 | Fitting SlidingEstimator : 1/31 [00:00<00:02, 14.82it/s]
6%|6 | Fitting SlidingEstimator : 2/31 [00:00<00:01, 14.82it/s]
10%|9 | Fitting SlidingEstimator : 3/31 [00:00<00:01, 17.96it/s]
13%|#2 | Fitting SlidingEstimator : 4/31 [00:00<00:01, 16.99it/s]
16%|#6 | Fitting SlidingEstimator : 5/31 [00:00<00:01, 18.74it/s]
19%|#9 | Fitting SlidingEstimator : 6/31 [00:00<00:01, 17.85it/s]
23%|##2 | Fitting SlidingEstimator : 7/31 [00:00<00:01, 19.10it/s]
26%|##5 | Fitting SlidingEstimator : 8/31 [00:00<00:01, 18.32it/s]
29%|##9 | Fitting SlidingEstimator : 9/31 [00:00<00:01, 19.31it/s]
32%|###2 | Fitting SlidingEstimator : 10/31 [00:00<00:01, 20.18it/s]
35%|###5 | Fitting SlidingEstimator : 11/31 [00:00<00:01, 19.37it/s]
39%|###8 | Fitting SlidingEstimator : 12/31 [00:00<00:00, 20.13it/s]
42%|####1 | Fitting SlidingEstimator : 13/31 [00:00<00:00, 20.69it/s]
45%|####5 | Fitting SlidingEstimator : 14/31 [00:00<00:00, 21.31it/s]
48%|####8 | Fitting SlidingEstimator : 15/31 [00:00<00:00, 20.48it/s]
52%|#####1 | Fitting SlidingEstimator : 16/31 [00:00<00:00, 21.06it/s]
55%|#####4 | Fitting SlidingEstimator : 17/31 [00:00<00:00, 20.32it/s]
58%|#####8 | Fitting SlidingEstimator : 18/31 [00:00<00:00, 20.86it/s]
61%|######1 | Fitting SlidingEstimator : 19/31 [00:00<00:00, 21.37it/s]
65%|######4 | Fitting SlidingEstimator : 20/31 [00:00<00:00, 20.66it/s]
68%|######7 | Fitting SlidingEstimator : 21/31 [00:01<00:00, 20.06it/s]
71%|####### | Fitting SlidingEstimator : 22/31 [00:01<00:00, 20.55it/s]
74%|#######4 | Fitting SlidingEstimator : 23/31 [00:01<00:00, 21.01it/s]
77%|#######7 | Fitting SlidingEstimator : 24/31 [00:01<00:00, 20.41it/s]
81%|######## | Fitting SlidingEstimator : 25/31 [00:01<00:00, 20.85it/s]
84%|########3 | Fitting SlidingEstimator : 26/31 [00:01<00:00, 20.30it/s]
87%|########7 | Fitting SlidingEstimator : 27/31 [00:01<00:00, 19.81it/s]
90%|######### | Fitting SlidingEstimator : 28/31 [00:01<00:00, 20.25it/s]
94%|#########3| Fitting SlidingEstimator : 29/31 [00:01<00:00, 19.78it/s]
97%|#########6| Fitting SlidingEstimator : 30/31 [00:01<00:00, 20.20it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 20.31it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 20.20it/s]
[Parallel(n_jobs=1)]: Using backend SequentialBackend with 1 concurrent workers.
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 0.1s remaining: 0.0s
[Parallel(n_jobs=1)]: Done 1 out of 1 | elapsed: 0.1s finished
[Parallel(n_jobs=1)]: Done 5 out of 5 | elapsed: 8.3s finished
To investigate weights, we need to retrieve the patterns of a fitted model
# The fitting needs not be cross validated because the weights are based on
# the training sets
time_decod.fit(X, y)
# Retrieve patterns after inversing the z-score normalization step:
patterns = get_coef(time_decod, 'patterns_', inverse_transform=True)
stc = stcs[0] # for convenience, lookup parameters from first stc
vertices = [stc.lh_vertno, np.array([], int)] # empty array for right hemi
stc_feat = mne.SourceEstimate(np.abs(patterns), vertices=vertices,
tmin=stc.tmin, tstep=stc.tstep, subject='sample')
brain = stc_feat.plot(views=['lat'], transparent=True,
initial_time=0.1, time_unit='s',
subjects_dir=subjects_dir)
0%| | Fitting SlidingEstimator : 0/31 [00:00<?, ?it/s]
3%|3 | Fitting SlidingEstimator : 1/31 [00:00<00:03, 9.85it/s]
6%|6 | Fitting SlidingEstimator : 2/31 [00:00<00:02, 11.89it/s]
10%|9 | Fitting SlidingEstimator : 3/31 [00:00<00:02, 12.77it/s]
13%|#2 | Fitting SlidingEstimator : 4/31 [00:00<00:02, 13.27it/s]
16%|#6 | Fitting SlidingEstimator : 5/31 [00:00<00:01, 13.59it/s]
19%|#9 | Fitting SlidingEstimator : 6/31 [00:00<00:01, 13.80it/s]
23%|##2 | Fitting SlidingEstimator : 7/31 [00:00<00:01, 13.96it/s]
26%|##5 | Fitting SlidingEstimator : 8/31 [00:00<00:01, 14.09it/s]
29%|##9 | Fitting SlidingEstimator : 9/31 [00:00<00:01, 14.19it/s]
32%|###2 | Fitting SlidingEstimator : 10/31 [00:00<00:01, 15.17it/s]
35%|###5 | Fitting SlidingEstimator : 11/31 [00:00<00:01, 15.12it/s]
39%|###8 | Fitting SlidingEstimator : 12/31 [00:00<00:01, 15.09it/s]
42%|####1 | Fitting SlidingEstimator : 13/31 [00:00<00:01, 15.07it/s]
45%|####5 | Fitting SlidingEstimator : 14/31 [00:00<00:01, 15.83it/s]
48%|####8 | Fitting SlidingEstimator : 15/31 [00:00<00:01, 15.73it/s]
52%|#####1 | Fitting SlidingEstimator : 16/31 [00:01<00:00, 15.64it/s]
55%|#####4 | Fitting SlidingEstimator : 17/31 [00:01<00:00, 15.56it/s]
58%|#####8 | Fitting SlidingEstimator : 18/31 [00:01<00:00, 15.47it/s]
61%|######1 | Fitting SlidingEstimator : 19/31 [00:01<00:00, 16.08it/s]
65%|######4 | Fitting SlidingEstimator : 20/31 [00:01<00:00, 15.98it/s]
68%|######7 | Fitting SlidingEstimator : 21/31 [00:01<00:00, 15.88it/s]
71%|####### | Fitting SlidingEstimator : 22/31 [00:01<00:00, 15.80it/s]
74%|#######4 | Fitting SlidingEstimator : 23/31 [00:01<00:00, 16.35it/s]
77%|#######7 | Fitting SlidingEstimator : 24/31 [00:01<00:00, 16.23it/s]
81%|######## | Fitting SlidingEstimator : 25/31 [00:01<00:00, 16.12it/s]
84%|########3 | Fitting SlidingEstimator : 26/31 [00:01<00:00, 16.64it/s]
87%|########7 | Fitting SlidingEstimator : 27/31 [00:01<00:00, 16.50it/s]
90%|######### | Fitting SlidingEstimator : 28/31 [00:01<00:00, 16.37it/s]
94%|#########3| Fitting SlidingEstimator : 29/31 [00:01<00:00, 16.26it/s]
97%|#########6| Fitting SlidingEstimator : 30/31 [00:01<00:00, 15.62it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 15.87it/s]
100%|##########| Fitting SlidingEstimator : 31/31 [00:01<00:00, 15.70it/s]
Using control points [ 3.50724321 4.37721471 15.59502306]
Total running time of the script: ( 0 minutes 19.864 seconds)
Estimated memory usage: 384 MB