.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "auto_examples/biosignal-erp/plot_embedding_MEG.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. .. rst-class:: sphx-glr-example-title .. _sphx_glr_auto_examples_biosignal-erp_plot_embedding_MEG.py: ============================================ 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]_. .. GENERATED FROM PYTHON SOURCE LINES 19-35 .. code-block:: Python # Authors: Pedro Rodrigues , # Gabriel Wagner vom Berg # 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__) .. GENERATED FROM PYTHON SOURCE LINES 36-38 Set parameters and read data ---------------------------- .. GENERATED FROM PYTHON SOURCE LINES 38-61 .. code-block:: Python 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] .. rst-class:: sphx-glr-script-out .. code-block:: none Using default location ~/mne_data for sample... 0%| | 0.00/1.65G [00:00`_ A. Goh and R Vidal, in 2008 IEEE Conference on Computer Vision and Pattern Recognition. .. rst-class:: sphx-glr-timing **Total running time of the script:** (1 minutes 13.315 seconds) .. _sphx_glr_download_auto_examples_biosignal-erp_plot_embedding_MEG.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_embedding_MEG.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_embedding_MEG.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: plot_embedding_MEG.zip ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_