Next Article in Journal
An Electrochemical Nucleic Acid Biosensor for Triple-Negative Breast Cancer Biomarker Detection
Previous Article in Journal
Research on Trajectory Planning of Autonomous Vehicles in Constrained Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variational Mode Decomposition Analysis of Electroencephalograms during General Anesthesia: Using the Grey Wolf Optimizer to Determine Hyperparameters

by
Kosuke Kushimoto
1,
Yurie Obata
2,
Tomomi Yamada
1,
Mao Kinoshita
1,
Koichi Akiyama
3 and
Teiji Sawa
4,*
1
Department of Anesthesiology, Graduate School of Medical Science, Kyoto Prefectural University of Medicine, Kyoto 602-8566, Japan
2
Department of Anesthesiology, Yodogawa Christian Hospital, Osaka 533-0024, Japan
3
Department of Anesthesiology, Kindai University, Higashiosaka 577-8502, Japan
4
Hospital of Kyoto Prefectural University of Medicine, Kyoto 602-8566, Japan
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(17), 5749; https://doi.org/10.3390/s24175749
Submission received: 6 July 2024 / Revised: 28 August 2024 / Accepted: 3 September 2024 / Published: 4 September 2024
(This article belongs to the Section Biomedical Sensors)

Abstract

:
Frequency analysis via electroencephalography (EEG) during general anesthesia is used to develop techniques for measuring anesthesia depth. Variational mode decomposition (VMD) enables mathematical optimization methods to decompose EEG signals into natural number intrinsic mode functions with distinct narrow bands. However, the analysis requires the a priori determination of hyperparameters, including the decomposition number (K) and the penalty factor (PF). In the VMD analysis of EEGs derived from a noninterventional and noninvasive retrospective observational study, we adapted the grey wolf optimizer (GWO) to determine the K and PF hyperparameters of the VMD. As a metric for optimization, we calculated the envelope function of the IMF decomposed via the VMD method and used its envelope entropy as the fitness function. The K and PF values varied in each epoch, with one epoch being the analytical unit of EEG; however, the fitness values showed convergence at an early stage in the GWO algorithm. The K value was set to 2 to capture the α wave enhancement observed during the maintenance phase of general anesthesia in intrinsic mode function 2 (IMF-2). This study suggests that using the GWO to optimize VMD hyperparameters enables the construction of a robust analytical model for examining the EEG frequency characteristics involved in the effects of general anesthesia.

1. Introduction

The decoding technology used to monitor electroencephalograms (EEGs) during general anesthesia (GA) and extract the information required to understand the depth of anesthesia or the state of consciousness is crucial not only for comprehending brain function during GA but also for assessing the appropriateness of the anesthesia level [1,2,3]. There have been numerous reports on new technologies for EEG signal processing that can significantly broaden the scope of EEG monitoring for assessing patients’ anesthetic sedation states during GA [4,5]. In addition, there are growing expectations for advancements in the intelligent sensors that are involved in detecting biological signals, spurred by the rapid growth of flexible electronics [6,7]. The potential benefits of these technologies, when combined with fast AI diagnostic algorithms, are significant. This combination holds promise for more reliable clinical diagnostics, especially in the analysis of EEGs. However, the analysis parameters of commercially available EEG monitors have only been optimized for certain inhaled anesthetics and propofol, and other anesthetics or the effects of age have not been accounted for [8,9,10]. Additionally, each algorithm is proprietary and not disclosed, which limits comparative studies and further developments. Therefore, establishing new technologies for effective feature extraction of EEGs during GA using open algorithms is crucial for anesthesiologists who aim to provide safe and effective anesthesia.
Mode decomposition methods have emerged as a new approach to feature extraction, with research on these methods beginning at the end of the 20th century [11]. We have examined and assessed the usefulness of various mode decomposition methods for extracting features from EEG signals under GA [12,13]. The origin of mode decomposition is the empirical mode decomposition (EMD) method of Huang et al. [11], which analyzes nonstationary signals, and the use of EMD to analyze EEG waves during GA has been reported. However, EMD has the drawback that it decomposes data into intrinsic mode functions (IMFs) without band limitations, leading to frequency mixing issues and difficulties in result interpretation. In 2013, Dragomiretsky and Zosso proposed variational mode decomposition (VMD) [14], a method for overcoming the shortcomings of EMD. VMD views the decomposition of signals as a variational problem, solving a functional equation by applying mathematical optimization algorithms. This procedure continuously updates the central frequency of each mode online in the frequency domain, decomposing the signals into multiple narrowband IMFs. However, the VMD method heavily depends on its hyperparameters, namely, the mode decomposition number (K) and the penalty factor (PF), which affect the accuracy of capturing useful modality information features. Determination of the K value requires some prior knowledge of the original signal, and the bandwidth of each modality obtained through VMD is determined by the PF, with a larger PF narrowing the bandwidth. Implementing the VMD method for EEG feature extraction necessitates the careful selection of the acceptable values for these hyperparameters.
Recently, the grey wolf optimizer (GWO) algorithm, a new heuristic algorithm for optimization problems, was proposed to solve this issue [15,16,17]. The GWO algorithm, developed by Mirjalili [18], is a swarm intelligence optimization method that mimics the hunting behavior and leadership hierarchy of grey wolves, and it has been successful in many real-world applications [19]. In this study, we used envelope entropy as the fitness function for the GWO algorithm, which we used to optimize the VMD method by determining the number of K modes and the PF in EEG signals during GA.
The signals analyzed are complex and volatile, making it challenging to determine the two parameters influencing them. Therefore, selecting appropriate parameters is key to analyzing signal data through VMD. Attempts to optimize VMD using intelligent algorithms have been reported in various fields. It is essential to use an algorithm that is easy to implement, has an adaptive convergence coefficient and information feedback mechanisms, and balances local optimization with global search. The GWO algorithm, which mimics the hunting behavior and leadership hierarchy of grey wolves, has been successful in many real-world applications. In the same number of iterations, the GWO algorithm has been reported to converge significantly better than other algorithms. Therefore, in this study, we implemented the GWO algorithm in the EEG analysis process with the aim of optimizing VMD. To do this, we used envelope entropy as the fitness function of the GWO algorithm to optimize the number of K modes and the PF in EEG signals during GA. As the effects of GA drugs diminish, there are changes in the frequency components of EEG signals. In such dynamically changing EEG signals, it is assumed that the hyperparameters of VMD also need to change dynamically. Based on this hypothesis, we conducted an experiment where the hyperparameters of VMD were calculated in real time for each EEG analysis epoch to optimize them. Specifically, during the approximately 10 min from the maintenance of GA to awakening, we observed how VMD’s hyperparameters were determined via the GWO algorithm and how the associated VMD process was optimized. We also examined whether the extraction of EEG features using VMD could lead to the estimation of anesthesia depth through the acquisition of IMF components and the analysis of their characteristic frequency components. In this study, we report on the use of an optimized VMD method for decomposing EEG data recorded under GA, demonstrating how it can be used to acquire a series of IMF components and examining its utility.

2. Materials and Methods

2.1. Anesthesia Management and Data Acquisition

All experimental protocols involving human participants adhered to the principles outlined in the Declaration of Helsinki. The current study (No. ERB-C-1074) received approval from the Institutional Review Board (IRB) for human experiments at the Kyoto Prefectural University of Medicine (KPUM). For this noninterventional and noninvasive retrospective observational study, the need for informed patient consent was waived by the KPUM IRB. Nevertheless, patients were informed of their opportunity to opt out and were provided this opportunity during their preoperative anesthesia clinic visit. The EEG data recorded from three patients who were treated with total anesthesia using propofol (for 10 min before waking from GA) were used for analysis (Table S1). The EEG dataset used is available on the author’s GitHub site (https://github.com/teijisw/EEG_DataSet, accessed on 4 July 2024) [20].
Anesthesia was induced with small doses of fentanyl (1 µg kg−1 per dose) and a continuous intravenous infusion of remifentanil (0.125–0.25 µg kg−1 minute−1), followed by a single intravenous injection of propofol (2 mg kg−1) and rocuronium (0.8–1.0 mg kg−1). Tracheal intubation was performed, and anesthesia was maintained with propofol using a target-controlled infusion pump (Terufusion™ TIC pump TE-371, Terumo, Tokyo), aiming for a blood concentration of 3 µg/mL. Additional maintenance doses of fentanyl (1 µg/kg per dose) and rocuronium (0.2 mg/kg at 20–30-minute intervals), along with a continuous intravenous infusion of remifentanil (0.125–0.25 µg/kg per minute), were administered. To acquire EEG data, we utilized the EEG Analyzer software (ver. 54_GP; available at http://anesth-kpum.org/blog_ts/?p=3169 4 July 2024) [13,21]. Raw EEG signals were recorded as text files on a personal computer via the RS-232 interface of a VISTA A-3000 BIS monitor (32-bit raw EEG signals at a sampling frequency of 128 Hz, eight packets/s, VISTA Application revision 3.22, Platform revision 2.03, Medtronic, Minneapolis, MN, USA) with a BIS Quatro sensor mounted on the frontal regions, as previously reported [13].
VMD and GWO programming codes were constructed in the computer programming languages Python (ver.3.8) and Processing (ver.4.3), with the Apache Common Mathematics Library (https://commons.apache.org/proper/commons-math/, 4 July 2024, ver.3.6.1), Digital Signal Processing in Java (JDSP ver.0.5.0, https://jdsp.dev 4 July 2024), and FloatTable.pde from Example for Visualizing Data by Ben Fry (https://benfry.com/writing/archives/3/ 4 July 2024), as shown in Supplementary Program Codes #1 and #2.

2.2. Algorithm for VMD and the Hilbert Transform

VMD divides an input signal f into k narrowband IMFs uk characterized by a discrete central frequency ωk with a particular sparsity property [14]. The following three mathematical steps are involved in the VMD algorithm [14]:
1.
The first step computes the analytic signal associated with each IMF according to the Hilbert transform (Formula (1)).
    δ t + j π t u k t
where * is the convolution, uk(t) is the kth IMF, and δ(t) is the pulse signal.
2.
The second step involves demodulating the analytic signal to the baseband by multiplying it with an exponential function, which is adjusted to the respective estimated central frequency ω k .
t δ t + j π t × u k t e j ω k t
3.
The third step estimates the bandwidth by using the H 1 Gaussian smoothness of the demodulated signal, which involves calculating the squared L2 norm of the slope of the demodulated analytic signal. It then computes the associated analytic signal using the Hilbert transform to obtain the one-sided frequency spectrum of each intrinsic mode uk.
The constrained variational evaluation of a given signal x(t) is described as
min u k , ω k k t δ t + j π t × u k t e j ω k t 2 2 ,   s . t . k u k = f ( t ) ,
where uk and ωk are shorthand notations for the sets {u1, ……, uk} and {ω1, ……, ωk}, representing all modes and their central frequencies, respectively. The sum k : = k = 1 K of all mode functions uk is equal to the specified time series signal f(t).
This constrained variational problem is addressed by incorporating a quadratic penalty term and introducing the Lagrange multiplier, λ, into the mathematical optimization method. This methodology effectively transforms the problem into one that involves an unconstrained extended Lagrangian function, denoted as L in Formula (4).
L u k ,   ω k ,   λ   α   k t δ t + j π t × u k t e j ω k t 2 2 + f t k u k t 2 2 + λ t ,   f t k u k t
The solution to this minimization problem is found as a saddle point of the augmented Lagrangian L in a series of iterative partial optimizations called the alternating direction method of multipliers (ADMM, Algorithm 1).
Algorithm 1. Optimization Concept for VMD [14].
initialization u ^ k 1 ,   ω ^ k 1 ,   λ ^ 1 ,   n   0
  repeat
    n n + 1
   for k = 1 : K   d o
    Update u ^ k for all ω   0 :  
u ^ k n + 1 ( ω )     f ^ ω i < k u ^ i n + 1 ω i > k u ^ i n ω + λ ^ ( ω ) 2 1 + 2 α ( ω ω k n ) 2 (5)
    Update ω ^ k :
ω k n + 1 = 0 ω u ^ k k n + 1 ω 2 d ω 0 u ^ k k n + 1 ω 2 d ω (6)
  end for
  Dual ascent for all ω   0 :
λ ^ n + 1     λ ^ n + τ f ^ ( ω ) k u ^ k n + 1   ( ω )         (7)
until convergence:
k u ^ k n + 1 u ^ k n   2 2 u ^ k n   2 2 < ϵ (8)
For an EEG signal x(t), VMD decomposes it into a series of IMFs, Cn (where n = 1, 2, …, N), with N being the total number of IMFs. After VMD, the signal x(t) can be expressed as
x t = i = 1 N i m f ( t ) u
Applying the Hilbert transform to the IMF components, it follows that
Z t = i m f t + i H i m f t = a t e i ω t d t
in which
a t = i m f 2 t + H 2 [ imf t ]
ω t = d d t a r t a n H i m f t i m f t
h ω = H ω ,   t d t
where a(t) is the amplitude of the IMF and ω (t) is the instantaneous frequency for obtaining a time–frequency distribution for signal x(t) and the Hilbert amplitude spectrum H(x, t).

2.3. Grey Wolf Optimizer

In the original VMD method, users need to provide the number of K modes and the PF in advance, although there is no guarantee whether these parameters are optimal or whether they will achieve satisfactory results [21,22]. The smaller the PF, the larger the bandwidth of each IMF component. Furthermore, improper selection of the number of K modes can directly lead to signaling component distortion and unrealistic denoising effects. In practice, the two influencing parameters are usually difficult to determine because the actual signal to be analyzed is complex and changeable. Thus, selecting suitable parameters is key to the VMD method. Next, we introduce the principle of the GWO algorithm (Figure 1) and its convergence and the steps of the parameter-optimized VMD method.
The GWO is an optimization algorithm simulating the social hierarchical relationships and hunting behaviors of grey wolves in nature [18]. To mathematically model the social hierarchy of wolves, the population (a candidate group of hyperparameter values to be optimized) corresponding to a pack of wolves is divided into four social classes, and the optimal solution is considered to be an alpha (α) wolf. The second and third best solutions are treated as the beta (β) and delta (δ) wolves, respectively, and the remaining candidate solutions in the swarm are allocated to the omega (ω) wolves. In the GWO algorithm, hunting (optimization) is guided by the top three solutions—α, β, and δ—with the ω wolf following these three wolves. Although it is unclear during the analysis where the optimal solution will be, it is assumed that a possible optimal solution is surrounded by three optimal solutions and that the movement of the ω wolf population approaching the optimal solutions (α, β, and δ) is randomly controlled. We adopted a method that approaches the optimal solution by incorporating and repeating variables, introducing a fitness function, and sequentially replacing the α, β, and δ wolves in the population.
Convergence of the GWO algorithm: Grey wolves surround their prey when hunting. To mathematically model this enveloping behavior, the following equation is proposed (Figure 2).
D = C · X p ( i ) X ( i )
  X ( i + 1 ) = X p ( i ) A · D
where i indicates the current iteration, A and C are the coefficient vectors, X p is the position vector of the prey, and X indicates the position vector of a grey wolf.
The vectors A and C are calculated as follows:
A = a · 2 · r 1 1
C = 2 · r 2
Here, the components of a in Formula (16) decrease linearly from 2 to 0 throughout the iterations, and r 1 and r 2 in Formulas (16) and (17) are random vectors in [0, 1]. The coefficient vector A of Formula (15) is related to the distance D between each wolf X ( i ) and the prey X p . However, the position of the prey X p itself is also multiplied by a coefficient vector C within a random range of [0, 2] in Formula (14), providing random weight to the prey by either emphasizing ( C > 1) or de-emphasizing ( C < 1) it. This allows GWO to exhibit more random behavior throughout the optimization, facilitating exploration and avoidance of local optima. The random variation in A , which initially ranges from −2 to +2, decreases gradually with each loop, forcing the search agents also to move away from the prey. This emphasis on exploration enables the GWO algorithm to search globally. To mathematically simulate the hunting behavior of grey wolves, we need to know the location of the prey. However, we cannot actually know it directly. Therefore, GWO assumes that α, β, and δ are related to the potential location of the prey. In other words, it is assumed that the prey position (optimal solution) is surrounded by three wolves, α, β, and δ. The first three best solutions obtained thus far are preserved, and it is mandated that the other search agents (including ω) update their positions according to the positions of the best search agents. In this regard, the following formula is proposed. In this case, the positions of wolves α, β, and δ are used instead of the prey position ( X p ) in Formulas (14) and (15), and the ω wolf position is updated in three ways using the following formula.
D α = C 1 · X α X ,   D β = C 2 · X β X ,   D δ = C 3 · X δ X
X 1 = X α A 1 · D α ,   X 2 = X β A 2 · D β ,   X 1 = X α A 1 · D α
The center of the three possible updated positions is selected as a new ω wolf update position candidate, as follows:
X ( i + 1 ) = X 1 + X 2 + X 3 3
However, the position will be updated only when the fitness function shows better fitting, and after updating the positions of the wolves in the entire pack using this method, the top three wolves showing the best values will be updated using the fitness function. By repeating the task of specifying the α, β, and δ wolves a specified number of times, we hope that most of the wolves will converge to the optimal solution.
When using these equations for VMD, the search agent updates its position in the two-dimensional search space for K and PF according to α, β, and δ (Figure 3). Furthermore, the final position is a random location within a circle defined by the positions of α, β, and δ in the search space. In other words, α, β, and δ estimate the position of the prey, and other wolves randomly update their positions around the prey. Through the fitness function, the positions of the wolves ranked again are used to update α, β, and δ, and by repeating this loop, the solution offering an optimized fitness function (a combination of optimal solutions) is progressively refined.

2.4. The Fitness Function: Envelope Entropy

In a situation where only one of several events occurs, and each can occur with equal probability P, the information I gained upon learning which event actually occurred can be expressed as
I = l o g 2 1 P = l o g 2 P
Therefore, when this process is repeated i times and the probability of an event occurring is P i , the information obtained when it occurs is l o g 2 P i . From this, the average amount of information obtained is
I = i P i × l o g 2 P i  
The information entropy S before the outcome is known can also be determined using the same formula.
S = i P i × l o g 2 P i  
Here, if B i   is the envelope signal of x(i), then B i can be normalized to obtain the probability P i as follows:
P i = B i i = 1 N B ( i )
Then, by applying the formula for information entropy to the envelope, we obtain the envelope entropy as follows:
E p = i = 1 N P i × l o g 2 P i  
In the GWO algorithm, the envelope entropy of the IMF decomposed via VMD is determined, and K and PF are optimized such that the envelope entropy values of the separated IMFs are reduced (indicating less information). Here, by calculating the envelope entropy of each of the multiple IMFs obtained, the IMF with the highest envelope entropy (the most information, using the maximum value method) is identified. Then, by ensuring that these values become lower (more realistic information) through new combinations of K and PF, it is possible to optimize the separated IMFs overall.
An example using synthetic cosine waves: The following is a synthetic function represented by a combination of three cosine waves with different known frequencies (2 Hz, 12 Hz, and 36 Hz).
y = c o s ( 2 × 2 π t ) + 1 2 c o s ( 12 × 2 π t ) + 1 4 c o s ( 36 × 2 π t )
where t = [1/1024, 2/1024, …, 1024/1024] (Figure 3A, left). Here, we used this synthetic function to demonstrate an example of the GWO algorithm being used to optimize K and PF. The envelope entropy of this synthetic function was −9.9767506 (Figure 3A, right). If this synthetic function is decomposed with maximum efficiency, it should result in decomposition into cosine waves y = c o s ( 2 × 2 π t ) , y = 1 2 c o s ( 12 × 2 π t ) , and y = 1 4 c o s ( 36 × 2 π t ) (Figure 3B, left). Each envelope of the efficiently decomposed components will become a flat line, and the envelope entropy will be −10.0 (Figure 3B, right). This trigonometric function’s composite waveform was analyzed using the GWO algorithm with 20 wolves and 20 iterations, setting K to range from 2 to 6 and PF from 10 to 5000 (randomization in increments of 1 for K and 1 for PF). The envelope entropy of each IMF was calculated, and K and PF were optimized using both the mean value method and the maximum value method. As a result, the GWO algorithm determined the optimal values of K = 3 and PF = 4984. The envelope entropies of the three IMFs were −9.9999944, −9.9999929 (max), and −9.9999958 (Figure 3C). With settings of K = 3 and PF = 10 (Figure 3D) or K = 3 and PF = 70 (Figure 3E), the IMFs were superimposed on the other cosine waves, indicating that complete frequency separation was not achieved.
Then, with PF fixed at the optimized value of 4984 and K set to 2, the application resulted in two frequency components being present in IMF-1, and the envelopes of the IMFs were not flat (Figure 3F). Finally, with PF fixed at the optimized value of 4984 and K set to 4, when VMD was applied, the envelopes of all the IMFs were flat (Figure 3G). However, although the highest envelope entropy was obtained with IMF-4 (−9.9997846), this value was larger than that of IMF-1 (−9.9999944) obtained with GWO, suggesting that GWO with K = 3 produced a better model than VMD with K = 4.

3. Results

3.1. Convergence Status of K, PF, and Fitness Values in GWO

First, the GWO algorithm optimization process (20 wolves, 20 optimization loops) for the VMD hyperparameters K and PF was performed using one epoch (8 s, 1024 data points) during the maintenance phase of the GA (induced by continuous intravenous propofol administration) in three patients (Table S1). This involved monitoring the position of each wolf during two-dimensional hunting.
In Patient #1, with 20 optimization loops, convergence to K = 2 occurred at the ninth loop, except for three ω wolves (Figure 4(1)). However, the PF values of the leader wolves, α, β, and δ, converged to 3640, 3620, and 3470, respectively, while the convergence for the ω wolves was unstable. Nonetheless, the fitness function converged to approximately −9.87 to −9.88, indicating that the PF values did not significantly impact the fitness values. In the case of Patient #2 (Figure 4(2)), except for the six ω wolves, convergence to K = 2 occurred in the sixth loop. Similar to Patient #1, the PF values for the α, β, and δ leader wolves were 4900, 4800, and 4680, respectively. The fitness values for these leader wolves had already converged to approximately −9.85 by the second loop, but the fitness function for the ω wolves had not fully converged. In Patient #3 (Figure 4(3)), except for the five ω wolves, there was convergence with K = 3. The PF values of the α, β, and δ leader wolves were 1580, 1830, and 1850, respectively, whereas those of the ω wolves tended to converge to values below 2500. By approximately the fourth loop, the fitness values converged to approximately −9.72 to −9.73. Table 1 shows the convergence status for K, PF, and fitness, with the final convergence value set as 100%. For K, convergence occurred between 6 and 13 loops; for PF and fitness, it was between 9 and 18 loops. Although 20 loops were conducted in this case, sufficient convergence had been achieved by then.
The results show that convergence of the K values was achieved by the α, β, and δ leader wolves at an early stage. However, since the impact on the fitness function was much greater for the K values than for the PF values, it was observed that the convergence of the PF values tended to be delayed. Additionally, when the fitness function converged, the impact of the PF values was considered low.

3.2. Temporal Changes in the K, PF, and Fitness Values According to the GWO

Next, for the three cases of GA induced by continuous intravenous propofol administration, the optimal solutions for K and PF were determined through fitness function optimization using the GWO algorithm (20 wolves, 20 optimization loops) for all 73 epochs during the 10 min before and after awakening from GA. Each epoch of EEG data lasted 8 s (128 Hz, 1024 data points).
K was set in the range of 2 to 6 (in increments of 1), and the PF values ranged from 10 to 5000 (in increments of 10). K and PF fluctuated significantly, with K oscillating from a base of 2 up to 3, 4, 5, and 6, while the PF values varied greatly between 10 and 5000 for each epoch (Figure 5A). The Fourier spectrogram for the 10-minute interval was obtained using the multitaper method and is displayed in Figure 5B. Figure 5C shows a Hilbert spectrogram of the IMF of the EEG from approximately the last 10 min before awakening from GA, for K = 2, PF = 2000, and K = 3, PF = 2000. The envelope entropy, excluding the variations after awakening, ranged from approximately −9.8 (Figure 5A). The fitness values derived from the envelope entropy were relatively stable and consistent within a single epoch before awakening. When the fitness values showed a convergence trend, the impact of the PF values on the fitness function appeared low, suggesting that prioritizing the convergence of K would be advisable. After awakening, there were significant increases in fitness, likely influenced by the incorporation of electromyography during the postemergence period, such as orbicularis oculi muscle activity.

3.3. VMD Analysis under GWO Support during the Three Phases of General Anesthesia

Finally, we applied the analysis to 8 s EEG segments during the (1) maintenance, (2) transition, and (3) emergence phases of GA induced by propofol (Figure 6A). The power spectrum obtained through Fourier analysis is shown in Figure 6B. The analysis was conducted with 20 wolves, 20 repetitions, and settings of K = 2 to 6 (in increments of 1) and PF = 10 to 5000 (in increments of 10). The results show a convergence to K = 2 for the 8 s EEG in all three phases. Figure 6C presents the Hilbert spectrum, and Figure 6D shows the IMFs of the EEG after applying VMD with the optimized K and PF values. The Fourier spectrogram (multitaper method) for 64 s, including the aforementioned 8 s EEG segment, is shown in Figure 6E. While IMF-1 was almost fixed in the delta wave region, the peak values of the Hilbert spectrum for IMF-2 showed variation across the phases of GA. IMF-2 showed enhanced α waves at approximately 10–12 Hz under propofol-induced GA, as shown in Figure 6E, and shifted forward and backward with the transition to emergence from anesthesia.
Table 2 statistically examines the fluctuations in K, PF, and fitness for each patient during the phase from the anesthetized to awake states, as shown in Figure 5-A, by comparing the mean values before and after awakening. The results indicate significant differences in the K values for Patient #2 and the combined data from three patients but not for the other K or PF values. However, significant differences were detected in the fitness values for all three patients between the means before and after awakening. This suggests that while the K and PF values optimized by GWO are less likely to serve as indicators of the difference between the anesthetized and awake states, the fitness values, which are outputs of the envelope entropy function, may significantly increase during the transition from anesthesia to emergence.
The reasoning here is that if there is no large variation in the envelope entropy as a fitness function, it is acceptable to decompose using a lower K value. A K value of approximately 2 to 3 should suffice for the VMD of EEG signals during GA. Therefore, in the three case studies, K was fixed at 2, PF was fixed at 2000, and PF was optimized to obtain the 10-minute Hilbert spectrum and spectrogram (Figure 6C). It was observed that the α wave enhancement due to propofol-induced GA, as shown in the power spectrogram, corresponded with the separation of the narrow band of α waves represented by IMF-2. We understood that this characteristic α wave enhancement, a feature of propofol-induced GA, could be distinctively extracted in the IMF-2 of VMD.

4. Discussion

GA involves three effects: hypnosis, analgesia, and immobility. As anesthesiologists, we want to quantify the hypnotic state of patients; however, there are no other biomarker signals that can quantitatively capture the state of hypnosis apart from EEG, including some evoked potentials. Therefore, there is a focus on quantifying the hypnotic state induced by anesthetic drugs through frequency and amplitude changes in EEGs. Apart from the standard Fourier transform-based power spectral analysis for analyzing the frequency components of EEGs, methods combining mode decomposition with the Hilbert transform, such as the Hilbert–Huang transform method and its variants, have been reported. A disadvantage of Fourier transform-based methods is the inversely proportional relationship between the frequency resolution and temporal resolution, which is affected by the analysis window but is not concerned with the Hilbert transform. As an alternative approach, frequency analysis of EEGs using the continuous wavelet transform has been reported [23]. However, the wavelet transform has drawbacks, such as variations in time and frequency resolution depending on the scale of the basis function and the need to appropriately select the wavelet type and parameters (such as the number of scales and the range of analysis), which requires extensive experience and trial and error.
In this respect, the Hilbert–Huang transform method has the advantage of extracting narrowband signals through mode decomposition. Mode decomposition originates from the Hilbert–Huang transform method combined with the EMD method reported by Huang [11]. There have been several attempts to apply the EMD method to the frequency analysis of EEGs during GA [12,24,25]. Subsequently, methods such as VMD [14] and the empirical wavelet transform (EWT) have been reported [26], along with their combination with the Hilbert transform in variations of the Hilbert–Huang transform. We reported using VMD to analyze EEG signals during GA in 2023 [13]. The EMD method, based on spline interpolation for envelope calculation, is mathematically less robust and is known to be incapable of separating signals into narrow bands. VMD, supported by a more robust mathematical theory based on mathematical optimization, is more effective at separating EEGs into narrower bands and extracting features of their frequency components. However, the necessity for user-defined hyperparameters in VMD requires some form of guidance, which was the theme of this study. It would be interesting to conduct comparative studies to determine which mode decomposition method—EMD, VMD, or even EWT—is superior for analyzing EEGs during GA. However, it is important to emphasize that VMD is dependent on hyperparameters, including the number of decompositions, which is an integral part of the algorithm. Therefore, before comparing different mode decomposition methods, it was necessary for us to introduce the optimization of hyperparameter determination in the current study.
In this analysis, the number of wolves was set to 20 because this number is generally considered appropriate. Increasing the number of wolves might make finding an optimized solution more straightforward, but it can also delay the wolves’ overall convergence (hunting). When we used 20 loops, it was uncertain whether increasing the loop count would lead to overall wolf convergence (hunting). In our situation, convergence of the fitness function is important. Regarding the VMD hyperparameters, the optimization of the K value is crucial, while the PF value, if it has little impact on fitness, may not converge across all wolves and may show variability. However, if the PF value does not affect the fitness value, it might not be necessary to focus on it to a great extent.
VMD itself involves computer analysis based on mathematical optimization theory. Additionally, the K and PF values are selected randomly to calculate the envelope entropy of each IMF obtained through VMD as a fitness value. A single optimization requires the VMD algorithm to be run 800 to 2400 times for K values ranging from 2 to 6, with 20 wolves and 20 loops. This presents a substantial challenge to the speed of the programming language and the computational speed of personal computers, with the analysis speed of the Processing application in Java being insufficient to implement the GWO algorithm in real time. It is, therefore, necessary to consider application development using programming languages with default GPU computation support, such as Julia or Mojo, to achieve more efficient processing.
In this study, we applied the GWO algorithm to optimize the hyperparameters K and PF in VMD, using the envelope entropies of the decomposed IMFs as fitness functions to analyze EEG signals during GA. Our results show that although K and PF varied significantly for each epoch, the envelope entropies remained stable. This implies that in terms of the K value, when the main frequency components for clinical observation are considered, two bands—δ to θ waves and the alpha frequency band, including sleep spindles—are critical during the maintenance phase of GA, suggesting that a K value of 2 or 3 might be sufficient. Subsequently, we considered optimizing the PF value using the GWO algorithm as needed. This study demonstrated the effectiveness of GWO for optimizing hyperparameters in the VMD analysis of biological signals and highlights its promise for future applications.

5. Conclusions

In this study, we used the GWO to determine hyperparameters for the VMD method in order to analyze changes in the frequency components of EEG signals during GA. We used the envelope entropy of the IMF decomposed via VMD for the fitness function for optimizing the GWO algorithm. Regarding the decomposition number K for VMD, when a limit of 2 to 6 was set for clinical significance, K varied from 2 to 6 for each epoch of the EEG being analyzed, but, in general, K = 2 or 3, and we believe that the characteristics of the EEG during GA were captured. Regarding the PF value, if the fitness function is optimized to a low value by determining an appropriate K value, the low convergence of the PF value means that different PF values within an allowable range of 10–5000 will not have a significant impact. The determination of appropriate values of VMD hyperparameters by adapting the GWO algorithm provides important findings for the analysis of EEG signals using VMD during GA.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/s24175749/s1, Table S1: Patient characteristics and surgical and anesthesia times.; Program Code 1: Python for GWO; Program Code 2: Processing for GWO.

Author Contributions

Conceptualization, T.S.; methodology, T.S.; software, T.S.; validation, M.K., Y.O. and K.A.; formal analysis, T.S.; investigation, K.K., Y.O. and T.Y.; data curation, K.K., Y.O. and T.Y.; writing—original draft preparation, T.S.; writing—review and editing, Y.O., M.K. and K.A.; supervision, T.S.; project administration, T.S. All authors have read and agreed to the published version of the manuscript.

Funding

The authors performed this study with the support of internal funding (2023 Fiscal Year) from the Department of Anesthesiology, Kyoto Prefectural University of Medicine.

Institutional Review Board Statement

The EEG data used for the spectral analyses in this study were obtained from anesthetized patients with ethical approval (No. ERB-C-1074-2) from the Institutional Review Board for Human Experiments at the Kyoto Prefectural University of Medicine (IRB of KPUM).

Informed Consent Statement

For this noninterventional and noninvasive retrospective observational study including minors, the requirement for informed patient consent was waived by the IRB of KPUM; however, patients were provided with an opt-out option, of which they were notified in the preoperative anesthesia clinic.

Data Availability Statement

The EEG dataset used is available on the author’s GitHub site (https://github.com/teijisw/EEG_DataSet 4 July 2024) [18]. The programming codes used in the analysis in this paper are available as Supplementary Data.

Acknowledgments

The authors thank Karl Embleton and F Rosas, Edanz (https://jp.edanz.com/ac, 22 March 2024) for editing a draft of this manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Laferriere-Langlois, P.; Morisson, L.; Jeffries, S.; Duclos, C.; Espitalier, F.; Richebe, P. Depth of anesthesia and nociception monitoring: Current state and vision for 2050. Anesth. Analg. 2024, 138, 295–307. [Google Scholar] [CrossRef]
  2. Roche, D.; Mahon, P. Depth of anesthesia monitoring. Anesth. Clin. 2021, 39, 477–492. [Google Scholar] [CrossRef] [PubMed]
  3. Shander, A.; Lobel, G.P.; Mathews, D.M. Brain monitoring and the depth of anesthesia: Another goldilocks dilemma. Anesth. Analg. 2018, 126, 705–709. [Google Scholar] [CrossRef]
  4. Chaddad, A.; Wu, Y.; Kateb, R.; Bourid, A. Electroencephalography signal processing: A comprehensive review and analysis of methods and techniques. Sensors 2023, 23, 6434. [Google Scholar] [CrossRef] [PubMed]
  5. Sharma, R.; Meena, H.K. Emerging trends in EEG signal processing: A systematic review. SN Comput. Sci. 2024, 5, 415. [Google Scholar] [CrossRef]
  6. Ma, Z.; Hua, H.; You, C.; Ma, Z.; Guo, W.; Yang, X.; Qiu, S.; Zhao, N.; Zhang, Y.; Ho, D.; et al. FlexiPulse: A machine-learning-enabled flexible pulse sensor for cardiovascular disease diagnostics. Cell Rep. Phys. Sci. 2023, 4, 101690. [Google Scholar] [CrossRef]
  7. Velcescu, A.; Lindley, A.; Cursio, C.; Krachunov, S.; Beach, C.; Brown, C.A.; Jones, A.K.P.; Casson, A.J. Flexible 3D-printed EEG electrodes. Sensors 2019, 19, 1650. [Google Scholar] [CrossRef]
  8. Sun, Y.; Wei, C.; Cui, V.; Xiu, M.; Wu, A. Electroencephalography: Clinical applications during the perioperative period. Front. Med. 2020, 7, 251. [Google Scholar] [CrossRef]
  9. Yuan, I.; Xu, T.; Kurth, C.D. Using electroencephalography (EEG) to guide propofol and sevoflurane dosing in pediatric anesthesia. Anesth. Clin. 2020, 38, 709–725. [Google Scholar] [CrossRef]
  10. Schultz, B.; Schultz, M.; Boehne, M.; Dennhardt, N. EEG monitoring during anesthesia in children aged 0 to 18 months: Amplitude-integrated EEG and age effects. BMC Pediatr. 2022, 22, 156. [Google Scholar] [CrossRef]
  11. Huang, N.; Shen, Z.; Long, S.; Wu, M.; Shih, H.; Zheng, Q.; Yen, N.-C.; Tung, C.; Liu, H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  12. Obata, Y.; Yamada, T.; Akiyama, K.; Sawa, T. 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. BMC Anesth. 2023, 23, 125. [Google Scholar] [CrossRef]
  13. Yamada, T.; Obata, Y.; Sudo, K.; Kinoshita, M.; Naito, Y.; Sawa, T. Changes in EEG frequency characteristics during sevoflurane general anesthesia: Feature extraction by variational mode decomposition. J. Clin. Monit. Comput. 2023, 37, 1179–1192. [Google Scholar] [CrossRef]
  14. Dragomiretskiy, K.; Zosso, D. Variational mode decomposition. IEEE Trans. Signal Process. 2014, 62, 531–544. [Google Scholar] [CrossRef]
  15. Wang, Y. An adaptive variational mode decomposition technique with differential evolution algorithm and its application analysis. Shock Vib. 2021, 2030, 128. [Google Scholar] [CrossRef]
  16. Hou, Y.; Gao, H.; Wang, J.; Du, C. Improved grey wolf optimization algorithm and application. Sensors 2022, 22, 318. [Google Scholar] [CrossRef]
  17. Zang, X.; Li, D.; Li, J.; Li, Y. Grey wolf optimization-based variational mode decomposition for magnetotelluric data combined with detrended fluctuation analysis. Acta Genophys. 2022, 2022, 111–120. [Google Scholar] [CrossRef]
  18. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey wolf optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef]
  19. Nadweh, S.; Khaddam, O.; Hayeh, G.; Atieh, B.; Haes Alhelou, H. Steady state analysis of modern industrial variable speed drive systems using controllers adjusted via grey wolf algorithm & particle swarm optimization. Heliyon 2020, 6, e05438. [Google Scholar] [CrossRef] [PubMed]
  20. Sawa, T. teijisw/EEG_DataSet. Available online: https://github.com/teijisw/EEG_DataSet/tree/master/general_anesth_propofol (accessed on 6 July 2024).
  21. Sivavaraprasad, G.; Padmaja, R.S.; Ratnam, D.V. Mitigation of ionospheric scintillation effects on GNSS signals using variational mode decomposition. IEEE Geosci. Remote Sens. Lett. 2017, 14, 389–392. [Google Scholar] [CrossRef]
  22. Xu, Z.; Qin, C.; Tang, G. A novel deconvolution cascaded variational mode decomposition for weak bearing fault detection with unknown signal transmission path. IEEE Sens. J. 2021, 21, 1746–1755. [Google Scholar] [CrossRef]
  23. Mousavi, S.M.; Adamoğlu, A.; Demiralp, T.; Shayesteh, M.G. A wavelet transform based method to determine depth of anesthesia to prevent awareness during general anesthesia. Comput. Math. Methods Med. 2014, 2014, 354739. [Google Scholar] [CrossRef] [PubMed]
  24. Li, X.; Li, D.; Liang, Z.; Voss, L.J.; Sleigh, J.W. Analysis of depth of anesthesia with Hilbert-Huang spectral entropy. Clin. Neurophysiol. 2008, 119, 2465–2475. [Google Scholar] [CrossRef]
  25. Liu, Q.; Ma, L.; Fan, S.Z.; Abbod, M.F.; Ai, Q.; Chen, K.; Shieh, J.S. Frontal EEG temporal and spectral dynamics similarity analysis between propofol and desflurane induced anesthesia using Hilbert-Huang Ttansform. Biomed. Res. Int. 2018, 2018, 4939480. [Google Scholar] [CrossRef]
  26. Gilles, J. Empirical wavelet transform. IEEE Trans. Signal Process. 2013, 61, 3999–4010. [Google Scholar] [CrossRef]
Figure 1. Flowchart of grey wolf optimization (GWO) for variational mode decomposition (VMD).
Figure 1. Flowchart of grey wolf optimization (GWO) for variational mode decomposition (VMD).
Sensors 24 05749 g001
Figure 2. Position updating algorithm for the wolves in the grey wolf optimizer. Three leader/subleader wolves, α, β, and δ, surround the prey, and the position of the prey is inferred from the positions of these three leader wolves, assuming that they surround it. Vectors A and C are coefficient vectors and are calculated for each coordinate. The other wolves in the pack, ω, locate the leader wolves’ positions, which are initially adjusted with a coefficient C and then gradually adjusted in each loop to better approximate the leader wolves’ positions. A random coefficient D is then applied to their distance, allowing them to approach the leader wolves within a range of −1 to +1. As a result, wolves that are closer to the prey than the leader wolves may replace them as the new leaders, enabling the pack to surround the prey more closely. The algorithm specifies three leader wolves and divides the obtained average position by the number of leaders. If the positions of α, β, and δ are 6, 8, and 3, respectively, then the position of the prey would be the midpoint, calculated as (6 + 8 + 3)/3 = 5.7, and the integer 5 becomes the updated position for the prey. The new position of the ω1 wolf is adjusted according to the leaders’ positions, taking into account its current position.
Figure 2. Position updating algorithm for the wolves in the grey wolf optimizer. Three leader/subleader wolves, α, β, and δ, surround the prey, and the position of the prey is inferred from the positions of these three leader wolves, assuming that they surround it. Vectors A and C are coefficient vectors and are calculated for each coordinate. The other wolves in the pack, ω, locate the leader wolves’ positions, which are initially adjusted with a coefficient C and then gradually adjusted in each loop to better approximate the leader wolves’ positions. A random coefficient D is then applied to their distance, allowing them to approach the leader wolves within a range of −1 to +1. As a result, wolves that are closer to the prey than the leader wolves may replace them as the new leaders, enabling the pack to surround the prey more closely. The algorithm specifies three leader wolves and divides the obtained average position by the number of leaders. If the positions of α, β, and δ are 6, 8, and 3, respectively, then the position of the prey would be the midpoint, calculated as (6 + 8 + 3)/3 = 5.7, and the integer 5 becomes the updated position for the prey. The new position of the ω1 wolf is adjusted according to the leaders’ positions, taking into account its current position.
Sensors 24 05749 g002
Figure 3. The relationship between the hyperparameters K and PF of VMD and the envelope entropy of each IMF. VMD is applied to a synthetic cosine wave composed of three known frequency components (2, 12, and 36 Hz). y = c o s ( 2 × 2 π t ) + 1 2 c o s ( 12 × 2 π t ) + 1 4 c o s ( 36 × 2 π t ) (1024 data points). (A) The synthetic cosine wave, its envelope function, and the envelope entropy. (B) The three cosine waves, their envelope functions, and envelope entropies. (C) Decomposition of the synthetic cosine wave into three IMFs using the GWO algorithm with optimized values of K = 3 and PF = 4984, and calculation of each IMF’s envelope function and envelope entropy. (D) Application of VMD with K = 3 and PF = 10, and determination of IMFs, their envelope functions, and envelope entropies. (E) Application of VMD with K = 3 and PF = 70, and analysis of IMFs and their envelope functions, and envelope entropies. (F) Application of VMD with K = 2 and PF = 4984, and evaluation of the IMFs and their envelope functions, and envelope entropy. (G) Application of VMD with K = 4 and PF = 4984, and analysis of the IMFs and their envelope functions, and envelope entropies. K: decomposition number; PF: penalty factor; fitness: fitness value of VMD = the maximum envelope entropy of the IMF.
Figure 3. The relationship between the hyperparameters K and PF of VMD and the envelope entropy of each IMF. VMD is applied to a synthetic cosine wave composed of three known frequency components (2, 12, and 36 Hz). y = c o s ( 2 × 2 π t ) + 1 2 c o s ( 12 × 2 π t ) + 1 4 c o s ( 36 × 2 π t ) (1024 data points). (A) The synthetic cosine wave, its envelope function, and the envelope entropy. (B) The three cosine waves, their envelope functions, and envelope entropies. (C) Decomposition of the synthetic cosine wave into three IMFs using the GWO algorithm with optimized values of K = 3 and PF = 4984, and calculation of each IMF’s envelope function and envelope entropy. (D) Application of VMD with K = 3 and PF = 10, and determination of IMFs, their envelope functions, and envelope entropies. (E) Application of VMD with K = 3 and PF = 70, and analysis of IMFs and their envelope functions, and envelope entropies. (F) Application of VMD with K = 2 and PF = 4984, and evaluation of the IMFs and their envelope functions, and envelope entropy. (G) Application of VMD with K = 4 and PF = 4984, and analysis of the IMFs and their envelope functions, and envelope entropies. K: decomposition number; PF: penalty factor; fitness: fitness value of VMD = the maximum envelope entropy of the IMF.
Sensors 24 05749 g003
Figure 4. Optimization process for the hyperparameters K and PF of VMD using the GWO algorithm. Twenty wolves and 20 optimization loops, observed over one epoch (8 s, 1024 data points) of EEG during the maintenance phase of GA (obtained through continuous intravenous propofol administration) in three patients. This process involved monitoring the position of each wolf during two-dimensional wolf hunting. (A) EEG (8 s). (B) The optimization of K. (C) The optimization of PF. (D) The optimization of fitness. α wolf: red line; β wolf: orange line; δ wolf: yellow line; ω wolves: black lines. K: decomposition number; PF: penalty factor; fitness: fitness value of VMD = the maximum envelope entropy of the IMF.
Figure 4. Optimization process for the hyperparameters K and PF of VMD using the GWO algorithm. Twenty wolves and 20 optimization loops, observed over one epoch (8 s, 1024 data points) of EEG during the maintenance phase of GA (obtained through continuous intravenous propofol administration) in three patients. This process involved monitoring the position of each wolf during two-dimensional wolf hunting. (A) EEG (8 s). (B) The optimization of K. (C) The optimization of PF. (D) The optimization of fitness. α wolf: red line; β wolf: orange line; δ wolf: yellow line; ω wolves: black lines. K: decomposition number; PF: penalty factor; fitness: fitness value of VMD = the maximum envelope entropy of the IMF.
Sensors 24 05749 g004
Figure 5. Grey wolf optimization in the VMD of EEGs obtained from three patients of GA induced by continuous intravenous propofol. (A) The optimal solutions for K and PF were determined through fitness function optimization using the GWO algorithm. Twenty wolves and 20 optimization loops for all 73 epochs. Each epoch lasted for 8 s (128 Hz, 1024 data points) of the 10-minute period before and after awakening from GA. Additionally, K was set in the range of 2 to 6 (in increments of 1), and the PF values ranged from 10 to 5000 (in increments of 10). (B) The spectrogram for the 10-minute interval was obtained using the multitaper method and is displayed in Figure 4B. (C) Using VMD with K = 2 or 3 and PF = 2000, the signal was decomposed into IMFs and a Hilbert spectrogram was obtained. The gradient color bar shows the relative power value of the signal. K: decomposition number; PF: penalty factor; fitness: fitness value of VMD = the maximum envelope entropy of the IMF.
Figure 5. Grey wolf optimization in the VMD of EEGs obtained from three patients of GA induced by continuous intravenous propofol. (A) The optimal solutions for K and PF were determined through fitness function optimization using the GWO algorithm. Twenty wolves and 20 optimization loops for all 73 epochs. Each epoch lasted for 8 s (128 Hz, 1024 data points) of the 10-minute period before and after awakening from GA. Additionally, K was set in the range of 2 to 6 (in increments of 1), and the PF values ranged from 10 to 5000 (in increments of 10). (B) The spectrogram for the 10-minute interval was obtained using the multitaper method and is displayed in Figure 4B. (C) Using VMD with K = 2 or 3 and PF = 2000, the signal was decomposed into IMFs and a Hilbert spectrogram was obtained. The gradient color bar shows the relative power value of the signal. K: decomposition number; PF: penalty factor; fitness: fitness value of VMD = the maximum envelope entropy of the IMF.
Sensors 24 05749 g005
Figure 6. Analysis of 8 s EEG segments during the (1) maintenance, (2) transition, and (3) emergence phases of general anesthesia induced by propofol. (A) Original EEG (8 s, 128 Hz, 1024 data points). (B) The power spectrum obtained through Fourier analysis. (C) Hilbert spectrum. To optimize K and PF, GWO analysis was conducted with 20 wolves, 20 repetitions, and settings of K = 2 to 6 (in increments of 1) and PF = 10 to 5000 (in increments of 10). (D) Decomposed intrinsic mode functions (IMFs). (E) The power spectrogram for a 64 s period was determined using the multitaper method. The first 8 s were the subject of the VMD analysis. The gradient color bar shows the relative power value of the signal. K: decomposition number; PF: penalty factor; fitness: fitness value of VMD = the maximum envelope entropy of the IMF.
Figure 6. Analysis of 8 s EEG segments during the (1) maintenance, (2) transition, and (3) emergence phases of general anesthesia induced by propofol. (A) Original EEG (8 s, 128 Hz, 1024 data points). (B) The power spectrum obtained through Fourier analysis. (C) Hilbert spectrum. To optimize K and PF, GWO analysis was conducted with 20 wolves, 20 repetitions, and settings of K = 2 to 6 (in increments of 1) and PF = 10 to 5000 (in increments of 10). (D) Decomposed intrinsic mode functions (IMFs). (E) The power spectrogram for a 64 s period was determined using the multitaper method. The first 8 s were the subject of the VMD analysis. The gradient color bar shows the relative power value of the signal. K: decomposition number; PF: penalty factor; fitness: fitness value of VMD = the maximum envelope entropy of the IMF.
Sensors 24 05749 g006
Table 1. Convergence rate in the optimization process for the hyperparameters K and PF of VMD using the GWO algorithm.
Table 1. Convergence rate in the optimization process for the hyperparameters K and PF of VMD using the GWO algorithm.
KPFFitness
Patient Patient Patient
#1#2#3Mean ± SD#1#2#3Mean ± SD#1#2#3Mean ± SD
Converge value2.352.72.82.6 ± 0.23178291823152804 ± 443−9.8769−9.8156−9.7289−9.8071 ± 0.07
Iteration%%%
1143117133131 ± 13787310886 ± 1999.924899.994099.760699.8361 ± 0.0829
2149113120127 ± 19828012094 ± 2399.832699.832799.760799.8562 ± 0.0662
3143111104119 ± 219682145108 ± 3399.832599.834699.892799.8946 ± 0.0610
4149111100120 ± 269682138106 ± 2999.834699.834699.937199.9141 ± 0.0610
514311196117 ± 2410090144111 ± 2999.972999.834699.937099.9200 ± 0.0708
612110096106 ± 149189133104 ± 2499.990299.988199.937099.9802 ± 0.0678
712110093105 ± 159189130104 ± 2399.990299.988199.967799.9820 ± 0.0155
810610095100 ± 69795125106 ± 1799.993999.990999.973399.9860 ± 0.0125
91001009398 ± 4100100118106 ± 10100.000100.00099.973399.9945 ± 0.0111
101001009398 ± 4100100112104 ± 7100.000100.00099.991099.9970 ± 0.0096
111001009899 ± 1100100103101 ± 2100.000100.00099.997099.9990 ± 0.0052
121001009899 ± 1100100103101 ± 2100.000100.00099.997699.9992 ± 0.0017
13100100100100 ± 0100100101100 ± 1100.000100.00099.998699.9995 ± 0.0014
14100100100100 ± 0100100101100 ± 1100.000100.00099.998899.9996 ± 0.0008
15100100100100 ± 0100100101100 ± 1100.000100.00099.998899.9996 ± 0.0007
16100100100100 ± 0100100101100 ± 0100.000100.00099.998999.9996 ± 0.0007
17100100100100 ± 0100100101100 ± 0100.000100.00099.999099.9997 ± 0.0006
18100100100100 ± 0100100100100 ± 0100.000100.000100.000100.0000 ± 0.000
19100100100100 ± 0100100100100 ± 0100.000100.000100.000100.0000 ± 0.000
20100100100100 ± 0100100100100 ± 0100.000100.000100.000100.0000 ± 0.000
Bold numbers indicate reaching the convergence value (=100%). K: decomposition number; PF: penalty factor; fitness: fitness value of VMD = the maximum envelope entropy of the IMF; SD: standard deviation.
Table 2. Grey wolf optimization in the VMD of EEGs obtained from three patients of general anesthesia induced by continuous intravenous propofol.
Table 2. Grey wolf optimization in the VMD of EEGs obtained from three patients of general anesthesia induced by continuous intravenous propofol.
KPFFitness
Before EmergenceAfter Emergence Before EmergenceAfter Emergence Before EmergenceAfter Emergence
PatientMean ± SDMean ± SDp-ValueMean ± SDMean ± SDp-ValueMean ± SDMean ± SDp-Value
#12.3 ± 0.72.5 ± 1.00.5252722 ± 17923190 ± 19140.309−5.7779 ± 0.0898−9.4826 ± 0.4004<0.001 *
#22.4 ± 0.83.0 ± 1.40.014 *2881 ± 16783071 ± 17720.675−9.7762 ± 0.0728−9.6799± 0.23410.016 *
#32.6 ± 1.02.8 ± 1.40.6242271 ± 18342271 ± 19220.107−9.7607 ± 0.0714−9.6886 ± 0.15630.008 *
Mean2.4 ± 0.82.7 ± 1.30.016 *2625 ± 17682844 ± 18690.989−8.4383 ± 0.0780−9.6170± 0.2636<0.001 *
SD0.2 ± 0.10.3 ± 0.2 316 ± 81500 ± 84 2.3039 ± 0.01030.1165 ± 0.1247
The optimal solutions for the K, PF, and fitness values were determined through fitness function optimization using the GWO algorithm (20 wolves and 20 optimization loops for all 73 epochs). The values between before-emergence and after-emergence were statistically compared. K: decomposition number; PF: penalty factor; fitness: fitness value of VMD = the maximum envelope entropy of the IMF; SD: standard deviation. * p < 0.05 between before emergence and after emergence.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kushimoto, K.; Obata, Y.; Yamada, T.; Kinoshita, M.; Akiyama, K.; Sawa, T. Variational Mode Decomposition Analysis of Electroencephalograms during General Anesthesia: Using the Grey Wolf Optimizer to Determine Hyperparameters. Sensors 2024, 24, 5749. https://doi.org/10.3390/s24175749

AMA Style

Kushimoto K, Obata Y, Yamada T, Kinoshita M, Akiyama K, Sawa T. Variational Mode Decomposition Analysis of Electroencephalograms during General Anesthesia: Using the Grey Wolf Optimizer to Determine Hyperparameters. Sensors. 2024; 24(17):5749. https://doi.org/10.3390/s24175749

Chicago/Turabian Style

Kushimoto, Kosuke, Yurie Obata, Tomomi Yamada, Mao Kinoshita, Koichi Akiyama, and Teiji Sawa. 2024. "Variational Mode Decomposition Analysis of Electroencephalograms during General Anesthesia: Using the Grey Wolf Optimizer to Determine Hyperparameters" Sensors 24, no. 17: 5749. https://doi.org/10.3390/s24175749

APA Style

Kushimoto, K., Obata, Y., Yamada, T., Kinoshita, M., Akiyama, K., & Sawa, T. (2024). Variational Mode Decomposition Analysis of Electroencephalograms during General Anesthesia: Using the Grey Wolf Optimizer to Determine Hyperparameters. Sensors, 24(17), 5749. https://doi.org/10.3390/s24175749

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop