HF-SEF dataset

This example looks at high frequency SEF responses.

# Author: Jussi Nurminen (jnu@iki.fi)
# License: BSD (3-clause)

import mne
import os
from mne.datasets import hf_sef

fname_evoked = os.path.join(hf_sef.data_path(),


Read evoked data

evoked = mne.Evoked(fname_evoked)


Found the data of interest:
        t =     -50.00 ...     250.00 ms (SEF)
        0 CTF compensation matrices available
        nave = 2790 - aspect type = 100
No projector specified for this dataset. Please consider the method self.add_proj.

Create a highpass filtered version

evoked_hp = evoked.copy()
evoked_hp.filter(l_freq=300, h_freq=None, fir_design='firwin')


Setting up high-pass filter at 3e+02 Hz
l_trans_bandwidth chosen to be 75.0 Hz
Filter length of 133 samples (0.044 sec) selected

Compare high-pass filtered and unfiltered data on a single channel

ch = 'MEG0443'
pick = evoked.ch_names.index(ch)
edi = {'HF': evoked_hp, 'Regular': evoked}
mne.viz.plot_compare_evokeds(edi, picks=pick)

Total running time of the script: ( 0 minutes 1.836 seconds)

Estimated memory usage: 10 MB

Gallery generated by Sphinx-Gallery