Note
Go to the end to download the full example code.
Embedding ERP MEG data in 2D Euclidean space¶
Riemannian embeddings via Laplacian Eigenmaps (LE) and Locally Linear Embedding (LLE) of a set of ERP data. Embedding via Laplacian Eigenmaps is referred to as Spectral Embedding (SE).
Locally Linear Embedding (LLE) assumes that the local neighborhood of a point on the manifold can be well approximated by the affine subspace spanned by the k-nearest neighbors of the point and finds a low-dimensional embedding of the data based on these affine approximations.
Laplacian Eigenmaps (LE) are based on computing the low dimensional representation that best preserves locality instead of local linearity in LLE [1].
# Authors: Pedro Rodrigues <pedro.rodrigues01@gmail.com>,
# Gabriel Wagner vom Berg <gabriel@bccn-berlin.de>
# License: BSD (3-clause)
import matplotlib.pyplot as plt
import mne
from mne import io
from mne.datasets import sample
from sklearn.model_selection import train_test_split
from pyriemann.estimation import XdawnCovariances
from pyriemann.utils.viz import plot_embedding
print(__doc__)
Set parameters and read data¶
data_path = str(sample.data_path())
raw_fname = data_path + "/MEG/sample/sample_audvis_filt-0-40_raw.fif"
event_fname = data_path + "/MEG/sample/sample_audvis_filt-0-40_raw-eve.fif"
tmin, tmax = -0., 1
event_id = dict(aud_l=1, aud_r=2, vis_l=3, vis_r=4)
# Setup for reading the raw data
raw = io.Raw(raw_fname, preload=True, verbose=False)
raw.filter(2, None, method="iir") # replace baselining with high-pass
events = mne.read_events(event_fname)
raw.info["bads"] = ["MEG 2443"] # set bad channels
picks = mne.pick_types(raw.info, meg=True, eeg=False, stim=False, eog=False,
exclude="bads")
# Read epochs
epochs = mne.Epochs(raw, events, event_id, tmin, tmax, proj=False,
picks=picks, baseline=None, preload=True, verbose=False)
X = epochs.get_data(copy=False)
y = epochs.events[:, -1]
Using default location ~/mne_data for sample...
0%| | 0.00/1.65G [00:00<?, ?B/s]
0%| | 4.03M/1.65G [00:00<00:40, 40.3MB/s]
1%|▏ | 9.71M/1.65G [00:00<00:32, 49.9MB/s]
1%|▎ | 14.7M/1.65G [00:00<00:44, 36.5MB/s]
1%|▍ | 20.9M/1.65G [00:00<00:36, 44.6MB/s]
2%|▌ | 27.6M/1.65G [00:00<00:31, 51.7MB/s]
2%|▊ | 34.4M/1.65G [00:00<00:28, 56.7MB/s]
2%|▉ | 41.2M/1.65G [00:00<00:26, 60.2MB/s]
3%|█ | 47.9M/1.65G [00:00<00:25, 62.4MB/s]
3%|█▏ | 54.8M/1.65G [00:00<00:24, 64.3MB/s]
4%|█▍ | 61.6M/1.65G [00:01<00:24, 65.6MB/s]
4%|█▌ | 68.5M/1.65G [00:01<00:23, 66.5MB/s]
5%|█▋ | 75.2M/1.65G [00:01<00:37, 42.0MB/s]
5%|█▊ | 80.6M/1.65G [00:01<00:36, 42.7MB/s]
5%|█▉ | 85.9M/1.65G [00:01<00:34, 44.9MB/s]
6%|██ | 91.1M/1.65G [00:01<00:34, 45.6MB/s]
6%|██▏ | 97.9M/1.65G [00:01<00:30, 51.2MB/s]
6%|██▍ | 103M/1.65G [00:02<00:46, 33.2MB/s]
7%|██▍ | 108M/1.65G [00:02<00:42, 36.3MB/s]
7%|██▋ | 115M/1.65G [00:02<00:36, 42.5MB/s]
7%|██▊ | 121M/1.65G [00:02<00:32, 46.7MB/s]
8%|██▉ | 127M/1.65G [00:02<00:29, 51.6MB/s]
8%|███ | 134M/1.65G [00:02<00:40, 37.2MB/s]
8%|███▏ | 139M/1.65G [00:03<01:15, 20.0MB/s]
9%|███▎ | 143M/1.65G [00:03<01:09, 21.7MB/s]
9%|███▎ | 146M/1.65G [00:03<01:07, 22.5MB/s]
9%|███▍ | 149M/1.65G [00:03<01:13, 20.4MB/s]
9%|███▍ | 152M/1.65G [00:04<01:59, 12.6MB/s]
9%|███▌ | 154M/1.65G [00:04<02:08, 11.6MB/s]
9%|███▌ | 156M/1.65G [00:04<02:26, 10.2MB/s]
9%|███▌ | 157M/1.65G [00:05<02:44, 9.11MB/s]
10%|███▋ | 158M/1.65G [00:05<02:43, 9.12MB/s]
10%|███▋ | 161M/1.65G [00:05<02:03, 12.0MB/s]
10%|███▋ | 162M/1.65G [00:05<02:13, 11.2MB/s]
10%|███▊ | 164M/1.65G [00:05<02:03, 12.1MB/s]
10%|███▊ | 167M/1.65G [00:05<01:28, 16.7MB/s]
10%|███▉ | 170M/1.65G [00:05<01:18, 19.0MB/s]
10%|███▉ | 173M/1.65G [00:05<01:12, 20.3MB/s]
11%|████ | 175M/1.65G [00:06<01:33, 15.8MB/s]
11%|████ | 178M/1.65G [00:06<01:18, 18.7MB/s]
11%|████▏ | 181M/1.65G [00:06<01:08, 21.4MB/s]
11%|████▏ | 183M/1.65G [00:06<01:53, 13.0MB/s]
11%|████▎ | 185M/1.65G [00:07<02:10, 11.2MB/s]
11%|████▎ | 187M/1.65G [00:07<01:53, 12.9MB/s]
11%|████▎ | 189M/1.65G [00:07<02:35, 9.40MB/s]
12%|████▍ | 191M/1.65G [00:07<02:21, 10.3MB/s]
12%|████▍ | 194M/1.65G [00:07<01:37, 14.9MB/s]
12%|████▌ | 197M/1.65G [00:07<01:22, 17.7MB/s]
12%|████▌ | 200M/1.65G [00:08<01:56, 12.4MB/s]
12%|████▋ | 202M/1.65G [00:08<01:36, 15.0MB/s]
12%|████▋ | 206M/1.65G [00:08<01:17, 18.7MB/s]
13%|████▊ | 209M/1.65G [00:08<01:05, 22.1MB/s]
13%|████▉ | 213M/1.65G [00:08<00:58, 24.8MB/s]
13%|████▉ | 215M/1.65G [00:08<01:41, 14.2MB/s]
13%|█████ | 218M/1.65G [00:09<02:02, 11.7MB/s]
13%|█████ | 220M/1.65G [00:09<01:48, 13.2MB/s]
13%|█████ | 223M/1.65G [00:09<01:31, 15.7MB/s]
14%|█████▏ | 225M/1.65G [00:09<02:00, 11.8MB/s]
14%|█████▏ | 226M/1.65G [00:09<02:07, 11.2MB/s]
14%|█████▎ | 229M/1.65G [00:10<01:45, 13.5MB/s]
14%|█████▎ | 231M/1.65G [00:10<01:37, 14.5MB/s]
14%|█████▎ | 233M/1.65G [00:10<01:25, 16.6MB/s]
14%|█████▍ | 236M/1.65G [00:10<01:18, 17.9MB/s]
15%|█████▌ | 240M/1.65G [00:10<00:55, 25.3MB/s]
15%|█████▌ | 243M/1.65G [00:10<01:21, 17.4MB/s]
15%|█████▋ | 245M/1.65G [00:11<01:50, 12.7MB/s]
15%|█████▋ | 250M/1.65G [00:11<01:17, 18.0MB/s]
15%|█████▊ | 253M/1.65G [00:11<01:25, 16.3MB/s]
15%|█████▊ | 255M/1.65G [00:11<01:50, 12.6MB/s]
16%|█████▉ | 256M/1.65G [00:11<02:06, 11.0MB/s]
16%|██████ | 261M/1.65G [00:12<01:24, 16.5MB/s]
16%|██████ | 265M/1.65G [00:12<01:05, 21.3MB/s]
16%|██████▏ | 269M/1.65G [00:12<00:55, 24.8MB/s]
16%|██████▎ | 273M/1.65G [00:12<00:51, 26.7MB/s]
17%|██████▍ | 278M/1.65G [00:12<00:42, 32.5MB/s]
17%|██████▍ | 281M/1.65G [00:12<00:49, 27.6MB/s]
17%|██████▌ | 287M/1.65G [00:12<00:40, 34.0MB/s]
18%|██████▋ | 291M/1.65G [00:12<00:39, 34.3MB/s]
18%|██████▊ | 296M/1.65G [00:12<00:37, 36.1MB/s]
18%|██████▉ | 301M/1.65G [00:13<00:32, 41.6MB/s]
19%|███████ | 308M/1.65G [00:13<00:27, 48.4MB/s]
19%|███████▏ | 315M/1.65G [00:13<00:24, 53.8MB/s]
19%|███████▍ | 322M/1.65G [00:13<00:22, 58.0MB/s]
20%|███████▌ | 329M/1.65G [00:13<00:21, 61.5MB/s]
20%|███████▋ | 336M/1.65G [00:13<00:20, 64.0MB/s]
21%|███████▉ | 343M/1.65G [00:13<00:19, 65.8MB/s]
21%|████████ | 350M/1.65G [00:13<00:19, 67.1MB/s]
22%|████████▏ | 357M/1.65G [00:13<00:19, 68.2MB/s]
22%|████████▎ | 364M/1.65G [00:13<00:18, 68.9MB/s]
22%|████████▌ | 371M/1.65G [00:14<00:18, 69.4MB/s]
23%|████████▋ | 378M/1.65G [00:14<00:18, 69.4MB/s]
23%|████████▊ | 385M/1.65G [00:14<00:18, 69.1MB/s]
24%|█████████ | 392M/1.65G [00:14<00:18, 68.9MB/s]
24%|█████████▏ | 399M/1.65G [00:14<00:18, 68.5MB/s]
25%|█████████▎ | 406M/1.65G [00:14<00:18, 68.3MB/s]
25%|█████████▍ | 412M/1.65G [00:14<00:18, 68.1MB/s]
25%|█████████▋ | 419M/1.65G [00:14<00:18, 67.8MB/s]
26%|█████████▊ | 426M/1.65G [00:14<00:18, 67.5MB/s]
26%|█████████▉ | 433M/1.65G [00:14<00:18, 67.7MB/s]
27%|██████████ | 440M/1.65G [00:15<00:17, 67.8MB/s]
27%|██████████▎ | 446M/1.65G [00:15<00:17, 68.0MB/s]
27%|██████████▍ | 453M/1.65G [00:15<00:17, 68.1MB/s]
28%|██████████▌ | 460M/1.65G [00:15<00:17, 68.5MB/s]
28%|██████████▋ | 467M/1.65G [00:15<00:17, 69.1MB/s]
29%|██████████▉ | 474M/1.65G [00:15<00:16, 69.7MB/s]
29%|███████████ | 481M/1.65G [00:15<00:16, 70.0MB/s]
30%|███████████▏ | 489M/1.65G [00:15<00:16, 70.3MB/s]
30%|███████████▍ | 496M/1.65G [00:15<00:16, 70.3MB/s]
30%|███████████▌ | 503M/1.65G [00:15<00:16, 70.3MB/s]
31%|███████████▋ | 510M/1.65G [00:16<00:16, 70.4MB/s]
31%|███████████▉ | 517M/1.65G [00:16<00:16, 69.9MB/s]
32%|████████████ | 524M/1.65G [00:16<00:16, 69.5MB/s]
32%|████████████▏ | 531M/1.65G [00:16<00:16, 69.1MB/s]
33%|████████████▎ | 538M/1.65G [00:16<00:16, 68.6MB/s]
33%|████████████▌ | 544M/1.65G [00:16<00:16, 68.2MB/s]
33%|████████████▋ | 551M/1.65G [00:16<00:16, 68.1MB/s]
34%|████████████▊ | 558M/1.65G [00:16<00:16, 67.8MB/s]
34%|████████████▉ | 565M/1.65G [00:16<00:16, 67.5MB/s]
35%|█████████████▏ | 572M/1.65G [00:16<00:15, 67.8MB/s]
35%|█████████████▎ | 579M/1.65G [00:17<00:15, 68.0MB/s]
35%|█████████████▍ | 585M/1.65G [00:17<00:15, 68.1MB/s]
36%|█████████████▌ | 592M/1.65G [00:17<00:15, 68.3MB/s]
36%|█████████████▊ | 599M/1.65G [00:17<00:15, 68.0MB/s]
37%|█████████████▉ | 606M/1.65G [00:17<00:15, 68.1MB/s]
37%|██████████████ | 613M/1.65G [00:17<00:15, 68.2MB/s]
37%|██████████████▏ | 620M/1.65G [00:17<00:15, 68.1MB/s]
38%|██████████████▍ | 626M/1.65G [00:17<00:15, 68.1MB/s]
38%|██████████████▌ | 633M/1.65G [00:17<00:14, 68.1MB/s]
39%|██████████████▋ | 640M/1.65G [00:18<00:16, 62.7MB/s]
39%|██████████████▊ | 646M/1.65G [00:18<00:17, 57.8MB/s]
40%|███████████████ | 653M/1.65G [00:18<00:16, 59.9MB/s]
40%|███████████████▏ | 660M/1.65G [00:18<00:16, 61.6MB/s]
40%|███████████████▎ | 666M/1.65G [00:18<00:15, 63.3MB/s]
41%|███████████████▍ | 673M/1.65G [00:18<00:15, 64.6MB/s]
41%|███████████████▋ | 680M/1.65G [00:18<00:14, 65.3MB/s]
42%|███████████████▊ | 686M/1.65G [00:18<00:14, 65.2MB/s]
42%|███████████████▉ | 693M/1.65G [00:18<00:14, 65.6MB/s]
42%|████████████████ | 700M/1.65G [00:18<00:14, 66.2MB/s]
43%|████████████████▏ | 707M/1.65G [00:19<00:14, 66.8MB/s]
43%|████████████████▍ | 713M/1.65G [00:19<00:14, 63.9MB/s]
44%|████████████████▌ | 720M/1.65G [00:19<00:14, 64.2MB/s]
44%|████████████████▋ | 726M/1.65G [00:19<00:26, 35.2MB/s]
44%|████████████████▊ | 733M/1.65G [00:19<00:22, 41.3MB/s]
45%|█████████████████ | 740M/1.65G [00:19<00:19, 46.9MB/s]
45%|█████████████████▏ | 747M/1.65G [00:19<00:17, 51.9MB/s]
46%|█████████████████▎ | 754M/1.65G [00:20<00:15, 56.3MB/s]
46%|█████████████████▍ | 761M/1.65G [00:20<00:14, 59.6MB/s]
46%|█████████████████▋ | 767M/1.65G [00:20<00:14, 61.6MB/s]
47%|█████████████████▊ | 774M/1.65G [00:20<00:13, 63.3MB/s]
47%|█████████████████▉ | 781M/1.65G [00:20<00:13, 64.2MB/s]
48%|██████████████████ | 787M/1.65G [00:20<00:13, 64.9MB/s]
48%|██████████████████▎ | 794M/1.65G [00:20<00:13, 65.4MB/s]
48%|██████████████████▍ | 801M/1.65G [00:20<00:12, 65.6MB/s]
49%|██████████████████▌ | 807M/1.65G [00:20<00:12, 65.2MB/s]
49%|██████████████████▋ | 814M/1.65G [00:20<00:12, 65.6MB/s]
50%|██████████████████▊ | 821M/1.65G [00:21<00:12, 64.8MB/s]
50%|███████████████████ | 828M/1.65G [00:21<00:12, 66.1MB/s]
50%|███████████████████▏ | 835M/1.65G [00:21<00:12, 67.5MB/s]
51%|███████████████████▎ | 841M/1.65G [00:21<00:12, 67.4MB/s]
51%|███████████████████▍ | 848M/1.65G [00:21<00:13, 57.9MB/s]
52%|███████████████████▋ | 854M/1.65G [00:21<00:15, 50.4MB/s]
52%|███████████████████▊ | 860M/1.65G [00:21<00:15, 51.6MB/s]
52%|███████████████████▉ | 865M/1.65G [00:21<00:16, 48.1MB/s]
53%|████████████████████ | 872M/1.65G [00:21<00:14, 52.8MB/s]
53%|████████████████████▏ | 879M/1.65G [00:22<00:13, 57.3MB/s]
54%|████████████████████▎ | 886M/1.65G [00:22<00:12, 60.8MB/s]
54%|████████████████████▌ | 893M/1.65G [00:22<00:12, 62.9MB/s]
54%|████████████████████▋ | 899M/1.65G [00:22<00:11, 64.8MB/s]
55%|████████████████████▊ | 906M/1.65G [00:22<00:11, 66.0MB/s]
55%|████████████████████▉ | 913M/1.65G [00:22<00:11, 65.2MB/s]
56%|█████████████████████▏ | 920M/1.65G [00:22<00:11, 65.3MB/s]
56%|█████████████████████▎ | 926M/1.65G [00:22<00:11, 65.5MB/s]
56%|█████████████████████▍ | 933M/1.65G [00:22<00:10, 65.5MB/s]
57%|█████████████████████▌ | 940M/1.65G [00:22<00:10, 66.3MB/s]
57%|█████████████████████▊ | 946M/1.65G [00:23<00:10, 67.1MB/s]
58%|█████████████████████▉ | 953M/1.65G [00:23<00:10, 67.0MB/s]
58%|██████████████████████ | 960M/1.65G [00:23<00:10, 67.7MB/s]
59%|██████████████████████▏ | 967M/1.65G [00:23<00:10, 67.4MB/s]
59%|██████████████████████▍ | 974M/1.65G [00:23<00:10, 67.6MB/s]
59%|██████████████████████▌ | 980M/1.65G [00:23<00:09, 67.6MB/s]
60%|██████████████████████▋ | 987M/1.65G [00:23<00:09, 67.8MB/s]
60%|██████████████████████▊ | 994M/1.65G [00:23<00:14, 46.7MB/s]
60%|██████████████████████▍ | 1.00G/1.65G [00:24<00:19, 33.3MB/s]
61%|██████████████████████▌ | 1.01G/1.65G [00:24<00:17, 37.9MB/s]
61%|██████████████████████▋ | 1.01G/1.65G [00:24<00:15, 42.0MB/s]
62%|██████████████████████▊ | 1.02G/1.65G [00:24<00:13, 47.6MB/s]
62%|██████████████████████▉ | 1.02G/1.65G [00:24<00:11, 52.7MB/s]
62%|███████████████████████ | 1.03G/1.65G [00:24<00:11, 53.0MB/s]
63%|███████████████████████▏ | 1.04G/1.65G [00:24<00:11, 54.1MB/s]
63%|███████████████████████▎ | 1.04G/1.65G [00:24<00:10, 56.8MB/s]
63%|███████████████████████▍ | 1.05G/1.65G [00:25<00:10, 59.0MB/s]
64%|███████████████████████▋ | 1.06G/1.65G [00:25<00:10, 56.7MB/s]
64%|███████████████████████▊ | 1.06G/1.65G [00:25<00:10, 58.9MB/s]
65%|███████████████████████▉ | 1.07G/1.65G [00:25<00:09, 58.9MB/s]
65%|████████████████████████ | 1.07G/1.65G [00:25<00:11, 49.5MB/s]
65%|████████████████████████▏ | 1.08G/1.65G [00:25<00:10, 53.4MB/s]
66%|████████████████████████▎ | 1.09G/1.65G [00:25<00:09, 57.0MB/s]
66%|████████████████████████▍ | 1.09G/1.65G [00:25<00:09, 59.8MB/s]
67%|████████████████████████▋ | 1.10G/1.65G [00:25<00:08, 62.2MB/s]
67%|████████████████████████▊ | 1.11G/1.65G [00:26<00:08, 61.3MB/s]
67%|████████████████████████▉ | 1.11G/1.65G [00:26<00:09, 58.2MB/s]
68%|█████████████████████████ | 1.12G/1.65G [00:26<00:08, 60.3MB/s]
68%|█████████████████████████▏ | 1.13G/1.65G [00:26<00:08, 62.0MB/s]
69%|█████████████████████████▎ | 1.13G/1.65G [00:26<00:08, 61.2MB/s]
69%|█████████████████████████▌ | 1.14G/1.65G [00:26<00:08, 62.1MB/s]
69%|█████████████████████████▋ | 1.15G/1.65G [00:26<00:08, 62.6MB/s]
70%|█████████████████████████▊ | 1.15G/1.65G [00:26<00:07, 63.2MB/s]
70%|█████████████████████████▉ | 1.16G/1.65G [00:26<00:07, 63.9MB/s]
71%|██████████████████████████ | 1.17G/1.65G [00:26<00:07, 65.0MB/s]
71%|██████████████████████████▏ | 1.17G/1.65G [00:27<00:07, 65.9MB/s]
71%|██████████████████████████▍ | 1.18G/1.65G [00:27<00:07, 66.6MB/s]
72%|██████████████████████████▌ | 1.19G/1.65G [00:27<00:06, 67.2MB/s]
72%|██████████████████████████▋ | 1.19G/1.65G [00:27<00:06, 67.4MB/s]
73%|██████████████████████████▊ | 1.20G/1.65G [00:27<00:06, 67.6MB/s]
73%|███████████████████████████ | 1.21G/1.65G [00:27<00:06, 67.8MB/s]
73%|███████████████████████████▏ | 1.21G/1.65G [00:27<00:06, 68.1MB/s]
74%|███████████████████████████▎ | 1.22G/1.65G [00:27<00:06, 68.3MB/s]
74%|███████████████████████████▍ | 1.23G/1.65G [00:27<00:06, 68.4MB/s]
75%|███████████████████████████▌ | 1.23G/1.65G [00:27<00:06, 68.4MB/s]
75%|███████████████████████████▊ | 1.24G/1.65G [00:28<00:06, 68.4MB/s]
75%|███████████████████████████▉ | 1.25G/1.65G [00:28<00:05, 68.2MB/s]
76%|████████████████████████████ | 1.25G/1.65G [00:28<00:05, 68.3MB/s]
76%|████████████████████████████▏ | 1.26G/1.65G [00:28<00:05, 68.2MB/s]
77%|████████████████████████████▍ | 1.27G/1.65G [00:28<00:05, 68.1MB/s]
77%|████████████████████████████▌ | 1.27G/1.65G [00:28<00:05, 68.0MB/s]
78%|████████████████████████████▋ | 1.28G/1.65G [00:28<00:05, 67.8MB/s]
78%|████████████████████████████▊ | 1.29G/1.65G [00:28<00:05, 67.7MB/s]
78%|████████████████████████████▉ | 1.30G/1.65G [00:28<00:05, 67.6MB/s]
79%|█████████████████████████████▏ | 1.30G/1.65G [00:28<00:05, 67.5MB/s]
79%|█████████████████████████████▎ | 1.31G/1.65G [00:29<00:05, 67.8MB/s]
80%|█████████████████████████████▍ | 1.32G/1.65G [00:29<00:04, 68.0MB/s]
80%|█████████████████████████████▌ | 1.32G/1.65G [00:29<00:04, 68.1MB/s]
80%|█████████████████████████████▊ | 1.33G/1.65G [00:29<00:04, 68.1MB/s]
81%|█████████████████████████████▉ | 1.34G/1.65G [00:29<00:04, 68.1MB/s]
81%|██████████████████████████████ | 1.34G/1.65G [00:29<00:04, 68.3MB/s]
82%|██████████████████████████████▏ | 1.35G/1.65G [00:29<00:04, 68.3MB/s]
82%|██████████████████████████████▎ | 1.36G/1.65G [00:29<00:04, 68.4MB/s]
82%|██████████████████████████████▌ | 1.36G/1.65G [00:29<00:04, 68.5MB/s]
83%|██████████████████████████████▋ | 1.37G/1.65G [00:29<00:04, 68.3MB/s]
83%|██████████████████████████████▊ | 1.38G/1.65G [00:30<00:04, 68.3MB/s]
84%|██████████████████████████████▉ | 1.38G/1.65G [00:30<00:03, 68.6MB/s]
84%|███████████████████████████████▏ | 1.39G/1.65G [00:30<00:03, 68.9MB/s]
85%|███████████████████████████████▎ | 1.40G/1.65G [00:30<00:03, 68.9MB/s]
85%|███████████████████████████████▍ | 1.40G/1.65G [00:30<00:03, 68.9MB/s]
85%|███████████████████████████████▌ | 1.41G/1.65G [00:30<00:03, 65.8MB/s]
86%|███████████████████████████████▊ | 1.42G/1.65G [00:30<00:03, 61.4MB/s]
86%|███████████████████████████████▉ | 1.43G/1.65G [00:30<00:03, 63.0MB/s]
87%|████████████████████████████████ | 1.43G/1.65G [00:30<00:03, 64.7MB/s]
87%|████████████████████████████████▏ | 1.44G/1.65G [00:31<00:03, 66.0MB/s]
87%|████████████████████████████████▎ | 1.45G/1.65G [00:31<00:03, 66.7MB/s]
88%|████████████████████████████████▌ | 1.45G/1.65G [00:31<00:02, 67.6MB/s]
88%|████████████████████████████████▋ | 1.46G/1.65G [00:31<00:02, 68.1MB/s]
89%|████████████████████████████████▊ | 1.47G/1.65G [00:31<00:02, 68.5MB/s]
89%|████████████████████████████████▉ | 1.47G/1.65G [00:31<00:02, 68.9MB/s]
90%|█████████████████████████████████▏ | 1.48G/1.65G [00:31<00:02, 69.2MB/s]
90%|█████████████████████████████████▎ | 1.49G/1.65G [00:31<00:02, 69.5MB/s]
90%|█████████████████████████████████▍ | 1.49G/1.65G [00:31<00:02, 69.3MB/s]
91%|█████████████████████████████████▌ | 1.50G/1.65G [00:31<00:02, 69.5MB/s]
91%|█████████████████████████████████▊ | 1.51G/1.65G [00:32<00:02, 67.3MB/s]
92%|█████████████████████████████████▉ | 1.52G/1.65G [00:32<00:02, 54.7MB/s]
92%|██████████████████████████████████ | 1.52G/1.65G [00:32<00:02, 56.1MB/s]
92%|██████████████████████████████████▏ | 1.53G/1.65G [00:32<00:02, 50.5MB/s]
93%|██████████████████████████████████▎ | 1.53G/1.65G [00:32<00:02, 54.0MB/s]
93%|██████████████████████████████████▍ | 1.54G/1.65G [00:32<00:03, 33.1MB/s]
93%|██████████████████████████████████▌ | 1.54G/1.65G [00:33<00:03, 31.6MB/s]
94%|██████████████████████████████████▋ | 1.55G/1.65G [00:33<00:02, 38.1MB/s]
94%|██████████████████████████████████▊ | 1.56G/1.65G [00:33<00:02, 44.6MB/s]
95%|███████████████████████████████████ | 1.56G/1.65G [00:33<00:01, 50.4MB/s]
95%|███████████████████████████████████▏ | 1.57G/1.65G [00:33<00:01, 55.4MB/s]
95%|███████████████████████████████████▎ | 1.58G/1.65G [00:33<00:01, 59.3MB/s]
96%|███████████████████████████████████▍ | 1.59G/1.65G [00:33<00:01, 62.3MB/s]
96%|███████████████████████████████████▋ | 1.59G/1.65G [00:33<00:00, 64.5MB/s]
97%|███████████████████████████████████▊ | 1.60G/1.65G [00:33<00:00, 66.2MB/s]
97%|███████████████████████████████████▉ | 1.61G/1.65G [00:33<00:00, 67.3MB/s]
98%|████████████████████████████████████ | 1.61G/1.65G [00:34<00:00, 68.2MB/s]
98%|████████████████████████████████████▎| 1.62G/1.65G [00:34<00:00, 68.6MB/s]
98%|████████████████████████████████████▍| 1.63G/1.65G [00:34<00:00, 69.0MB/s]
99%|████████████████████████████████████▌| 1.63G/1.65G [00:34<00:00, 69.2MB/s]
99%|████████████████████████████████████▋| 1.64G/1.65G [00:34<00:00, 69.0MB/s]
100%|████████████████████████████████████▉| 1.65G/1.65G [00:34<00:00, 69.0MB/s]
0%| | 0.00/1.65G [00:00<?, ?B/s]
100%|█████████████████████████████████████| 1.65G/1.65G [00:00<00:00, 5.57TB/s]
Download complete in 01m08s (1576.2 MB)
Filtering raw data in 1 contiguous segment
Setting up high-pass filter at 2 Hz
IIR filter parameters
---------------------
Butterworth highpass zero-phase (two-pass forward and reverse) non-causal filter:
- Filter order 8 (effective, after forward-backward)
- Cutoff at 2.00 Hz: -6.02 dB
Embedding of Xdawn covariance matrices¶
nfilter = 4
xdwn = XdawnCovariances(estimator="scm", nfilter=nfilter)
split = train_test_split(X, y, train_size=0.25, random_state=42)
Xtrain, Xtest, ytrain, ytest = split
covs = xdwn.fit(Xtrain, ytrain).transform(Xtest)
Laplacian Eigenmaps (LE), also called Spectral Embedding (SE)¶
plot_embedding(covs, ytest, metric="riemann", embd_type="Spectral",
normalize=True)
plt.show()
Locally Linear Embedding (LLE)¶
plot_embedding(covs, ytest, metric="riemann", embd_type="LocallyLinear",
normalize=False)
plt.show()
References¶
Total running time of the script: (1 minutes 13.315 seconds)