Determination of autoregressive model orders for seizure detection

In the present study, a step-wise least square estimation algorithm (SLSA), implemented in a Matlab package called as ARfit, has been newly applied to clinical data for estimation of the accurate Auto-Regressive (AR) model orders of both normal and ictal EEG series where the power spectral density (PSD) estimations are provided by the Burg Method. The ARfit module is found to be usefull in comparison to a large variety of traditional methods such as Forward Prediction Error (FPE), Akaike's Information Criteria (AIC), Minimum Description Lenght (MDL), and Criterion of Autoregressive Transfer function (CAT) for EEG discrimination. According to tests, the FPE, AIC and CAT give the identical orders for both normal and epileptic series whereas the MDL produces lower orders. Considering the resulting PSD estimations, it can be said that the most descriptive orders are provided by the SLSA. In conclusion, the SLSA can mark the seizure, since the estimated AR model orders meet the EEG complexity/regularity such that the low orders indicate an increase of EEG regularity in seizure. Then, the SLSA is proposed to select the accurate AR orders of long EEG series in diagnose for many possible future applications. The SLSA implemented by ARfit module is found to be superior to traditional methods since it is not heuristic and it is less computational complex. In addition, the more reasonable orders can be provided by the SLSA.

Determination of autoregressive model orders for seizure detection

In the present study, a step-wise least square estimation algorithm (SLSA), implemented in a Matlab package called as ARfit, has been newly applied to clinical data for estimation of the accurate Auto-Regressive (AR) model orders of both normal and ictal EEG series where the power spectral density (PSD) estimations are provided by the Burg Method. The ARfit module is found to be usefull in comparison to a large variety of traditional methods such as Forward Prediction Error (FPE), Akaike's Information Criteria (AIC), Minimum Description Lenght (MDL), and Criterion of Autoregressive Transfer function (CAT) for EEG discrimination. According to tests, the FPE, AIC and CAT give the identical orders for both normal and epileptic series whereas the MDL produces lower orders. Considering the resulting PSD estimations, it can be said that the most descriptive orders are provided by the SLSA. In conclusion, the SLSA can mark the seizure, since the estimated AR model orders meet the EEG complexity/regularity such that the low orders indicate an increase of EEG regularity in seizure. Then, the SLSA is proposed to select the accurate AR orders of long EEG series in diagnose for many possible future applications. The SLSA implemented by ARfit module is found to be superior to traditional methods since it is not heuristic and it is less computational complex. In addition, the more reasonable orders can be provided by the SLSA.

___

  • Guyton A.C: Text Book of Medical Physiology, Saunders, sec.2, 968, (1986).
  • John M., Peter J.: Epilepsy as a Dynamic Disease, Bio and Med Physics, ISBN: 978-3-540-42762-9, (2003).
  • Kannathala N., Acharyab U.R., et al: Characterization of EEG-A comparative study, Comp Meth and Prog in Biomed, 80, 17- 23 (2005).
  • Kannathala N., Choob M.L., et al: Entropies for detection of epilepsy in EEG, Comp Meth Prog in Biomed, 80, –194 (2005).
  • Cromwell D.H., Jones R.H., Kapuniai L.E: Autoregressive representation of infant EEG for the purpose of hypothesis testing and classiŞcation, EEG and Clin. Neurophysiology, 43, 317–324 (1977).
  • Jansen B.H., Boume J.R., Ward J.W: Autoregressive estimation of short segment spectra for computerized EEG analysis, IEEE Transaction on BME, 28(9), 630–637 (1981).
  • Wright J., Kydd R., Sergejew A: Autoregression models of EEG, Biol. Cybern. 62, 201–210 (1990).
  • Wang T.W., Guohui W., Huanqing F: Autoregressive model order property for sleep EEG, J of Biomedical Eng, (3), 394-6 (2004).
  • Mohammadi G., Shoushtari P. et al: Person identiŞcation by using AR model for EEG signals, Proceed of World Academy of Science, Eng and Techn, 11, 281–285 (2006).
  • Proakis J.G., Manolakis D.G: Digital Signal Processing, Prentice Hall, 3rd Ed., Chp.12, p.925–956 (1996).
  • Manolakis D.G, Ingle V.K., Kogon S.M: Statistical and Adaptive Signal Processing, Artech House Press, Ch.1, Ch.9 (2005).
  • Andrzejak R.G., Lehnertz K., Mormann F., et al: Indications of nonlinear deterministic and Şnite dimensional structures in time seires of brain electrical activity: Dependence on recording region and brain state. Physical Review E, 64, 061907 (2001).
  • Aydın S: Comparison of power spectrum predictors in computing coherence functions for intracortical EEG signals, Annals of Biomedical Engineering, 37(1), 192–200 (2009).
  • Schneider T., Neumaier A: Algorithm 808: ARŞt-A Matlab package for the estimation of parameters and eigenmodes of multivariate autoregressive models. ACM Trans Math Softw, 27, 58–65 (2001).
  • Kay S.M: Modern Spectral Estimation, Ch.7, Prentice-Hall, Englewood-Cliffs N.J.(1988).
  • Neumaier A., Schneider T.: Estimation of parameters and eigenmodes of multivariate autoregressive models, ACM Trans on Math Softw, 27(1), 27–57 (2001).
  • ChatŞeld C.: The Analysis of Time Series, Chapman and Hall (1996).
  • L¨utkepohl H.: Introduction to Multiple Time Series Analysis, 2nd ed, Cp.3, Springer-Verlag, Germany (1993).
  • Hansen P.C.: Rank-DeŞcient and Discrete Ill-Posed Problems: Numerical Aspects of Linear Inversion. SIAM, Philadelphia (1997).
  • Aydın S.: Tikhonov regularized solutions for improvement of signal-to-noise-ratio in case of auditory-evoked potentials, Medical and Biological Engineering and Computing, 46(10), 1051–1056 (2008).
  • Kashyap R.L: Inconsistency of the AIC rule for estimating the order of autoregressive Models, IEEE Trans on Auto. Cont., AC-25, 996–998 (1980).
  • Schwarz G., Estimation of the dimension of a model, Annals of Statistics, 6, 461–464 (1978).
  • Bj¨orck, A.: Numerical Methods for Least Squares Problems, SIAM, Philadelphia (1996).