Scribd Logo
Upload

Welcome to Scribd!

  • Páginas de Pyriemann Documentação

    Uploaded by

    Lucas Custodio
    18 views2 pages
    This document summarizes an EEG decoding method that uses spatial filtering with Xdawn to estimate covariance matrices from EEG epochs. The covariances are then projected into tangent space and classified with logistic regression. It presents an example analysis pipeline applying this method to classify auditory and visual ERP responses from MEG data, obtaining a cross-validated classification accuracy.

    Original Description:

    Original Title

    Páginas de Pyriemann documentação

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    Report this Document
    This document summarizes an EEG decoding method that uses spatial filtering with Xdawn to estimate covariance matrices from EEG epochs. The covariances are then projected into tangent space and classified with logistic regression. It presents an example analysis pipeline applying this method to classify auditory and visual ERP responses from MEG data, obtaining a cross-validated classification accuracy.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Flag for inappropriate content
    Download as pdf or txt
    18 views2 pages

    Páginas de Pyriemann Documentação

    Uploaded by

    Lucas Custodio
    This document summarizes an EEG decoding method that uses spatial filtering with Xdawn to estimate covariance matrices from EEG epochs. The covariances are then projected into tangent space and classified with logistic regression. It presents an example analysis pipeline applying this method to classify auditory and visual ERP responses from MEG data, obtaining a cross-validated classification accuracy.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Flag for inappropriate content
    Download as pdf or txt
    of 2

    pyRiemann Documentation, Release 0.2.

    4.1.2 ERP EEG decoding in Tangent space.

    Decoding applied to EEG data in sensor space decomposed using Xdawn. After spatial filtering, covariances matrices
    are estimated, then projected in the tangent space and classified with a logistic regression.
    # Authors: Alexandre Barachant <[email protected]>
    #
    # License: BSD (3-clause)

    import numpy as np

    from pyriemann.estimation import XdawnCovariances


    from pyriemann.tangentspace import TangentSpace
    from pyriemann.utils.viz import plot_confusion_matrix

    import mne
    from mne import io
    from mne.datasets import sample

    from sklearn.pipeline import make_pipeline


    from sklearn.model_selection import KFold
    from sklearn.linear_model import LogisticRegression

    from matplotlib import pyplot as plt

    print(__doc__)

    data_path = sample.data_path()

    Set parameters and read data


    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=False, eeg=True, 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)

    labels = epochs.events[:, -1]


    evoked = epochs.average()

    Decoding in tangent space with a logistic regression


    n_components = 2 # pick some components

    # Define a monte-carlo cross-validation generator (reduce variance):


    cv = KFold(n_splits=10, shuffle=True, random_state=42)
    (continues on next page)

    12 Chapter 4. Examples Gallery


    pyRiemann Documentation, Release 0.2.6

    (continued from previous page)


    epochs_data = epochs.get_data()

    clf = make_pipeline(XdawnCovariances(n_components),
    TangentSpace(metric='riemann'),
    LogisticRegression())

    preds = np.zeros(len(labels))

    for train_idx, test_idx in cv.split(epochs_data):


    y_train, y_test = labels[train_idx], labels[test_idx]

    clf.fit(epochs_data[train_idx], y_train)
    preds[test_idx] = clf.predict(epochs_data[test_idx])

    # Printing the results


    acc = np.mean(preds == labels)
    print("Classification accuracy: %f " % (acc))

    names = ['audio left', 'audio right', 'vis left', 'vis right']


    plot_confusion_matrix(preds, labels, names)
    plt.show()

    Out:

    4.1. Classification of ERP 13

    You might also like