Epileptic seizure prediction by non-linear methods

Methods and apparatus for automatically predicting epileptic seizures monitor and analyze brain wave (EEG or MEG) signals. Steps include: acquiring the brain wave data from the patient; digitizing the data; obtaining nonlinear measures of the data via chaotic time series analysis tools; obtaining time serial trends in the nonlinear measures; comparison of the trend to known seizure predictors; and providing notification that a seizure is forthcoming.

Skip to:  ·  Claims  ·  References Cited  · Patent History  ·  Patent History

Claims

1. A method for automatically predicting an epileptic seizure in a patient comprising the steps of:

(a) providing at least one channel of a patient's raw brain wave data, called e-data, selected from the group consisting of electroencephalogram data and magnetoencephalogram data;
(b) separating the e-data into artifact data, called f-data, and artifact-free data, called g-data, while preventing phase distortions in the data;
(c) processing g-data through a low-pass filter to produce a low-pass-filtered version of g-data, called h-data;
(d) applying at least one measure selected from the group consisting of linear statistical measures minimum and maximum, standard deviation, absolute minimum deviation, skewedness, and kurtosis, and nonlinear measures time steps per cycle, Kolmogorov entropy, first minimum in mutual information function, and correlation dimension to at least one type of data selected from the group consisting of e-data, f-data, g-data, and h-data to provide at least one time serial sequence of nonlinear measures, from which at least one indicative trend selected from the group consisting of abrupt increases, abrupt decreases, peaks, valleys, and combinations thereof is determined;
(e) comparing at least one indicative trend with at least one known seizure predictor; and
(f) determining from said comparison whether an epileptic seizure is oncoming in the patient.

2. The method as described in claim 1 wherein said at least one time serial sequence of linear and nonlinear measures is selected from the group consisting of: the standard deviation of the correlation dimension for e-data; the standard deviation of the correlation dimension for f-data; the standard deviation of the correlation dimension for g-data; the standard deviation of the correlation dimension for h-data; the Kolmogorov entropy for f-data; the Kolmogorov entropy for h-data; the first minimum in the Mutual Information Function for e-data; the first minimum in the Mutual Information Function for g-data; the average Kolmogorov entropy for e-data; the skewedness for f-data; the kurtosis for f-data, and combinations thereof.

3. The method as described in claim 1 wherein said at least one seizure predictor is selected from the group consisting of:

(a) a long-time, large value of the standard deviation of the skewedness of e-data;
(b) a long-time, large value of the standard deviation of the skewedness of f-data;
(c) quasi-periodic maxima and minima in the average Kolmogorov entropy of e-data;
(d) quasi-periodic maxima and minima in the average Kolmogorov entropy of g-data;
(e) quasi-periodic maxima and minima in the average Kolmogorov entropy of h-data;
(f) quasi-periodic maxima and minima in the average correlation dimension of e-data;
(g) quasi-periodic maxima and minima in the average correlation dimension of g-data;
(h) quasi-periodic maxima and minima in the average correlation dimension of h-data;
(i) quasi-periodic maxima and minima in the average number of time steps per cycle of e-data;
(j) quasi-periodic maxima and minima in the average number of time steps per cycle of g-data;
(k) quasi-periodic maxima and minima in the average number of time steps per cycle of h-data;
(l) coincidence of the minima and maxima selected from the group consisting of items c through k of this list and combinations thereof;
(m) a valley in the average of the minimum in the mutual information function in g-data;
(n) a peak in the standard deviation of the minimum in the mutual information function in g-data;
(o) coincidence of the extrema selected from the group consisting of items m and n of this list and combinations thereof with at least one of extrema times selected from the group consisting of items c through k of this list and combinations thereof;
(p) large decreases and increases in the standard deviation of the correlation dimension of e-data;
(q) large decreases and increases in the standard deviation of the correlation dimension of f-data;
(r) large decreases and increases in the standard deviation of the correlation dimension of g-data;
(s) large decreases and increases in the standard deviation of the correlation dimension of h-data;
(t) large, aperiodic variations in Kolmogorov entropy of f-data, including some values close to zero for >10 seconds;
(u) large, aperiodic variations in Kolmogorov entropy of h-data, including some values close to zero for >10 seconds;
(v) the average of the first minimum of the mutual information function of e-data gradually decreases below a critical value;
(w) the average of the first minimum of the mutual information function of g-data gradually decreases below a critical value;
(x) quasi-periodic maxima and minima in the average of the first minimum of the mutual information function of f-data;
(y) aperiodic maxima and minima in the average e Kolmogorov entropy of e-data;
(z) oscillations in Kolmogorov entropy of e-data of increasing or decreasing magnitude, about the average Kolmogorov entropy of item y;
(aa) small variation ratio in average skewedness (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(ab) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of f-data;
(ac) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of g-data;
(ad) small variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of e-data;
(ae) small variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(af) small variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(ag) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of e-data;
(ah) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of g-data;
(ai) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of h-data;
(aj) small variation ratio in average number of timesteps per cycle (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(ak) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of number of timesteps per cycle of g-data;
(al) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of number of timesteps per cycle of h-data;
(am) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of Kolmogorov entropy of f-data;
(an) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of Kolmogorov entropy of h-data;
(ao) small variation ratio in the average of the first minimum in mutual information function (.DELTA..sub.t /.DELTA..sub.n) of f-data;
(ap) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of first minimum in mutual information function of f-data;
(aq) small variation ratio in average correlation dimension (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(ar) large variation ratio in average (.DELTA..sub.t /.DELTA..sub.n) of the standard deviation of e-data;
(as) large variation ratio in average (.DELTA..sub.t /.DELTA..sub.n) of the standard deviation of f-data;
(at) large variation ratio in average (.DELTA..sub.t /.DELTA..sub.n) of the standard deviation of g-data;
(au) large variation ratio in average (.DELTA..sub.t /.DELTA..sub.n) of the standard deviation of h-data;
(av) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of the standard deviation of e-data;
(aw) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of the standard deviation of f-data;
(ax) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of the standard deviation of h-data;
(ay) large variation ratio in the average of the absolute average deviation (.DELTA..sub.t /.DELTA..sub.n) of e-data;
(az) large variation ratio in the average of the absolute average deviation (.DELTA..sub.t /.DELTA..sub.n) of f-data;
(ba) large variation ratio in the average of the absolute average deviation (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(bb) large variation ratio in the average of the absolute average deviation (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(bc) large variation ratio in the standard deviation of the absolute average deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of e-data;
(bd) large variation ratio in the standard deviation of the absolute average deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of f-data;
(be) large variation ratio in the standard deviation of the absolute average deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of h-data;
(bf) large variation ratio in average skewedness (.DELTA..sub.t /.DELTA..sub.n) of e-data;
(bg) large variation ratio in average skewedness (.DELTA..sub.t /.DELTA..sub.n) of f-data;
(bh) large variation ratio in average skewedness (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(bi) large variation ratio in average skewedness (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(bj) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of e-data;
(bk) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of f-data;
(bl) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of g-data;
(bm) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of h-data;
(bn) large variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of e-data;
(bo) large variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of f-data;
(bp) large variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(bq) large variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(br) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of e-data;
(bs) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of f-data;
(bt) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of g-data;
(bu) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of h-data;
(bv) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of Kolmogorov entropy of g-data;
(bw) large variation ratio in the average of the first minimum in mutual information function (.DELTA..sub.t /.DELTA..sub.n) of f-data;
(bx) large variation ratio in the average of the first minimum in mutual information function (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(by) large variation ratio in the average of the first minimum in mutual information function (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(bz) large variation ratio in average correlation dimension (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(ca) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of correlation dimension of e-data;
(cb) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of correlation dimension of h-data; and
(cc) combinations thereof.

4. The method as described in claim 1 wherein the artifact data is separated from the raw data by use of a zero-phase filter.

5. The method as described in claim 1 wherein said low-pass filter comprises a standard low-pass filter selected from the group consisting of second-order, third-order and fourth-order low-pass filters at frequencies between about 35 Hz and about 60 hz.

6. The method as described in claim 5 wherein said low-pass filter comprises a standard fourth-order low-pass filter at 50 Hz.

7. Apparatus for automatically predicting an epileptic seizure in a patient comprising:

(a) data provision means for providing at least one channel of raw brain wave data, called e-data, selected from the group consisting of electroencephalogram data and magnetoencephalogram data;
(b) separation means for separating e-data into artifact data, called f-data, and artifact-free data, called g-data, while preventing phase distortions in the data, communicably connected to said data provision means;
(c) low-pass filter means for filtering g-data to produce a low-pass filtered version of g-data, called h-data, communicably connected to said separation means;
(d) application means for applying at least one measure selected from the group of consisting of linear statistical measures minimum and maximum, standard deviation, absolute minimum deviation, skewedness, and kurtosis, and nonlinear measures time steps per cycle, Kolmogorov entropy, first minimum in mutual information function, and correlation dimension to at least one type of data selected from the group consisting of e-data, f-data, g-data, and h-data to provide at least one time serial sequence of nonlinear measures, from which at least one indicative trend selected from the group consisting of abrupt increases, abrupt decreases, peaks, valleys, and combinations thereof is determined, communicably connected to said low-pass filter means;
(e) comparison means for comparing at least one indicative trend with known seizure predictors, connected to said application means; and,
(f) determination means for determining from the comparison whether an epileptic seizure is oncoming in the patient, communicably connected to said comparison means.

8. The apparatus as described in claim 7 wherein said at least one time serial sequence of linear and nonlinear measures is selected from the group consisting of: the standard deviation of the correlation dimension for e-data; the standard deviation of the correlation dimension for f-data; the standard deviation of the correlation dimension for g-data; the standard deviation of the correlation dimension for h-data; the Kolmogorov entropy for f-data; the Kolmogorov entropy for h-data; the first minimum in the Mutual Information Function for e-data; the first minimum in the Mutual Information Function for g-data; the average Kolmogorov entropy for e-data; the skewedness for f-data; the kurtosis for f-data, and combinations thereof.

9. The apparatus as described in claim 7 wherein said at least one seizure predictor is selected from the group consisting of:

(a) a long-time, large value of the standard deviation of the skewedness of e-data;
(b) a long-time, large value of the standard deviation of the skewedness of f-data;
(c) quasi-periodic maxima and minima in the average Kolmogorov entropy of e-data;
(d) quasi-periodic maxima and minima in the average Kolmogorov entropy of g-data;
(e) quasi-periodic maxima and minima in the average Kolmogorov entropy of h-data;
(f) quasi-periodic maxima and minima in the average correlation dimension of e-data;
(g) quasi-periodic maxima and minima in the average correlation dimension of g-data;
(h) quasi-periodic maxima and minima in the average correlation dimension of h-data;
(i) quasi-periodic maxima and minima in the average number of time steps per cycle of e-data;
(j) quasi-periodic maxima and minima in the average number of time steps per cycle of g-data;
(k) quasi-periodic maxima and minima in the average number of time steps per cycle of h-data;
(l) coincidence of the minima and maxima selected from the group consisting of items c through k of this list and combinations thereof;
(m) a valley in the average of the minimum in the mutual information function in g-data;
(n) a peak in the standard deviation of the minimum in the mutual information function in g-data;
(o) coincidence of the extrema selected from the group consisting of items m and n of this list and combinations thereof with at least one of extrema times selected from the group consisting of items c through k of this list and combinations thereof;
(p) large decreases and increases in the standard deviation of the correlation dimension of e-data;
(q) large decreases and increases in the standard deviation of the correlation dimension of f-data;
(r) large decreases and increases in the standard deviation of the correlation dimension of g-data;
(s) large decreases and increases in the standard deviation of the correlation dimension of h-data;
(t) large, aperiodic variations in Kolmogorov entropy of f-data, including some values close to zero for >10 seconds;
(u) large, aperiodic variations in Kolmogorov entropy of h-data, including some values close to zero for >10 seconds;
(v) the average of the first minimum of the mutual information function of e-data gradually decreases below a critical value;
(w) the average of the first minimum of the mutual information function of g-data gradually decreases below a critical value;
(x) quasi-periodic maxima and minima in the average of the first minimum of the mutual information function of f-data;
(y) aperiodic maxima and minima in the average Kolmogorov entropy of e-data;
(z) oscillations in Kolmogorov entropy of e-data of increasing or decreasing magnitude, about the average Kolmogorov entropy of item y;
(aa) small variation ratio in average skewedness (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(ab) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of f-data;
(ac) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of g-data;
(ad) small variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of e-data;
(ae) small variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(af) small variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(ag) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of e-data;
(ah) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of g-data;
(ai) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of h-data;
(aj) small variation ratio in average number of timesteps per cycle (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(ak) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of number of timesteps per cycle of g-data;
(al) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of number of timesteps per cycle of h-data;
(am) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of Kolmogorov entropy of f-data;
(an) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of Kolmogorov entropy of h-data;
(ao) small variation ratio in the average of the first minimum in mutual information function (.DELTA..sub.t /.DELTA..sub.n) of f-data;
(ap) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of first minimum in mutual information function of f-data;
(aq) small variation ratio in average correlation dimension (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(ar) large variation ratio in average (.DELTA..sub.t /.DELTA..sub.n) of the standard deviation of e-data;
(as) large variation ratio in average (.DELTA..sub.t /.DELTA..sub.n) of the standard deviation of f-data;
(at) large variation ratio in average (.DELTA..sub.t /.DELTA..sub.n) of the standard deviation of g-data;
(au) large variation ratio in average (.DELTA..sub.t /.DELTA..sub.n) of the standard deviation of h-data;
(av) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of the standard deviation of e-data;
(aw) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of the standard deviation of f-data;
(ax) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of the standard deviation of h-data;
(ay) large variation ratio in the average of the absolute average deviation (.DELTA..sub.t /.DELTA..sub.n) of e-data;
(az) large variation ratio in the average of the absolute average deviation (.DELTA..sub.t /.DELTA..sub.n) of f-data;
(ba) large variation ratio in the average of the absolute average deviation (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(bb) large variation ratio in the average of the absolute average deviation (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(bc) large variation ratio in the standard deviation of the absolute average deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of e-data;
(bd) large variation ratio in the standard deviation of the absolute average deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of f-data;
(be) large variation ratio in the standard deviation of the absolute average deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of h-data;
(bf) large variation ratio in average skewedness (.DELTA..sub.t /.DELTA..sub.n) of e-data;
(bg) large variation ratio in average skewedness (.DELTA..sub.t /.DELTA..sub.n) of f-data;
(bh) large variation ratio in average skewedness (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(bi) large variation ratio in average skewedness (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(bj) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of e-data;
(bk) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of f-data;
(bl) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of g-data;
(bm) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of skewedness of h-data;
(bn) large variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of e-data;
(bo) large variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of f-data;
(bp) large variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(bq) large variation ratio in average kurtosis (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(br) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of e-data;
(bs) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of f-data;
(bt) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of g-data;
(bu) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of kurtosis of h-data;
(bv) small variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of Kolmogorov entropy of g-data;
(bw) large variation ratio in the average of the first minimum in mutual information function (.DELTA..sub.t /.DELTA..sub.n) of f-data;
(bx) large variation ratio in the average of the first minimum in mutual information function (.DELTA..sub.t /.DELTA..sub.n) of g-data;
(by) large variation ratio in the average of the first minimum in mutual information function (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(bz) large variation ratio in average correlation dimension (.DELTA..sub.t /.DELTA..sub.n) of h-data;
(ca) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of correlation dimension of e-data;
(cb) large variation ratio in standard deviation (.DELTA..sigma..sub.t /.DELTA..sigma..sub.n) of correlation dimension of h-data; and
(cc) combinations thereof.

10. The apparatus as described in claim 7 wherein said separation means comprises a zero-phase filter.

11. The apparatus as described in claim 10 wherein said zero-phase filter is embodied in a programmed integrated-circuit semiconductor chip.

12. The apparatus as described in claim 7 wherein said low-pass filter means comprises a standard low-pass filter selected from the group consisting of second-order, third-order, and fourth-order low-pass filters at frequencies between about 35 Hz and about 60 Hz.

13. The apparatus as described in claim 7 wherein said low-pass filter comprises a standard fourth-order low-pass filter at about 50 Hz.

14. The apparatus as described in claim 7 further comprising notification means for providing notification that a seizure is oncoming in the patient, the notification means being communicably connected to said determination means.

Referenced Cited
U.S. Patent Documents
5311876 May 17, 1994 Olsen et al.
5349962 September 27, 1994 Lockard et al.
5392788 February 28, 1995 Hudspeth
5630425 May 20, 1997 Panescu et al.
Other references
  • A. Babloyantz and A. Destexhe, "Low-Dimensional Chaos in an Instance of Epilepsy," Proc. Natl. Acad. Sci. USA 83, 3513-3517 (1986). D.F. Elliott and K. R. Rao, Fast Transforms, Analyses, Applications, Academic Press, 1982. A. M. Fraser and H. L. Swinney, "Independent Coordinates for Strange Attractors from Mutual Information," Phys. Rev A 33, 1134-1140 (1986). J. Glanz, "Do Chaos-Control Techniques Offer Hope for Epilepsy?" Science 265, 1174 (1994). J. Gotman, "Seizure Recognition and Analysis," Long-term Monitoring in Epilepsy (EEG Suppl. No. 37), pp. 133-145 (1985). G. A. Korn and T. M. Korn, Mathematical Handbook for Scientists and Engineers, McGraw-Hill Book Company, New York, 1968. D. E. Olsen, J. A. Cristion, and C. W. Spaur, "Automatic Detection of Epileptic Seizures Using Electroencephalographic Signals," Johns Hopkins APL Techn. Digest 12, 182-191 (1991). J.P.M. Pijn, "Quantitative Evaluation of EEG Signals in Epilepsy--Nonlinear Associations, Time Delays, and Nonlinear Dynamics," Ph.D. Thesis, University of Amsterdam, 1990. L. R. Rabiner and B. Gold, Theory and Application of Digital Signal Processing, Prentice Hall Inc., 1975. S.M. Selby and B. Girling, Standard Mathematical Tables, p. 390, The Chemical Rubber Co. (14th edition, 1965). S. J. Schiff, K. Jerger, D. H. Duong, T. Chang, M. L. Spano, and W.L. Ditto, "Controlling Chaos in the Brain," Nature 370, 615-620 (1994). J. C. Schouten, F. Takens, and C. M. van den Bleek, "Maximum-Likelihood Estimation of the Entropy of an Attractor," Phys. Rev. E 49, 126-129 (1994). J. C. Schouten, F. Takens, and C. M. van den Bleek, "Estimation of the Dimension of a Noisy Attractor," Phys. Rev. E 50, 1851-1861 (1994). F. Takens, "On the Numerical Determination of the Dimension of an Attractor," in Dynamical Systems and Bifurcations, ed. B.L.J. Braaksma, H.W. Broer, and F. Takens, Lecture Notes in Mathematics 1125, 99-106 (1984) Springer-Verlag, Berlin. F. J. Taylor, Digital Filter Design Handbook, Marcel Dekker, Inc., Publ., New York and Basel, 1983. J. Theiler, "On the Evidence for Low-Dimensional Chaos in an Epileptic Electroencephalogram," Phys. Lett. A 196, 335-341 (1995). R. C. Watt and S. R. Hameroff, "Phase Space Analysis of Human EEG during General Anesthesia," Ann. N.Y. Acad. Sci. 504, 286-288 (1987). T. Elbert, W. J. Ray, Z. J. Kowalik, J. E. Skinner, K. E. Graf, N. Birbaumer, "Chaos and Physiology: Deterministic Chaos in Excitable Cell Assemblies," Physiological Reviews 74 (1994) 1-47. S. Blanco, R. Quian Quiroga, O. A. Rosso, and S. Kochen, "Time-frequency analysis of electroencephalogram series," Physical Review E, The American Physical Society, vol. 51, No. 3, Mar. 1995. Weiming Yang and Mingzhou Ding, "Preserving Chaos; Contrl strategies to preserve complex dynamics with potential relevance to biological disorders," Review E, The American Physical Society, vol. 51, No. 1 Jan. 1995.
Patent History
Patent number: 5857978
Type: Grant
Filed: Mar 20, 1996
Date of Patent: Jan 12, 1999
Assignee: Lockheed Martin Energy Systems, Inc. (Oak Ridge, TN)
Inventors: Lee M. Hively (Knoxville, TN), Ned E. Clapp (Knoxville, TN), C. Stuart Daw (Knoxville, TN), William F. Lawkins (Knoxville, TN)
Primary Examiner: Jennifer Bahr
Assistant Examiner: Ryan Carter
Attorney: J. Kenneth Davis
Application Number: 8/619,030
Classifications
Current U.S. Class: Detecting Brain Electric Signal (600/544)
International Classification: A61B 504;