Free
Clinical Science  |   September 2007
Median Frequency Revisited: An Approach to Improve a Classic Spectral Electroencephalographic Parameter for the Separation of Consciousness from Unconsciousness
Author Affiliations & Notes
  • Denis Jordan, M.Sc.
    *
  • Gudrun Stockmanns, Ph.D.
  • Eberhard F. Kochs, M.D.
  • Gerhard Schneider, M.D.
    §
  • * Research Fellow, ‡ Professor, Director and Chair, § Associate Professor, Department of Anesthesiology, Klinikum rechts der Isar, Technische Universität München. † Assistant Professor, Department of Computer Science and Applied Science, University Duisburg-Essen, Campus Duisburg, Germany.
Article Information
Clinical Science / Central and Peripheral Nervous Systems
Clinical Science   |   September 2007
Median Frequency Revisited: An Approach to Improve a Classic Spectral Electroencephalographic Parameter for the Separation of Consciousness from Unconsciousness
Anesthesiology 9 2007, Vol.107, 397-405. doi:10.1097/01.anes.0000278904.63884.4c
Anesthesiology 9 2007, Vol.107, 397-405. doi:10.1097/01.anes.0000278904.63884.4c
ELECTROENCEPHALOGRAM-based monitoring of the hypnotic component of anesthesia has gained popularity because it may help to complete the information given by vital signs during surgery. “Depth-of-anesthesia” monitors use algorithms to compute suitable values from the electroencephalographic signal (called parameters  ). The parameters are combined to calculate an index, which represents the hypnotic component of anesthesia. Currently, these monitors are of limited value in discerning consciousness from unconsciousness,1 and the question of whether the used algorithms can be improved should be considered. In the current investigation, signal analysis methods based on the spectrum of the electroencephalogram are used, and a generalization of the classic electroencephalographic parameters median frequency (MF)2 and spectral edge frequency (SEF95)3,4 is introduced. This new definition is adapted to the special characteristics of electroencephalographic signals, and the resulting parameters are tested for their ability to discern consciousness from unconsciousness, reflected by the ability to follow command. The parameters are assessed using data at the transition from consciousness to unconsciousness (or vice versa  ). Electroencephalographic data were from two different patient studies, one for parameter development, the other for parameter evaluation.
Materials and Methods
Protocol Design and Data Collection
The current analysis is based on data from two clinical studies, denoted as study A and study B, with similar study design approved by the ethics committee of the Technische Universität München, Faculty of Medicine, Munich, Germany. In each of the studies, 40 consenting adult patients undergoing general anesthesia were enrolled.1,5 In both studies, the study period was from induction of anesthesia until patients followed command after surgery and included a reduction of the hypnotic agent after tracheal intubation until patients followed command. Patients with contraindications to the study drugs, a history of psychiatric or neurologic disease, drug abuse or medication known to affect the central nervous system, pregnancy, or indication for rapid-sequence induction were excluded from the study. In study A, patients were randomly assigned to an anesthetic regimen with remifentanil (minimum infusion rate 0.2 μg · kg−1· min−1) and sevoflurane (n = 20) or remifentanil and propofol (n = 20). In study B, the sevoflurane and propofol groups were divided into two subgroups with either “low infusion rate” of remifentanil (0.1 μg · kg−1· min−1) or “high infusion rate” (0.2 μg · kg−1· min−1). Without premedication, remifentanil infusion was started via  a cannula in the cubital vein. In 30-s intervals, patients were asked to squeeze the investigator's hand. A response was verified by an immediate repetition of the command that also required a response. This prevents a misinterpretation of involuntary movement as a response. Anesthesia was slowly induced with sevoflurane inhalation or propofol injection (0.7 mg/kg, followed by 20 mg every 30 s). The first time when the patient did not squeeze the investigator's hand to command was labeled as loss of consciousness (LOC1). Additional propofol or sevoflurane was given to increase depth of anesthesia, and a tourniquet was used to occlude the circulation of the right forearm for 5 min to retain the ability to move the hand to command, before succinylcholine (1.0 mg/kg) was given (isolated forearm technique of Tunstall).6 After intubation, sevoflurane or propofol was stopped until patients followed command (return of consciousness [ROC1]). Thereafter, sevoflurane (5 vol%) or propofol (20-mg boluses) was readministered to induce anesthesia again slowly. The time when patients stopped responding to command again was defined as loss of consciousness 2 (LOC2), and requests to squeeze the hand were stopped. Anesthetic drugs were administered according to clinical practice, and surgery was performed. At the end of surgery, requests to squeeze the hand were recommenced, and sevoflurane, propofol, and remifentanil were discontinued. Return of consciousness 2 (ROC2) was observed at the first verified response to command. Recovered from anesthesia, patients were asked for signs of recall in the recovery room. This interview7 was repeated within 48 h in the ward. Standard monitoring parameters were measured with a Datex AS/3 (Datex-Ohmeda Division Instrumentation Corp., Helsinki, Finland) compact anesthesia monitor. For data transfer and storage, a personal computer with NeuMonD (Department of Anesthesiology, Klinikum rechts der Isar, Technische Universität München) was used. NeuMonD is a software program developed by members of the research group allowing the recording of monitoring data and the electronic storage of events and comments during the study.8 In addition to standard monitoring, electroencephalogram and auditory evoked potentials were measured (study A). Two channel electroencephalographic signals at electrode positions AT1, M2, Fpz (reference), and F7 (ground) were recorded on a second personal computer with synchronized system time. This auditory evoked potential/electroencephalographic device has been designed specifically for intraoperative use as described previously.9 Electroencephalogram and concomitant trigger information of auditory evoked potentials were stored with a sampling rate of f  s= 1,000 Hz and 12-bit amplitude resolution. In study B, electroencephalogram was recorded using the Aspect A-1000 electroencephalographic monitor (BIS® version 3.3; Aspect Medical Systems Inc., Newton, MA). A two-channel referential electroencephalogram was obtained with ZipPrep Ag/AgCl electrodes in positions AT1, AT2, Fz (reference), and Fp1 (ground, electrode positions according to the international 10–20 system). The high pass was set at 0.25 Hz, no low pass was used, and the notch Filter (50 Hz) was enabled. Electroencephalogram was continuously digitized at 256 Hz per channel and simultaneously recorded with standard monitoring parameters. For data analysis in both studies, a time window duration T  = 8 s was used.
Signal processing and statistical analysis were performed using LabVIEW 6.0 (National Instruments, Austin, TX; 2000), MATLAB 6.0 Release 12 (The MathWorks, Inc., Natick, MA; 2000), Microsoft Access 2002 (Microsoft Corporation, Redmond, WA), and R 2.4.0 (R Foundation for Statistical Computing, Vienna, Austria; 2006) on personal computers with Windows NT/XP (Microsoft Corporation; 2002).
Median Frequency and Spectral Edge Frequency
The electroencephalogram is obtained by recording the electrical activity of the brain. Scalp electrodes are used to detect electrical activity on the surface of the head. This reflects the activity of action potentials of cortical neurons, in particular superposition of postsynaptic potentials mainly generated in layer V of the cerebral cortex.10 Different rhythms are observed, and in the past, the following frequency bands related to brain activity were introduced for classification: the δ band, including a frequency range up to 4 Hz; the θ band [4 Hz, 8 Hz], the α band [8 Hz, 12 Hz], and the β band [12 Hz, 30 Hz]. Spectral analysis—as calculated by (discrete) Fourier transform—seems to be an appropriate method to analyze electroencephalographic time series, because it follows this approach of signal analysis. Although the spectrum including amplitude and phase contains not more information as the time series, some signal characteristics are more evident by observing the spectrum of the signal. To reduce dimension of the spectrum, suitable parameters can be extracted and assessed, such that they correspond to the observed patient states.
Parameters derived from the electroencephalographic power spectrum have been shown to indicate—at least to some extent—anesthetic effects. Two of those parameters are the median frequency (MF or SEF50)2 and spectral edge frequency (SEF95)3,4 between low and high cutoff frequency f  lowand f  high. MF provides an overall average frequency and roughly decreases with increasing “depth of anesthesia.” This is consistent with the expectation that in “deeper” anesthesia, high-frequency components of the superficial neuronal areas are replaced by low-frequency components.10 Unfortunately, some effects limit the usefulness of this general observation. These may in part be related to the presence of muscle activity, i.e.  , electromyogram superposed to the electroencephalographic signal. Electromyogram may change mainly the spectrum of the electroencephalographic time signal such that most of the electroencephalographic parameters are influenced by the presence of eye blinks and electromyogram-related high-frequency components. Changes of the electroencephalographic time signal may not primarily reflect the level of general anesthesia, but a specific drug effect (which subsequently influences the level of anesthesia). As an example, this occurs on the lower end of the frequency spectrum, where opioids induce increasing amplitudes of the electroencephalographic δ band (without primarily effecting the level of hypnosis). On the other hand, relatively deep propofol anesthesia is also reflected by increasing δ activity, but this time the level of hypnosis is increased. As this example illustrates, changes of an electroencephalographic parameter may be drug specific and can not be used equally for all types of general anesthesia. On the higher end of the frequency spectrum, the behavior of most parameters is affected by a characteristic biphasic reaction, which occurs during induction of anesthesia. This was shown for electroencephalographic amplitude and SEF95 during the transition from consciousness to unconsciousness and is reflected by an initial increase of the effect variable followed by a decrease at higher concentrations.11 
WSMF, a Generalization of MF
In general, a parameter is characterized by numerous variables, e.g.  , its low and high cutoff frequency f  lowand f  high, and each of them augments the parameter dimension by one. The setting of these variables may considerably affect the behavior of the parameter itself, e.g.  , its capability to distinguish consciousness from unconsciousness. Therefore, the influence of variable settings must be examined. This was mostly neglected in the past.
The parameter weighted spectral median frequency (WSMF) can be seen as a generalization of MF: Given an electroencephalographic signal s  , the amplitude spectrum A  (f  ) of s  is computed by a discrete Fourier transform. The amplitudes are then weighted by the p  thpower, where the influence of the value p  is explained below. The “area” of the weighted spectrum A  p  (f  ) delimited by the frequency axis between the cutoff frequencies f  lowand f  highis divided into two parts at the frequency called WSMF. The frequency WSMF is such that the complete “area” between f  lowand f  highmultiplied by r  is equal to the “area” between f  lowand WSMF.
A value p  = 1 defines the computation of WSMF based on the amplitude spectrum. p  = 2 computes WSMF using the power spectrum as for MF and SEF95. Particularly p  < 1 allows a non linear attenuation of high amplitudes and leads to an equalization of amplitude differences. It may provide a more stable behavior of WSMF when artifacts containing high amplitudes occur. The variable r  can be denoted as splitting ratio and defines how much the spectrum up to WSMF is weighted in relation to the total spectrum between f  lowand f  high. r  = 0.5 as defined for MF provides a symmetrical weight of the spectrum below and above the frequency WSMF, where r  = 0.95 as for SEF95 means that the spectrum above WSMF obtain a main emphasis. A mathematical description of WSMF can be found in 1.
A generalization WSMF of MF with four variables f  low, f  high, p  , r  for parameter configuration is obtained. Their values with suitable bounds are denoted by the four-dimensional set I  (see equation (4)in 1for an exact definition). Therefore, WSMF defines an own parameter for every setting f  low, f  high, p  , r  . To simplify notation, a particular parameter WSMF with settings f  low, f  high, p  , r  is denoted by WSMF f  lowf  high, if there is no risk of confusion by different settings of p  and r  . For a general configuration, the notation WSMF is used without indication of its variables.
Settings of f  low, f  high, p  , r  for the Development of WSMF
To obtain best possible classification results for separating consciousness and unconsciousness of the parameters WSMF with settings of f  low, f  high, p  , r  , a discrete and finite subset of I  is determined. This allows the computation of WSMF for a finite number of settings in I  .
The technical equipment for data acquisition described above determines the lower bound of 0.5 Hz for f  low. The upper bound for f  highis set to approximately 50 Hz. Therefore, the signal does not exclusively reflect electroencephalogram, but may also include electromyogram in further analysis. The exponent p  indicates how much WSMF is affected by high amplitudes in the spectrum A  . For calculation of the classic MF, p  = 2, i.e.  , the power spectrum is chosen, such that high frequencies with small amplitudes, which may be typical for consciousness, are not sufficiently weighted when electromyographic components are absent. Therefore, it seems very promising to consider also smaller values for the variable p  , and its range is defined between 0.1 and 2.4. The splitting ratio r  is generally defined in the open interval with 0 < r  < 1, and therefore settings between 0.05 and 0.95 are considered, such that MF and SEF95 are included in the analysis.3,4 It allows a crossover between MF and SEF95. The low-frequency resolution of the discrete Fourier transform and the time stability of WSMF are mainly affected by the setting of T  and increase monotonically with increasing values of T  . As stated above, in the current analysis, T  = 8 s.
Summarizing, the following settings for the variables f  low, f  high, p  , r  are used (subset of I  ; see equation (4)in 1):
The resulting configurations of WSMF are included for parameter development.
Statistical Analysis
Prediction probability (PK) has been established in anesthesia12 as a statistical method to asses the capability of a parameter to discern different levels of anesthesia, e.g.  , consciousness and unconsciousness. The aim of the PKanalysis is to quantify the association between the (clinically) observed anesthetic level and the parameter values, in this case of WSMF. A short description can be found in 2.
To give an estimation of suitable configurations of WSMF, the extremes of the PKvalues with respect to the settings of f  low, f  high, p  , r  in I  are computed using study A. PKanalysis is then applied on study B to provide a further assessment of the capability of WSMF indicating the changes from consciousness to unconsciousness (or vice versa  ).
The PKanalysis used below is defined on the basis of a set of electroencephalographic signals (samples) of T  = 8 s length directly preceding or after the transition between consciousness and unconsciousness. Study A and study B involved 40 patients each. Because two signals were analyzed at every clinical event LOC1, LOC2, ROC1, and ROC2, a maximum of 320 signals was available from study A and study B (see Protocol Design and Data Collection section). The signals are situated either completely in a phase of consciousness or in a phase of unconsciousness. For study A (development study), 306 signals (152 assigned to consciousness and 154 assigned to unconsciousness) and for study B (evaluation study), 286 signals (135 assigned to consciousness and 151 assigned to unconsciousness) were used for PKanalysis. Fourteen (study A) and 34 (study B) signals containing artifacts expressed by signals of constant amplitude (flat line) or values exceeding the measuring range of 250 μV were manually (study A) or automatically (study B) excluded from analysis.
Ninety-five percent confidence intervals based on the number of samples including Bonferroni correction were used to compare PKvalues of selected parameters of study B, where PKvalues are significantly different if the corresponding confidence intervals are disjoint.
Results
Particular Cases of WSMF
To provide an intuitive idea of the behavior of WSMF, both the stability of the parameter and its ability to distinguish consciousness from unconsciousness are examined using an exemplarily electroencephalographic signal of study A for some particular settings of f  low, f  high, p  , r  .
Figure 1shows time series of two categories of WSMF, one computed with the full “classic” electroencephalographic frequency range from 0.5 to 30 Hz, represented by MF (f  low= 0.5 Hz, f  high= 30 Hz, p  = 2, r  = 0.5)2 and SEF95 (f  low= 0.5 Hz, f  high= 30 Hz, p  = 2, r  = 0.95),3,4 and a second with an higher low cutoff frequency, represented by WSMF8–49(f  low= 8 Hz, f  high= 49 Hz, p  = 1, r  = 0.5) and WSMF8–30(f  low= 8 Hz, f  high= 30 Hz, p  = 0.4, r  = 0.5). WSMF8–49and WSMF8–30allow a much better appreciation of the hypnotic component of anesthesia than MF and SEF95; in particular, they are more stable and show improved monotonic behavior, i.e.  , their values are decreasing during increasing levels of hypnosis and vice versa  . Note that all four parameters show a certain amount of inappropriately low values during consciousness. Furthermore, all parameter time series are affected by the characteristic biphasic reaction during induction of anesthesia,11 expressed by an increase of the parameter values before loss of consciousness (LOC1 and LOC2) after propofol administration.
Fig. 1. Induction of anesthesia in a randomly selected patient (propofol): Parameter time series based on study A of weighted spectral median frequency (WSMF8–49  ,A  ; WSMF8–30  ,B  ), median frequency (MF;  C  ), and spectral edge frequency (SEF95;  D  ). Parameter values are normalized to the interval [0, 1]. The  gray lines  indicate changes of the level of consciousness: loss of consciousness at induction (LOC1: 214 s), return of consciousness after intubation (ROC1: 784 s), followed by second loss of consciousness (LOC2: 846 s). WSMF8–49and WSMF8–30indicate better the state of consciousness than MF and SEF95. Parameter calculation was performed in time steps of 4 s using a signal of 8 s length (overlapping factor 2). 
Image Not Available
Fig. 1. Induction of anesthesia in a randomly selected patient (propofol): Parameter time series based on study A of weighted spectral median frequency (WSMF8–49  ,A  ; WSMF8–30  ,B  ), median frequency (MF;  C  ), and spectral edge frequency (SEF95;  D  ). Parameter values are normalized to the interval [0, 1]. The  gray lines  indicate changes of the level of consciousness: loss of consciousness at induction (LOC1: 214 s), return of consciousness after intubation (ROC1: 784 s), followed by second loss of consciousness (LOC2: 846 s). WSMF8–49and WSMF8–30indicate better the state of consciousness than MF and SEF95. Parameter calculation was performed in time steps of 4 s using a signal of 8 s length (overlapping factor 2). 
×
The aim is to find a parameter WSMF with appropriate settings of f  low, f  high, p  , r  that reflects the hypnotic component of anesthesia, such that it is as robust as possible against electromyographic artifacts and that it operates in the “classic” electroencephalographic spectrum up to 30 Hz, because detection of consciousness should not be based on electromyographic activity.
Influence of the Exponent p  and the Splitting Ratio r  on WSMF
Two observations were made when WSMF was applied on development study A using the discrete settings of the variables as in equation (1).
First, PK≥ 0.82 for WSMF can be achieved at every level of the splitting ratio r  including frequencies up to 30 Hz only, i.e.  , with f  high≤ 30 Hz. PK≥ 0.82 is considered to be of “similar grade” as the maximal PKvalue of 0.824.
Second, for each chosen frequency range from f  lowto f  high, i.e.  , [f  low, f  high], it can be observed that at constant levels of PK, the splitting ratio r  decreases with decreasing exponent p  , as shown in figure 2. This implies that a high value of p  combined with a high value of r  leads to a better PKvalue than a high value of p  combined with a low value of r  . For example, there is PK= 0.73 in case of SEF95 (p  = 2, r  = 0.95) and PK= 0.55 in case of MF (p  = 2, r  = 0.5). Both observations suggest that the variables p  and r  are in a way related, and one of them can be chosen to some extent arbitrarily.
Fig. 2. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A at fixed cutoff frequencies  f  low= 5 Hz  ,f  high= 50 Hz (  A  ) and  f  low= 9 Hz  , f  high  = 30 Hz (  B  ). For the exponent  p  (  x-axis  ), settings between 0.2 and 2.0 are displayed; for the splitting ratio  r  (  y-axis  ), settings between 0.1 and 0.95 are displayed. The resulting PKvalues of WSMF with corresponding settings of  p  and  r  are plotted as third dimension (  z-axis  ,left  ) and represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . A relation between  p  and  r  can be observed if values of PK≥ 0.80 are assumed (  red area  ). The PKvalues are interpolated giving a continuous approximation of the discrete settings to obtain a better presentation (  left  and  right  ).
Image Not Available
Fig. 2. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A at fixed cutoff frequencies  f  low= 5 Hz  ,f  high= 50 Hz (  A  ) and  f  low= 9 Hz  , f  high  = 30 Hz (  B  ). For the exponent  p  (  x-axis  ), settings between 0.2 and 2.0 are displayed; for the splitting ratio  r  (  y-axis  ), settings between 0.1 and 0.95 are displayed. The resulting PKvalues of WSMF with corresponding settings of  p  and  r  are plotted as third dimension (  z-axis  ,left  ) and represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . A relation between  p  and  r  can be observed if values of PK≥ 0.80 are assumed (  red area  ). The PKvalues are interpolated giving a continuous approximation of the discrete settings to obtain a better presentation (  left  and  right  ).
×
Further ranges [f  low, f  high] show that the splitting ratio r  < 0.6 is favorable if f  high≤ 30 Hz is assumed. r  = 0.5 is selected, giving a direct relation of WSMF to the classic parameter MF. This leads to a reduced number of settings with
containing three variables f  low, f  high, p  only and 5,760 configurations of WSMF. The totality of settings of the variables is denoted by I  discrete(three-dimensional finite set).
Optimization of the Configuration of WSMF Using Study A
Based on data from development study A, figures 3A and Bshow the behavior of WSMF at the transition between consciousness and unconsciousness with settings of f  low, f  high, p  in the three-dimensional set I  discreteas defined in (2) and visualized in drawing a four-dimensional cube. The fourth dimension indicates the results of PKanalysis and is superposed as colors at every point of I  discrete.
Fig. 3. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A. Settings of the exponent  p  (  x-axis  ) are displayed between 0.4 and 2.4; settings of the cutoff frequencies  f  high(  y-axis  ) and  f  low(  z  -axis) are displayed between 25 Hz and 53 Hz and between 0.5 Hz and 15 Hz, respectively. In  A  and  B  , the PKvalues at the corresponding setting of  f  low  ,f  high  ,p  are represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . A maximal PKvalue of 0.824 results for WSMF with  f  low= 8 Hz  ,f  high= 30 Hz, and  p  = 0.4 that is situated near the frontal edge of the cube (  dark red area  ). The cube in  B  is interpolated, giving a continuous approximation of the discrete settings to obtain a better presentation. The cube in  C  shows only values of PK≥ 0.82, which are represented by  black dots.  They occur as two “clusters,” one as “high-frequency cluster” with  f  high≥ 40 Hz and one as “low-frequency cluster” with  f  high≤ 30 Hz. 
Image Not Available
Fig. 3. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A. Settings of the exponent  p  (  x-axis  ) are displayed between 0.4 and 2.4; settings of the cutoff frequencies  f  high(  y-axis  ) and  f  low(  z  -axis) are displayed between 25 Hz and 53 Hz and between 0.5 Hz and 15 Hz, respectively. In  A  and  B  , the PKvalues at the corresponding setting of  f  low  ,f  high  ,p  are represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . A maximal PKvalue of 0.824 results for WSMF with  f  low= 8 Hz  ,f  high= 30 Hz, and  p  = 0.4 that is situated near the frontal edge of the cube (  dark red area  ). The cube in  B  is interpolated, giving a continuous approximation of the discrete settings to obtain a better presentation. The cube in  C  shows only values of PK≥ 0.82, which are represented by  black dots.  They occur as two “clusters,” one as “high-frequency cluster” with  f  high≥ 40 Hz and one as “low-frequency cluster” with  f  high≤ 30 Hz. 
×
Maximum PKvalues are obtained for settings of f  low, f  high, p  within the cube, i.e.  , inside I  discrete, where two “clusters” can be observed (fig. 3C). The first “high-frequency cluster” provides parameters WSMF operating at frequencies with f  highabove 40 Hz and exponents p  ∈[1.0, 1.5]. The frequency range is bounded at the lower limit around f  low≅ 8 Hz. The second “low-frequency cluster” is much more interesting, because it reflects analysis of the “classic” electroencephalographic range, i.e.  , frequencies up to 30 Hz, whereas in high-frequency ranges, electroencephalographic activity overlaps with electromyographic activity. The frequency f  lowremains the same as in the “high-frequency cluster,”i.e.  , f  low≅ 8 Hz. To illustrate the behavior of WSMF, some slices through I  discreteare selected at fixed f  low(fig. 4). A comparison of the three slices confirms favorable results for WSMF defining f  lowin the range of 8 Hz. A change of f  lowin the range between 7 and 10 Hz does not considerably change the characteristics of WSMF, and PK= 0.82 while f  high≤ 30 Hz is maintained. This behavior can be explained in part by the reduced influence of opioid-induced δ-band activation and eye blinks during induction of anesthesia, which are mainly present in the lower-frequency range.
Fig. 4. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A at three fixed levels of the low cutoff frequency  f  lowwith variable settings of the exponent  p  (  x-axis  ) between 0.4 and 2.4 and high cutoff frequency  f  high(  y-axis  ) between 25 Hz and 53 Hz. In  A  , f  low= 0.5 Hz; in  B, f  low= 8 Hz; and in  C, f  low= 13 Hz. The resulting PKvalues of WSMF with the corresponding settings of  p  and  f  highare represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . PKvalues are interpolated.  A  contains the configuration of median frequency (MF) with  f  low= 0.5 Hz,  f  high= 30 Hz,  p  = 2,  r  = 0.5 leading to a PKof 0.55 (all PKvalues smaller than 0.66 are represented by  dark blue  ). B  contains the configurations of WSMF that provide maximal PKvalues in the range of 0.82 within both “clusters” shown in  figure 3C(  dark red areas  ):  f  low= 8 Hz,  f  high= 49 Hz,  p  = 1,  r  = 0.5 (denoted by WSMF8–49) and  f  low= 8 Hz,  f  high= 30 Hz  , p  = 0.4  , r  = 0.5 (denoted by WSMF8–30). 
Image Not Available
Fig. 4. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A at three fixed levels of the low cutoff frequency  f  lowwith variable settings of the exponent  p  (  x-axis  ) between 0.4 and 2.4 and high cutoff frequency  f  high(  y-axis  ) between 25 Hz and 53 Hz. In  A  , f  low= 0.5 Hz; in  B, f  low= 8 Hz; and in  C, f  low= 13 Hz. The resulting PKvalues of WSMF with the corresponding settings of  p  and  f  highare represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . PKvalues are interpolated.  A  contains the configuration of median frequency (MF) with  f  low= 0.5 Hz,  f  high= 30 Hz,  p  = 2,  r  = 0.5 leading to a PKof 0.55 (all PKvalues smaller than 0.66 are represented by  dark blue  ). B  contains the configurations of WSMF that provide maximal PKvalues in the range of 0.82 within both “clusters” shown in  figure 3C(  dark red areas  ):  f  low= 8 Hz,  f  high= 49 Hz,  p  = 1,  r  = 0.5 (denoted by WSMF8–49) and  f  low= 8 Hz,  f  high= 30 Hz  , p  = 0.4  , r  = 0.5 (denoted by WSMF8–30). 
×
In contrast, the classic MF with a frequency range defined by f  low≤ 1 Hz provides PK≤ 0.60 for any chosen frequency f  highand seems to be less adequate in separating consciousness from unconsciousness (fig. 4A). Finally, choosing f  low> 10 Hz, essential information in the α or even β band is ignored and the performance of WSMF decreases, as shown in figure 4C. However, this decrease is not as dramatic as with f  low< 5 Hz. PKanalysis shows that the correlation to the anesthetic state is improved when p  decreases.
As a result of these observations, two configurations of WSMF within both clusters of I  discreteand leading to maximal PKvalues are selected: WSMF8–49(f  low= 8 Hz, f  high= 49 Hz, p  = 1, r  = 0.5) and WSMF8–30(f  low= 8 Hz, f  high= 30 Hz, p  = 0.4, r  = 0.5). MF and SEF95 are outside both clusters, and their computation is based on the power spectrum (p  = 2) using the full “classic” electroencephalographic frequency range from 0.5 to 30 Hz.
Evaluation of the Configuration of WSMF Using Study B
Study B was used to evaluate the results obtained by study A. The main focus of this analysis was to evaluate whether the selected parameters WSMF8–49and WSMF8–30within both “clusters” in figure 3Cmaintain their ability to separate consciousness from unconsciousness when applied to data of study B.
Table 1shows the PKstatistics of both studies at the transition from consciousness to unconsciousness. Although the corresponding PKvalues are different for both studies, WSMF8–49and WSMF8–30obtain significantly better results (PK= 0.82 and PK= 0.79) in comparison with MF and SEF95 (PK= 0.58 and PK= 0.59), where confidence intervals considering the number of samples were performed for the involved parameters.
Table 1. Prediction Probability Values for Different Settings 
Image Not Available
Table 1. Prediction Probability Values for Different Settings 
×
If frequencies higher than 30 Hz are included in the analysis, the parameter WSMF8–49in the “high-frequency cluster” of I  discretein figure 3Cprovide PK= 0.82 using study B, which is similar to WSMF8–30with PK= 0.79.
Discussion
Spectral analysis based on the linear Fourier transform is still a powerful mathematical tool, which can be used effectively to design electroencephalographic parameters indicating the hypnotic component of anesthesia. Apparently, despite nonlinear characteristics of the electroencephalogram,4,13–15 the signal properties described by harmonic signal decomposition allow the development of electroencephalographic parameters which produce a good resolution during the transition from consciousness to unconsciousness (or vice versa  ).
So far, different settings and their interaction have mostly been disregarded in attempts to improve the quality of such a parameter. The parameter should be adjusted carefully. For this purpose, it must be examined not only at a few discrete settings, but as far as possible with respect to a wide range of configurations.
Generalization of MF and SEF95 gives a class of new parameters denoted by WSMF. To provide appropriate configurations with respect to the requirements above, it is crucial to find suitable frequency ranges in which the parameter leads to reliable results. The choice of a low cutoff frequency f  low∈[7 Hz, 10 Hz] and a high cutoff frequency f  high∈[25 Hz, 30 Hz] in combination with a real-valued exponent of the amplitude p  ≤ 1 leads to improved results, in particular when compared with MF and SEF95 based on the power spectrum (p  = 2) of the full “classic” electroencephalographic frequency range from 0.5 to 30 Hz. This can be explained by the low-frequency spectrum signals generated by muscle activity of eye blinks. MF and SEF95 are in particular affected, because a choice of f  low= 0.5 Hz may even compromise a differentiation between eye blinks and δ activity (fig. 1). With respect to separation of consciousness to unconsciousness, a classification rate of PK≥ 0.79 for WSMF8–30(f  low= 8 Hz, f  high= 30 Hz, p  = 0.4, r  = 0.5) is reached both in study A (development study) and in study B (evaluation study). As a consequence of the limited frequency range, the influence of electromyogram or eye blinks is avoided, and the influence of drug specific effects on the δ band is reduced. While a low cutoff frequency of around 8 Hz leads to an improved behavior of WSMF at the transition between consciousness and unconsciousness, it may also limit the capability of WSMF to indicate “deeper” anesthesia that may be characterized by δ activity. The statistical comparison of WSMF8–30and WSMF8–49(f  low= 8 Hz, f  high= 49 Hz, p  = 1, r  = 0.5) shows a significant improvement in comparison with MF and SEF95.
On the data basis of study B, WSMF8–49provides similar PK= 0.82 in comparison with WSMF8–30(PK= 0.79). In general, inclusion of high-frequency components (> 30 Hz) may be useful in the detection of awareness. But, in particular if measured on the forehead, electroencephalographic signals may overlap with electromyographic signals, i.e.  , activity of forehead muscles, which are characterized by high frequency activity. As a consequence, increase of high-frequency activity during awareness may reflect increasing muscle activity on the patient's forehead. Therefore, inclusion of frequencies above 30 Hz bears the risk that a parameter does not primarily reflect activity of the main target organ of anesthesia, the brain, but may also rely on increasing muscle activity. Because the frequency range of electromyogram overlaps with the high-frequency electroencephalographic range, the presence of (high-frequency) electromyographic activity may lead to an overestimation of high-frequency electroencephalographic activity, and an increase of spectral energy at high frequencies, which will increase calculated parameter values. This increase, however, is based on electromyographic rather than electroencephalographic activation. Therefore, a parameter that includes high-frequency components may also be a surrogate measure (muscle activity) of the hypnotic component of anesthesia. The use of this surrogate parameter bears the risk that a patient who is fully awake but paralyzed by neuromuscular blockers may not be detected as “awake,” because no electromyogram is detected.
Using recorded signals of both studies, the analysis of BIS®, a commercially available monitor based on electroencephalogram, produces PK= 0.74 (study A) and PK= 0.68 (study B),1,5 which is significantly lower than the results of WSMF8–30or even WSMF8–49. These results are consistent with previous studies, where BIS values show a high correlation with propofol target concentrations, whereas the discrimination between consciousness and unconsciousness is less than ideal.16 
The current investigation shows that a more general approach of a parameter is essential to understand the interaction between the configuration and the behavior of the parameter applied on electroencephalographic signals. Suitable settings lead to a significant improvement in representing the hypnotic component of anesthesia. In a further step, analysis must be performed not only for the two states consciousness and unconsciousness, but for the entire range from light sedation to general anesthesia. It should be evaluated whether the good performance of WSMF can be further improved by combining the parameter with other electroencephalographic parameters.
The authors thank the members of the research group EEG and Neurophysiology of the Department of Anaesthesiology, Technische Universität München, Munich, Germany, and the research group Knowledge-based Signal Processing for Medical Applications of the Department of Computer Science and Applied Science, University Duisburg-Essen, Campus Duisburg, Germany, for their help.
Definition of WSMF
In the continuous case, i.e.  , assuming a time-continuous electroencephalographic signal s  , the complex spectrum is computed by the Fourier transform
and A  (f  ) =|c  (f  )| is the amplitude spectrum of the electroencephalographic signal s  (s  contains in the presented case T  ·f  s= 8,000 samples, T  = 8 s, f  s= 1 kHz, where s  (t  ) = 0 if t  ∉[0, T  ]). Then the weighted spectral median frequency (WSMF) is defined as
with high and low cutoff frequencies f  low, f  high, and p  as the exponent of the amplitude A  . The variable r  can be denoted as splitting ratio and defines how much the spectrum up to WSMF is weighted in relation to the total spectrum between f  lowand f  high. A value p  = 1 defines the computation of WSMF based on the amplitude spectrum, and p  = 2 defines a computation using the power spectrum. A different notation of equation (2)is obtained if the total band is split according [f  low, f  high]=[f  low, WSMF][WSMF, f  high], such that
with q  =r  /(1 −r  ). In case of r  = 0.5, i.e.  , q  = 1 as given for MF, WSMF is the frequency that splits [f  low, f  high] in two sub-bands of equal weights.
WSMF is a generalization of MF using 4 variables f  low, f  high, p  , r  for parameter configuration. Their values are in the set denoted by
for suitable bounds of I  (ă  denotes the lower and â  the upper bound of a variable a  ).
Because of the discrete recording of electroencephalogram in the current investigation, all computations were performed using a discrete Fourier transform. To reduce the leakage effect due to the boundary of the analyzed electroencephalographic signal s  , a Hamming window was used for the calculation of WSMF.
PKStatistics
PKstatistics are closely related to the definition of Kim's dy*x  12 : Let x  = (x  1, x  2, …, x  n  ) be the vector of obtained parameter values computed from n  signal intervals s  1, …, s  n  and y  = (y  1, y  2, …, y  n  ) the corresponding observed anesthetic state, in our case y  i  ∈{0 (unconscious), 1 (conscious)}. Then dy*x  differs three cases, (A) the probability pc  , that a pair (x  i  , x  j  ) (i  , j  ∈{1, …, n  }, i  ≠j  ) of arbitrary drawn parameter values is in concordance with the anesthetic level (y  i  , y  j  ), i.e.  , y  i  > y  j  ⇒x  i  > x  j  , (B) the probability pd  that it is in discordance, i.e.  , y  i  > y  j  ⇒x  i  < x  j  and (C) the probability pt  that a parameter value is assigned to more than one anesthetic level, i.e.  , y  i  ≠y  j  ⇒x  i  =x  j  . The definition
implies dy*x  ↑ with pc  ↑, pd  ↓, pt  ↓, and dy*x  ∈[−1, 1]. By the affine transformation PK= (dy*x  + 1)/2, the index dy*x  is mirrored at 0.5, such that PK= 1 gives a completely concordant relation between x  and y  , PK= 0 implies completely discordant relation, and PK= 0.5 means that there is no relation between x  and y  .
References
Schneider G, Hollweck R, Ningler M, Stockmanns G, Kochs E: Detection of consciousness by electroencephalogram and auditory evoked potentials. Anesthesiology 2005; 103:934–43Schneider, G Hollweck, R Ningler, M Stockmanns, G Kochs, E
Schwilden H, Stoeckel H: Investigations on several EEG-parameters as indicators of the state of anaesthesia the median: A quantitative measure of the depth of anaesthesia. Anasth Intensivther Notfallmed 1980; 15:279–86Schwilden, H Stoeckel, H
Rampil IJ: A primer for EEG signal processing in anesthesia. Anesthesiology 1998; 89:980–1002Rampil, IJ
Thakor NV, Tong S: Advances in quantitative electroencephalogram analysis methods. Annu Rev Biomed Eng 2004; 6:453–95Thakor, NV Tong, S
Schneider G, Gelb AW, Schmeller B, Tschakert R, Kochs E: Detection of awareness in surgical patients with EEG-based indices: Bispectral index and patient state index. Br J Anaesth 2003; 91:329–35Schneider, G Gelb, AW Schmeller, B Tschakert, R Kochs, E
Tunstall ME: Detecting wakefulness during general anaesthesia for caesarean section. BMJ 1977; 1:1321Tunstall, ME
Brice DD, Hetherington RR, Utting JE: A simple study of awareness and dreaming during anaesthesia. Br J Anaesth 1970; 42:535–42Brice, DD Hetherington, RR Utting, JE
Stockmanns G, Ningler M, Omerovich A, Kochs E, Schneider G: NeuMonD: A tool for the development of new indicators of anaesthetic effect. Biomed Tech 2007; 52:96–101Stockmanns, G Ningler, M Omerovich, A Kochs, E Schneider, G
Jordan C, Weller C, Thornton C, Newton DEF: Monitoring evoked potentials during surgery to assess the level of anaesthesia. J Med Eng Technol 1995; 19:77–9Jordan, C Weller, C Thornton, C Newton, DEF
Widman G, Schreiber T, Rehberg B, Hoeft A, Elger CE: Quantification of depth of anesthesia by nonlinear time series analysis of brain electrical activity. Phys Rev E 2000; 62:4898–903Widman, G Schreiber, T Rehberg, B Hoeft, A Elger, CE
Kuizenga K, Wierda JM, Kalkman CJ: Biphasic EEG changes in relation to loss of consciousness during induction with thiopental, propofol, etomidate, midazolam or sevoflurane. Br J Anaesth 2001; 86:354–60Kuizenga, K Wierda, JM Kalkman, CJ
Smith WD, Dutton RC, Smith N: Measuring the performance of anesthetic depth indicators. Anesthesiology 1996; 84:38–51Smith, WD Dutton, RC Smith, N
Bruhn J, Röpcke H, Hoeft A: Approximate entropy as an electroencephalographic measure of anesthetic drug effect during desflurane anesthesia. Anesthesiology 2000; 92:715–26Bruhn, J Röpcke, H Hoeft, A
Zhang XS, Roy RJ, Jensen EW: EEG complexity as a measure of depth of anesthesia for patients. IEEE Trans Biomed Eng 2001; 48:1424–33Zhang, XS Roy, RJ Jensen, EW
Pincus S, Singer BH: Randomness and degree of irregularity. Proc Natl Acad Sci U S A 1996; 93:2083–8Pincus, S Singer, BH
Gajraj RJ, Doi M, Mantzaridis H, Kenny GN: Comparison of bispectral EEG analysis and auditory evoked potentials for monitoring depth of anaesthesia during propofol anaesthesia. Br J Anaesth 1999; 82:672–8Gajraj, RJ Doi, M Mantzaridis, H Kenny, GN
Fig. 1. Induction of anesthesia in a randomly selected patient (propofol): Parameter time series based on study A of weighted spectral median frequency (WSMF8–49  ,A  ; WSMF8–30  ,B  ), median frequency (MF;  C  ), and spectral edge frequency (SEF95;  D  ). Parameter values are normalized to the interval [0, 1]. The  gray lines  indicate changes of the level of consciousness: loss of consciousness at induction (LOC1: 214 s), return of consciousness after intubation (ROC1: 784 s), followed by second loss of consciousness (LOC2: 846 s). WSMF8–49and WSMF8–30indicate better the state of consciousness than MF and SEF95. Parameter calculation was performed in time steps of 4 s using a signal of 8 s length (overlapping factor 2). 
Image Not Available
Fig. 1. Induction of anesthesia in a randomly selected patient (propofol): Parameter time series based on study A of weighted spectral median frequency (WSMF8–49  ,A  ; WSMF8–30  ,B  ), median frequency (MF;  C  ), and spectral edge frequency (SEF95;  D  ). Parameter values are normalized to the interval [0, 1]. The  gray lines  indicate changes of the level of consciousness: loss of consciousness at induction (LOC1: 214 s), return of consciousness after intubation (ROC1: 784 s), followed by second loss of consciousness (LOC2: 846 s). WSMF8–49and WSMF8–30indicate better the state of consciousness than MF and SEF95. Parameter calculation was performed in time steps of 4 s using a signal of 8 s length (overlapping factor 2). 
×
Fig. 2. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A at fixed cutoff frequencies  f  low= 5 Hz  ,f  high= 50 Hz (  A  ) and  f  low= 9 Hz  , f  high  = 30 Hz (  B  ). For the exponent  p  (  x-axis  ), settings between 0.2 and 2.0 are displayed; for the splitting ratio  r  (  y-axis  ), settings between 0.1 and 0.95 are displayed. The resulting PKvalues of WSMF with corresponding settings of  p  and  r  are plotted as third dimension (  z-axis  ,left  ) and represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . A relation between  p  and  r  can be observed if values of PK≥ 0.80 are assumed (  red area  ). The PKvalues are interpolated giving a continuous approximation of the discrete settings to obtain a better presentation (  left  and  right  ).
Image Not Available
Fig. 2. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A at fixed cutoff frequencies  f  low= 5 Hz  ,f  high= 50 Hz (  A  ) and  f  low= 9 Hz  , f  high  = 30 Hz (  B  ). For the exponent  p  (  x-axis  ), settings between 0.2 and 2.0 are displayed; for the splitting ratio  r  (  y-axis  ), settings between 0.1 and 0.95 are displayed. The resulting PKvalues of WSMF with corresponding settings of  p  and  r  are plotted as third dimension (  z-axis  ,left  ) and represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . A relation between  p  and  r  can be observed if values of PK≥ 0.80 are assumed (  red area  ). The PKvalues are interpolated giving a continuous approximation of the discrete settings to obtain a better presentation (  left  and  right  ).
×
Fig. 3. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A. Settings of the exponent  p  (  x-axis  ) are displayed between 0.4 and 2.4; settings of the cutoff frequencies  f  high(  y-axis  ) and  f  low(  z  -axis) are displayed between 25 Hz and 53 Hz and between 0.5 Hz and 15 Hz, respectively. In  A  and  B  , the PKvalues at the corresponding setting of  f  low  ,f  high  ,p  are represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . A maximal PKvalue of 0.824 results for WSMF with  f  low= 8 Hz  ,f  high= 30 Hz, and  p  = 0.4 that is situated near the frontal edge of the cube (  dark red area  ). The cube in  B  is interpolated, giving a continuous approximation of the discrete settings to obtain a better presentation. The cube in  C  shows only values of PK≥ 0.82, which are represented by  black dots.  They occur as two “clusters,” one as “high-frequency cluster” with  f  high≥ 40 Hz and one as “low-frequency cluster” with  f  high≤ 30 Hz. 
Image Not Available
Fig. 3. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A. Settings of the exponent  p  (  x-axis  ) are displayed between 0.4 and 2.4; settings of the cutoff frequencies  f  high(  y-axis  ) and  f  low(  z  -axis) are displayed between 25 Hz and 53 Hz and between 0.5 Hz and 15 Hz, respectively. In  A  and  B  , the PKvalues at the corresponding setting of  f  low  ,f  high  ,p  are represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . A maximal PKvalue of 0.824 results for WSMF with  f  low= 8 Hz  ,f  high= 30 Hz, and  p  = 0.4 that is situated near the frontal edge of the cube (  dark red area  ). The cube in  B  is interpolated, giving a continuous approximation of the discrete settings to obtain a better presentation. The cube in  C  shows only values of PK≥ 0.82, which are represented by  black dots.  They occur as two “clusters,” one as “high-frequency cluster” with  f  high≥ 40 Hz and one as “low-frequency cluster” with  f  high≤ 30 Hz. 
×
Fig. 4. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A at three fixed levels of the low cutoff frequency  f  lowwith variable settings of the exponent  p  (  x-axis  ) between 0.4 and 2.4 and high cutoff frequency  f  high(  y-axis  ) between 25 Hz and 53 Hz. In  A  , f  low= 0.5 Hz; in  B, f  low= 8 Hz; and in  C, f  low= 13 Hz. The resulting PKvalues of WSMF with the corresponding settings of  p  and  f  highare represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . PKvalues are interpolated.  A  contains the configuration of median frequency (MF) with  f  low= 0.5 Hz,  f  high= 30 Hz,  p  = 2,  r  = 0.5 leading to a PKof 0.55 (all PKvalues smaller than 0.66 are represented by  dark blue  ). B  contains the configurations of WSMF that provide maximal PKvalues in the range of 0.82 within both “clusters” shown in  figure 3C(  dark red areas  ):  f  low= 8 Hz,  f  high= 49 Hz,  p  = 1,  r  = 0.5 (denoted by WSMF8–49) and  f  low= 8 Hz,  f  high= 30 Hz  , p  = 0.4  , r  = 0.5 (denoted by WSMF8–30). 
Image Not Available
Fig. 4. Prediction probability (PK) of weighted spectral median frequency (WSMF) based on development study A at three fixed levels of the low cutoff frequency  f  lowwith variable settings of the exponent  p  (  x-axis  ) between 0.4 and 2.4 and high cutoff frequency  f  high(  y-axis  ) between 25 Hz and 53 Hz. In  A  , f  low= 0.5 Hz; in  B, f  low= 8 Hz; and in  C, f  low= 13 Hz. The resulting PKvalues of WSMF with the corresponding settings of  p  and  f  highare represented as  superposed colors  , where the scaling of the PKvalues is indicated by  gradient bars  . PKvalues are interpolated.  A  contains the configuration of median frequency (MF) with  f  low= 0.5 Hz,  f  high= 30 Hz,  p  = 2,  r  = 0.5 leading to a PKof 0.55 (all PKvalues smaller than 0.66 are represented by  dark blue  ). B  contains the configurations of WSMF that provide maximal PKvalues in the range of 0.82 within both “clusters” shown in  figure 3C(  dark red areas  ):  f  low= 8 Hz,  f  high= 49 Hz,  p  = 1,  r  = 0.5 (denoted by WSMF8–49) and  f  low= 8 Hz,  f  high= 30 Hz  , p  = 0.4  , r  = 0.5 (denoted by WSMF8–30). 
×
Table 1. Prediction Probability Values for Different Settings 
Image Not Available
Table 1. Prediction Probability Values for Different Settings 
×