Time-trend analysis of the center frequency of the intrinsic mode function from the Hilbert–Huang transform of electroencephalography during general anesthesia: a retrospective observational study

BACKGROUND: Anesthesiologists are required to maintain an optimal depth of anesthesia during general anesthesia, and several electroencephalogram (EEG) processing methods have been developed and approved for clinical use to evaluate anesthesia depth. Recently, the Hilbert–Huang transform (HHT) was introduced to analyze nonlinear and nonstationary data. In this study, we assessed whether the changes in EEG characteristics during general anesthesia that are analyzed by the HHT are useful for monitoring the depth of anesthesia. METHODS: This retrospective observational study enrolled patients who underwent propofol anesthesia. Raw EEG signals were obtained from a monitor through a previously developed software application. We developed an HHT analyzer to decompose the EEG signal into six intrinsic mode functions (IMFs) and estimated the instantaneous frequencies (HHT_IF) for each IMF. Changes over time in the raw EEG waves and parameters such as HHT_IF, BIS, spectral edge frequency 95 (SEF95), and electromyogram parameter (EMGlow) were assessed, and a Gaussian process regression model was created to assess the association between BIS and HHT_IF. RESULTS: We analyzed EEG signals from 30 patients. The beta oscillation frequency range (13–25 Hz) was detected in IMF1 and IMF2 during the awake state, then after loss of consciousness, the frequency decreased and alpha oscillation (8–12 Hz) was detected in IMF2. At the emergence phase, the frequency increased and beta oscillations were detected in IMF1, IMF2, and IMF3. BIS and EMGlow changed significantly during the induction and emergence phases, whereas SEF95 showed a wide variability and no significant changes during the induction phase. The root mean square error between the observed BIS values and the values predicted by a Gaussian process regression model ranged from 4.69 to 9.68. CONCLUSIONS: We applied the HHT to EEG analyses during propofol anesthesia. The instantaneous frequency in IMF1 and IMF2 identified changes in EEG characteristics during induction and emergence from general anesthesia. Moreover, the HHT_IF in IMF2 showed strong associations with BIS and was suitable for depicting the alpha oscillation. Our study suggests that the HHT is useful for monitoring the depth of anesthesia. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12871-023-02082-4.