This example looks at high-frequency SEF responses.
# Author: Jussi Nurminen (firstname.lastname@example.org) # # License: BSD-3-Clause
Read evoked data
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
Setting up high-pass filter at 3e+02 Hz FIR filter parameters --------------------- Designing a one-pass, zero-phase, non-causal highpass filter: - Windowed time-domain design (firwin) method - Hamming window with 0.0194 passband ripple and 53 dB stopband attenuation - Lower passband edge: 300.00 - Lower transition bandwidth: 75.00 Hz (-6 dB cutoff frequency: 262.50 Hz) - Filter length: 133 samples (0.044 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 306 out of 306 | elapsed: 0.1s finished
|Timepoints||901 samples||Channels||306 channels|
|Number of averaged epochs||2790|
|Time range (secs)||-0.05 – 0.25|
Compare high-pass filtered and unfiltered data on a single channel
Total running time of the script: ( 0 minutes 6.481 seconds)
Estimated memory usage: 9 MB