 Review
 Open access
 Published:
The applied principles of EEG analysis methods in neuroscience and clinical neurology
Military Medical Research volume 10, Article number: 67 (2023)
Abstract
Electroencephalography (EEG) is a noninvasive measurement method for brain activity. Due to its safety, high resolution, and hypersensitivity to dynamic changes in brain neural signals, EEG has aroused much interest in scientific research and medical fields. This article reviews the types of EEG signals, multiple EEG signal analysis methods, and the application of relevant methods in the neuroscience field and for diagnosing neurological diseases. First, three types of EEG signals, including timeinvariant EEG, accurate eventrelated EEG, and random eventrelated EEG, are introduced. Second, five main directions for the methods of EEG analysis, including power spectrum analysis, time–frequency analysis, connectivity analysis, source localization methods, and machine learning methods, are described in the main section, along with different submethods and effect evaluations for solving the same problem. Finally, the application scenarios of different EEG analysis methods are emphasized, and the advantages and disadvantages of similar methods are distinguished. This article is expected to assist researchers in selecting suitable EEG analysis methods based on their research objectives, provide references for subsequent research, and summarize current issues and prospects for the future.
Introduction
Neuroscience, also called brain science, is a discipline that explores the structure and function of the brain [1]. It is a typical interdiscipline that involves multiple disciplines, such as biology, psychology, information science, medicine, engineering, and artificial intelligence. Neuroscience has been developed for about 100 years and extensively applied to diagnose neurological disorders. With the development of research methods, the focus of neuroscience has gradually transitioned from brain structure to brain function in the past decade [2,3,4]. The brain nerve response is known as the core of cognitive production. Accurate identification of the brain nerve response can contribute to identifying important human cognitive functions, developing intelligent algorithms, and advancing medical developments regarding neurological diseases [3].
With the development of research tools for neuroscience, multiple neuroimaging tools are available for exploring brain function, including electroencephalography (EEG)/intracranial electroencephalography (iEEG), functional magnetic resonance imaging (fMRI), magnetoencephalography (MEG), positron emission tomography, and optogenetic techniques [4,5,6,7]. Among these methods, EEG/iEEG has been the most widely used tool for functional brain imaging due to its excellent temporal resolution and low equipment cost [8]. From the perspective of neurophysiology, EEG/iEEG reflects postsynaptic potential, which is generated when neurotransmitters bind to receptors on the postsynaptic membrane [9]. These postsynaptic potentials generate electric fields around neurons. Once sufficient neurons are activated, electroneurographic signals with specific patterns can be captured through a voltage amplifier. Owing to the shorter spatial distance between iEEG and neuronal groups, iEEG has higher accuracy and signaltonoise ratio compared to EEG [9]. The electric signals captured by EEG have poor spatial resolution and signaltonoise ratio since they are passed through the skull. However, EEG is a noninvasive technology, so it can be applied in a wider range of scenarios [8]. The information captured by both methods is the discharge of neuronal groups, so the capture equipment in both cases is a voltage amplifier, and the captured signals have basically the same manifestation. Hitherto, EEG/iEEG has been extensively applied to research diverse aspects of brain function, including attention [10], memory [11], language [12], emotions [13], and brain function disorders [14].
Although EEG/iEEG has good practicality, its application requires a certain foundation in signal processing technologies due to the complex representation of EEG signals; this leads to a problem in that some researchers lack clarity in selecting and applying analytical methods for EEG/iEEG. Therefore, we attempted to provide a brief introduction to commonly used EEG signal processing methods in this article. In this review, we first comprehensively classified EEG signals based on their characteristics. Next, we introduced commonly used analytical methods for EEG in terms of characteristics such as power spectrum and connectivity, and presented their adaptability to various types of EEG to assist researchers in method selection. We also summarized current issues and prospects for the future, which is expected to expedite the application of EEG/iEEG in brain science and neurological disease research.
Types of EEG signals
Generally, in research articles, especially those on neurological disease, EEG is classified based on the research subjects. For instance, in sleep study, EEG is classified into EEG during wakefulness and sleep EEG [15]; in epilepsyrelated study, EEG is subdivided into interictal EEG, preictal EEG, ictal EEG, and postictal EEG [16]; and in research on eventrelated potentials, EEG is categorized into restingstate EEG and taskstate EEG [17]. In addition, it can be classified according to the shape of the EEG itself. For example, it can be divided into delta, theta, and alpha based on frequency [18, 19] or slow wave, fast wave, sharp wave, and spike wave based on shape.
However, from the perspective of EEG analysis, we believe that EEG can be classified into the following three categories. (1) An EEG in which the functional state of the brain remains unchanged over time is called a timeinvariant EEG for short [20]. In this type of EEG, the state of the brain does not show significant changes during the capture process, for example, a restingstate EEG without psychological activity [16]. Alternatively, some changes in brain characteristics are not included among the main features to be studied. For example, in epilepsy research, researchers pay more attention to the pathological EEG; in this case, the interictal period without epileptic discharge can also be considered a timeinvariant EEG [16]. Figure 1 is an example of sleep EEG. During sleep, the EEG is in a stable state for a long time. Figure 1a shows an EEG of a 5s period of sleep, while Fig. 1b shows a 150s sleep EEG. Although the EEG is unstable, relatively stable data segments can be found within the unstable EEG by analyzing these two segments. (2) Accurate eventrelated EEG can be regarded as an extension of eventrelated potentials; it refers to the EEG induced by a certain event where the induction time of the event is definite. EEG with timevarying characteristics caused by stimuli with a definite time, such as visuoauditory, transcranial magnetic stimulation, and electrical cortical stimulation [17]. Figure 2 is an example of eventrelated potential. Events will appear at a clear time, so accurate brain electrical responses can be obtained through time. Figure 2a is a section of an EEG containing eventrelated potentials, which marks the exact moment when the event occurred. Figure 2b shows the EEG response after superimposing multiple eventrelated potentials with the event as the 0 moment. Figure 2c shows the EEG response after the superposition of multiple eventrelated potentials, with the reaction time as the 0 moment. (3) Random eventrelated EEG refers to an EEG induced by a certain event in which the induction time of the event is random and cannot be determined. In research on diseases such as epilepsy or Parkinson’s disease, pathological EEG is triggered by abnormal neural activity in the lesion area, but the timing of pathological induction is difficult to determine, resulting in a timevarying EEG [19]. Figure 3 displays the EEG signals of an epilepsy patient. Figure 3a shows the EEG signal during the interictal period, while Fig. 3b, c show the EEG signal in the early and late stages of the seizure, respectively. Figure 3d shows the EEG signal of the entire seizure process. Epilepsy is a random event, so the time of occurrence of the event needs to be retrospectively located after event onset, which poses challenges to the realtime analysis of epilepsy EEG data. However, Fig. 3a, c show that epileptic EEG is still relatively stable within a period. In Fig. 3d, the data observed during the interictal and postictal periods, which represent two stable stages, reveals a significant difference.
It should also be noted that the main classification criteria for these three types of EEG were based on the EEG features to be analyzed, and specific analysis is required according to the features of interest. In sleep disorder research, we can consider stage N1 EEG as a timeinvariant EEG. However, if the study target is a dream or memory, the main characteristics to be studied may also change during the N1 phase, and this EEG should be classified as a random eventrelated EEG. Therefore, researchers should accurately identify the target features to be studied before selecting an analytical method.
Common EEG analysis methods
In this section, the extraction methods of common brain functional features based on the characteristics of EEG signals are introduced.
Power spectrum analyses
The power spectrum is a very common analytical method in EEG analysis that can analyze the energy changes of various frequency components in EEG signals. This method can be applied in studies on brain science and neurological diseases that can trigger EEG energy changes upon state changes, such as sleep stage changes, seizures, and emotional changes. Multiple power analysis methods are available for selection, such as the fast Fourier transform (FFT), Welch, and autoregressive (AR) model, with different characteristics and usage limitations. The articles in Table 1 covered a range of topics related to EEG analysis, including sleep onset, transitions between sleep stages, classification of neurological disorders, detection of poststroke EEG signals, analysis of EEG background activity in autism and dyslexia, and the impact of various factors such as focused ultrasound stimulation, cognitive impairment in diabetes, and neurofeedback training in autism [20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41].
FFT
FFT is a fast algorithm for computing discrete Fourier transform (DFT) [42]. DFT:
Here, \(X\left(k\right)\) denotes the DFT, \(N\) represents the length of the available data, \(x\left(n\right)\) refers to the input signal in the time domain, \(e\) signifies the exponential operation, \(i\) denotes the imaginary part, and \(k\) represents the sampling frequency. The calculation process of Equation (1) is known as the FFT algorithm. Using the symmetric and periodic nature of the exponential factor in the DFT calculation equation, FFT can reduce repetitive calculations [42]. FFT calculations have a highfrequency resolution but are also easily affected by EEG signal noise, so an average period method has been proposed for improvement.
The average period method splits the original signal into N nonoverlapping consecutive segments and then calculates the periodogram of each segment individually, all of which are finally averaged [43]. Through the window technique, averaging N periodograms can reduce the variance of the power spectral density estimation, but spectral leakage easily occurs because signal segmentation leads to increased boundaries of the data, while Fourier transform has a poor ability to process these data boundaries [44].
Welch
The Welch method has two improvements that enhance resolution and reduce errors in results compared to the average period method. First, this method allows overlap between data segments. Second, the Hamming window function is used for each segment instead of the rectangular window function, which ameliorates the potential distortion caused by too many rectangular windows [45]. \({\{x}_{l}\left(n\right)\},l=1,\dots ,S\) refer to data segments, and \(M\) represents the length of each segment. The overlapping coefficient is usually set as 50% (\(M\)/2). The Welch spectrum estimate is given by the following equation:
Here, \({\widehat{P}}_{l}\left(f\right)\) represents the periodogram estimate of segment \(l\) and \(v\left(n\right)\) denotes the window function. \(P\) refers to the general average of \(\left{v}_{\left(n\right)}^{2}\right\): \(P=1/M{\Sigma }_{n=1}^{M}{\left{v}_{\left(n\right)}^{2}\right}^{2}\), \(e\) signifies the exponential operation, \(i\) denotes the imaginary part, \(f\) represents the sampling frequency, \({\widehat{P}}_{w}\left(f\right)\) refers to the Welch power spectral density estimate, and \(S\) signifies the number of segments.
Currently, Welch’s method is one of the most widely used power spectrum analysis methods because it reduces the influence of boundary effects on the power spectrum, providing more stable power spectral results than FFT/shorttime Fourier transform (STFT) methods.
Multitaper
Multitaper can solve the bias and variance problems of nonparametric spectral estimation simultaneously in an optimal manner [46]. Windowing the signal using different tapers allows multiple independent estimates to be derived from the same signal since the different windows are uncorrelated with each other.
Assuming that \({X}_{k}={\{{x}_{1,k},{x}_{2,k},\dots ,{x}_{p,k}\}}^{T}\) is the signal sequence, where p denotes the number of channels and k signifies the length of the sequence, the multitaper of the channel \(l,m\) is calculated as follows:
Here, K refers to the crossspectrum estimate, N represents the sequence length, \({\widehat{S}}_{l,m}^{\left(i\right)}\left(f\right)\) signifies the \(k\)th direct crossspectrum estimate of the channel \(l,m\), \(\Delta\) indicates the sampling interval, \(f\) represents the sampling frequency, \(e\) signifies the exponential operation, \(i\) denotes the imaginary part, and \({h}_{k}^{\left(i\right)}\) represents taper. Since the final result is obtained by processing multiple tapers, the problem of information loss caused by singlescale analysis can be reduced.
The multitaper method is a modified Welch’s method that provides features similar to those of the STFT and Welch’s methods, but its stability is improved and the number of parameters to be determined is reduced because it uses multiple tapers for superposition. For instance, in the wellknown tool Fieldtrip, multitaper is employed as the major power spectrum analysis method [47].
AR model
The AR model can achieve the linear prediction modeling of the original signal as with the signal \(x\left(n\right),0\le n\le (N1)\) deemed as white noise with a mean value of 0 and a variance equal to \({\sigma }^{2}\). The signal’s amplitude during a specific period is determined by aggregating the various amplitudes from preceding signals and incorporating the estimation error. The model’s order, or filter, is contingent upon the quantity of AR coefficients employed.
Here, \(a\left(k\right)\) represents the coefficients of the AR model, \(w\left(n\right)\) signifies white noise with a variance equal to \({\sigma }^{2}\), and \(p\) refers to the model order. The AR(P) model can be characterized in terms of parameters \(\left\{a\left[1\right],a\left[2\right],\dots ,a\left[p\right],{\sigma }^{2}\right\}\). Power spectral density:
\(A\left(f\right)=1+{a}_{1}\mathrm{exp}\left(j2\pi f\right)+{\dots +a}_{p}\mathrm{exp}\left(j2\pi fp\right).\)AR model parameters can be derived using Burg or least squares [48]. The AR model can process the shortterm signal in contrast to the FFT. The AR model is less used because of its high coefficient requirement, but its scalability is so excellent that it has been used in many researches [37,38,39,40,41].
Time–frequency analyses
Since the EEG of human beings in a task state shows timedependent changes, a time–frequency analysis is quite suitable for EEG analysis. For shortterm signals, time–frequency analysis can replace power spectrum analysis to characterize the signals in two dimensions. The commonly used methods for time–frequency analysis include STFT, the wavelet transform (WT), empirical mode decomposition (EMD), and the WignerVille distribution (WVD). The articles in Table 2 investigated various methods for EEG signal analysis, including the use of rational discrete STFT and deep learning for epileptic seizure classification, a hybrid approach for alcohol and control EEG signal classification, connectivity analysis in autism disorders using STFT and coherence values, drowsiness detection based on relative band power and STFT, automatic sleep stage classification using time–frequency images, and detection of deception using smoothed pseudo WVD, among other techniques [49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68].
STFT
The STFT is a technique that divides longterm signals into shorter segments of uniform length. It then computes the Fourier transform separately for each of these shorter segments. The Fourier transform is defined as follows:
\(f\left(t\right)\) refers to the original signal, \(\tau\) signifies the translation parameter, \({\psi }^{*}\left(t\tau \right)\) denotes the window function (usually Hamming window), and when the window function uses a Gaussian function, the STFT is called a Gabor transform. Moreover, \(e\) signifies the exponential operation, \(i\) denotes the imaginary part, \(t\) refers to time, and \(\omega\) represents frequency. The STFT has a limitation in that its fixed time window results in a fixed time–frequency resolution [69].
WVD
The WVD is a classical method for time–frequency analysis that excels in handling nonstationary signals. Unlike STFT, the WVD remains unaffected by leakage effects. This distribution represents a secondary energy density, derived by correlating the signal with time and frequency translations along with their complex conjugates. The instantaneous autocorrelation function of the signal \(f(t)\) is as follows:
\(t\) refers to time, \(\tau\) signifies the timelag correlation coefficient, \(*\) denotes the complex conjugate, and the WVD of \(f(t)\) refers to the Fourier transform of \({R}_{f}\left(t,\tau \right)\) about \(\tau\).
\(e\) signifies the exponential operation, \(j\) denotes the imaginary part, and \(\omega\) represents frequency.
The WVD has a series of good properties such as conjugate symmetry, timemarginal properties, frequencymarginal properties, and energy distribution properties. However, when the signals have multiple frequency components, the WVD is affected by cross terms, that is, it is easily affected by noise [70].
WT
The WT overcomes the time–frequency resolution limitation observed in the STFT algorithm. This is achieved by introducing varying time–frequency resolutions in the outcomes through the translation and scaling of wavelets. Wavelets:
\(\psi \left(t\right)\) refers to the mother wavelet, \({\psi }_{a,b}\left(t\right)\) signifies the subwavelet, \(a,b\) refer to the modulation and translation parameters, respectively, and the WT of the signal \(f(t)\) is as follows:
Now, there are many optional mother wavelet functions (such as morse, morlet, db, and Harr), which include discrete and continuous wavelets. For EEG signal analysis, discrete wavelets are commonly used for signal decomposition, and continuous wavelets are commonly used for signal presentation. Thus far, although continuous wavelets cause massive data redundancy, they have been the preferred time–frequency analysis method, with more accurate and smooth time–frequency representation [69].
EMD
EMD is a selfadaptive multiresolution technique that decomposes the original signals into components of different resolutions and can analyze nonlinear and nonsmooth signals [71]. EMD decomposes the input signals into several intrinsic mode functions (IMFs) and a residual:
\(I\left(n\right)\) refers to a multicomponent signal, \({IMF}_{M}\left(n\right)\) signifies the Mth IMF, and \({Res}_{M}\left(n\right)\) denotes the corresponding residual intrinsic modes. The IMF components usually extract time–frequency features using the Hilbert spectral analysis. EMD is characterized by selfadaptability and high efficiency. However, it may exhibit aliasing effects due to the presence of IMFs that contain significantly different characteristic time scales or when similar characteristic time scales are dispersed across different IMFs. EMD does not rely on the primary function of the fixed frequency, so the time–frequency results obtained by EMD are affected by Gibs; however, its positioning performance to the frequency is poor [72].
Connectivity analyses
Connectivity analysis of EEG is an analytical method that has gained much attention in recent years and is fundamental to research on brain networks and connectivity. Connectivity analysis includes multiple types, such as signal morphologybased, signal phasebased, statisticsbased, and informationbased analyses. Different correlation analysis methods are oriented to different principles and the obtained results express different characteristics. Hence, it is important to accurately choose suitable connectivity analysis methods during brain network research. We introduce several common connectivity analysis methods here. Table 3 summarized some articles that explored various aspects of EEG signal analysis, including analysis for different severities of obstructive sleep apnea, synchrony measures for early Alzheimer’s disease diagnosis, correlation between EEG abnormalities and symptoms of autism spectrum disorder, EEG channel correlation for emotion recognition, quantitative EEG in ischemic stroke correlation with functional status, variability of EEG functional connectivity in young attention deficit hyperactivity disorder subjects, and the identification of causal relationships between EEG activity and intracranial pressure changes in neurocritical care patients, among other topics [73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106].
Correlation (CORR)
CORR measures the similarity between two signals by calculating the variance of signals [107]. The CORR for each given frequency is as follows:
\({C}_{AB}\left(x\right)\) represents the crosscovariance between signal \(A\) and signal \(B\), and \({C}_{AA}\left(x\right)\) and \({C}_{BB}\left(x\right)\) refer to the autocovariance of signal \(A\) and signal \(B\), respectively. CORR is sensitive to both phase and polarity.
Coherence (COH)
COH measures the similarity between two signals by calculating the power spectral density [108]. The COH for each given frequency is:
\({S}_{AB}\left(x\right)\) represents the crossspectrum between signal \(A\) and signal \(B\), and \({S}_{AA}\left(x\right)\) and \({S}_{BB}\left(x\right)\) refer to the autocovariance of signal \(A\) and signal \(B\), respectively. Because COH is calculated through crossspectrum and autospectrum, it is very sensitive to the phase changes of the signal but is little affected by energy changes.
Wavelet coherence (WTC)
WTC can represent the timevarying relationships between different signals in the time–frequency domain by producing different time–frequency resolutions through wavelet translation and dilation [109].
The WT of signal \(x\) is a function of time and frequency, defined as the convolution of an input with a wavelet family \(\theta (u)\):
With given input signals \(x\) and \(y\), the wavelet crossspectrum around time \(t\) and frequency \(f\) can be derived through the WT of \(x\) and \(y\):
Here, \(*\) represents the complex conjugate, and \(\delta\) denotes the scalar dependent on the frequency. The WTC of time \(t\) and frequency \(f\) is represented by \(C{W}_{xx}\left(t,f\right)\), and \(C{W}_{yy}\left(t,f\right)\) refers to the Fourier transform of the autocorrelation function of signal \(x\) and signal \(y\). WTC can view the phase correlation between signals on the time spectrum and reduce the interference of energy.
Phase lag value/index (PLV/PLI)
PLV and PLI are commonly applied to acquire the strength of phase synchronization [110]. The instantaneous phase of signal \(x(t)\) is generated using the following formula:
Here, \(\widetilde{x}\left(t\right)\) signifies the Hilbert transform of \(x\left(t\right)\), defined as follows:
PV refers to the Cauchy principal value. The PLV of two signals is defined as follows:
Here, \(\Delta t\) denotes the sampling period, \(N\) represents the number of samples per signal, \(j\) refers to the imaginary part, and \(e\) signifies the exponent. PLV signifies phase synchronization, with values ranging from 0 to 1. A value of 0 indicates a lack of synchronization, while 1 represents strict phase synchronization. On the other hand, PLI characterizes the asymmetry in the phase difference distribution between two signals. It is computed based on the relative phase difference between the two signals:
\(E\) represents expectation, the result value is located within the interval [0, 1], and a higher value indicates a higher phase synchronization.
Mutual information (MI)
MI is designed based on information theory, which presents how one signal provides information for another signal [109]. \(P({x}_{j})\) and \(P({y}_{j})\) are the probability distributions of signal \(X=\{{x}_{j}\}\) and signal \(Y=\{{y}_{j}\}\), respectively. The entropy of \(X\) and \(Y\) is defined as follows:
N signifies the window length. \(H\left(\left.Y\rightX\right)\) and \(H\left(X,Y\right)\) refer to the conditional entropy and joint entropy between \(X\) and \(Y\), respectively, which are defined as:
Here, E denotes the expected value function. The MI of signal X and signal Y is calculated as follows:
MI can simultaneously detect the linear and nonlinear correlations between two signals, but it requires mass data.
Granger causality (GC)
GC is a linear vector AR model based on random timeseries data, which can estimate effective interactions from timeseries data [111]. For this method, if the past value of the signal \({X}_{1}(t)\) contains information that contributes to the prediction of \({X}_{2}(t)\), which exceeds the information contained only in the past value of Y, the signal \({X}_{1}(t)\) “Granger causes” the signal \({X}_{2}(t)\). Therefore, the bivariate AR model is as follows:
\(p\) refers to the maximum number of delayed observations, \(j\) denotes the lag coefficient, the matrix \(A\) represents the contribution of each delayed observation to the predicted signal value, and \({E}_{1}\left(t\right)\), \({E}_{2}\left(t\right)\) signify the residual of each time series. GC can only provide information on the linear characteristics of the signal and cannot analyze nonlinear situations.
Crossfrequency analysis (CFA)
CFA is a kind of rapidly developing connectivity analysis method, which mainly includes crossfrequency coupling (CFC) and crossfrequency directionality (CFD).
CFC describes the interaction of brain oscillations across different frequency bands and manifests in four modes: phasetoamplitude, powertopower, phasetophase, and phasetofrequency interactions. The Kullback–Leibler distance serves as an effective metric for quantifying CFC [112]. Notably, CFC holds significance in working memory processes [113]. According to the theta/gamma neural code hypothesis, conserved memory items are encoded through thetanested gamma cycles in sensory regions, facilitating communication between different brain cortices during memory and sensory processes [114]. A study leveraging iEEG data in epilepsy patients, coupled with behavioral outcomes, underscore the association between theta/gamma CFC across diverse brain regions and working memory performance [115]. Key findings reveal the widespread distribution of theta/gamma phase amplitude coupling across the cortex, with increased coupling strength observed in more cognitively demanding working memory tasks [116].
CFD, measuring information flow direction between brain regions, involves the modulation of highfrequency signal amplitude by the phase of a lowfrequency signal [117]. It relies on the phase slope index, quantifying the phase slope in the crossspectrum of two signals [117]. CFD has proven valuable for inferring causal relationships and estimating signal delays [118]. Additionally, it has been employed in exploring information flow directions between distinct brain regions during various cognitive tasks [117].
Source localization analysis
With the continuous development of EEG and MEG devices, the number of channels in scalp EEG or MEG has increased to over 100. Multichannel and multilocation EEG/MEG signals have accelerated the development of EEG source localization.
First, structural MRI is often used as a prior in source localization analysis because it provides a highresolution threedimensional (3D) image of the brain’s anatomy. This image can be used to create a head model that accurately represents the geometry and conductivity of the brain and skull [119]. The head model is then used to calculate the forward solution, which describes how electrical activity generated by the brain is measured at the scalp [120]. By incorporating structural MRI information into the forward solution, the accuracy of the source localization can be improved.
Moreover, structural MRI has the potential to generate an accurate boundary element model of the head, facilitating the computation of the lead field matrix [121]. This matrix characterizes the propagation of electrical activity generated by the brain to the scalp electrodes [122]. The utilization of a realistic boundary element model enhances the precision of lead field matrix calculations, thereby improving the accuracy of source localization [123].
The source localization method can infer the intracranial discharge status of the brain through multichannel signals from the scalp, human brain physical models, and finite element calculations. Source localization methods are commonly employed to localize functional areas and lesion areas, among others, under noninvasive conditions. Common source localization methods are introduced below. The articles in Table 4 explored EEG source localization techniques, including dipole analysis, beamforming approaches, and methods like lowresolution electromagnetic tomography (LORETA) and standardized LORETA (sLORETA), to study various conditions such as epilepsy, visual working memory tasks, auditory attention, depression, obsessive–compulsive disorder, pain perception, agerelated hearing loss, and different neurological disorders, providing insights into the localization of brain activity in these contexts [124,125,126,127,128,129,130,131,132,133,134,135,136,137,138,139,140,141,142,143,144,145,146,147,148,149,150,151].
Minimum norm estimation
The minimum norm estimation method uses MEG for analysis and solves the current distribution by estimating the linear combination of the magnetometer lead field. \({L}_{i}\) signifies the vector field at the position \(i\), so the output of the magnetometer is defined as follows:
\(J\left(r\right)\) denotes the conversion of various energy types into electrical energy, and the linear relationship among the magnetometer reading, current distribution, and lead field is expressed as:
Consequently, the shortest current vector needed to explain the magnetometer output is defined by multiplying the output vector \(B\) by the pseudoinverse of \(L\):
Here, \({L}^{+}={L}^{T}{(L{L}^{T})}^{+}\) represents the Moore–Penrose generalized inverse, predicting minimum norm solutions for pure signals, noisecontaminated signals, and smoothed noise signals. Due to the harmonic nature of the minimum norm solution, the method faces challenges in resolving deep source localization within the outermost cortex, leading to localization errors [152].
Focal underdetermined system solution (FOCUSS)
FOCUSS is a highresolution nonparametric technique that allocates current to each element within a predetermined reconstruction region using a forward model [153]. The weighted minimum norm method is used to perform mathematical calculations in the recursive steps in focusing. The calculation formula for the unknown current element \(I\) is as follows:
Here, \(W\) is an \(n \times n\) matrix that refers to a constraint on the results to strengthen some elements in \(I\), \(B\) denotes the measured value of the radial magnetic field, and \(G\) signifies the spatial weight of the element:
\({I}_{ik1}\) represents the \(i\)th element of vector \(I\) in the \((k1)\) iteration, and \(k\) signifies the index of the iteration step. By continuously constructing \(W\) and calculating the weighted minimum norm, the model results are converged, but the computation time of FOCUSS is longer than that of other algorithms.
LORETA
LORETA is an innovative method in the high temporal resolution neuroimaging field that allows for the 3D reconstruction of the EEG activity distribution [154]. A head model is used for LORETA, and the intensity and direction of electrical activity at each point determine the electromagnetic field measured on the scalp. It is defined as:
with
In the above equation, \(\Phi\) represents a vector of potential difference, \(K\) denotes the lead field matrix of the volume, \(J\) signifies the current density, \(W\) denotes the discrete Laplace operator in the square space, and \(\alpha\) refers to the Tikhonov regularization parameter.
The \({T}_{W}\) value can be calculated using the following formula:
\(H\) denotes the mean reference operator, which is realized using the discrete spatial Laplacian operator, so the spatial resolution of LORETA is relatively low.
sLORETA is also a common and popular source localization method. sLORETA incorporates additional assumptions regarding the smoothing and weighting of the values [155]. An advantage of sLORETA is that it has “guaranteed accuracy” in the presence of a single dipole, while LORETA does not [155]. sLORETA has been used in various studies to estimate the sources of EEG signals in the brain [156]. For example, sLORETA has been used to study the neural correlates of cognitive processes such as attention, memory, and language [156]. sLORETA has also been used to study the neural correlates of various disorders such as depression, schizophrenia, and Alzheimer’s disease [156].
Dipole
The dipole method can predict the electric field generated by a theoretical dipole in the brain using dipole propertyrelated principles [157]. Location and orientation are two parameters of the dipole model, with the location indicating the position of active region within the brain in this model and the orientation indicating the arrangement of brain cells in the active region.
The six parameters of a dipole source consist of three coordinates in \({r}_{d}\in {R}^{3\times 1}\) and three dipole components in \(d=({d}_{x},{d}_{y},{d}_{x})\in {R}^{3\times 1}\) (equivalently two orientation angles and an intensity parameter). For each dipole position \({r}_{d}\) within the head, the relation between \(d\) and the potential measured at the \(m\)th electrode \({V}_{mod}\in {R}^{m\times 1}\) can be written as:
The matrix \(L\in {R}^{m\times 3}\) is a lead field matrix, determined by dipole position, electrode position, and head geometry.
A more realistically shaped head model is often required for patient EEG data analysis, but in this case, boundary element methods or numerical methods such as the finitedifference method are needed to compute the lead field matrix.
Beamforming
Beamforming is a spatial filtering technique for signals measured by discrete sensors [158]. Beamforming refocuses the signals captured on the scalp to their original locations by finding the weights of each location in the source space, thus minimizing the variance of the current dipole at each location. It is often desirable to extract signals from a small region of the brain that is modeled by dipoles at the location \({r}_{d}\) with a specific orientation. With a given dipole and its components, the potential distribution is defined as follows:
\({r}_{d}\) denotes the dipole coordinate, \(d\) represents the dipole component, \(L\) signifies the lead field matrix, and \(w\) refers to the weight vector.
The output variance or output power of a beamformer is calculated as follows:
\(y(k)\) represents the output and \(\varepsilon \left\{{\cdot }^{2}\right\}\) denotes the expected value of its parameter. The results are constrained with different restrictions.
Current source density (CSD)
CSD calculates an estimate of the current projected radially from the underlying neuronal tissue at a given surface location to the skull and scalp and calculates a spatially weighted sum of the potential gradients pointing to that location from some or all of the recorded locations.
CSD estimates:
Here, \(C\left(E\right)\) denotes the current density value at any point \(E\) on the sphere surface, \({c}_{i}\) refers to a constant to express an \(i\) surface potential set, and \(\mathrm{cos}\left(E,{E}_{i}\right)\) refers to the cosine of the angle between the surface point \(E\) and the electrode projection \({E}_{i}\). The function \(h\left(x\right)\) is defined as the sum of the grades:
Here, \(m\) is a constant greater than 1 and \({P}_{n}\) is the \(n\)th Legendre polynomial, defined as follows:
CSD does not necessitate reference information but is susceptible to noise. CSD is a source localization method designed for scalp EEG, treating the entire head as a conductor with equal conductivity. It concentrates signals from multiple EEG channels to their respective channels by adjusting the parameters. This method cannot focus EEG signals to the intracranial region. However, its arithmetic is simple and fast, and it is still used in some EEG analyses for scalp localization without the assistance of brain models [159].
Machine learning
Machine learning is a very popular class of signal processing methods currently applied in the medical field [160], and with the rapid development of deep learning, machine learning methods in EEG analyses have gained attention. Machine learning methods are commonly used for classification and regression problems in EEG analyses and have yielded substantial results in disease research [161]. The studies in Table 5 utilized various machine learning and signal processing techniques, including common spatial pattern (CSP), deep learning, wavelet analysis, support vector machine (SVM), convolutional neural network (CNN), recurrent neural network (RNN), and long shortterm memory (LSTM) network, to address diverse applications such as seizure detection, diagnosis of neurological disorders (autism, schizophrenia, Parkinson’s disease), mental fatigue measurement, and emotion recognition using EEG signals [162,163,164,165,166,167,168,169,170,171,172,173,174,175,176,177,178,179,180,181,182,183,184,185].
CSP
The CSP algorithm uses a linear transformation to maximize the variance ratio of two signals after mapping, which is a common spatialfiltering algorithm used for multichannel EEG analysis [186].
\({X}_{1},{X}_{2}\) refer to the signal data of \(\left(n,{T}_{1}\right),(n,{T}_{2})\) size, where \(n\) is the number of channels and \({T}_{1},{T}_{2}\) are the length of the respective signal:
\(w\) denotes the projection matrix, which can be solved using matrix diagonalization.
In contrast to other spatial feature extraction methods, the CSP method is simple and efficient, but it is only suitable for processing two categories of signal data.
Linear discriminant analysis (LDA)
LDA, a classical linear method, is mainly used to find features that characterize or separate two classes and is also applicable for the dimensionality reduction of data [186]. Regarding projection, the projected data have high cohesion and low coupling characteristics.
\({S}_{1},{S}_{2}\) are intraclass scatters, \({\left({M}_{1}{M}_{2}\right)}^{2}\) refers to the interclass scatter, and \(\widehat{w}\) represents the mapping matrix. Because LDA assumes that the data obey the Gaussian distribution, it does not perform satisfactorily in processing data with nonGaussian distributions.
SVM
SVM is a class of generalized linear classifiers for the classification of binary data in a supervised learning manner [176]. SVM constructs a hyperplane in highdimensional space to distinguish between two classes of data. Assuming that the dataset is \([\left({x}_{1},{y}_{1}\right),\left({x}_{2},{y}_{2}\right),\dots \left({x}_{n},{y}_{n})\right]\), wherein \({y}_{i}\in [1, 1]\), the hyperplane is defined as:
The plane separating the two classes of data is as follows:
Here, \({w}^{T}\) represents the normal vector and \(b\) denotes the offset, so the data interval is \(2/\Vert w\Vert\). This method maximizes \(2/\Vert w\Vert\) while ensuring that all data satisfy the conditions. Methods such as Lagrangian duals can be used to solve such constrained optimization problems, computing the hyperplane of the separated data. SVM performs poorly in resolving multiclassification problems.
CNN
The unique convolutional layer of CNN can effectively extract EEG signals and structural information in the spatiotemporal frequency domains [187,188,189].
For feature extraction, the dot product is completed using the input data with the filterbank regionbyregion, and each kernel is scanned using step length, sharing equal weight. The resulting output is a set of Kdimensional feature maps.
Here, \({B}_{j}^{l}\) signifies the \(j\)th deviation in the layer \(l\), \({W}_{ j,i}^{l}\) refers to the weight matrix connecting with the feature map in the neighboring layer (\({Z}_{j}^{l},{Z}_{j}^{l1}\)), \(*\) represents the convolution operator, and \(\sigma (\cdot )\) denotes the nonlinear activation function.
The extracted feature maps are recognized by a classifier, which often uses the crossentropy loss function:
Here, \({y}_{i}\) represents the sample’s true value and \({\widehat{y}}_{i}\) signifies the modelpredicted value. CNN requires less preprocessing than other algorithms but also has the risk of overfitting.
RNN
The RNN model performs well for temporal signals, wherein connections between nodes generate directed or undirected graphs along the time series, effectively extracting feature information in the time dimension [190]. However, gradient explosion and gradient vanishing problems are present due to RNN’s structure of backpropagation through time. Later, LSTM was developed, which has broader applications than RNN.
LSTM
LSTM has improved on the problems of the RNN network [191] and selectively transmits data utilizing forgetting gates, input gates, and output gates.
The forgetting gate determines which information to remove from the state of the unit:
The input gate determines which values will be updated:
Then, the unit value state is updated based on the equations above:
Finally, the output gate determines which parts of the unit state will be the final output:
wherein \(\sigma\) signifies the sigmoid activation function that compresses numbers to the range 0, 1, \(tanh\) denotes the hyperbolic tangent activation function that compresses numbers to the range − 1, 1, \({W}_{f}\), \({W}_{i}\), \({W}_{c}\), and \({W}_{o}\) are the weight matrixes, \({x}_{t}\) represents the input vector, \({h}_{t1}\) represents past hidden states, and \({b}_{f}\), \({b}_{i}\), \({b}_{c}\), and \({b}_{o}\) are deviation vectors. LSTM has a slow training speed due to its performing and processing difficulties [192].
Joint application of EEG analysis and machine learning methods
The joint application of EEG analysis and machine learning methods has been an active area of research in neuroscience and disease diagnosis. EEG is a noninvasive method of measuring the electrical activity of the brain, and machine learning algorithms can be used to extract information from EEG signals to help diagnose various disorders and identify different brain states [193]. Machine learning algorithms have been developed to extract features from EEG signals, such as frequency bands, time–frequency representations, and connectivity measures [193]. These features can then be used to train machine learning models to classify different brain states or diagnose various disorders [193]. Machine learning algorithms have been developed to detect seizures in EEG signals with high accuracy [194] and classify EEG signals from patients with Alzheimer’s disease and healthy controls [195].
The joint application of EEG analysis and machine learning methods has several advantages in neuroscience and disease diagnosis. It allows for the identification of patterns in EEG signals that are difficult to detect using traditional methods. Machine learning algorithms can be used to extract features from EEG signals that are not easily visible to the human eye, such as subtle changes in frequency or amplitude [193]. These features can then be used to train machine learning models to classify different brain states or diagnose various disorders. Another advantage of this combination is that it can help reduce the subjectivity of EEG analysis. Traditional EEG analysis methods rely on visual inspection of the EEG signal by a trained expert, which can be timeconsuming and subjective [193]. Machine learning algorithms can be used to automate the process of EEG analysis, reducing the time and subjectivity involved in the analysis [194].
Moreover, different models may be better suited to different aspects of the data. One model may be better at detecting certain types of patterns in the data, while another model may be better at classifying the data into different categories. By combining the strengths of different models, it is possible to create a more accurate and robust analysis [193] to reduce the risk of overfitting, and create a more generalizable analysis [196].
In summary, the joint application of EEG analysis and machine learning methods has great potential for the diagnosis of various disorders and the identification of different brain states. It has several advantages, including the ability to identify patterns in EEG signals that are difficult to detect using traditional methods, and the ability to reduce the subjectivity of EEG analysis.
Discussion
This article reviewed several commonly used EEG/iEEG analysis methods in neuroscience and introduced the applied principles based on the data generation characteristics. Due to different EEG data generation approaches, there are fundamental differences in data time, events, and variability. Therefore, the methods for data analysis should be selected based on these characteristics to ensure theoretical accuracy. Figure 4 presents a summary diagram of method selection. The required method can be selected based on the characteristics of the EEG signal and application requirements. The various methods are discussed below.
The power spectrum analysis method is used to reflect the energy changes in various brain regions. FFT has a highfrequency resolution and accuracy but is easily affected by noise and requires a large amount of data [197]. Therefore, FFT is more suitable for analyzing timeinvariant EEG signals of good quality [44]. Welch and multitaper can suppress noise in EEG signals using the window averaging method, but the frequency resolution is decreased and the requirement for data length is increased [198]. Hence, for timeinvariant EEG, if the signals contain slight burst noise, Welch or multitaper is a good choice [45]. Additionally, for accurate eventrelated EEG, EEG can be segmented into epochs based on the onset time of events. Windowed superposition analysis of EEG in the same state, such as baseline EEG before stimulation, can also yield good power spectrum results. If the signal can be directly spliced, signal jumps easily occur at the splicing site, which can cause severe Gibs. Therefore, it is not recommended to use signal splicing before FFT calculation or Welch power spectrum analysis. The AR model can calculate the power spectrum after signal prediction and modeling, which is very suitable for analyzing shortterm signal power spectra [58], such as the power spectrum of a small segment of rapidly changing EEG signals in eventrelated EEG. However, it should be noted that the AR model method is based on the model, so the selected model may not fit the signal well. Selection of the wrong model will result in large deviations [197]. The application of AR models requires a more accurate evaluation of the signals [199]. Additionally, for signals with rapid changes in signal amplitude, commonly used models cannot fit effectively, so AR models cannot effectively analyze the power spectrum. In summary, if the data obtained are of good quality, such as highquality sleep data, FFT is recommended to obtain more accurate results. If the data contain random noise, such as in the longterm monitoring of epilepsy, Welch or multitaper is recommended. If the data are short, such as EEG after physical instant stimulation, the AR model is recommended for fitting. Power spectrum analysis is the basis of EEG signal analysis and one of the most important analysis methods. With upgrades in EEG acquisition equipment, the signaltonoise ratio of the obtained EEG signal is also increasing. At this time, the transient power spectrum will play an important role, and related research will enhance the development of braincomputer interfaces (BCIs) and machine learning fields.
Time–frequency analysis is a powerful tool for analyzing eventrelated EEG signals because it can describe the changes in timevarying EEG from two dimensions based on the changes in time and frequency. Among the various time–frequency analysis methods, the continuous WT is the most commonly used method due to its good performance in balancing time and frequency. However, its high computational and spatial complexity makes it unsuitable for longterm data analysis [200]. Correspondingly, DWT shows a good performance in signal decomposition but a poor visualization effect and unsatisfactory frequency resolution [201]. For timeinvariant EEG, although the changes in EEG are insignificant, the energy of diverse frequency bands will certainly change in response to longterm changes. STFT can present this response well, having good frequency resolution [3] and low time resolution [69]. If WT is used, although it is more informative, the resultant massive data redundancy is not conducive to observing the main features. WVD has a very high time–frequency resolution for shortterm signals but is highly susceptible to noise. Therefore, WVD is more suitable for analyzing shortterm signals with less noise [70]. EMD is very suitable for analyzing signals with many abrupt amplitudes [72] and will not be affected by Gibs, but WVD has poor frequency localization of lownoise signals [201]. With advances in computer performance, time–frequency analysis methods have gradually replaced power spectrum methods as the first choice for observing spectrum changes. The integration of time–frequency graphs and imagebased deep learning has also produced many highquality applications. However, it should be noted that time–frequency analysis is still not detailed enough for the calculation of instantaneous changes. For example, it is difficult to use time–frequency diagram analysis for the EEG signal during epileptic discharge. In this case, it should be combined with the time domain method to improve the description accuracy.
Connectivity analysis is an important part of neural signal analysis. The selection of relevant analysis methods is crucial to the main features that need to be analyzed. This article introduced several classical correlation analysis methods from the perspectives of time, frequency, and nonlinearity. However, many other methods have not been introduced, such as event statistics, CFC coupling, and amplitude frequency coupling. These methods are similar to the ones in this article and can be selected based on the data characteristics. Several correlation analysis methods applicable to longterm signals can be used to observe the correlations between different brain regions in timeinvariant EEG due to insignificant changes in the characteristics of timeinvariant EEG. CORR focuses more on the timescale similarity of signals, while COH focuses more on their frequencyscale similarity [107]. The PLV/PLI method refers to phasebased statistics, with PLV being stricter. Such a method can be easily converted into correlation statistics of crossfrequency or other events [202]. These three methods can be used to measure the correlation between channels in terms of signal similarity and signal phase synchronization [112]. MI and GC can analyze the driving force between signals from the perspective of information transmission [109]. These two methods can be used to study the relationship between signals. MI observes the relationship from the perspective of information transmission, while GC infers the relationship from the perspective of regression. Both methods are nonlinear but differ in their approach. MI uses information entropy, while GC uses regression. As a result, the characteristics of the observations differ between the two methods [111]. However, WTC requires clear events to produce more accurate results, making it suitable for accurate eventrelated EEG but not for timeinvariant EEG. All the methods mentioned above are also suitable for accurate eventrelated EEG. However, it is important to note that the relevant changes presented by eventrelated EEG during the event are accurate and real. Therefore, correlation analysis of shortterm signals after the event is recommended to obtain related results with a higher signaltonoise ratio. For random eventrelated EEG, due to the inaccuracy of its events, WTC methods are also unable to provide accurate results. However, due to the timevarying characteristics, other correlation analysis methods are required to segment and classify the signals according to the main analysis features to improve the signaltonoise ratio of the results [203]. Connectivity analysis methods have emerged as a powerful tool for studying brain networks, which are important components of brain cognition. Moreover, these methods can be combined with neural networks to develop new bionic operations. Different connection patterns can represent different ways of connecting neurons. The existing neural networks usually use direct signal connections, but other connection methods can be used to produce more intelligent network interfaces. Spiking neural networks are a representative example of this type of research.
Source localization analysis can convert multichannel EEG from the scalp to deeper brain regions, which can more clearly localize the position of signal generation. A personalized brain model is generated through structured MRI, following which intracranial nerve activity is inferred from extracranial nerve activity signals to estimate the discharge location of the personalized brain. This process can locate key intracranial locations without surgery and has been widely used in the fields of epilepsy focus location and functional area location. Existing source localization methods have different implementation principles. The minimum norm estimation can partly localize the source into the intracranial region but has low accuracy in deep source localization and a high requirement for signal quality due to the use of MEG [139]. The dipole and beamforming methods can effectively localize the source into the intracranial region and have good localization accuracy, but they require more accurate parameters and an accurate head model [142, 163]. Beamforming balances speed and accuracy and is a rapidly developing method [204]. The FOCUSS method has good resolution but low computational efficiency [140]. LORETA is currently the most frequently used source localization method [205]. The spatial resolution of this method is not high, but its temporal resolution is good; thus, it has received considerable attention in EEG analysis [139]. CSD does not require a brain model and its calculation process is simple and fast, but it cannot localize the source in the intracranial cavity [144]. Signal quality plays a crucial role in the accuracy of multichannel brain source localization methods. Therefore, the source localization method best suits shortterm EEG localization analysis in eventrelated EEG. By leveraging the high signaltonoise ratio features at the onset time of events, more accurate localization results can be obtained. For the source localization analysis of timeinvariant EEG, the discharge characteristics of the EEG should be converted before localization, or the superposition method should be used to increase the signaltonoise ratio. For random eventrelated EEG, it is advisable to choose EEG signals with a high signaltonoise ratio for localization as far as possible. Source location methods have significant restrictions, including signal quality, number of channels, and brain model accuracy, which can affect positioning accuracy. However, targeted selection can improve the accuracy of these methods. With the continuous development of computing power and artificial intelligence technology, the quality requirements of source localization methods will gradually decrease, allowing them to be widely used in clinical practice.
Machine learning methods are rapidly developing methods that learn features, which can be the original EEG signals or converted features of EEG [148], such as power and connectivity [149]. Hence, although machine learning has a wide range of adaptability due to its selflearning nature, the generation and selection of features remain points of discussion. For timeinvariant EEG, the features of EEG can be the average of longterm features, such as the power spectrum and connectivity of eachchannel EEG. For accurate eventrelated EEG, eventrelated features can be learned. For random eventrelated EEG, machine learning methods with clustering properties are more suitable for semisupervised or unsupervised feature learning [150]. Deep learning has been successfully applied in multiple EEG signal tasks, such as motor imagery, epilepsy detection, severe depression detection, sleep stage scoring, and eventrelated potential tasks [206]. There are differences between the data for different tasks, such as signal window length and channel count. Given these differences, selecting the suitable type of deep learning network can achieve better classification performance. Using CNN to classify spectrograms can also produce good results [207], but CNN models are suitable for data classification without time information, while LSTM models are suitable for regression analysis with time information [190]. However, compared with traditional single models such as CNN and LSTM, a mixed model is recommended [208]. If a model fails to achieve the expected results, researchers can opt for a fusion of multiple models to improve accuracy. For skilled machine learning researchers, this is a simple and fast way to build applications. The mixed model is expected to perform well in classification accuracy and transfer learning. Nevertheless, deep learning models also have shortcomings. For example, due to the small size of disease datasets, the model is prone to overfitting. Therefore, methods such as introducing regularization terms into the model should be considered to minimize the impact of overfitting [209]. Moreover, the introduction of model interpretability can aid in understanding the feature selection for model classification. There may be causal relationships within the EEG features, and the introduction of causal algorithms can be considered to further optimize the models [210]. At present, most deep learning models are designed based on images and cannot adapt well to EEG signal data. The transformer, as a new type of neural network model, is being used in the diagnosis and treatment of brain diseases, but its application in EEG needs to be studied further [211]. Therefore, it is necessary to consider encoding and decoding the data based on EEG signal data features and developing new model structures. Machine learning methods are currently the fastestgrowing neural signal processing methods, and many researchers have proposed new processing ideas in EEG analysis. With the rapid development of brainlike intelligence, a largescale model may emerge to perform bionic simulations of human brain functions. This kind of research will help develop the field of artificial intelligence to a higher level.
This article discussed some commonly used EEG analysis methods. However, in practice, a combination of methods is frequently used. For instance, BCI, a frontier field of neuroscience and neurological diseases, requires the complex processing of brain electrical signals using multiple methods. A neural device known as a BCI translates the neural activity of an individual into external responses or directives. These interfaces find applications in restoring functionality for conditions such as epilepsy, stroke, spinal cord injuries, ALS, cerebral palsy, narcolepsy, Parkinson’s disease, and neuromuscular disorders [212]. In the realm of mental health, BCIs are under investigation as potential treatments for conditions like depression, anxiety, obsessive–compulsive disorder, and other neuropsychiatric disorders [213]. BCIs are versatile in acquiring a diverse array of signals, each associated with different objects. Despite this variability, the decoding of brain signals generally follows a fivestage process: signal acquisition, preprocessing, feature extraction, classification, and control interface. These stages involve the integration of various methodologies. Ongoing research in BCI analysis methods aims to enhance accuracy and reliability. Notably, the application of deep learning algorithms for EEG data analysis is a promising avenue. Another focus is on leveraging explainable artificial intelligence techniques to gain insights into BCI analysis outcomes. Like BCI, numerous studies necessitate the integration of traditional and innovative technologies to continually enhance the efficacy of EEG analysis methods and establish a foundation for further research.
Conclusions
EEG/iEEG is commonly applied in functional neuroimaging and is one of the leading tools in neuroscience. Clinical medicine, BCIs, and psychological research all require EEG/iEEG analysis. In recent decades, a variety of analysis methods have emerged for researchers to choose from, and interest in such techniques is high. However, the abundance of analysis methods has led researchers to question their applicability.
This review categorizes representative research methods based on the characteristics of EEG/iEEG signals. The methods are classified into power spectrum analysis, time–frequency analysis, connectivity analysis, source localization analysis, and machine learning. Other methods with wide application scenarios, such as nonlinear analysis, predictive analysis, and graph theory analysis, are not introduced in this review. These methods are considered to have certain similarities to or be the extension of the classical methods in this review from the perspective of analysis purpose. The methods introduced in this review are only a subset of common methods, and users need to make choices based on the characteristics of the data and methods.
Availability of data and materials
Not applicable.
Abbreviations
 AR:

Autoregressive
 BCI:

Braincomputer interface
 CFA:

Crossfrequency analysis
 CFC:

Crossfrequency coupling
 CFD:

Crossfrequency directionality
 CNN:

Convolutional neural network
 COH:

Coherence
 CORR:

Correlation
 CSD:

Current source density
 CSP:

Common spatial pattern
 EEG:

Electroencephalography
 EMD:

Empirical mode decomposition
 FFT:

Fast Fourier transform
 FOCUSS:

Focal underdetermined system solution
 GC:

Granger causality
 iEEG:

Intracranial electroencephalography
 IMF:

Intrinsic mode function
 LDA:

Linear discriminant analysis
 LORETA:

Lowresolution electromagnetic tomography
 LSTM:

Long shortterm memory
 MEG:

Magnetoencephalography
 MI:

Mutual information
 MRI:

Magnetic resonance imaging
 PLV/PLI:

Phase lag value/index
 RNN:

Recurrent neural network
 sLORETA:

Standardized lowresolution electromagnetic tomography
 STFT:

Shorttime Fourier transform
 SVM:

Support vector machine
 WT:

Wavelet transform
 WTC:

Wavelet coherence
 WVD:

WignerVille distribution
References
Boone W, Piccinini G. The cognitive neuroscience revolution. Synthese. 2016;193(5):1509–34.
Mukamel EA, Ngai J. Perspectives on defining cell types in the brain. Curr Opin Neurobiol. 2019;56:61–8.
Bullmore ET, Sporns O. Complex brain networks: graph theoretical analysis of structural and functional systems. Nat Rev Neurosci. 2009;10(3):186–98.
Fingelkurts AA, Fingelkurts AA. Timing in cognition and EEG brain dynamics: discreteness versus continuity. Cogn Process. 2006;7(3):135–62.
Michel CM, Murray MM. Towards the utilization of EEG as a brain imaging tool. Neuroimage. 2012;61(2):371–85.
Caria A, Sitaram R, Birbaumer N. Realtime fMRI: a tool for local brain regulation. Neuroscientist. 2012;18(5):487–501.
Brookes MJ, Leggett J, Rea M, Hill RM, Holmes N, Boto E, et al. Magnetoencephalography with optically pumped magnetometers (OPMMEG): the next generation of functional neuroimaging. Trends Neurosci. 2022;45(8):621–34.
Bowman FD, Guo Y, Derado G. Statistical approaches to functional neuroimaging data. Neuroimag Clin N Am. 2007;17(4):441.
Olejniczak P. Neurophysiologic basis of EEG. J Clin Neurophysiol. 2006;23(3):186–9.
Sauseng P, Klimesch W, Stadler W, Schabus M, Doppelmayr M, Hanslmayr S, et al. A shift of visual spatial attention is selectively associated with human EEG alpha activity. Eur J Neurosci. 2005;22(11):2917–26.
van Driel J, Gunseli E, Meeter M, Olivers CNL. Local and interregional alpha EEG dynamics dissociate between memory for search and memory for recognition. Neuroimage. 2017;149:114–28.
Alday PM. M/EEG analysis of naturalistic stories: a review from speech to language processing. Lang Cogn Neurosci. 2019;34(4):457–73.
Hu W, Huang G, Li L, Li Z, Zhang Z, Liang Z. Videotriggered EEGemotion public databases and current methods: a survey. Brain Sci Adv. 2020;6(3):255–87.
Ponciano V, Pires IM, Ribeiro FR, Villasana MV, Teixeira MC, Zdravevski E. Experimental study for determining the parameters required for detecting ECG and EEG related diseases during the timedup and go test. Computers. 2020;9(3):67.
Miskovic V, MacDonald KJ, Rhodes LJ, Cote KA. Changes in EEG multiscale entropy and powerlaw frequency scaling during the human sleep cycle. Hum Brain Mapp. 2019;40(2):538–51.
Tsipouras M. Spectral information of EEG signals with respect to epilepsy classification. EURASIP J Adv Signal Process. 2019;2019:10.
Kappenman ES, Farrens JL, Zhang W, Stewart AX, Luck SJ. ERP CORE: an open resource for human eventrelated potential research. Neuroimage. 2021;225:117465.
Mahato S, Paul S. Classification of depression patients and normal subjects based on electroencephalogram (EEG) signal using alpha power and theta asymmetry. J Med Syst. 2020;44(1):28.
Suárez AAP, Batista SB, Ibáñez IP, Fernández EC, Campos MF, Chacón LM. EEGderived functional connectivity patterns associated with mild cognitive impairment in Parkinson’s disease. Behav Sci. 2021;11(3):40.
Ogilvie RD, Simons IA, Kuderian RH, Macdonald T, Rustenburg J. Behavioral, eventrelated potential, and EEG/FFT changes at sleep onset. Psychophysiology. 1991;28(1):54–64.
Hadjiyannakis K, Ogilvie RD, Alloway CED, Shapiro C. FFT analysis of EEG during stage 2toREM transitions in narcoleptic patients and normal sleepers. Electroen Clin Neuro. 1997;103(5):543–53.
Sun J, Cao R, Zhou MN, Hussain W, Wang B, Xue JY, et al. A hybrid deep neural network for classification of schizophrenia using EEG data. Sci Rep. 2021;11(1):4706.
Behnam H, Sheikhani A, Mohammadi MR, Noroozian M, Golabi P. Analyses of EEG background activity in autism disorders with fast Fourier transform and short time Fourier measure. In: ICIAS 2007: International Conference on Intelligent & Advanced Systems. Kuala Lumpur, Malaysia; 2007. p. 12404.
Djamal EC, Furi WI, Nugraha F. Detection of EEG signal poststroke using FFT and convolutional neural network. In: 6th International Conference on Electrical Engineering, Computer Science and Informatics (EECSI). Bandung, Indonesia; 2019. p. 1823.
Farihah SN, Lee KY, Mansor W, Mohamad NB, Mahmoodin Z, Saidi SA. EEG average FFT index for dyslexic children with writing disorder. In: IEEE Student Symposium in Biomedical Engineering & Sciences (ISSBES). Shah Alam, Malaysia; 2015. p. 118–21
Melinda M, Enriko IK, Furqan M, Irhamsyah M, Yunidar Y, Basir N. The effect of power spectral density on the electroencephalography of autistic children based on the welch periodogram method. J Infotel. 2023;15:111–20.
Bian ZJ, Li QL, Wang L, Lu CB, Yin SM, Li XL. Relative power and coherence of EEG series are related to amnestic mild cognitive impairment in diabetes. Front Aging Neurosci. 2014;6:11.
Yuan Y, Lu CB, Li XL. Effect of focused ultrasound stimulation at different ultrasonic power levels on the local field potential power spectrum. Chinese Phys B. 2015;24(8):088704.
Wang Y, Sokhadze EM, ElBaz AS, Li XL, Sears L, Casanova MF, et al. Relative power of specific eeg bands and their ratios during neurofeedback training in children with autism spectrum disorder. Front Hum Neurosci. 2016;9:723.
Göker H. Automatic detection of migraine disease from EEG signals using bidirectional longshort term memory deep learning model. Signal Image Video P. 2023;17(4):1255–63.
Hu M, Li JJ, Li G, Tang XW, Ding QP. Classification of normal and hypoxia EEG based on approximate entropy and Welch powerspectraldensity. The 2006 IEEE International Joint Conference on Neural Network Proceedings. Vancouver, BC; 2006. p. 321822.
Wijaya SK, Badri C, Misbach J, Soemardi TP, Sutarmo V. Electroencephalography (EEG) for detecting acute ischemic stroke. In: 4th International Conference on Instrumentation, Communications, Information Technology, and Biomedical Engineering (ICICIBME). Bandung, Indonesia; 2015. p. 42–8.
Cornelissen L, Kim SE, Purdon PL, Brown EN, Berde CB. Agedependent electroencephalogram (EEG) patterns during sevoflurane general anesthesia in infants. Elife. 2015;4:e06513.
Yang C, Han X, Wang YJ, Gao XR. A frequency recognition method based on Multitaper Spectral Analysis and SNR estimation for SSVEPbased braincomputer interface. IEEE Eng Med Biol. 2017:1930–3.
Oliva JT, Rosa JLG. Binary and multiclass classifiers based on multitaper spectral features for epilepsy detection. Biomed Signal Proces. 2021;66:102469.
Oliveira GHBS, Coutinho LR, da Silva JC, Pinto IJP, Ferreira JMS, Silva FJS, et al. Multitaperbased method for automatic kcomplex detection in human sleep EEG. Expert Syst Appl. 2020;151:113331.
Moharnmadi G, Shoushtari P, Ardekani BM, Shamsollahi MB. Person identification by using AR model for EEG signals. Proc Wrld Acad Sci E. 2006;11:281–5.
Perumalsamy V, Sankaranarayanan S, Rajamony S. Sleep spindles detection from human sleep EEG signals using autoregressive (AR) model: a surrogate data approach. J Biomed Sci Eng. 2009;02:294–303.
Saidatul A, Paulraj MP, Yaacob SB, Yusnita MA. Analysis of EEG signals during relaxation and mental stress condition using AR modeling techniques. In: IEEE International Conference on Control System, Computing and Engineering. Penang, Malaysia; 2011. p. 477–81.
Lawhern V, Hairston WD, McDowell K, Westerfield M, Robbins K. Detection and classification of subjectgenerated artifacts in EEG signals using autoregressive models. J Neurosci Meth. 2012;208(2):181–9.
Mousavi SR, Niknazar M, Vahdat BV. Epileptic seizure detection using AR model on EEG signals. In: Cairo International Biomedical Engineering Conference. Cairo, Egypt; 2008. p. 14,
Akin M. Comparison of wavelet transform and FFT methods in the analysis of EEG signals. J Med Syst. 2002;26(3):241–7.
Evans DH. Doppler signal analysis. Ultrasound Med Biol. 2000;26:S13–5.
SpyersAshby JM, Bain PG, Roberts SJ. A comparison of fast Fourier transform (FFT) and autoregressive (AR) spectral estimation techniques for the analysis of tremor data. J Neurosci Meth. 1998;83(1):35–43.
Günes S, Polat K, Yosunkaya S. Efficient sleep stage recognition system based on EEG signal using kmeans clustering based feature weighting. Expert Syst Appl. 2010;37(12):7922–8.
Prerau MJ, Brown RE, Bianchi MT, Ellenbogen JM, Purdon PL. Sleep neurophysiological dynamics through the lens of multitaper spectral analysis. Physiology. 2017;32(1):60–92.
Oostenveld R, Fries P, Maris E, Schoffelen JM. FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput Intel Neurosc. 2011;2011:156869.
Zhang Y, Liu B, Ji XM, Huang D. Classification of EEG signals based on autoregressive model and wavelet packet decomposition. Neural Process Lett. 2017;45(2):365–78.
Samiee K, Kovács P, Gabbouj M. Epileptic seizure classification of EEG timeseries using rational discrete shorttime Fourier transform. IEEE Trans BioMed Eng. 2015;62(2):541–52.
Beeraka SM, Kumar A, Sameer M, Ghosh S, Gupta B. Accuracy enhancement of epileptic seizure detection: a deep learning approach with hardware realization of SIFT. Circuits Syst Signal Process. 2022;41(1):461–84.
Bajaj V, Guo YH, Sengur A, Siuly S, Alcin OF. A hybrid method based on timefrequency images for classification of alcohol and control EEG signals. Neural Comput Appl. 2017;28(12):3717–23.
Sheikhani A, Behnam H, Mohammadi MR, Noroozian M, Golabi P. Connectivity analysis of quantitative electroencephalogram background activity in autism disorders with short time Fourier transform and coherence values. In: Congress on Image and Signal Processing. Sanya, China; 2008. p. 207–12.
Krishnan P, Yaacob S, Krishnan AP, Rizon M, Ang CK. EEG based drowsiness detection using relative band power and shorttime Fourier transform. J Robot Netw Artif L. 2020;7(3):147–51.
Bajaj V, Pachori RB. Automatic classification of sleep stages based on the timefrequency image of EEG signals. Comput Meth Prog Bio. 2013;112(3):320–8.
Yan AY, Zhou WD, Yuan Q, Yuan SS, Wu Q, Zhao XH, et al. Automatic seizure detection using Stockwell transform and boosting algorithm for longterm EEG. Epilepsy Behav. 2015;45:8–14.
Ebrahimzadeh E, Alavi SM, Bijar A, Pakkhesal A. A novel approach for detection of deception using smoothed Pseudo Wigner–Ville distribution (SPWVD). J Biomed Sci Eng. 2013;2013:8–18.
Khare SK, Bajaj V, Acharya UR. SPWVDCNN for automated detection of schizophrenia patients using EEG signals. IEEE T Instrum Meas. 2021;70:1–9.
Faust O, Acharya UR, Adeli H, Adeli A. Waveletbased EEG processing for computeraided seizure detection and epilepsy diagnosis. SeizureEur J Epilep. 2015;26:56–64.
Gandhi T, Panigrahi BK, Bhatia M, Anand S. Expert model for detection of epileptic activity in EEG signature. Expert Syst Appl. 2010;37(4):3513–20.
Anuragi A, Sisodia DS. Alcohol use disorder detection using EEG signal features and flexible analytical wavelet transform. Biomed Signal Proces. 2019;52:384–93.
Murugappan M, Ramachandran N, Sazali Y. Classification of human emotion from EEG using discrete wavelet transform. J Biomed Sci Eng. 2010;3:390–6.
Adeli H, Zhou Z, Dadmehr N. Analysis of EEG records in an epileptic patient using wavelet transform. J Neurosci Meth. 2003;123(1):69–87.
Hamad A, Houssein EH, Hassanien AE, Fahmy AA. Feature extraction of epilepsy EEG using discrete wavelet transform. In: 12th International Computer Engineering Conference (ICENCO). Cairo, Egypt; 2016. p. 1905.
Li SF, Zhou WD, Yuan Q, Geng SJ, Cai DM. Feature extraction and recognition of ictal EEG using EMD and SVM. Comput Biol Med. 2013;43(7):807–16.
Siuly S, Khare SK, Bajaj V, Wang H, Zhang YC. A computerized method for automatic detection of schizophrenia using EEG signals. IEEE T Neur Sys Reh. 2020;28(11):2390–400.
Chen SJ, Peng CJ, Chen YC, Hwang YR, Lai YS, Fan SZ, et al. Comparison of FFT and marginal spectra of EEG using empirical mode decomposition to monitor anesthesia. Comput Meth Prog Biol. 2016;137:77–85.
Babiker A, Faye I, Mumtaz W, Malik AS, Sato H. EEG in classroom: EMD features to detect situational interest of students during learning. Multimed Tools Appl. 2019;78(12):16261–81.
Priya A, Yadav P, Jain S, Bajaj V. Efficient method for classification of alcoholic and normal EEG signals using EMD. J Eng. 2018;3:166–72.
Kiymik MK, Güler I, Dizibüyük A, Akin M. Comparison of STFT and wavelet transform methods in determining epileptic seizure activity in EEG signals for realtime application. Comput Biol Med. 2005;35(7):603–16.
Mert A, Akan A. Emotion recognition based on timefrequency distribution of EEG signals using multivariate synchrosqueezing transform. Digit Signal Process. 2018;81:106–15.
Riaz F, Hassan A, Rehman S, Niazi IK, Dremstrup K. EMDbased temporal and spectral features for the classification of EEG signals using supervised learning. IEEE Trans Neural Syst Rehabilit. 2016;24(1):28–35.
Zhuang N, Zeng Y, Tong L, Zhang C, Zhang HM, Yan B. Emotion recognition from EEG signals using multidimensional information in EMD domain. Biomed Res Int. 2017;2017:8317357.
Kang JM, Cho SE, Na KS, Kang SG. Spectral power analysis of sleep electroencephalography in subjects with different severities of obstructive sleep apnea and healthy controls. Nat Sci Sleep. 2021;13:477–86.
Dauwels J, Vialatte F, Musha T, Cichocki A. A comparative study of synchrony measures for the early diagnosis of Alzheimer’s disease based on EEG. Neuroimage. 2010;49(1):668–93.
Yasuhara A. Correlation between EEG abnormalities and symptoms of autism spectrum disorder (ASD). Brain DevJpn. 2010;32(10):791–8.
Islam MR, Islam MM, Rahman MM, Mondal C, Singha SK, Ahmad M, et al. EEG channel correlation based model for emotion recognition. Comput Biol Med. 2021;136:104757.
Sheorajpanday RVA, Nagels G, Weeren AJTM, van Putten MJAM, De Deyn PP. Quantitative EEG in ischemic stroke: correlation with functional status after 6 months. Clin Neurophysiol. 2011;122(5):874–83.
Alba G, Pereda E, Mañas S, Méndez LD, Duque MAR, González A, et al. The variability of EEG functional connectivity of young ADHD subjects in different resting states. Clin Neurophysiol. 2016;127(2):1321–30.
CarrascoGómez M, Keijzer HM, Ruijter BJ, Bruña R, TjepkemaCloostermans MC, Hofmeijer J, et al. EEG functional connectivity contributes to outcome prediction of postanoxic coma. Clin Neurophysiol. 2021;132(6):1312–20.
Barry RJ, Clarke AR, McCarthy R, Selikowitz M, Johnstone SJ, Rushby JA. Age and gender effects in EEG coherence: I. Development trends in normal children. Clin Neurophysiol. 2004;115(10):2252–8.
Locatelli T, Cursi M, Liberati D, Franceschi M, Comi G. EEG coherence in Alzheimer’s disease. Electroen Clin Neuro. 1998;106(3):229–37.
Coben R, Clarke AR, Hudspeth W, Barry RJ. EEG power and coherence in autistic spectrum disorder. Clin Neurophysiol. 2008;119(5):1002–9.
Catarino A, Andrade A, Churches O, Wagner AP, BaronCohen S, Ring H. Taskrelated functional connectivity in autism spectrum conditions: an EEG study using wavelet transform coherence. Mol Autism. 2013;4(1):1.
Omidvarnia A, Pedersen M, Vaughan DN, Walz JM, Abbott DF, Zalesky A, et al. Dynamic coupling between fMRI local connectivity and interictal EEG in focal epilepsy: a wavelet analysis approach. Hum Brain Mapp. 2017;38(11):5356–74.
Khan DM, Masroor K, Jailani MFM, Yahya N, Yusoff MZ, Khan SM. Development of wavelet coherence EEG as a biomarker for diagnosis of major depressive disorder. IEEE Sens J. 2022;22(5):4315–25.
Sankari Z, Adeli H, Adeli A. Wavelet coherence model for diagnosis of Alzheimer disease. Clin Eeg Neurosci. 2012;43(4):268–78.
Briels CT, Schoonhoven DN, Stam CJ, de Waal H, Scheltens P, Gouw AA. Reproducibility of EEG functional connectivity in Alzheimer’s disease. Alzheimers Res Ther. 2020;12(1):68.
Olejarczyk E, Jernajczyk W. Graphbased analysis of brain connectivity in schizophrenia. PLoS One. 2017;12(11):e0188629.
Wang ZM, Chen ZY, Zhang J. EEG emotion recognition based on PLVrichclub dynamic brain function network. Appl Intell. 2023;53(14):17327–45.
Huang H, Zhang JH, Zhu L, Tang JJ, Lin G, Kong WZ, et al. EEGbased sleep staging analysis with functional connectivity. Sensors. 2021;21(6):1988.
Zuchowicz U, WozniakKwasniewska A, Szekely D, Olejarczyk E, David O. EEG phase synchronization in persons with depression subjected to transcranial magnetic stimulation. Front Neurosci. 2019;12:1037.
Chen H, Song Y, Li XL. A deep learning framework for identifying children with ADHD using an EEGbased brain network. Neurocomputing. 2019;356:83–96.
Aydin S, Tunga MA, Yetkin S. Mutual information analysis of sleep EEG in detecting psychophysiological insomnia. J Med Syst. 2015;39(5):43.
Piho L, Tjahjadi T. A mutual information based adaptive windowing of informative EEG for emotion recognition. IEEE T Affect Comput. 2020;11(4):722–35.
Hassan KM, Islam MR, Nguyen TT, Molla MKI. Epileptic seizure detection in EEG using mutual informationbased best individual feature selection. Expert Syst Appl. 2022;193:116414.
Yin ZL, Li J, Zhang Y, Ren AF, Von Meneen KM, Huang LY. Functional brain network analysis of schizophrenic patients with positive and negative syndrome based on mutual information of EEG time series. Biomed Signal Proces. 2017;31:331–8.
SanzGarcía A, PérezRomero M, Pastor J, Sola RG, VegaZelaya L, Monasterio F, et al. Identifying causal relationships between EEG activity and intracranial pressure changes in neurocritical care patients. J Neural Eng. 2018;15(6):066029.
de Tommaso M, Trotta G, Vecchio E, Ricci K, Van de Steen F, Montemurno A, et al. Functional connectivity of EEG signals under laser stimulation in migraine. Front Hum Neurosci. 2015;9:640.
Nicolaou N, Hourris S, Alexandrou P, Georgiou J. EEGbased automatic classification of “awake” versus “anesthetized” state in general anesthesia using granger causality. PLoS ONE. 2012;7(3):e33869.
Nicolaou N, Georgiou J. Neural networkbased classification of anesthesia/awareness using granger causality features. Clin Eeg Neurosci. 2014;45(2):77–88.
Barrett AB, Murphy M, Bruno MA, Noirhomme Q, Boly M, Laureys S, et al. Granger causality analysis of steadystate electroencephalographic signals during propofolinduced anaesthesia. PLoS ONE. 2012;7(1):e29072.
Coben R, MohammadRezazadeh I. Neural connectivity in epilepsy as measured by granger causality. Front Hum Neurosci. 2015;9:194.
Chen X, Li YJ, Li RR, Yuan X, Liu M, Zhang W, et al. Multiple crossfrequency coupling analysis of restingstate EEG in patients with mild cognitive impairment and Alzheimer’s disease. Front Aging Neurosci. 2023;15:1142085.
Lynn PA, Sponheim SR. Disturbed theta and gamma coupling as a potential mechanism for visuospatial working memory dysfunction in people with schizophrenia. Neuropsychiatr Electrophysiol. 2016;2:1–30.
Papadaniil CD, Kosmidou VE, Tsolaki A, Tsolaki M, Kompatsiaris I, Hadjileontiadis LJ, et al., editors. Phaseamplitude crossfrequency coupling in EEGderived cortical time series upon an auditory perception task. In: 37th Annual International Conference of the IEEEEngineeringinMedicineandBiologySociety (EMBC). Milan, Italy; 2015. p. 41503.
Park H, Lee DS, Kang E, Kang H, Hahm J, Kim JS, et al. Formation of visual memories controlled by gamma power phaselocked to alpha oscillations. Sci Rep. 2016;6:28092.
Li XW, La R, Wang Y, Hu B, Zhang XM. A deep learning approach for mild depression recognition based on functional connectivity using electroencephalography. Front Neurosci. 2020;14:192.
Thatcher RW, North D, Biver C. EEG and intelligence: relations between EEG coherence, EEG phase delay and power. Clin Neurophysiol. 2005;116(9):2129–41.
Tafreshi TF, Daliri MR, Ghodousi M. Functional and effective connectivity based features of EEG signals for object recognition. Cogn Neurodyn. 2019;13(6):555–66.
van Mierlo P, Papadopoulou M, Carrette E, Boon P, Vandenberghe S, Vonck K, et al. Functional brain connectivity from EEG in epilepsy: seizure prediction and epileptogenic focus localization. Prog Neurobiol. 2014;121:19–35.
Seth AK, Barrett AB, Barnett L. Granger causality analysis in neuroscience and neuroimaging. J Neurosci. 2015;35(8):3293–7.
Tort ABL, Komorowski R, Eichenbaum H, Kopell N. Measuring phaseamplitude coupling between neuronal oscillations of different frequencies. J Neurophysiol. 2010;104(2):1195–210.
Heusser AC, Poeppel D, Ezzyat Y, Davachi L. Episodic sequence memory is supported by a thetagamma phase code. Nat Neurosci. 2016;19(10):1374–80.
Lisman JE, Jensen O. The thetagamma neural code. Neuron. 2013;77(6):1002–16.
Roux F, Uhlhaas PJ. Working memory and neural oscillations: Alphagamma versus thetagamma codes for distinct WM information? Trends Cogn Sci. 2014;18(1):16–25.
Abubaker M, Al Qasem W, Kvasnák E. Working memory and crossfrequency coupling of neuronal oscillations. Front Psychol. 2021;12:756661.
LozanoSoldevilla D, ter Huurne N, Oostenveld R. Neuronal oscillations with nonsinusoidal morphology produce spurious phasetoamplitude coupling and directionality. Front Comput Neurosci. 2016;10:87.
la Tour TD, Tallot L, Grabot L, Doyère V, van Wassenhove V, Grenier Y, et al. Nonlinear autoregressive models for crossfrequency coupling in neural time series. PLoS Comput Biol. 2017;13(12):e1005893.
He B, Sohrabpour A, Brown E, Liu ZM. Electrophysiological source imaging: a noninvasive window to brain dynamics. Annu Rev Biomed Eng. 2018;20:171–96.
Lei X, Wu TY, ValdesSosa PA. Incorporating priors for EEG source imaging and connectivity analysis. Front Neurosci. 2015;9:284.
Hinkley LBN, Dale CL, Cai C, Zumer J, Dalal S, Findlay A, et al. NUTMEG: open source software for M/EEG source reconstruction. Front Neurosci. 2020;14:710.
Rosenow F, Lüders H. Presurgical evaluation of epilepsy. Brain. 2001;124:1683–700.
Symms M, Jäger HR, Schmierer K, Yousry TA. A review of structural magnetic resonance neuroimaging. J Neurol Neurosurg Psychiatry. 2004;75(9):1235–44.
Toole C, MartinezJuárez IE, Gaitanis JN, Sunderam S, Ding L, DiCecco J, et al. Source localization of highfrequency activity in tripolar electroencephalography of patients with epilepsy. Epilepsy Behav. 2019;101(Pt A):106519.
Galaris E, Gallos I, Myatchin I, Lagae L, Siettos C. Electroencephalography source localization analysis in epileptic children during a visual workingmemory task. Int J Numer Method Biomed Eng. 2020;36(12):e3404.
Lee AKC, Larson E, Maddox RK. Mapping cortical dynamics using simultaneous MEG/EEG and anatomicallyconstrained minimumnorm estimates: an auditory attention example. J Vis Exp. 2012;68:e4262.
Sperli F, Spinelli L, Seeck M, Kurian M, Michel CM, Lantz G. EEG source imaging in pediatric epilepsy surgery: a new perspective in presurgical workup. Epilepsia. 2006;47(6):981–90.
Plummer C, Wagner M, Fuchs M, Vogrin S, Litewka L, Farish S, et al. Clinical utility of distributed source modelling of interictal scalp EEG in focal epilepsy. Clin Neurophysiol. 2010;121(10):1726–39.
Wu W, Jia WY, Liu HS, Gao XR, Zhang GJ, Wang YP, editors. Localization of epileptic foci from preictal EEG data using standardized shrinking LORETAFOCUSS algorithm. In: 1st international conference on neural interface and control (CNIC). Wuhan, China; 2005. p. 20912.
Ye JC, Tak S, Han Y, Park HW. Projection reconstruction MR imaging using FOCUSS. Magn Reson Med. 2007;57(4):764–75.
Saletu B, Anderer P, SaletuZyhlarz GM. EEG topography and tomography (LORETA) in diagnosis and pharmacotherapy of depression. Clin EEG Neurosci. 2010;41(4):203–10.
Clemens B, Bessenyei M, Fekete I, Puskás S, Kondákor I, Tóth M, et al. Theta EEG source localization using LORETA in partial epilepsy patients with and without medication. Clin Neurophysiol. 2010;121(6):848–58.
Kopřivová J, Congedo M, Horácek J, Prasko J, Raszka M, Brunovsky M, et al. EEG source analysis in obsessivecompulsive disorder. Clin Neurophysiol. 2011;122(9):1735–43.
Shao SY, Shen KQ, Yu K, WilderSmith EPV, Li XP. Frequencydomain EEG source analysis for acute tonic cold pain perception. Clin Neurophysiol. 2012;123(10):2042–9.
Loughrey DG, Jordan C, Ibanez A, Parra MA, Lawlor BA, Reilly RB. Agerelated hearing loss associated with differences in the neural correlates of feature binding in visual working memory. Neurobiol Aging. 2023;132:233–45.
Dubová D, Dvořácková D, Pavlů D, Pánek D. Cerebral projection of mirrored touch via sLORETA imaging. Life. 2023;13(5):1201.
Liu JH, Zhang Q, Liang MJ, Wang YJ, Chen YB, Wang JB, et al. Altered processing of visual stimuli in vestibular migraine patients between attacks: a combined VEP and sLORETA study. Front Hum Neurosci. 2021;15:762970.
Yoshinaga H, Koutroumanidis M, Shirasawa A, Kikumoto K, Ohtsuka Y, Oka E. Dipole analysis in panayiotopoulos syndrome. Brain Dev. 2005;27(1):46–52.
Ebersole JS. Noninvasive localization of the epileptogenic focus by EEG dipole modeling. Acta Neurol Scand. 1994;89:20–8.
Nakajima Y, Homma S, Musha T, Okamoto Y, Ackerman RH, Correia JA, et al. Dipoletracing of abnormal slow brain potentials after cerebral stroke—EEG, PET, MRI correlations. Neurosci Lett. 1990;112(1):59–64.
Verhellen E, Boon P. EEG source localization of the epileptogenic focus in patients with refractory temporal lobe epilepsy, dipole modelling revisited. Acta Neurol Belg. 2007;107(3):71–7.
Ntolkeras G, Tamilia E, AlHilani M, Bolton J, Grant PE, Prabhu SP, et al. Presurgical accuracy of dipole clustering in MRInegative pediatric patients with epilepsy: validation against intracranial EEG and resection. Clin Neurophysiol. 2022;141:126–38.
Knyazev GG, Savostyanov AN, Bocharov AV, Tamozhnikov SS, Saprigyn AE. Taskpositive and tasknegative networks and their relation to depression: EEG beamformer analysis. Behav Brain Res. 2016;306:160–9.
Neugebauer F, Antonakakis M, Unnwongse K, Parpaley Y, Wellmer J, Rampp S, et al. Validating EEG, MEG and combined MEG and EEG Beamforming for an estimation of the epileptogenic zone in focal cortical dysplasia. Brain Sci. 2022;12(1):114.
Ward DM, Jones RD, Bones PJ, Carroll GJ. Enhancement of deep epileptiform activity in the EEG via 3D adaptive spatial filtering. IEEE Trans Biomed Eng. 1999;46(6):707–16.
Kouchaki S, Boostani R, Razavipour F. Fatigue detection using the strength of dominant EEG source: a beamforming approach. Biomed Eng Appl Basis Commun. 2018;30(3):1850023.
Vergallo P, LayEkuakille A, Giannoccaro NI, Caratelli D, Massaro A, Trabacca A, et al, editors. Processing EEG signals through Beamforming techniques for seizure diagnosis. In: 6th International Conference on Sensing Technology (ICST). Kolkata, India; 2012. p. 497501.
Ponomarev VA, Mueller A, Candrian G, GrinYatsenko VA, Kropotov JD. Group independent component analysis (gICA) and current source density (CSD) in the study of EEG in ADHD adults. Clin Neurophysiol. 2014;125(1):83–97.
Stewart JL, Coan JA, Towers DN, Allen JJB. Resting and taskelicited prefrontal EEG alpha asymmetry in depression: support for the capability model. Psychophysiology. 2014;51(5):446–55.
GrinYatsenko VA, Ponomarev VA, Pronina MV, Poliakov YI, Plotnikova IV, Kropotov JD. Local and widely distributed EEG Activity in schizophrenia with prevalence of negative symptoms. Clin EEG Neurosci. 2017;48(5):307–15.
Kamarajan C, Rangaswamy M, Manz N, Chorlian DB, Pandey AK, Roopesh BN, et al. Topography, power, and current source density of theta oscillations during reward processing as markers for alcohol dependence. Hum Brain Mapp. 2012;33(5):1019–39.
Jatoi MA, Kamel N, Malik AS, Faye I, Begum T. A survey of methods used for source localization using EEG signals. Biomed Signal Process Control. 2014;11:42–52.
Asadzadeh S, Rezaii TY, Beheshti S, Delpak A, Meshgini S. A systematic review of EEG source localization techniques and their applications on diagnosis of brain abnormalities. J Neurosci Methods. 2020;339:108740.
Dattola S, Morabito FC, Mammone N, La Foresta F. Findings about LORETA applied to highdensity EEG—a review. Electronics. 2020;9(4):660.
SadatNejad Y, Beheshti S, IEEE, editors. Higher resolution sLORETA (HRsLORETA) in EEG source Imaging. In: 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC). Berlin, Germany; 2019. p. 16903.
Jatoi MA, Kamel N, Malik AS, Faye I. EEG based brain source localization comparison of sLORETA and eLORETA. Australas Phys Eng Sci Med. 2014;37(4):713–21.
Chikara RK, Ko LW. Modulation of the visual to auditory human inhibitory brain network: an EEG dipole source localization study. Brain Sci. 2019;9(9):216.
Dippel G, Chmielewski W, Mückschel M, Beste C. Response modedependent differences in neurofunctional networks during response inhibition: an EEGbeamforming study. Brain Struct Funct. 2016;221(8):4091–101.
Kayser J, Tenke CE. On the benefits of using surface Laplacian (current source density) methodology in electrophysiology. Int J Psychophysiol. 2015;97(3):171–3.
Peng H, Siddiqui MM, Rhind SG, Zhang J, da Luz LT, Beckett A. Artificial intelligence and machine learning for hemorrhagic trauma care. Mil Med Res. 2023;10(1):6.
Craik A, He YT, ContrerasVidal JL. Deep learning for electroencephalogram (EEG) classification tasks: a review. J Neural Eng. 2019;16(3):031001.
Qaraqe M, Ismail M, Serpedin E. Bandsensitive seizure onset detection via CSPenhanced EEG features. Epilepsy Behav. 2015;50:77–87.
Dissanayake T, Fernando T, Denman S, Sridharan S, Fookes C. Deep learning for patientindependent epileptic seizure prediction using scalp EEG signals. IEEE J Biomed Health Inform. 2022;26(2):527–38.
Liu Y, Zhang H, Chen M, Zhang L. A boostingbased spatialspectral model for stroke patients’ EEG analysis in rehabilitation training. IEEE Trans Neural Syst Rehabil Eng. 2016;24(1):169–79.
Alturki FA, Aljalal M, Abdurraqeeb AM, AlSharabi K, AlShamma’a AA. Common spatial pattern technique with EEG signals for diagnosis of autism and epilepsy disorders. IEEE Access. 2021;9:24334–49.
Jamal W, Das S, Oprescu IA, Maharatna K, Apicella F, Sicca F. Classification of autism spectrum disorder using supervised learning of brain connectivity measures extracted from synchrostates. J Neural Eng. 2014;11(4):046019.
Jeong DH, Kim YD, Song IU, Chung YA, Jeong J. Wavelet energy and wavelet coherence as EEG biomarkers for the diagnosis of Parkinson’s diseaserelated dementia and Alzheimer’s disease. Entropy. 2016;18(1):8.
Boostani R, Sadatnezhad K, Sabeti M. An efficient classifier to diagnose of schizophrenia based on the EEG signals. Expert Syst Appl. 2009;36(3):6492–9.
Rajaguru H, Kumar Prabhakar S. Time frequency analysis (dB2 and dB4) for epilepsy classification with LDA classifier. In: 2017 2nd international conference on communication and electronics systems (ICCES). Coimbatore, India; 2017. p. 708–11.
Kang J, Han X, Song J, Niu Z, Li X. The identification of children with autism spectrum disorder by SVM approach on EEG and eyetracking data. Comput Biol Med. 2020;120:103722.
Fu K, Qu J, Chai Y, Dong Y. Classification of seizure based on the timefrequency image of EEG signals using HHT and SVM. Biomed Signal Process Control. 2014;13:15–22.
Shen KQ, Li XP, Ong CJ, Shao SY, WilderSmith EPV. EEGbased mental fatigue measurement using multiclass support vector machines with confidence estimate. Clin Neurophysiol. 2008;119(7):1524–33.
Liu YX, Zhou WD, Yuan Q, Chen SS. Automatic seizure detection using wavelet transform and SVM in longterm intracranial EEG. IEEE Trans Neural Syst Rehabil Eng. 2012;20(6):749–55.
Zhou MN, Tian C, Cao R, Wang B, Niu Y, Hu T, et al. Epileptic seizure detection based on EEG signals and CNN. Front Neuroinform. 2018;12:95.
Hassan F, Hussain SF, Qaisar SM. Epileptic seizure detection using a hybrid 1D CNNmachine learning approach from EEG data. J Healthc Eng. 2022;2022:9579422.
Hassan F, Hussain SF, Qaisar SM. Fusion of multivariate EEG signals for schizophrenia detection using CNN and machine learning techniques. Inf Fus. 2023;92:466–78.
Dong HY, Chen D, Zhang L, Ke HJ, Li XL. Subject sensitive EEG discrimination with fast reconstructable CNN driven by reinforcement learning: a case study of ASD evaluation. Neurocomputing. 2021;449:136–45.
Aliyu I, Lim YB, Lim CG, Acm, editors. Epilepsy detection in EEG signal using recurrent neural network. In: 3rd International Conference on Intelligent Systems, Metaheuristics and Swarm Intelligence (ISMSI). Male, Maldives; 2019. p. 503.
Lee S, Hussein R, Ward R, Wang ZJ, McKeown MJ. A convolutionalrecurrent neural network approach to restingstate EEG classification in Parkinson’s disease. J Neurosci Methods. 2021;361:109282.
Sarkar AR, Singh A, Chakraborty R. A deep learningbased comparative study to track mental depression from EEG data. Neurosci Inform. 2022;2:100039.
Mishra S, Birok R, editors. Sleep classification using CNN and RNN on Raw EEG singlechannel. In: International Conference on Computational Performance Evaluation (ComPE). Shillong, India; 2020. p. 2327.
Michielli N, Acharya UR, Molinari F. Cascaded LSTM recurrent neural network for automated sleep stage classification using singlechannel EEG signals. Comput Biol Med. 2019;106:71–81.
Hu XM, Yuan SS, Xu FZ, Leng Y, Yuan KJ, Yuan Q. Scalp EEG classification using deep BiLSTM network for seizure detection. Comput Biol Med. 2020;124:103919.
Koya JR, Rao S, Pothunoori SK, Malyala S, IEEE, editors. EEGbased emotion recognition using LSTMRNN machine learning algorithm. In: International conference on innovation in information communication and technology (ICIICT). Chennai, India; 2019. p. 14.
Lee C, An J. LSTMCNN model of drowsiness detection from multiple consciousness states acquired by EEG. Expert Syst Appl. 2023;213:119032.
Alotaiby TN, Alshebeili SA, Alotaibi FM, Alrshoud SR. Epileptic seizure prediction using CSP and LDA for scalp EEG signals. Comput Intell Neurosci. 2017;2017:1240323.
Kaur T, Gandhi TK. Deep convolutional neural networks with transfer learning for automated brain image classification. Mach Vis Appl. 2020;31(3):20.
Chen PL, Chen D, Zhang L, Tang YB, Li XL. Automated sleep spindle detection with mixed EEG features. Biomed Signal Process Control. 2021;70:103026.
Wen D, Li P, Zhou YH, Sun YB, Xu J, Liu YJ, et al. Feature classification method of restingstate EEG signals from amnestic mild cognitive impairment with type 2 diabetes mellitus based on multiview convolutional neural network. IEEE Trans Neural Syst Rehabil Eng. 2020;28(8):1702–9.
Roy Y, Banville H, Albuquerque I, Gramfort A, Falk TH, Faubert J. Deep learningbased electroencephalography analysis: a systematic review. J Neural Eng. 2019;16(5):051001.
Zeng H, Yang C, Dai GJ, Qin FW, Zhang JH, Kong WZ. EEG classification of driver mental states by deep learning. Cogn Neurodyn. 2018;12(6):597–606.
Zhang YQ, Chen JL, Tan JH, Chen YX, Chen YY, Li DH, et al. An investigation of deep learning models for EEGbased emotion recognition. Front Neurosci. 2020;14:622759.
Shoorangiz R, Weddell SJ, Jones RD. EEGbased machine learning: theory and applications. In: Thakor NV, editors. Handbook of Neuroengineering. Springer, Singapore; 2021.
Si Y. Machine learning applications for electroencephalograph signals in epilepsy: a quick review. Acta Epileptol. 2020;2:1–7.
Gao ZK, Dang WD, Wang XM, Hong XL, Hou LH, Ma K, et al. Complex networks and deep learning for EEG signal analysis. Cogn Neurodyn. 2021;15(3):369–88.
Huang XY, Zhou N, Choi KS. A generalizable and discriminative learning method for deep EEGbased motor imagery classification. Front Neurosci. 2021;15:760979.
Güler I, Kiymik MK, Akin M, Alkan A. AR spectral analysis of EEG signals by using maximum likelihood estimation. Comput Biol Med. 2001;31(6):441–50.
Fdez J, Guttenberg N, Witkowski O, Pasquali A. Crosssubject EEGbased emotion recognition through neural networks with stratified normalization. Front Neurosci. 2021;15:626277.
Ding ZH, Ouyang GX, Chen H, Li XL. Closedloop transcranial magnetic stimulation of realtime EEG based on the AR mode method. Biomed Phys Eng Express. 2020;6(3):035010.
Ieracitano C, Mammone N, Hussain A, Morabito FC. A novel multimodal machine learning based approach for automatic classification of EEG recordings in dementia. Neural Netw. 2020;123:176–90.
Ji N, Ma L, Dong H, Zhang XJ. EEG signals feature extraction based on DWT and EMD combined with approximate entropy. Brain Sci. 2019;9(8):201.
RuizGómez SJ, Hornero R, Poza J, MaturanaCandelas A, Pinto N, Gómez C. Computational modeling of the effects of EEG volume conduction on functional connectivity metrics. Application to Alzheimer’s disease continuum. J Neural Eng. 2019;16(6):066019.
Luo CH, Li FL, Li PY, Yi CL, Li CB, Tao Q, et al. A survey of brain network analysis by electroencephalographic signals. Cogn Neurodyn. 2022;16(1):17–41.
Westner BU, Dalal SS, Gramfort A, Litvak V, Mosher JC, Oostenveld R, et al. A unified view on beamformers for M/EEG source reconstruction. Neuroimage. 2022;246:118789.
Wang J, Yi M, Zhang C, Bian ZJ, Wan Y, Chen RX, et al. Cortical activities of heatsensitization responses in suspended moxibustion: an EEG source analysis with sLORETA. Cogn Neurodyn. 2015;9(6):581–8.
Fogel SM, Smith CT. The function of the sleep spindle: a physiological index of intelligence and a mechanism for sleepdependent memory consolidation. Neurosci Biobehav Rev. 2011;35(5):1154–65.
Huang DM, Huang J, Qiao K, Zhong NS, Lu HZ, Wang WJ. Deep learningbased lung sound analysis for intelligent stethoscope. Mil Med Res. 2023;10(1):44.
Thoduparambil PP, Dominic A, Varghese SM. EEGbased deep learning model for the automatic detection of clinical depression. Phys Eng Sci Med. 2020;43(4):1349–60.
Lashgari E, Liang DH, Maoz U. Data augmentation for deeplearningbased electroencephalography. J Neurosci Methods. 2020;346:108885.
Gao YY, Wang XK, Potter T, Zhang JH, Zhang YC. Singletrial EEG emotion recognition using Granger Causality/Transfer Entropy analysis. J Neurosci Methods. 2020;346:108904.
Chen C, Wang H, Chen Y, Yin Z, Yang X, Ning H, et al. Understanding the brain with attention: a survey of transformers in brain sciences. BrainX. 2023;1(3):e29.
Qi Y, Ma FQ, Ge TT, Wang YM, Zhu JM, Zhang JM, et al. Abidirectional brain–computer interface for effective epilepsy control. J Zhejiang Univ Sci C Comput Electron. 2014;15(10):839–47.
Di Giamberardino P, Iacoviello D, Placidi G, Polsinelli M, Spezialetti M, editors. A brain computer interface by EEG signals from selfinduced emotions. In: 6th ECCOMAS thematic conference on computational vision and medical image processing (VipIMAGE). Porto, Portugal; 2018. p. 71321.
Acknowledgements
Not applicable.
Funding
This study was supported by STI2030 Major Projects (2021ZD0204300), and the National Natural Science Foundation of China (61803003, 62003228).
Author information
Authors and Affiliations
Contributions
HZ, QQZ and HC wrote the first draft with assistance from WGL, YB, JXH and YW. ZHL, DC, FYC, JQY, XLL edited the final draft. XQH participated in critically revising the article. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhang, H., Zhou, QQ., Chen, H. et al. The applied principles of EEG analysis methods in neuroscience and clinical neurology. Military Med Res 10, 67 (2023). https://doi.org/10.1186/s40779023005027
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40779023005027