Time-Frequency and Oscillations
Estimate time-varying spectral content, normalize it correctly, and decide whether a "rhythm" is real before reporting it.
You can now…
- Explain evoked vs induced
- Choose wavelet parameters and predict the resolution trade-off
- Use multitaper and filter-Hilbert appropriately
- Baseline-normalize TF power
- Compute ITC
- Detect bursts and parameterize spectra
- Run TF cluster statistics
Level 4 is about what the trial average cannot see. It opens by showing that a burst present on every single trial can be made to vanish from the ERP by letting its phase wander, which is the case for looking at power and phase separately; then it builds the estimators that do that — the short-time Fourier transform and the Morlet wavelet with their fixed uncertainty product, multitaper with a bandwidth you state, filter-Hilbert with a filter you designed — and shows the three disagreeing on the same three seconds of real data. From there: normalizing a map against a baseline, which is where most time-frequency results go wrong, and the event-related desynchronization that is the level’s worked example; the phase itself, intertrial coherence and its small-sample floor, and cross-frequency coupling with the waveform-shape confound that produces it out of nothing; and then the lesson the level exists for.
Two threads run through it. The first is that every number here is a ratio or a smoothing, and both have parameters — cycles, bandwidth, band edges, baseline window, normalization, averaging order, threshold — so a time-frequency result travels with its parameters or it is not a result. The second is stated as a question in L4.6’s title: band power is not the same claim as a rhythm. A band with no spectral peak in it has no rhythm to desynchronize, a narrow filter makes anything look rhythmic, a non-sinusoidal rhythm produces its own harmonics and its own cross-frequency coupling, and a burst count moves by a factor of eleven on a threshold nobody has agreed on. L4.7 is the third stop on the statistics thread that runs from L3.7 through this level to L6.1 and L6.4: the cluster-permutation test is unchanged by a frequency axis, and the discipline needed to use it honestly is not. The capstone puts it together on motor imagery across a documented subset of subjects, including the ones with no mu peak.
Lessons
-
L4.1 Why time-frequency
What trial averaging hides: evoked versus induced versus total power, and when a time-frequency analysis is warranted.
-
L4.2 STFT and Morlet wavelets
The STFT and its window trade-off, Morlet wavelets defined by frequency and cycles, the time-frequency uncertainty, and edge effects.
-
L4.3 Multitaper and filter-Hilbert
Multitaper spectral concentration and bandwidth, band-limited amplitude and phase from filter-Hilbert, and when to prefer each.
-
L4.4 Baseline normalization and ERD/ERS
Why raw TF power needs normalization, dB versus percent versus z baselines, choosing an edge-safe baseline, and computing ERD/ERS.
-
L4.5 Phase, ITC and cross-frequency coupling
Intertrial phase coherence, phase resetting versus additive evoked models as competing accounts, and phase-amplitude coupling with its confounds.
-
L4.6 Is it really an oscillation?
Test for a spectral peak before reporting band power, detect bursts, run cycle-by-cycle analysis, and recognize non-sinusoidal waveforms.
-
L4.7 Statistics for time-frequency
Cluster-permutation tests over channel × time × frequency, TFCE, limiting the search space with a priori ROIs and bands, and reporting.