ANNALS OF BIOMEDICAL ENGINEERING, cilt.37, ss.192-200, 2009 (SCI-Expanded)
https://link.springer.com/article/10.1007/s10439-008-9579-8
https://pubmed.ncbi.nlm.nih.gov/18941895/
ABSTRACT
The present study compares two Auto-Regressive (AR)
model based (Burg Method (BM) and Yule Walker Method) and two subspace based
(Eigen Method and Multiple Signal Classification Method) power spectral density
predictors in computing the Coherence Function (CF) to observe EEG
synchronization between right and left hemispheres. For this purpose, two
channels intracortical EEG series recorded from WAG/Rij rats (a genetic model
for human absence epilepsy) are analyzed. In tests, AR model-based predictors
result the close performance such that the CF estimations are sensitive to the
AR model order. Dealing with the subspace-based predictors; certain peaks in CF
estimations can also be detected in case of low noise subspace dimension. Besides,
they are more computational complexity. In conclusion, high order BM is
proposed in EEG synchronization. The results support that each EEG sequence
probably meets a high order AR model where the dimension of the related noise
subspace is relatively low in comparison to the model order.
DISCUSSION AND CONCLUSION
Overall results show that Set A
consists of pre-lesion EEG series, whereas the others are related to epileptic
lesions: There are more lesions in both Set D and Set E in comparison to Set B
and Set C. All data sets show sleep-wave like oscillations.
Parametric PSD predictors (BM, YWM)
are based on the AR model. The data is not windowed in performing these
methods. Therefore, application of them eliminates the assumption that the
autocorrelation sequence is zero outside the window. The BM and YWM are applied
to short time EEG series at low SNR. It is stated that these predictors exhibit
spectral line splitting in case of high SNR and sensitive to the initial phase
in case of sinusoidal noisy signals.10 In the present application, the
high order (p ≥ 20) parametric predictors provide us to useful
CF estimations in the present EE synchronization where the number of FFT is
128.
Subspace-based PSD predictors (EM,
MM) are based on an eigen-decomposition of the correlation matrix of the noisy
data. These methods assume that the data consists of some real sinusoids and
then the frequencies of sinusoidal components are estimated. For the data at
low SNR, it is difficult to determine the number of principal eigenvectors.10 In the present case study,
subspace-based predictors can also provide useful estimations when the
empirically selected subspace dimension is low (K = 10) such
that the EM is less sensitive to this parameter. Therefore, the low dimensional
EM can also be proposed in EEG synchronization for short time EEG series at low
SNR. However, adequate estimations are obtained when the number of FFT is 256
at least.
In conclusion, it can be said that
high order AR model-based predictors are more suitable than the subspace-based
predictors. The reasons of this superiority can be summarized as: When the
sampling frequency is low, short epochs are advised in EEG analysis as best
guarantee for wide sense stationary. Therefore, a stable AR model can be
yielded by the parametric estimations. High frequency resolution can be
obtained by using both BM and YWM without consuming a large memory.
Note that, the recent results support that the short time intracortical EEG series can be represented by a high order (p ≥ 20) AR model as indicated in past literature.2,6 In the former EEG study, p = 10 was provided for successful AR modeling with high accuracy. It was suggested in the latter study that p = 30 is the best order to verify normality in EEG series collected during anesthesia.