Given an observation (A, red), we estimate model parameters in an initial interval of data. At subsequent times, we use the observed data to (B) predict and update the state estimate (purple), and (C…
(A) Example 2 s of simulated observed data (thick curves) and the true phase (thin blue curves) for each scenario. (B) For each scenario, example spectra of (B.i) the signal, and (B.ii) the …
Circular standard deviation for all methods.
(A) Example observed signal (black) and phase for two simulated rhythms, the confounding rhythm at 5 Hz (blue) and the target rhythm at 6 Hz (red). (B) Example phase estimates for the target rhythm …
Circular standard deviation for all methods on two simultaneous oscillations simulation.
(A) Example phase of a 6 Hz sinusoid (blue) with four phase resets (red dashed lines). (B) Example phase estimates (thin curves) and the true phase (thick blue curve) at the indicated phase reset in …
(A) Example 2 s of simulated observed data (thick curves) and true phase (thin blue curves) for each SNR value. (B) For each SNR, (i) the average credible interval width (colored circles) of the …
Circular standard deviation for all methods on signal-to-noise ratio (SNR) simulation.
At each signal-to-noise ratio (SNR), for a smooth spectrum, we used 15 s of data and applied the multitaper method with 29 tapers and a bandwidth of 1 Hz. The SNR is controlled by shifting the power …
(A) The phase error for a single instance of the model fit (using data at times 10–20 s, blue curve), and the 90% interval for phase error (red bands) at time t derived using parameter estimates …
Error and spectrogram for in vivo local field potential (LFP) data.
The error (circular deviation) in performance of each algorithm across all LFP data. Each dot represents the error computed for 1 s of data. To compute the error, we compare the real-time phase …
(A) Example rodent local field potential (LFP) data (red, solid) with a consistent broadband peak in the theta band (4–8 Hz). The estimated state of the SSPE method (purple, dashed) tracks the …
Amplitude and credible interval widths for local field potential (LFP) data.
(A.i.) Example rodent local field potential (LFP) data. (A.ii) Phase estimates from each method (see legend) are consistent when the SNR of the theta rhythm is high.
(A) Phase error for a model fit using data from the first 10 s (blue curve), and the 95% confidence intervals (red) in the error for phase estimates at time t using models fit at times prior to time …
Error and spectrogram for in vivo electroencephalogram (EEG) analysis.
(A.i) Example mu rhythm clearly visible in the early portion of the electroencephalogram (EEG) trace. (A.ii) Phase estimates from each method (see legend) are consistent when the SNR of the mu …
The error (circular deviation) in performance of each algorithm across all EEG data. Each dot represents the error computed for 1 s of data. To compute the error, we compare the real-time phase …
(A.i) The observed data (orange) is well tracked by the (summed) state estimate (purple, dotted) of the 10 and 22 Hz oscillators. (A.ii) Phase estimates for the 10 Hz (blue, dotted) and 20 Hz …
(A) Log-transformed MEP amplitude versus phase with no threshold applied to the credible interval width. White dot indicates the median of the MEP amplitude. Yellow line is the best fit circular …
(A) Open Ephys GUI for using SSPE. The user specifies the number of frequencies to track, the center frequencies to track, the frequencies of interest for phase calculation and output (FOI), …
The amplitude of the 6 Hz sinusoid waxes and wanes in time.
We find that SSPE still performs well for two oscillations at nearby frequencies, with random initial phases; compare to Figure 3 of the main manuscript.
(A,B) Two example traces showing intervals in which the rhythm amplitude is small (below the 65th percentile, white intervals, also indicated by purple trace equaling 0), yet the credible intervals …
SSPE | Blackwood | Zrenner | Acausal FIR | |
---|---|---|---|---|
Error (std. dev.) | 2.85 (0.89)* | 62.08 (0.38) | 44.68 (0.38) | 15.04 (0.23) |
Time to convergence(std. dev.) | 34 ms (162) | 190 ms (33) | 747 ms (22) | 555 ms (206) |
p~0 compared to all other methods.
Circular standard deviation and convergence time information for phase reset simulation.