Processing for multi-channel signals

Information

  • Patent Grant
  • 8786624
  • Patent Number
    8,786,624
  • Date Filed
    Wednesday, June 2, 2010
    14 years ago
  • Date Issued
    Tuesday, July 22, 2014
    9 years ago
Abstract
Method and apparatus for improved processing for multi-channel signals. In an exemplary embodiment, an anomaly metric is computed for a multi-channel signal over a time window. The magnitude of the anomaly metric may be used to determine whether an anomaly is present in the multi-channel signal over the time window. In an exemplary embodiment, the anomaly metric may be a condition number associated with the singular values of the multi-channel signal over the time window, as further adjusted by the number of channels to produce a data condition number. Applications of the anomaly metric computation include the scrubbing of signal archives for epileptic seizure detection/prediction/counter-prediction algorithm training, pre-processing of multi-channel signals for real-time monitoring of bio-systems, and boot-up and/or adaptive self-checking of such systems during normal operation.
Description
TECHNICAL FIELD

The present disclosure relates generally to methods for processing multi-channel signals. In particular, the present disclosure relates to improved processing of multi-channel signals by detecting and/or treating possible anomalies in the multi-channel signals


INCORPORATION BY REFERENCE

All publications and patent applications mentioned in this specification are herein incorporated by reference to the same extent as if each individual publication or patent application was specifically and individually indicated to be incorporated by reference.


BACKGROUND OF THE INVENTION

A common task encountered in the field of signal processing is the sampling and processing of a physical state using multiple, ideally independent, signal sensors. The diversity of the resulting multi-sensor or multi-channel signal typically reveals more information about the underlying sampled state than can be obtained from employing a single sensor.


Multi-channel signal processing is utilized in biomedical applications. For example, in the field of neurological monitoring for epileptic seizure prediction, multiple electrodes may be implanted in diverse locations on or in a patient's brain to monitor the susceptibility of the patient to enter into an epileptic seizure. The multi-channel signals generated by the electrodes may be processed to, e.g., alert the patient and/or medical personnel of a high likelihood of imminent seizure. See, e.g., U.S. Patent Publication No. 2008/0183096 A1, “Systems and Methods for Identifying a Contra-ictal Condition in a Subject,” filed Jan. 25, 2008, assigned to the assignee of the present application, the contents of which are hereby incorporated by reference in their entirety. The signals may also be stored and processed offline to, e.g., train customized algorithms for estimating the likelihood that a patient will experience an imminent seizure. See, e.g., U.S. Pat. No. 6,678,548, “Unified probabilistic framework for predicting and detecting seizure onsets in the brain and multitherapeutic device,” the contents of which are hereby incorporated by reference in their entirety.


Other applications of multi-channel signal processing include the reception of wireless signals by a communications device using multiple antennas, geological monitoring of seismic activity for earthquake prediction, stereo imaging using multiple video cameras, etc.


When multi-channel signals are sampled over an extended period of time, artifacts or anomalies often appear in the signal. Such anomalies may be due to interference from external sources, disruptions to the power supply of the sensors, and/or other sources. Left untreated, such anomalies may degrade the quality of the measured signal and disrupt the accuracy of any subsequent processing of the multi-channel signal.


It would be desirable to have techniques to detect the presence of anomalies in a multi-channel signal, and to optimize the processing of a signal containing such anomalies.


SUMMARY OF THE INVENTION

An aspect of the present disclosure provides a method for displaying information associated with a multi-channel signal to a user, the method comprising: accepting input from the user selecting a metric to be displayed; displaying at least one time-series plot of the selected metric associated with the multi-channel signal; using a backdrop pattern, indicating on the at least one time-series plot portions of the plot having at least one identified characteristic; for each of the at least one time-series plot displayed, displaying a time event index wherein the corresponding metric meets a predetermined condition; accepting input from the user as to whether to display further information associated with the time event index; and displaying the further information associated with the time event index if the user so specifies.


Another aspect of the present disclosure provides a method for detecting anomalies in a multi-channel signal, the method comprising: sampling the multi-channel signal over a time window; computing an anomaly metric for the multi-channel signal over the time window; and identifying the presence of an anomaly based on the magnitude of the anomaly metric; the computing an anomaly metric comprising: computing a condition number of the multi-channel signal over the time window; and adjusting the condition number based on a parameter of the multi-channel signal to generate a data condition number (DCN); the identifying the presence of an anomaly comprising comparing the magnitude of the data condition number (DCN) to at least one threshold; the method further comprising generating one of the at least one threshold, the generating comprising: generating a histogram of number of anomalies detected for each of a plurality of candidate thresholds; generating a cumulative distribution function from the histogram, the cumulative distribution function mapping each candidate threshold to a percentage value; and determining the generated threshold as the candidate threshold mapped to a corresponding predetermined percentage value by the cumulative distribution function.


Yet another aspect of the present disclosure provides a method for detecting anomalies in a multi-channel signal, the method comprising: sampling the multi-channel signal over a time window; computing an anomaly metric for the multi-channel signal over the time window; and identifying the presence of an anomaly based on the magnitude of the anomaly metric; the computing an anomaly metric comprising: computing a condition number of the multi-channel signal over the time window; and adjusting the condition number based on a parameter of the multi-channel signal to generate a data condition number (DCN); the identifying the presence of an anomaly comprising comparing the magnitude of the rate of change of the data condition number (DCN) to at least one threshold.


Yet another aspect of the present disclosure provides a method for detecting anomalies in a multi-channel signal, the method comprising: sampling the multi-channel signal over a time window; computing an anomaly metric for the multi-channel signal over the time window; and identifying the presence of an anomaly based on the magnitude of the anomaly metric; the computing an anomaly metric comprising: computing a condition number of the multi-channel signal over the time window; and adjusting the condition number based on a parameter of the multi-channel signal to generate a data condition number (DCN); the multi-channel signal comprising a signal sampled from a plurality of electrodes implanted on or in a patient's brain, the method further comprising: generating a DCN time series corresponding to a plurality of time windows; generating an anomaly log based on the DCN time series; merging anomalies in the anomaly log separated by less than a minimum separation to generate a modified anomaly log; identifying segments of the multi-channel signal corresponding to anomalies in the modified anomaly log; and outputting time-expanded versions of the identified segments to a record; the identifying the presence of an anomaly comprising matching the DCN time series to at least one known pattern of DCN time series corresponding to an anomalous condition.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 depicts an exemplary embodiment of the present disclosure of a system for quality management of measured electrical signals generated by electrodes placed on or within the brain.



FIGS. 2A, 2B and 2C illustrate possible signal anomalies that may be detected by the anomaly detector/processor in the multi-channel signal.



FIG. 3A depicts an exemplary method according to the present disclosure, wherein a metric known as a “data condition number,” or DCN, is derived for a given multi-channel signal.



FIG. 3B shows a matrix A[k] derived from the multi-channel signal.



FIG. 4A illustrates a data condition number (DCN) time series obtained from a sample multi-channel signal using the steps of FIG. 3A and thresholding for event detection.



FIG. 4B illustrates a plot used in a procedure by which specific thresholds may be automatically chosen for any patient or group of patients.



FIGS. 4C and 4D illustrate an alternative exemplary embodiment of a technique for choosing an optimum threshold T1*.



FIG. 5 depicts an exemplary method according to the present disclosure for taking action in response to the detection of an anomaly.



FIG. 6 depicts an alternative exemplary embodiment of the present disclosure, wherein the techniques disclosed hereinabove are applied in the context of a real-time patient monitoring and neurological event detection system.



FIG. 7 depicts a generalized block diagram of the real-time analysis system featuring an anomaly pre-processing block according to the present disclosure.



FIG. 8 depicts an exemplary method according to the present disclosure that may be implemented, e.g., using the system shown in FIG. 7.



FIGS. 9A and 9B illustrate exemplary plots of data condition number (DCN) versus time, corresponding to certain specific operating scenarios.



FIGS. 10A and 10B illustrate an exemplary embodiment of a graphical display interface for conveniently displaying information associated with the techniques of the present disclosure.





DETAILED DESCRIPTION OF THE INVENTION

The detailed description set forth below in connection with the appended drawings is intended as a description of exemplary embodiments of the present invention and is not intended to represent the only exemplary embodiments in which the present invention can be practiced. The term “exemplary” used throughout this description means “serving as an example, instance, or illustration,” and should not necessarily be construed as preferred or advantageous over other exemplary embodiments. The detailed description includes specific details for the purpose of providing a thorough understanding of the exemplary embodiments of the invention. It will be apparent to those skilled in the art that the exemplary embodiments of the invention may be practiced without these specific details. In some instances, well known structures and devices are shown in block diagram form in order to avoid obscuring the novelty of the exemplary embodiments presented herein.


While the discussion below focuses on measuring electrical signals generated by electrodes placed near, on, or within the brain or nervous system (EEG signals) of subjects and subject populations for the determination of when an epileptic subject is in a condition susceptible to seizure, it should be appreciated that the techniques of the present disclosure are not limited to measuring EEG signals or to determining when the subject is susceptible to seizure. For example, the techniques could also be used in systems that measure one or more of a blood pressure, blood oxygenation (e.g., via pulse oximetry), temperature of the brain or of portions of the subject, blood flow measurements, ECG/EKG, heart rate signals, respiratory signals, chemical concentrations of neurotransmitters, chemical concentrations of medications, pH in the blood, or other physiological or biochemical parameters of a subject.


The present disclosure may also be applicable to monitoring other neurological or psychiatric disorders and identifying a condition or state for such disorders in which the subject is unlikely to experience some adverse effect. For example, the present disclosure may also be applicable to monitoring and management of sleep apnea, Parkinson's disease, essential tremor, Alzheimer's disease, migraine headaches, depression, eating disorders, cardiac arrhythmias, bipolar spectrum disorders, or the like.


Non-biomedical applications of the techniques described herein are also contemplated to be within the scope of the present disclosure.



FIG. 1 depicts an exemplary embodiment of the present disclosure of a system for quality management of measured electrical signals generated by electrodes implanted into the brain, subdurally, epidurally, partially or fully in the skull, between the skull and one or more layers of the patient's scalp, or on the exterior of the patient's head. In some embodiments, the electrodes are configured to be inserted through a single opening in a skull of a patient and to be dispersed on a surface of the brain. The surface is preferably the cortical surface of the brain, but may alternatively be any other suitable surface or layer of the brain, such as the dura mater. Note the exemplary embodiment depicted in FIG. 1 is shown for illustrative purposes only, and is not meant to limit the scope of the present disclosure to any particular embodiment shown.


In FIG. 1, electrodes 110 are implanted on or in the brain of a patient 100, e.g., underneath the dura mater and on a cortical surface of the patient's brain. Each of the electrodes generates a corresponding signal that is input to a data pre-processing module 120. The electrodes may be directly connected to the data pre-processing module 120 through a series of electrical leads 115a, or they may be wirelessly connected to the data pre-processing module 120 over a wireless link. In an exemplary embodiment, data pre-processing module 120 may, e.g., digitize the plurality of received electrical signals 115b to generate a multi-channel signal 120a for further processing.


Multi-channel signal 120a is input to an anomaly detector/processor 130. In an exemplary embodiment, the anomaly detector/processor 130 may utilize techniques further disclosed hereinbelow to identify the presence of signal anomalies in the multi-channel signal 120a. As further disclosed hereinbelow, anomaly detector/processor 130 may also take further action to address the anomalies detected, e.g., flagging the portions of the multi-channel signal corresponding to the detected anomalies in a log file 130a for optional manual review by a human operator.


In the exemplary embodiment shown, the log file 130a from the anomaly detection/processing module 130 is provided along with the multi-channel signal 120a to a data processing/adaptive algorithm training module 140. In an exemplary embodiment, module 140 may utilize the multi-channel signal 120a, coupled with information from the log file 130a about which portions of the multi-channel signal 120a contain anomalies, to train an adaptive algorithm to identify conditions under which patient 100 is susceptible to seizure. An exemplary system is described in Snyder, et al., “The statistics of a practical seizure warning system,” Journal of Neural Engineering, vol. 5, pp. 392-401 (2008), the contents of which are incorporated by referenced herein in its entirety. In an exemplary embodiment, module 140 may, e.g., automatically de-emphasize portions of multi-channel signal 120a corresponding to signal anomalies, and emphasize other portions of the signal 120a, to configure adaptive weights for a seizure prediction algorithm 140a. In alternative exemplary embodiments, a human operator may manually review portions of multi-channel signal 120a that have been flagged in the log file 130a, and decide whether such portions may be used for adaptive algorithm training.


In an exemplary embodiment, the log file 130a need not be limited to a single file residing in a single piece of storage hardware. For example, the log file 130a can be an extensive archive of intracranial EEG patterns that can be used to develop a predictive neurosensing device for managing seizures by mining the archive for signal patterns over the patient population. This archive may be stored for multi-user access and processing in, e.g., a server or cloud computing system.



FIGS. 2A, 2B, and 2C illustrate possible signal anomalies that may be detected by the anomaly detector/processor 130 in the multi-channel signal 120a. Note the anomaly patterns identified in FIGS. 2A, 2B, and 2C are for illustrative purposes only, and are not meant to limit the scope of the present disclosure to any particular anomaly patterns highlighted.


In FIG. 2A, reference numeral 210 points to an instance of a “spikes without phase reversal” anomaly in a referential multi-channel signal, or “spikes without double-phase reversal” anomaly in a bipolar multi-channel signal, possibly due to non-physiological interference sources. Reference numeral 220 points to an instance of a “flatline” anomaly in the multi-channel signal. In FIG. 2B, reference numeral 230 points to an instance of a “line noise” anomaly in the multi-channel signal. Such an anomaly may correspond to, e.g., powerline noise, i.e., 50 or 60 Hz noise and/or harmonics superimposed on the multi-channel signal. In FIG. 2C, reference numeral 240 points to an instance of a “saturation” anomaly in the multi-channel signal.


Other possible anomalies in a multi-channel signal (not shown) include episodic artifacts such as motion (large swings in the multi-channel signal), DC shifts (different DC levels between different channels or across a single channel), pops (exponential decay from amplifier highpass characteristic of a sudden change in the DC level of a signal), and glitches (e.g., 50 ms burst transients in the signal). Long-term anomalies may include deterioration trends in the system, and/or channels of persistently poor quality. Such anomalies and others not explicitly enumerated are contemplated to be within the scope of the present disclosure.



FIG. 3A depicts an exemplary method 300 according to the present disclosure, wherein a specific anomaly metric known as a “data condition number,” or DCN, is derived for a given multi-channel signal. In an exemplary embodiment, the DCN may be used to help identify the presence of anomalies in the multi-channel signal. The DCN as described with reference to FIG. 3A has been found to provide a reliable indicator for the presence of anomalies in a multi-channel signal, and is readily applicable to a wide variety of scenarios.


In FIG. 3A, at step 310, the number of dimensions of a multi-channel signal is defined as the variable N. The number of dimensions may correspond to the total number of sensors, e.g., the number of independent electrodes placed in a patient's brain in the epilepsy monitoring system of FIG. 1.


At step 320, the multi-channel signal is divided in time into windows of duration T, with each of the time windows being indexed by a counter k. According to the present disclosure, the multi-channel signal may be a discrete-time signal sampled at a rate of S Hz. In this case, each time window k may contain a total of T·S discrete-time samples multiplied by N channels (or sensors), which may be arranged to form a T·S matrix by N A[k] at step 330. The matrix A[k] is also shown in FIG. 3B. In an exemplary embodiment, the matrix A[k] is generally rectangular.


In an exemplary embodiment, T may be 5 seconds, S may be 400 Hz, and N may be 128 for an epilepsy monitoring unit such as depicted in FIG. 1. In an alternative exemplary embodiment, such as the real-time patient monitoring and neurological event detection system 60 of FIG. 6 utilizing implanted electrodes, N may be 16.


In an exemplary embodiment, the time windows k may be chosen to collectively span the entire duration of the multi-channel signal, i.e., the time windows are non-overlapping and contiguous in time. In alternative exemplary embodiments, the time windows need not be contiguous in time, and may be spread out over the duration of the multi-channel signal. In this way, the time windows effectively sub-sample the total duration of the multi-channel signal. This sub-sampling may result in fewer matrices A[k] to be processed as compared to using contiguous time windows, reducing the computational complexity. In yet other alternative exemplary embodiments, the time windows also may be made overlapping in time. In an exemplary embodiment, the signals for the channels may be represented using a referential montage, wherein each signal is measured as an electrical potential with respect to a common (“ground”) contact placed somewhere else in the body; e.g., an “earlobe” in a scalp electroencephalogram.


At step 340, a condition number C[k] is computed for each matrix A[k] as follows:











C


[
k
]


=


max


{

σ

i





k


}



min


{

σ

i





k


}




;




(

Equation





1

)








wherein {σik} is the set of singular values of each matrix A[k]. In an exemplary embodiment, the singular values of each matrix A[k] may be computed using a singular-value decomposition (SVD) well-known in the art:

A[k]=USVT;  (Equation 2)

wherein U and V are both square unitary matrices, and S contains the singular values of A[k]. One of ordinary skill in the art will appreciate that a variety of software tools are available to calculate the singular values, or approximate singular values, of a given matrix, for the purpose of deriving the matrix S. Such software tools, include, e.g., the publicly available software packages LAPACK or EISPACK.


At step 350, the condition number C[k] is further refined by conversion into a “data condition number” DCN[k]. In an exemplary embodiment, DCN[k] may be computed as:










DCN


[
k
]


=

1
+



(


C


[
k
]


-
1

)

N

.






(

Equation





3

)







The conversion from a condition number into a data condition number may be done to compensate for an expected quasi-linear increase in the condition number due to the number of channels N. Following the conversion, the data condition number may generally be processed independently of the number of channels. One of ordinary skill in the art will appreciate that alternative techniques for accounting for the number of channels may be employed. For example, the condition number need not be converted to a data condition number, and the thresholds used to determine the presence of anomalies may instead be adjusted by the number of channels. Such alternative exemplary embodiments are also contemplated to be within the scope of the present disclosure.


One of ordinary skill in the art will further appreciate that, in alternative exemplary embodiments, either the condition number C[k] or thresholds used may be alternatively, or further, normalized with respect to any source of variation in C[k] that is not of interest, i.e., not indicative of an anomaly. For example, variables such as the window size, signal measurement bandwidth, electrode montage, etc., may also be accounted for in converting the condition number C[k] to the data condition number DCN[k], and/or choosing the thresholds against which the DCN[k] is compared. Such additional coefficients and parameters may be calibrated empirically, and one of ordinary skill in the art may readily modify Equation 3 accordingly to account for such coefficients and parameters.


According to the present disclosure, the magnitude of the data condition number DCN[k] may serve as an indicator of whether a time window k of a multi-channel signal contains anomalies. For example, the DCN[k] may range from a value of 1, indicative of “healthy” data that is lacking in anomalies, to an arbitrarily large value ∞, indicative of “ill” data corresponding, e.g., to complete flatlining over the time window of interest. Evaluating the magnitude of each DCN[k] may provide an indication of whether a signal anomaly is present in the corresponding time window k.


In alternative exemplary embodiments (not shown), the DCN may be computed for a subset of the total number of signal sensors by constructing the matrix A[k] using such subset of signals, and employing the number of signals in the subset for the variable N. For example, instead of employing all N channels to construct the matrix A[k], only the signals from a subset N−1 of the channels may be used. This may be advantageous when, e.g., one of the signal sensors is known to be faulty. Such exemplary embodiments are contemplated to be within the scope of the present disclosure.


In alternative exemplary embodiments, the DCN may be computed from a subset N−1 of the channels to account for the effect of average montages, circular bipolar montages, or any other electrode montages, wherein one channel is known a priori to be a linear combination of all remaining channels, and therefore the matrix A[k] is rank-deficient.


One of ordinary skill in the art will appreciate that in light of the present disclosure, various alternative anomaly metrics to the condition number for detecting the presence of anomalies may be derived, based on, e.g., determining the independence or correlation between two or more of the channels. These alternative metrics may be derived based on the assumption that healthy data, especially in neurological multi-channel signals, is associated with independence among the channels, while data containing anomalies is associated with a lack of independence among the channels.


For example, for a square matrix A[k], a generalized condition number C′[k] may be derived as:

C′[k]=Norm(A[k])·Norm(A−1[k]);  (Equation 4)

wherein Norm(·) denotes a norm of the matrix in parentheses, and A[k] is the inverse of the square matrix A[k]. For example, norms such as the L1-norm and L-infinity norm are well-known in the art, and may be applied to compute an anomaly metric associated with a matrix A[k] in the following manner. One of ordinary skill in the art will appreciate that the L1-norm may be defined as the maximum of the column sums of A[k], and the L-infinity norm may be defined as the maximum of the row sums.


As the computation of the generalized condition number C′ [k] and other types of condition numbers may require that the matrix A[k] be, e.g., a square matrix, or have other pre-specified dimensions, the data from the multi-channel signals may be suitably modified to ensure that the matrix A[k] takes on the proper form. For example, to ensure that the matrix A[k] is a square matrix, the size of the time window may be chosen such that the number of discrete time samples for each channel is equal to the total number of channels. Alternatively, the discrete time samples over a given time window may be sub-sampled at regular intervals to arrive at a square matrix A[k]. Such exemplary embodiments are contemplated to be within the scope of the present disclosure.



FIG. 4A illustrates a data condition number time series 400 obtained from a sample multi-channel signal using the steps of FIG. 3A and thresholding for event detection. In FIG. 4A, the horizontal axis depicts the time window index k, while the vertical axis depicts the value of DCN[k] corresponding to k. Note the absolute values of DCN[k] shown in FIG. 4A and discussed below are given for illustration only, and are not meant to restrict the techniques of the present disclosure to any particular values of DCN[k] shown. In particular, the absolute values of DCN[k] may depend on various parameters of the application.


In FIG. 4A, DCN[k] is shown to take on values generally less than 500 prior to time window k=K1. Between k=K1 and k=K2, DCN[k] is approximately 1000, while at time k=K2, DCN[k] is greater than 2000. In an exemplary embodiment, one or more predetermined thresholds may be chosen to identify the presence of anomalies in the multi-channel signal. For example, in FIG. 4A, the DCN[k] may be compared to a single threshold T1, and values of DCN[k] that are greater than the threshold T1 may indicate the presence of an “event” in the corresponding time window k. In an exemplary embodiment, the presence of an event may be taken to indicate the presence of an anomaly.


Further shown in FIG. 4A is an anomaly plot 410 generated from the time series 400 by comparison with threshold T1. Note anomaly plot 410 identifies the presence of two events A and B, possibly corresponding to two anomalies.


In alternative exemplary embodiments, two or more thresholds may be chosen for more precise categorization of the DCN. For example, in an exemplary embodiment, two thresholds T1 and T2 may be chosen, wherein T1<T2. In this exemplary embodiment, if the DCN is less than T1, then a lack of anomaly in the time window may be declared. If the DCN is greater than T2, then an anomaly can be automatically declared. If the DCN is between T1 and T2, then further processing, such as manual inspection of the multi-channel signal, may be performed to determine whether an anomaly is actually present. In an exemplary embodiment, patient-specific thresholds may be chosen that are customized to an individual patient whose neurological or other biological state is being monitored by the multiple sensors. Thresholds may be set differently for different patients, to account for the unique characteristics of each patient's bio-signals.


According to an aspect of the present disclosure, specific thresholds may be automatically chosen for any patient or group of patients using a procedure as described with reference to the plot 400A shown in FIG. 4B. In FIG. 4B, the horizontal axis shows a candidate threshold T1, while the vertical axis tallies a number of events detected for the corresponding candidate threshold. Based on the plot 400A, an optimum threshold T1* may be chosen to, e.g., maximize the probability of detecting an anomalous event while minimizing the probability of “false” positives. In an exemplary embodiment, as shown in FIG. 4B, such an optimum threshold T1* may be chosen as, e.g., an inflection point of the plot 400A having the largest abscissa from amongst all inflection points, wherein an inflection is defined as a point where the curvature changes sign, e.g., from facing down to facing up. Alternatively, if such an inflection point does not exist (e.g., in exponential decay), the threshold T1* may be chosen as the minimum abscissa where the plot 400A drops to between 1% and 2% of its maximum value Nmax.



FIGS. 4C and 4D illustrate an alternative exemplary embodiment of a technique for choosing an optimum threshold T1*. FIG. 4C illustrates an exemplary histogram of DCN (horizontal axis), with the vertical axis showing the number of instances of the corresponding DCN in an arbitrary multi-channel signal (not shown). From the histogram of FIG. 4C, one of ordinary skill in the art may readily derive a corresponding cumulative distribution function (CDF), as shown in FIG. 4D, which indicates what percentage of events correspond to DCN's equal to or less than the shown DCN. One of ordinary skill in the art will appreciate that by selecting an appropriate percentage (e.g., 90% as illustrated in FIG. 4D), a corresponding threshold DCN of T1*, such as T1*90% in FIG. 4D, may be automatically determined using the functional mapping provided by the CDF. One of ordinary skill in the art will further appreciate that in an exemplary embodiment, a computer may be programmed to perform such threshold selection, such that the choice of threshold may be made fully automatic without the need for manual intervention.


While certain exemplary techniques for automatically a choosing a suitable threshold T1* have been disclosed hereinabove, one of ordinary skill in the art will appreciate that alternative techniques not explicitly described may be readily derived in light of the present disclosure. Such alternative exemplary embodiments are contemplated to be within the scope of the present disclosure.


In alternative exemplary embodiments, additional properties of the DCN may be analyzed to further aid in the detection of anomalies in the multi-channel signal. For example, the rate of change of the DCN over a predetermined interval of time may also be utilized to detect the presence of an anomaly. Such modifications to the DCN and others not explicitly described will be clear to one of ordinary skill in the art, and are contemplated to be within the scope of the present disclosure.



FIG. 5 depicts an exemplary method according to the present disclosure for taking action in response to the detection of an anomaly.


At step 510, a DCN time series such as 400 in FIG. 4A is generated from a multi-channel signal. In an exemplary embodiment, such a time series may be generated according to the method 300 depicted in FIG. 3A. However, one of ordinary skill in the art will appreciate that alternative methods may be employed to generate a suitable time series in light of the present disclosure.


At step 520, an anomaly log is generated based on the DCN time series. Such an anomaly log may identify, e.g., time window indices k in the DCN time series corresponding to detected anomalies. For example, in the exemplary embodiment wherein DCN[k] is compared to a single threshold T1 to determine the presence of an anomaly, the anomaly log may record all time window indices k wherein DCN[k] is larger than T1. As such, the anomaly log may effectively capture the relevant information from the anomaly plot 410.


One of ordinary skill in the art will appreciate that the information in an anomaly log may be recorded in several ways. For example, each line in the anomaly log may record the time index k associated with the beginning of a detected anomaly, and the corresponding time duration of the detected anomaly. Alternatively, the start and stop time indices associated with each detected anomaly may be recorded. Such exemplary embodiments are contemplated to be within the scope of the present disclosure.


At step 530, the complexity of the anomaly log may be reduced by merging separate anomalies that are separated by less than a minimum time separation. For example, assume two anomalies each of duration 10 are found in an anomaly log starting at time windows k=1 and k=12, i.e., the two anomalies are separated by a time duration of Δk=1. If a minimum time separation is defined as Δkmin=5, then the two anomalies may be merged to form a single anomaly, which can be recorded in a simplified anomaly log as a single merged anomaly of duration Δk=21 starting at k=1.


By performing the merging as described at step 530, the number of recorded anomalies and the size of the resulting anomaly log may be reduced to facilitate subsequent processing.


At step 540, the segments of the original multi-channel signal corresponding to the detected anomalies are identified.


At step 550, the identified segments of the multi-channel signal may be stored in an output record for post-processing. For example, the output record may be a computer file stored in a storage medium such as a computer hard drive, or it may be a paper print-out. In an exemplary embodiment, the identified segments output to the file may be expanded beyond those strictly associated with the anomalies. For example, fixed time segments of the multi-channel signal both immediately prior to and immediately subsequent to each identified data anomaly may also be output for each identified segment corresponding to an anomaly. The additional segments may further aid in the post-processing of the anomalies in the multi-channel signal, as further described hereinbelow.


In an exemplary embodiment (not shown), the output record generated by the method 500 may be manually reviewed, or “scrubbed,” by a human technician to verify the presence of anomalies in the identified multi-channel signal segments. If the identified segment is verified to contain an anomaly, the segment may be, e.g., omitted from further post-processing, or other measures may be taken.



FIG. 6 depicts an alternative exemplary embodiment of the present disclosure, wherein the techniques disclosed hereinabove are applied in the context of a real-time patient monitoring and neurological event detection system 60. For a more detailed description of the system in FIG. 6, see, e.g., “Minimally Invasive Monitoring Methods,” U.S. patent application Ser. No. 11/766,751, filed Jun. 21, 2007, published as U.S. Patent Publication No. 2008/0027347A1 and assigned to the assignee of the present application, the contents of which are hereby incorporated by reference in their entirety. Note that FIG. 6 is provided for illustrative purposes only, and is not meant to limit the scope of the present disclosure in any way.


In FIG. 6, system 60 includes one or more implantable sensors or devices 62 that are configured to sample electrical activity from the patient's brain (e.g., EEG signals). The implantable devices may be active (with internal power source), passive (no internal power source), or semi-passive (internal power source to power components, but not to transmit data signal). The implantable devices 62 may be implanted anywhere in the patient. In an exemplary embodiment, one or more of the devices 62 may be implanted adjacent to a previously identified epileptic focus or a portion of the brain where the focus is believed to be located. Alternatively, the devices 62 themselves may be used to help determine the location of an epileptic focus


In one aspect, the neural signals of the patient are sampled substantially continuously with the electrodes coupled to the electronic components of the implanted leadless device. A wireless signal is transmitted that is encoded with data that is indicative of the sampled neural signal from the implanted device to an external device. The wireless signal that is encoded with data that is indicative of the sampled neural signal is derived from the wireless signal received from the external device. The wireless signal can be any type of wireless signal—radiofrequency signal, magnetic signal, optical signal, acoustic signal, infrared signal, or the like.


The physician may implant any desired number of devices in the patient. As noted above, in addition to monitoring brain signals, one or more additional implanted devices 62 may be implanted to measure other physiological signals from the patient.


Implantable devices 62 may be configured to substantially continuously sample the brain activity of the groups of neurons in the immediate vicinity of the implanted device. The implantable devices 62 may be interrogated and powered by a signal from an external device 64 to facilitate the substantially continuous sampling of the brain activity signals. Sampling of the brain activity may be carried out at a rate above about 200 Hz, and preferably between about 200 Hz and about 1000 Hz, and most preferably at about 400 Hz, but it could be higher or lower, depending on the specific condition being monitored, the patient, and other factors. Each sample of the patient's brain activity may contain between about 8 bits per sample and about 32 bits per sample, and preferably between about 12 bits per sample and about 16 bits per sample.


In alternative embodiments, it may be desirable to have the implantable devices sample the brain activity of the patient on a non-continuous basis. In such embodiments, the implantable devices 62 may be configured to sample the brain activity signals periodically (e.g., once every 10 seconds) or aperiodically.


Implantable devices 16 may comprise a separate memory module for storing the recorded brain activity signals, a unique identification code for the device, algorithms, other programming, or the like.


A patient instrumented with the implanted devices 62 may carry a data collection device 64 that is external to the patient's body. The external device 64 would receive and store the signals from the implanted devices 62 with the encoded EEG data (or other physiological signals). The signals received from the plurality of implanted devices 62 may be represented as a multi-channel signal, and may be pre-processed according to the techniques of the present disclosure. The external device 64 is typically of a size so as to be portable and carried by the patient in a pocket or bag that is maintained in close proximity to the patient. In alternative embodiments, the device may be configured to be used in a hospital setting and placed alongside a patient's bed. Communication between the data collection device 64 and the implantable device 62 may take place through wireless communication. The wireless communication link between implantable device 62 and external device 64 may provide a communication link for transmitting data and/or power. External device 64 may include a control module 66 that communicates with the implanted device through an antenna 68. In the illustrated embodiment, antenna 68 is in the form of a necklace that is in communication range with the implantable devices 62.


Transmission of data and power between implantable device 62 and external device 64 may be carried out through a radiofrequency link, infrared link, magnetic induction, electromagnetic link, Bluetooth® link, Zigbee link, sonic link, optical link, other types of wireless links, or combinations thereof.


In an exemplary embodiment, the external device 64 may include software to pre-process the data according to the present disclosure and analyze the data in substantially real-time. For example, the received RF signal with the sampled EEG may be analyzed for the presence of anomalies according to the present disclosure, and further by EEG analysis algorithms to estimate the patient's brain state which is typically indicative of the patient's propensity for a neurological event. The neurological event may be a seizure, migraine headache, episode of depression, tremor, or the like. The estimation of the patient's brain state may cause generation of an output. The output may be in the form of a control signal to activate a therapeutic device (e.g., implanted in the patient, such as a vagus nerve stimulator, deep brain or cortical stimulator, implanted drug pump, etc.).


In an exemplary embodiment, the output may be used to activate a user interface on the external device to produce an output communication to the patient. For example, the external device may be used to provide a substantially continuous output or periodic output communication to the patient that indicates their brain state and/or propensity for the neurological event. Such a communication could allow the patient to manually initiate self-therapy (e.g., wave wand over implanted vagus nerve stimulator, cortical, or deep brain stimulator, take a fast acting anti-epileptic drug, etc.).


In an alternative exemplary embodiment, the external device 64 may further communicate with an auxiliary server (not shown) having more extensive computational and storage resources than can be supported in the form factor of the external device 64. In such an exemplary embodiment, the anomaly pre-processing and EEG analysis algorithms may be performed by an auxiliary server, or the computations of the external device 64 may be otherwise facilitated by the computational resources of the auxiliary server.



FIG. 7 depicts a generalized block diagram 700 of the real-time analysis system 60 featuring an anomaly pre-processing block 735 according to the present disclosure. One of ordinary skill in the art will appreciate that the block diagram 700 need not be limited to the exemplary system 60 shown in FIG. 6, but may also be broadly applicable to other types of multi-channel sensing and processing systems.


In FIG. 7, biosensor signals 710 form a multi-channel signal that is provided over a wireless link to wireless unit 720. Wireless unit 720 communicates the multi-channel signal to a processing module 730 that may be resident either on the wireless unit 720 itself, or separately from the wireless unit 720, as described with reference to FIG. 6. When residing separately from the wireless unit 720, the processing module 730 may be configured to run algorithms, perform computations, or perform anomaly checking that may be too complex or intensive for a low-power wireless unit 720 to implement. Such anomaly checking may correspond to the “self-checking” techniques as further described hereinbelow with reference to FIG. 8.


Processing module 730 includes a pre-processing block 735 that identifies and processes anomalies in the multi-channel signal. The output of pre-processing block 735 is provided to a data analysis block 737, which may output an event indicator 730a. The output event indicator 730a may correspond to the output of the estimation of the patient's brain state as described with reference to FIG. 6.


In the exemplary embodiment shown, the pre-processing block 735 communicates with an anomaly data service 740. The anomaly data service 740 may reside remotely from the processing module 730, and may provide the pre-processing block 735 with dynamically adjusted thresholds and/or other parameters to aid the pre-processing block 735 in identifying anomalies in the multi-channel signal. For example, the anomaly data service 740 may analyze anomalies from a plurality of multi-channel signals sampled over a population of seizure detection systems, seizure prediction systems, and/or seizure counter-prediction systems. The anomaly data service 740 may periodically derive preferred DCN comparison thresholds for use in the individual real-time analysis system 700. In an exemplary embodiment, the real-time analysis system 700 may also upload data samples to the anomaly data service 740 to aid the anomaly data service 740 in deriving preferred thresholds.


In an exemplary embodiment, the anomaly data service 740 may communicate with the processing module 730 wirelessly. Alternatively, the anomaly data service 740 may communicate with the processing module 730 over a wired connection. In yet another exemplary embodiment, the anomaly data service 740 may be omitted altogether, and the pre-processing block 735 may simply rely on pre-programmed threshold values. Such exemplary embodiments are contemplated to be within the scope of the present disclosure.



FIG. 8 depicts an exemplary method 800 according to the present disclosure that may be implemented, e.g., using the system 700. Note the exemplary method is shown for illustrative purposes only, and is not meant to limit the scope of the present disclosure to any particular method disclosed.


At step 810, the wireless unit 720 and processing module 730 are powered on.


At step 820, the wireless unit 720 receives the multi-channel signal from, e.g., a plurality of biosensors such as implanted devices 62 in FIG. 6.


At step 830, anomaly pre-processor 735 in processing module 730 identifies the presence of boot-up anomalies in the multi-channel signal received at step 820. This step may also be termed “self-checking,” or “self-test diagnostics.”


In an exemplary embodiment, “boot-up” anomalies may be any anomalies identified in the multi-channel signal during an initial boot-up phase. The boot-up phase may correspond to a time when software in the processing module 730 is initialized, and/or other parameters of the system 700 are initially configured. For example, the boot-up phase may last for a fixed amount of time after the wireless unit 720 and processing module 730 are powered on at step 800.


In an exemplary embodiment, the identification of anomalies in the multi-channel signal may be performed using the DCN computation techniques earlier described herein with reference to FIG. 3A. However, anomaly identification need not be performed using DCN computation. Anomaly identification in the method 800 may generally be performed using any suitable anomaly detection metric or metrics derivable by one of ordinary skill in the art in light of the present disclosure. Such alternative exemplary embodiments are contemplated to be within the scope of the present disclosure.


At step 840, an anomaly central data service may be continuously updated during operation of the method 800 with appropriate thresholds and/or algorithms for detecting the presence of anomalies in the multi-channel signal. In an exemplary embodiment, the anomaly data service may update a series of thresholds T1, T2, etc., against which the data condition number (DCN) is compared to detect the presence of anomalies in the multi-channel signal. The anomaly data service may vary the value of such thresholds over time, based on, e.g., offline analysis of anomalies and associated anomaly metrics as computed over an entire population of multi-channel signals.


At step 850, operation of the system 700 proceeds with the processing module 730 processing the multi-channel signal, taking into account the information in the anomaly data service.


At step 860, the anomaly processor 735 checks for anomalies in the multi-channel signal during normal operation of the system 700. The checking at step 860 may be termed “adaptive” anomaly identification and processing, as contrasted with the “boot-up” anomaly identification and processing described with reference to step 830. Information about anomalies identified during step 860 may be used to update the anomaly data service, as illustrated by the return arrow from step 860 to step 840, and as earlier described with reference to block 740 hereinabove. The steps 840, 850, 860 may be continuously repeated during normal operation of the system 700. An advantage of the adaptive anomaly identification and processing techniques described herein is that they may be varied over an extended temporal context used to monitor the multi-channel signal, as compared to the one-time self-checking diagnostics provided during a boot-up phase.


In an exemplary embodiment, entries from the anomaly data service may also be removed from the data service if anomaly processor 735 determines that such anomalies are no longer applicable. Such exemplary embodiments are contemplated to be within the scope of the present disclosure.



FIG. 9A illustrates an exemplary plot of data condition number (DCN) versus time, corresponding to an operating scenario wherein two of the sensors used to generate a multi-channel signal are gradually electrically shorted together during the interval starting from Time=5 minutes to Time=10 minutes. The scenario shown may arise, e.g., when two adjacent connectors at a distal end of a lead assembly implanted in a patient are gradually shorted together, effectively forming a path for interconnect ionic migration.


As seen in FIG. 9A, the DCN increases to extremely high levels in response to the gradual electrical shorting, and takes on values as high as 1015 when two of the leads are eventually completely shorted together, and converging to their common average, after Time=10 minutes. This expected behavior of the DCN in response to a physical condition such as electrical shorting allows the DCN to be used as an effective detection metric for indicating when such conditions might arise in an actual scenario. For example, by detecting when the DCN exceeds a suitably high threshold, the error condition wherein there is electrode shortage may be flagged to an automated software module, or to a manual user. One of ordinary skill in the art will appreciate that alternative techniques may also be employed to detect such an error condition, e.g., pattern matching the observed DCN to expected behavior of the DCN in the presence of the error condition, such expected behavior being obtained through, e.g., prior simulation or historical data. Furthermore, through further processing, e.g., by computing and evaluating the DCN over specific subsets of channels of the multi-channel signal, the specific leads causing the electrical shortage problem may be specifically identified and dealt with.



FIG. 9B illustrates an exemplary plot of data condition number (DCN) versus time, corresponding to an operating scenario wherein there is a gradual loss of contact of a reference electrode used to sample the multi-channel signal, over the span of about three days. As seen in FIG. 9B, there is a noticeable upward trend in the DCN over the time interval shown. In an exemplary embodiment, the presence of such an upward trend in the DCN may be used to detect the condition corresponding to gradual loss of contact of an electrode. Such a condition may then be flagged to an automated software module, or to a manual user. One of ordinary skill in the art will appreciate that various techniques may be employed to quantify any upward trend in the DCN. For example, a low-pass filter may be applied to the DCN time-series, and absolute increases in the DCN over a suitably long time period may be assessed and compared to some threshold. Alternatively, other pattern matching techniques may also be employed. Such alternative exemplary embodiments are contemplated to be within the scope of the present disclosure.



FIG. 10A illustrates an exemplary embodiment of a graphical display interface 1000 for simultaneously displaying various types of information associated with the techniques of the present disclosure. Note the exemplary embodiment of FIG. 10A is shown for illustrative purposes only, and is not meant to limit the scope of the present disclosure to any particular graphical display interface shown. One of ordinary skill in the art will appreciate that the graphical display interface 1000 may be a “screen shot,” i.e., an instantaneous sample of the image on a computer screen displayed when a user (e.g., a technician) is using the disclosed graphical interface to analyze data. The graphical display may be generated by, e.g., computer software running on a general-purpose computer, and may be readily designed by one of ordinary skill in the art in light of the further description given hereinbelow. For example, in an exemplary embodiment, such computer software may be conveniently accessed via a Web browser.


In FIG. 10A, graphical display interface 1000 includes a graphical representation of a catalog configuration block 1010, a graphical representation of a legend 1020, and a series of time-series displays 1040. The time-series displays 1040 may include time series 1041, 1042, 1043, each associated with a corresponding subject identification number as denoted in the subject ID column 1030. In the exemplary embodiment shown, there are three subject ID's 0001, 0002, and 0003. One of ordinary skill in the art will appreciate that the techniques disclosed herein may readily be applied to display information associated with any number of subjects, per the requirements of the application.


The catalog configuration block 1010 includes a series of drop-down menus 1011, 1012, and 1013.


Drop-down menu 1012 allows the user to select a scheme for backdrop patterns used to highlight certain information contained in the time-series displays 1040. In an exemplary embodiment, the mapping between a type of information and a corresponding backdrop pattern used to highlight such information type may be as shown according to the legend 1020. The legend 1020 indicates backdrop patterns used to show an ‘interictal’ segment 1021, a ‘preseizure’ segment 1022 (e.g., a fixed preictal time interval), other segment 1023, dropout 1024 (e.g., lapsed or invalid data), and seizure event 1025. As shown in the time-series 1041, 1042, and 1043, the backdrop patterns identified in legend 1020 may be used to annotate the time-series displays 1040 with information in a concise and easily presented form for convenient analysis by a user of the graphical display interface 1000.


One of ordinary skill in the art will appreciate that various modifications to the backdrop pattern scheme shown are readily derivable in light of the present disclosure. For example, in a color graphical display interface (not shown), backdrop colors may be used in place of, or in addition to, the backdrop patterns shown. A color scheme may assign a distinct color to each of the following types of segments: “interictal” segments, “preseizure” segments, “other” segments, “dropouts” (e.g., lapsed or invalid data), and “seizure” events. Alternative exemplary embodiments may also include other types of information not explicitly shown in FIG. 10A. Such alternative exemplary embodiments are contemplated to be within the scope of the present disclosure.


Drop-down menu 1013 allows the user to select the quantity to be displayed on the vertical axes of the time-series displays 1040. For example, one selectable quantity may be a data condition number (DCN) as calculated for each subject according to the present disclosure, plotted versus time. Alternative types of quantities include, e.g., the channel-sum of line-lengths (sum of absolute deviations in time) or its spatiotemporal and normalized variants, which would be appropriate for large-scale temporal localization of seizures rather than EEG anomalies.


As further shown in FIG. 10A, time series 1041 is accompanied by a time event index 1041.1. In the exemplary embodiment shown, the time event index 1041.1 is shown in bold and underlined, and has a numerical value of 20108. The time event index 1041.1 may indicate a time index in the time series 1041 wherein the corresponding DCN value exceeds some pre-specified threshold. See, e.g., the portion of time series 1041 indicated by 1041.1a. As more clearly shown in FIG. 10B, by rolling a cursor 1080A over the time event index 1041.1, an auxiliary window 1090A may “pop up” in the same display. Alternatively, the time event index 1041.1 may be “hyperlinked,” such that placing the cursor 1080A over the time event index 1041.1 and clicking the index 1041.1 causes the auxiliary window 1090A to simultaneously pop up.


The information content, or type of event, shown in the auxiliary window 1090A may be selected via the drop-down menu 1011. In the exemplary embodiment shown in FIG. 10B, the auxiliary window 1090A may show, e.g., the portion of the original multi-channel signal (such as EEG data) from which the DCN value corresponding to the time event index 1041.1 was derived. The auxiliary window 1090A advantageously allows the graphical interface user to visually inspect the portion of the multi-channel signal giving rise to the above-threshold DCN value alongside the DCN time-series, and manually determine whether the identified portion is in fact anomalous.


Other pop-up events to display may include context menus that direct the user to the raw EEG data where an anomaly can be further inspected in full detail.


As shown for indices 1043.1 and 1043.2 associated with time-series 1043, there may be multiple points in a time series wherein a DCN value exceeds a pre-specified threshold, and thus there may generally be multiple event indices associated with each time-series.


Based on the teachings described herein, it should be apparent that an aspect disclosed herein may be implemented independently of any other aspects and that two or more of these aspects may be combined in various ways. In one or more exemplary embodiments, the functions described may be implemented in hardware, software, firmware, or any combination thereof. If implemented in software, the functions may be stored on or transmitted over as one or more instructions or code on a computer-readable medium. Computer-readable media includes both computer storage media and communication media including any medium that facilitates transfer of a computer program from one place to another. A storage media may be any available media that can be accessed by a computer. By way of example, and not limitation, such computer-readable media can comprise RAM, ROM, EEPROM, CD/DVD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, solid-state flash cards or drives, or any other medium that can be used to carry or store desired program code in the form of instructions or data structures and that can be accessed by a computer. Also, any connection is properly termed a computer-readable medium. For example, if the software is transmitted from a website, server, or other remote source using a coaxial cable, fiber optic cable, twisted pair, digital subscriber line (DSL), or wireless technologies such as infrared, radio, and microwave, then the coaxial cable, fiber optic cable, twisted pair, DSL, or wireless technologies such as infrared, radio, and microwave are included in the definition of medium. Disk and disc, as used herein, includes compact disc (CD), laser disc, optical disc, digital versatile disc (DVD), floppy disk and Blu-Ray disc where disks usually reproduce data magnetically, while discs reproduce data optically with lasers. Combinations of the above should also be included within the scope of computer-readable media.


In this specification and in the claims, it will be understood that when an element is referred to as being “connected to” or “coupled to” another element, it can be directly connected or coupled to the other element or intervening elements may be present. In contrast, when an element is referred to as being “directly connected to” or “directly coupled to” another element, there are no intervening elements present.


A number of aspects and examples have been described. However, various modifications to these examples are possible, and the principles presented herein may be applied to other aspects as well. These and other aspects are within the scope of the following claims.

Claims
  • 1. A method for displaying information to a user, the method comprising: accepting input from the user selecting a metric to be displayed, said selected metric being associated with a multi-channel signal derived from outputs of a plurality of sensors positioned to acquire physiological signals, wherein each channel of the multi-channel signal corresponds to a different one of the plurality of sensors;displaying with a processor and a computer graphical display interface at least one time-series plot of the selected metric associated with the multi-channel signal;using a backdrop pattern, indicating on the at least one time-series plot portions of the plot having at least one identified characteristic;for each of the at least one time-series plot displayed, displaying with the processor and the computer graphical display interface a time event index wherein the corresponding metric meets a predetermined condition, wherein the predetermined condition comprises the selected metric being greater than a predetermined threshold;accepting input from the user as to whether to display further information associated with the time event index; anddisplaying the further information associated with the time event index when the user so specifies.
  • 2. The method of claim 1, the metric comprising a data condition number calculated for the multi-channel signal.
  • 3. The method of claim 1, the at least one time-series plot comprising a plurality of time-series plots, each time-series plot associated with a subject ID.
  • 4. The method of claim 1, the at least one identified characteristic including a portion of the time-series corresponding to a seizure segment.
  • 5. The method of claim 1, the backdrop pattern comprising a color scheme, the color scheme representing “interictal” segments, “preseizure” segments, “other” segments, “dropouts,” and “seizure” events.
  • 6. The method of claim 1, the input from the user as to whether to display further information associated with the time event index comprising positioning of a cursor over the time event index.
  • 7. The method of claim 1, further comprising: accepting input from the user indicating a color scheme for the backdrop pattern to be used.
  • 8. The method of claim 1, the displaying the further information associated with the time event index comprising: displaying a pop-up plot of a portion of the multi-channel signal corresponding to the time event index.
  • 9. The method of claim 8, the multi-channel signal comprising a multi-channel EEG signal.
  • 10. A method for detecting anomalies in a multi-channel signal, the method comprising: sampling the multi-channel signal over a time window;computing with a processor an anomaly metric for the multi-channel signal over the time window; andidentifying the presence of an anomaly based on the magnitude of the anomaly metric; the computing an anomaly metric comprising: computing a condition number of the multi-channel signal over the time window; andadjusting the condition number based on a parameter of the multi-channel signal to generate a data condition number (DCN); the multi-channel signal comprising a signal sampled from a plurality of electrodes positioned to acquire physiological signals from a body, the method further comprising:generating a DCN time series corresponding to a plurality of time windows;generating an anomaly log based on the DCN time series;merging anomalies in the anomaly log separated by less than a minimum separation to generate a modified anomaly log;identifying segments of the multi-channel signal corresponding to anomalies in the modified anomaly log; andoutputting time-expanded versions of the identified segments to a record; the identifying the presence of an anomaly comprising matching the DCN time series to at least one known pattern of DCN time series corresponding to an anomalous condition.
  • 11. A method for detecting anomalies in a multi-channel signal, the method comprising: sampling the multi-channel signal over a time window;computing with a processor an anomaly metric for the multi-channel signal over the time window; andidentifying the presence of an anomaly by comparing the magnitude of the anomaly metric to an optimum threshold; the computing an anomaly metric comprising:computing a condition number of the multi-channel signal over the time window; andadjusting the condition number based on a parameter of the multi-channel signal to generate a data condition number (DCN); the method further comprising:generating a DCN time series corresponding to a plurality of time windows;generating an event log based on the DCN time series, the generating an event log comprising identifying an event as a contiguous set of DCN values greater than a candidate threshold;merging events in the event log separated by less than a minimum separation to generate a modified event log;repeating the steps of generating a DCN time series, generating an event log, and merging events using a plurality of candidate thresholds;generating a plot of number of events identified in a modified event log versus candidate threshold used;attempting to identify at least one inflection point in the generated plot; andsetting the optimum threshold to be the candidate threshold corresponding to the inflection point with the largest abscissa in the plot.
  • 12. The method of claim 11, further comprising, when the attempt to identify at least one inflection point is unsuccessful: identifying a maximum number of events corresponding to a candidate threshold in the histogram; andsetting the optimum threshold to be a candidate threshold whose corresponding number of events is less than or equal to a fixed percentage of the maximum number of events.
CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of U.S. Provisional Application No. 61/183,449, filed Jun. 2, 2009, which is incorporated herein by reference in its entirety.

US Referenced Citations (523)
Number Name Date Kind
3218638 Honig Nov 1965 A
3498287 Ertl Mar 1970 A
3522811 Schwartz Aug 1970 A
3575162 Gaarder Apr 1971 A
3837331 Ross Sep 1974 A
3850161 Liss Nov 1974 A
3863625 Viglione et al. Feb 1975 A
3882850 Bailin et al. May 1975 A
3918461 Cooper Nov 1975 A
3967616 Ross Jul 1976 A
3993046 Fernandez Nov 1976 A
4201224 John May 1980 A
4214591 Sato et al. Jul 1980 A
4279258 John Jul 1981 A
4305402 Katims Dec 1981 A
4334545 Shiga Jun 1982 A
4407299 Culver Oct 1983 A
4408616 Duffy et al. Oct 1983 A
4421122 Duffy Dec 1983 A
4471786 Inagaki Sep 1984 A
4494950 Fischell Jan 1985 A
4505275 Chen Mar 1985 A
4545388 John Oct 1985 A
4556061 Barreras et al. Dec 1985 A
4566464 Piccone et al. Jan 1986 A
4573481 Bullara Mar 1986 A
4579125 Strobl et al. Apr 1986 A
4590946 Loeb May 1986 A
4612934 Borkan Sep 1986 A
4679144 Cox et al. Jul 1987 A
4686999 Snyder et al. Aug 1987 A
4702254 Zabara Oct 1987 A
4735208 Wyler et al. Apr 1988 A
4768176 Kehr et al. Aug 1988 A
4768177 Kehr et al. Aug 1988 A
4785827 Fischer Nov 1988 A
4793353 Borkam Dec 1988 A
4817628 Zealear Apr 1989 A
4838272 Lieber Jun 1989 A
4844075 Liss et al. Jul 1989 A
4852573 Kennedy Aug 1989 A
4867164 Zabara Sep 1989 A
4873981 Abrams et al. Oct 1989 A
4878498 Abrams et al. Nov 1989 A
4903702 Putz Feb 1990 A
4920979 Bullara May 1990 A
4926865 Oman May 1990 A
4955380 Edell Sep 1990 A
4978680 Sofia Dec 1990 A
4979511 Terry Dec 1990 A
4991582 Byers et al. Feb 1991 A
5010891 Chamoun Apr 1991 A
5016635 Graupe May 1991 A
5025807 Zabara Jun 1991 A
5031618 Mullett Jul 1991 A
5070873 Graupe et al. Dec 1991 A
5082861 Sofia Jan 1992 A
5097835 Putz Mar 1992 A
RE34015 Duffy Aug 1992 E
5154172 Terry Oct 1992 A
5167229 Peckham et al. Dec 1992 A
5179950 Stanislaw Jan 1993 A
5181520 Wertheim et al. Jan 1993 A
5186170 Varichio Feb 1993 A
5188104 Wernicke Feb 1993 A
5190029 Byron et al. Mar 1993 A
5193539 Schulman et al. Mar 1993 A
5193540 Schulman et al. Mar 1993 A
5205285 Baker, Jr. Apr 1993 A
5215086 Terry, Jr. et al. Jun 1993 A
5215088 Normann Jun 1993 A
5215089 Baker, Jr. Jun 1993 A
5222494 Baker, Jr. Jun 1993 A
5222503 Ives Jun 1993 A
5231988 Wernicke et al. Aug 1993 A
5235980 Varichio et al. Aug 1993 A
5237991 Baker, Jr. Aug 1993 A
5251634 Weinberg Oct 1993 A
5263480 Wernicke et al. Nov 1993 A
5265619 Comby et al. Nov 1993 A
5269302 Swartz et al. Dec 1993 A
5269303 Wernicke et al. Dec 1993 A
5269315 Leuchter et al. Dec 1993 A
5292772 Sofia Mar 1994 A
5293879 Vonk Mar 1994 A
5299118 Martens et al. Mar 1994 A
5299569 Wernicke et al. Apr 1994 A
5300094 Kallok et al. Apr 1994 A
5304206 Baker, Jr. et al. Apr 1994 A
5311876 Olsen et al. May 1994 A
5312439 Loeb May 1994 A
5314458 Najafi et al. May 1994 A
5324316 Schulman et al. Jun 1994 A
5330515 Rutecki et al. Jul 1994 A
5335657 Terry et al. Aug 1994 A
5342408 deCoriolis et al. Aug 1994 A
5342409 Mullett Aug 1994 A
5343064 Spangler et al. Aug 1994 A
5349962 Lockard et al. Sep 1994 A
5351394 Weinberg Oct 1994 A
5361760 Normann Nov 1994 A
5365939 Ochs Nov 1994 A
5376359 Johnson Dec 1994 A
5392788 Hudspeth Feb 1995 A
5405365 Hoegnelid et al. Apr 1995 A
5405367 Schulman et al. Apr 1995 A
5411540 Edell et al. May 1995 A
5458117 Chamoun Oct 1995 A
5474547 Aebischer et al. Dec 1995 A
5476494 Edell et al. Dec 1995 A
5486999 Mebane Jan 1996 A
5513649 Gevins May 1996 A
5517115 Prammer May 1996 A
5531778 Maschino et al. Jul 1996 A
5540730 Terry Jul 1996 A
5540734 Zabara Jul 1996 A
5549656 Reiss Aug 1996 A
5555191 Hripcsak Sep 1996 A
5571148 Loeb et al. Nov 1996 A
5571150 Wernicke Nov 1996 A
5575813 Edell et al. Nov 1996 A
5578036 Stone et al. Nov 1996 A
5611350 John Mar 1997 A
5626145 Clapp et al. May 1997 A
5626627 Krystal et al. May 1997 A
5638826 Wolpaw Jun 1997 A
5649068 Boser et al. Jul 1997 A
5672154 Sillen et al. Sep 1997 A
5683422 Rise Nov 1997 A
5683432 Goedeke et al. Nov 1997 A
5690681 Geddes et al. Nov 1997 A
5690691 Chen et al. Nov 1997 A
5697369 Long Dec 1997 A
5700282 Zabara Dec 1997 A
5704352 Tremblay et al. Jan 1998 A
5707400 Terry et al. Jan 1998 A
5711316 Elsberry et al. Jan 1998 A
5713923 Ward et al. Feb 1998 A
5715821 Faupel Feb 1998 A
5716377 Rise et al. Feb 1998 A
5720294 Skinner Feb 1998 A
5735814 Elsberry et al. Apr 1998 A
5743860 Hively et al. Apr 1998 A
5752979 Benabid May 1998 A
5769778 Abrams et al. Jun 1998 A
5776434 Purewal et al. Jul 1998 A
5782798 Rise Jul 1998 A
5782874 Loos Jul 1998 A
5792186 Rise Aug 1998 A
5800474 Benabid et al. Sep 1998 A
5813993 Kaplan Sep 1998 A
5814014 Elsberry et al. Sep 1998 A
5815413 Hively et al. Sep 1998 A
5816247 Maynard Oct 1998 A
5824021 Rise Oct 1998 A
5832932 Elsberry et al. Nov 1998 A
5833709 Rise et al. Nov 1998 A
5857978 Hively et al. Jan 1999 A
5862803 Besson et al. Jan 1999 A
5876424 O'Phelan et al. Mar 1999 A
5899922 Loos May 1999 A
5913881 Benz et al. Jun 1999 A
5916239 Geddes et al. Jun 1999 A
5917429 Otis, Jr. et al. Jun 1999 A
5928272 Adkins Jul 1999 A
5931791 Saltzstein et al. Aug 1999 A
5938689 Fischell et al. Aug 1999 A
5941906 Barreras et al. Aug 1999 A
5950632 Reber et al. Sep 1999 A
5957861 Combs et al. Sep 1999 A
5971594 Sahai et al. Oct 1999 A
5975085 Rise Nov 1999 A
5978702 Ward et al. Nov 1999 A
5978710 Prutchi et al. Nov 1999 A
5995868 Dorfmeister et al. Nov 1999 A
6006124 Fischell et al. Dec 1999 A
6016449 Fischell et al. Jan 2000 A
6018682 Rise Jan 2000 A
6042548 Giuffre Mar 2000 A
6042579 Elsberry et al. Mar 2000 A
6051017 Loeb et al. Apr 2000 A
6052619 John Apr 2000 A
6061593 Fischell et al. May 2000 A
6066163 John May 2000 A
6081744 Loos Jun 2000 A
6094598 Elsberry et al. Jul 2000 A
6109269 Rise et al. Aug 2000 A
6117066 Abrams et al. Sep 2000 A
6128537 Rise et al. Oct 2000 A
6128538 Fischell et al. Oct 2000 A
6134474 Fischell et al. Oct 2000 A
6161045 Fischell et al. Dec 2000 A
6167304 Loos Dec 2000 A
6171239 Humphrey Jan 2001 B1
6176242 Rise Jan 2001 B1
6205359 Boveja Mar 2001 B1
6208893 Hofmann Mar 2001 B1
6221011 Bardy Apr 2001 B1
6227203 Rise et al. May 2001 B1
6230049 Fischell et al. May 2001 B1
6248126 Lesser et al. Jun 2001 B1
6249703 Stanton Jun 2001 B1
6263237 Rise Jul 2001 B1
6280198 Calhoun et al. Aug 2001 B1
6304775 Iasemidis et al. Oct 2001 B1
6309406 Jones et al. Oct 2001 B1
6328699 Eigler Dec 2001 B1
6337997 Rise Jan 2002 B1
6339725 Naritoku Jan 2002 B1
6341236 Osorio et al. Jan 2002 B1
6343226 Sunde et al. Jan 2002 B1
6353754 Fischell et al. Mar 2002 B1
6354299 Fischell et al. Mar 2002 B1
6356784 Lozano et al. Mar 2002 B1
6356788 Boveja Mar 2002 B2
6358203 Bardy Mar 2002 B2
6358281 Berrang et al. Mar 2002 B1
6360122 Fischell Mar 2002 B1
6366813 DiLorenzo Apr 2002 B1
6366814 Boveja Apr 2002 B1
6374140 Rise Apr 2002 B1
6386882 Linberg May 2002 B1
6402678 Fischell et al. Jun 2002 B1
6411854 Tziviskos et al. Jun 2002 B1
6427086 Fischell et al. Jul 2002 B1
6434419 Gevins et al. Aug 2002 B1
6442421 Le Van Quyen et al. Aug 2002 B1
6443891 Grevious Sep 2002 B1
6453198 Torgerson Sep 2002 B1
6463328 John Oct 2002 B1
6466822 Pless Oct 2002 B1
6471645 Warkentin et al. Oct 2002 B1
6473639 Fischell et al. Oct 2002 B1
6473644 Terry et al. Oct 2002 B1
6480743 Kirkpatrick Nov 2002 B1
6484132 Hively et al. Nov 2002 B1
6488617 Katz Dec 2002 B1
6496724 Levendowski et al. Dec 2002 B1
6505077 Kast et al. Jan 2003 B1
6510340 Jordan Jan 2003 B1
6511424 Moore-Ede Jan 2003 B1
6529774 Greene Mar 2003 B1
6534693 Fischell et al. Mar 2003 B2
6547746 Marino Apr 2003 B1
6549804 Osorio et al. Apr 2003 B1
6553262 Lang et al. Apr 2003 B1
6560486 Osorio et al. May 2003 B1
6571123 Ives et al. May 2003 B2
6571125 Thompson May 2003 B2
6572528 Rohan et al. Jun 2003 B2
6587719 Barrett et al. Jul 2003 B1
6591132 Gotman et al. Jul 2003 B2
6591137 Fischell et al. Jul 2003 B1
6591138 Fischell et al. Jul 2003 B1
6597954 Pless et al. Jul 2003 B1
6600956 Maschino Jul 2003 B2
6609025 Barrett et al. Aug 2003 B2
6618623 Pless et al. Sep 2003 B1
6620415 Donovan Sep 2003 B2
6622036 Suffin Sep 2003 B1
6622038 Barrett et al. Sep 2003 B2
6622041 Terry et al. Sep 2003 B2
6622047 Barrett et al. Sep 2003 B2
6658287 Litt et al. Dec 2003 B1
6665562 Gluckman et al. Dec 2003 B2
6668191 Boveja Dec 2003 B1
6671555 Gielen Dec 2003 B2
6678548 Echauz et al. Jan 2004 B1
6684105 Cohen et al. Jan 2004 B2
6687538 Hrdlicka et al. Feb 2004 B1
6735467 Wilson May 2004 B2
6760626 Boveja Jul 2004 B1
6768969 Nikitin et al. Jul 2004 B1
6778854 Puskas Aug 2004 B2
6782292 Whitehurst Aug 2004 B2
6819956 DiLorenzo Nov 2004 B2
6879859 Boveja Apr 2005 B1
6893395 Kraus et al. May 2005 B1
6901294 Whitehurst et al. May 2005 B1
6901296 Whitehurst et al. May 2005 B1
6912419 Hill Jun 2005 B2
6921538 Donovan Jul 2005 B2
6921541 Chasin et al. Jul 2005 B2
6923784 Stein Aug 2005 B2
6931274 Williams Aug 2005 B2
6934580 Osorio Aug 2005 B1
6937891 Leinders et al. Aug 2005 B2
6944501 Pless Sep 2005 B1
6950706 Rodriguez Sep 2005 B2
6973342 Swanson Dec 2005 B1
6990372 Perron et al. Jan 2006 B2
7010351 Firlik et al. Mar 2006 B2
7089059 Pless Aug 2006 B1
7174212 Klehn et al. Feb 2007 B1
7177701 Pianca Feb 2007 B1
7185283 Takahashi Feb 2007 B1
7209787 DiLorenzo Apr 2007 B2
7212851 Donoghue et al. May 2007 B2
7231254 DiLorenzo Jun 2007 B2
7242984 DiLorenzo Jul 2007 B2
7277758 DiLorenzo Oct 2007 B2
7324851 DiLorenzo Jan 2008 B1
7373198 Bibian et al. May 2008 B2
7403820 DiLorenzo Jul 2008 B2
7463917 Martinez Dec 2008 B2
7623928 DiLorenzo Nov 2009 B2
7631015 Gupta et al. Dec 2009 B2
7747325 DiLorenzo Jun 2010 B2
7805196 Miesel et al. Sep 2010 B2
7853329 DiLorenzo Dec 2010 B2
7881798 Miesel et al. Feb 2011 B2
8036736 Snyder et al. Oct 2011 B2
8055348 Heruth et al. Nov 2011 B2
8125484 Gering Feb 2012 B2
20010051819 Fischell et al. Dec 2001 A1
20010056290 Fischell et al. Dec 2001 A1
20020002390 Fischell et al. Jan 2002 A1
20020035338 Dear et al. Mar 2002 A1
20020054694 Vachtsevanos et al. May 2002 A1
20020072770 Pless Jun 2002 A1
20020072776 Osorio et al. Jun 2002 A1
20020072782 Osorio et al. Jun 2002 A1
20020077670 Archer et al. Jun 2002 A1
20020095099 Quyen et al. Jul 2002 A1
20020099412 Fischell et al. Jul 2002 A1
20020103512 Echauz et al. Aug 2002 A1
20020109621 Khalr et al. Aug 2002 A1
20020111542 Warkentin et al. Aug 2002 A1
20020116042 Boling Aug 2002 A1
20020126036 Flaherty et al. Sep 2002 A1
20020147388 Mass et al. Oct 2002 A1
20020169485 Pless et al. Nov 2002 A1
20030004428 Pless Jan 2003 A1
20030009207 Paspa et al. Jan 2003 A1
20030013981 Gevins et al. Jan 2003 A1
20030018367 DiLorenzo Jan 2003 A1
20030028072 Fischell et al. Feb 2003 A1
20030050549 Sochor Mar 2003 A1
20030050730 Greeven et al. Mar 2003 A1
20030073917 Echauz et al. Apr 2003 A1
20030074033 Pless et al. Apr 2003 A1
20030083716 Nicolelis et al. May 2003 A1
20030114886 Gluckman et al. Jun 2003 A1
20030144709 Zabara et al. Jul 2003 A1
20030144711 Pless et al. Jul 2003 A1
20030144829 Geatz et al. Jul 2003 A1
20030149457 Tcheng et al. Aug 2003 A1
20030158587 Esteller et al. Aug 2003 A1
20030167078 Weisner et al. Sep 2003 A1
20030174554 Dunstone et al. Sep 2003 A1
20030176806 Pineda et al. Sep 2003 A1
20030187621 Nikitin et al. Oct 2003 A1
20030195574 Osorio et al. Oct 2003 A1
20030195588 Fischell et al. Oct 2003 A1
20030195602 Boling Oct 2003 A1
20040034368 Pless et al. Feb 2004 A1
20040039427 Barrett et al. Feb 2004 A1
20040039981 Riedl et al. Feb 2004 A1
20040054297 Wingeier et al. Mar 2004 A1
20040059761 Hively Mar 2004 A1
20040068199 Echauz et al. Apr 2004 A1
20040073273 Gluckman et al. Apr 2004 A1
20040077995 Ferek-Petric Apr 2004 A1
20040078160 Frei et al. Apr 2004 A1
20040082984 Osorio et al. Apr 2004 A1
20040087835 Hively May 2004 A1
20040097802 Cohen May 2004 A1
20040122281 Fischell et al. Jun 2004 A1
20040122335 Sackellares et al. Jun 2004 A1
20040127810 Sackellares et al. Jul 2004 A1
20040133119 Osorio et al. Jul 2004 A1
20040133248 Frei et al. Jul 2004 A1
20040133390 Osorio et al. Jul 2004 A1
20040138516 Osorio et al. Jul 2004 A1
20040138517 Osorio et al. Jul 2004 A1
20040138536 Frei et al. Jul 2004 A1
20040138578 Pineda et al. Jul 2004 A1
20040138579 Deadwyler et al. Jul 2004 A1
20040138580 Frei et al. Jul 2004 A1
20040138581 Frei et al. Jul 2004 A1
20040138647 Osorio et al. Jul 2004 A1
20040138711 Osorio et al. Jul 2004 A1
20040138721 Osorio et al. Jul 2004 A1
20040147969 Mann et al. Jul 2004 A1
20040152958 Frei et al. Aug 2004 A1
20040153129 Pless et al. Aug 2004 A1
20040158119 Osorio et al. Aug 2004 A1
20040169723 Kawasaki et al. Sep 2004 A1
20040172089 Whitehurst et al. Sep 2004 A1
20040176359 Wermeling Sep 2004 A1
20040181263 Balzer et al. Sep 2004 A1
20040199212 Fischell Oct 2004 A1
20040210269 Shalev et al. Oct 2004 A1
20040243146 Chesbrough et al. Dec 2004 A1
20040267152 Pineda et al. Dec 2004 A1
20050004621 Boveja et al. Jan 2005 A1
20050010113 Hall et al. Jan 2005 A1
20050010261 Luders et al. Jan 2005 A1
20050010864 Horikiri et al. Jan 2005 A1
20050015128 Rezai et al. Jan 2005 A1
20050015129 Mische Jan 2005 A1
20050021105 Firlik et al. Jan 2005 A1
20050021108 Klosterman et al. Jan 2005 A1
20050021313 Nikitin et al. Jan 2005 A1
20050027328 Greenstein Feb 2005 A1
20050033369 Badelt Feb 2005 A1
20050043772 Stahmann et al. Feb 2005 A1
20050043774 Devlin et al. Feb 2005 A1
20050049649 Luders et al. Mar 2005 A1
20050059867 Cheng Mar 2005 A1
20050070970 Knudson et al. Mar 2005 A1
20050075067 Lawson et al. Apr 2005 A1
20050096710 Kieval May 2005 A1
20050113885 Haubrich et al. May 2005 A1
20050124863 Cook Jun 2005 A1
20050131493 Boveja et al. Jun 2005 A1
20050137640 Freeberg et al. Jun 2005 A1
20050143786 Boveja Jun 2005 A1
20050143787 Boveja et al. Jun 2005 A1
20050149123 Lesser et al. Jul 2005 A1
20050182308 Bardy Aug 2005 A1
20050182464 Schulte et al. Aug 2005 A1
20050187789 Hatlestad Aug 2005 A1
20050197590 Osorio et al. Sep 2005 A1
20050203366 Donoghue et al. Sep 2005 A1
20050203584 Twetan et al. Sep 2005 A1
20050209218 Meyerson et al. Sep 2005 A1
20050222503 Dunlop et al. Oct 2005 A1
20050222641 Pless Oct 2005 A1
20050228249 Boling Oct 2005 A1
20050228461 Osorio et al. Oct 2005 A1
20050231374 Diem et al. Oct 2005 A1
20050234355 Rowlandson Oct 2005 A1
20050240245 Bange et al. Oct 2005 A1
20050245970 Erickson et al. Nov 2005 A1
20050245971 Brockway et al. Nov 2005 A1
20050245984 Singhal et al. Nov 2005 A1
20050266301 Smith et al. Dec 2005 A1
20050277844 Strother Dec 2005 A1
20060015034 Martinerie et al. Jan 2006 A1
20060015153 Gliner et al. Jan 2006 A1
20060094970 Drew May 2006 A1
20060111644 Guttag et al. May 2006 A1
20060122469 Martel Jun 2006 A1
20060129056 Leuthardt et al. Jun 2006 A1
20060142822 Tulgar Jun 2006 A1
20060173259 Flaherty et al. Aug 2006 A1
20060173510 Besio et al. Aug 2006 A1
20060200038 Savit et al. Sep 2006 A1
20060212092 Pless et al. Sep 2006 A1
20060212093 Pless et al. Sep 2006 A1
20060212096 Stevenson Sep 2006 A1
20060217792 Hussein et al. Sep 2006 A1
20060224191 Dilorenzo Oct 2006 A1
20060253096 Blakley et al. Nov 2006 A1
20060293578 Rennaker, II Dec 2006 A1
20060293720 DiLorenzo et al. Dec 2006 A1
20070027367 Oliver et al. Feb 2007 A1
20070027514 Gerber Feb 2007 A1
20070035910 Stevenson Feb 2007 A1
20070043459 Abbott, III et al. Feb 2007 A1
20070055320 Weinand Mar 2007 A1
20070060973 Ludvig et al. Mar 2007 A1
20070073355 DiLorenzo Mar 2007 A1
20070073357 Rooney et al. Mar 2007 A1
20070100398 Sloan May 2007 A1
20070149952 Bland et al. Jun 2007 A1
20070150024 Leyde et al. Jun 2007 A1
20070150025 DiLorenzo et al. Jun 2007 A1
20070161919 DiLorenzo Jul 2007 A1
20070167991 DiLorenzo Jul 2007 A1
20070185890 VanEpps et al. Aug 2007 A1
20070213629 Greene Sep 2007 A1
20070213785 Osorio et al. Sep 2007 A1
20070217121 Fu et al. Sep 2007 A1
20070238939 Giftakis et al. Oct 2007 A1
20070244407 Osorlo Oct 2007 A1
20070250077 Skakoon et al. Oct 2007 A1
20070250901 McIntire et al. Oct 2007 A1
20070287931 DiLorenzo Dec 2007 A1
20080021341 Harris et al. Jan 2008 A1
20080027347 Harris et al. Jan 2008 A1
20080027348 Harris et al. Jan 2008 A1
20080027515 Harris et al. Jan 2008 A1
20080033502 Harris et al. Feb 2008 A1
20080082019 Ludving et al. Apr 2008 A1
20080091090 Guillory et al. Apr 2008 A1
20080103556 Li et al. May 2008 A1
20080114417 Leyde May 2008 A1
20080119900 DiLorenzo May 2008 A1
20080161712 Leyde Jul 2008 A1
20080161713 Leyde et al. Jul 2008 A1
20080163085 Subbu et al. Jul 2008 A1
20080183096 Snyder et al. Jul 2008 A1
20080183097 Leyde et al. Jul 2008 A1
20080208074 Snyder et al. Aug 2008 A1
20080221876 Holdrich Sep 2008 A1
20080255582 Harris Oct 2008 A1
20080273287 Iyer et al. Nov 2008 A1
20080319281 Aarts Dec 2008 A1
20090018609 DiLorenzo Jan 2009 A1
20090062682 Bland et al. Mar 2009 A1
20090062696 Nathan et al. Mar 2009 A1
20090171168 Leyde et al. Jul 2009 A1
20090171420 Brown et al. Jul 2009 A1
20090264952 Jassemidis et al. Oct 2009 A1
20090290971 Shamseldin et al. Nov 2009 A1
20100023089 DiLorenzo Jan 2010 A1
20100125219 Harris et al. May 2010 A1
20100145176 Himes Jun 2010 A1
20100168603 Himes et al. Jul 2010 A1
20100168604 Echauz et al. Jul 2010 A1
20100179627 Floyd et al. Jul 2010 A1
20100217348 DiLorenzo Aug 2010 A1
20100265072 Goetz et al. Oct 2010 A1
20110166430 Harris et al. Jul 2011 A1
20110172554 Leyde et al. Jul 2011 A1
20110201944 Higgins et al. Aug 2011 A1
20110213222 Leyde et al. Sep 2011 A1
20110218820 Himes et al. Sep 2011 A1
20110219325 Himes et al. Sep 2011 A1
20110260855 John et al. Oct 2011 A1
20110319785 Snyder et al. Dec 2011 A1
Foreign Referenced Citations (84)
Number Date Country
2251852 Apr 1999 CA
2423840 Feb 2002 CA
2428116 May 2002 CA
2428383 May 2002 CA
2425122 Jun 2002 CA
2425004 Aug 2002 CA
2456443 Jan 2003 CA
2491987 Jan 2004 CA
69832022D Dec 2005 DE
0124663 Nov 1984 EP
0898460 Mar 1999 EP
1145735 Oct 2001 EP
1145736 Oct 2001 EP
1307260 May 2003 EP
1335668 Aug 2003 EP
1525551 Apr 2005 EP
0911061 Oct 2005 EP
1609414 Dec 2005 EP
24033673 Feb 2004 JP
1074484 Feb 1984 SU
WO 8501213 Mar 1985 WO
WO 9200119 Jan 1992 WO
WO 9726823 Jul 1997 WO
WO 9734522 Sep 1997 WO
WO 9734524 Sep 1997 WO
WO 9734525 Sep 1997 WO
WO 9739797 Oct 1997 WO
WO 9742990 Nov 1997 WO
WO 9745160 Dec 1997 WO
WO 9849935 Nov 1998 WO
WO 9920342 Apr 1999 WO
WO 9956821 Nov 1999 WO
WO 9956822 Nov 1999 WO
WO 0007494 Feb 2000 WO
WO 0010455 Mar 2000 WO
WO 0141867 Jun 2001 WO
WO 0148676 Jul 2001 WO
WO 0149364 Jul 2001 WO
WO 0167288 Sep 2001 WO
WO 0175660 Oct 2001 WO
WO 0209610 Feb 2002 WO
WO 0209811 Feb 2002 WO
WO 0236003 May 2002 WO
WO 0238031 May 2002 WO
WO 0238217 May 2002 WO
WO 0249500 Jun 2002 WO
WO 02058536 Aug 2002 WO
WO 02067122 Aug 2002 WO
WO 03001996 Jan 2003 WO
WO 03009207 Jan 2003 WO
WO 03030734 Apr 2003 WO
WO 03035165 May 2003 WO
WO 03084605 Oct 2003 WO
WO 2004008373 Jan 2004 WO
WO 2004032720 Apr 2004 WO
WO 2004034231 Apr 2004 WO
WO 2004034879 Apr 2004 WO
WO 2004034880 Apr 2004 WO
WO 2004034881 Apr 2004 WO
WO 2004034882 Apr 2004 WO
WO 2004034885 Apr 2004 WO
WO 2004034982 Apr 2004 WO
WO 2004034997 Apr 2004 WO
WO 2004034998 Apr 2004 WO
WO 2004035130 Apr 2004 WO
WO 2004036370 Apr 2004 WO
WO 2004036372 Apr 2004 WO
WO 2004036376 Apr 2004 WO
WO 2004036377 Apr 2004 WO
WO 20041034883 Apr 2004 WO
WO 2004037342 May 2004 WO
WO 2004043536 May 2004 WO
WO 2004091718 Oct 2004 WO
WO 2005007236 Jan 2005 WO
WO 2005028026 Mar 2005 WO
WO 2005028028 Mar 2005 WO
WO 2005031630 Apr 2005 WO
WO 2005051167 Jun 2005 WO
WO 2005051306 Jun 2005 WO
WO 2005117693 Dec 2005 WO
WO 2006014971 Feb 2006 WO
WO 2006014972 Feb 2006 WO
WO 2006020794 Feb 2006 WO
WO2006035392 Apr 2006 WO
Non-Patent Literature Citations (196)
Entry
Spector et al.; High and Low Perceived Self-Control of Epileptic Seizures; Epilepsia, vol. 42(4), Apr. 2001; pp. 556-564.
Adjouadi, et al. A new mathematical approach based on orthogonal operators for the detection of interictal spikes in epileptogenic data. Biomed. Sci. Instrum. 2004; 40: 175-80.
Adjouadi, et al. Detection of interictal spikes and artifactual data through orthogonal transformations. J. Clin. Neurophysiol. 2005; 22(1):53-64.
Adjouadi, et al. Interictal spike detection using the Walsh transform. IEEE Trans. Biomed. Eng. 2004; 51(5): 868-72.
Aksenova, et al. Nonparametric on-line detection of changes in signal spectral characteristics for early prediction of epilepsy seizure onset. J. Automation and Information Sciences. 2004; 36(8): 35-45.
Aksenova, et al. On-line disharmony detection for early prediction of epilepsy seizure onset. 5th International Workshop Neural Coding 2003. Aulla (Italy) Sep. 20-25, 2003. (Abstract).
Andrzejak, et al. Bivariate surrogate techniques: necessity, strengths, and caveats. Physical Review E. 2003; 68: 066202-1-066202-15.
Andrzejak, et al. Testing the null hypothesis of the nonexistence of a preseizure state. Physical Review E. 2003; 67: 010901-1-010901-4.
Aschenbrenner-Scheibe, et al. How well can epileptic seizures be predicted? An evaluation of a nonlinear method. Brain. 2003; 126: 2616-26.
Bangham et al. Diffusion of univalent ions across the lamellae of swollen phospholipids. 1965. J Mol. Biol. 13: 238-252.
Baruchi, et al. Functional holography of complex networks activity—From cultures to the human brain. Complexity. 2005; 10(3): 38 R 51.
Baruchi, et al. Functional holography of recorded neuronal networks activity. Neuroinformatics. 2004; 2(3): 333-51.
Ben-Hur, et al. Detecting stable clusters using principal component analysis. Methods Mol. Biol. 2003; 224: 159-82.
Bergey, et al. Epileptic seizures are characterized by changing signal complexity. Clin. Neurophysiol. 2001; 112(2): 241-9.
Betterton, et al. Determining State of Consciousness from the Intracranial Electroencephalogram (IEEG) for Seizure Prediction. From Proceeding (377) Modeling, Identification, and Control. 2003; 377-201: 313-317.
Bhattacharya, et al. Enhanced phase synchrony in the electroencephalograph gamma band for musicians while listening to music. Phys. Rev. E. 2001; 64:012902-1-4.
Boley, et al. Training Support Vector Machine using Adaptive Clustering. 2004 SIAM International Conference on Data Mining, Apr. 22-Apr. 24, 2004. Lake Buena Vista, FL, USA. 12 pages.
Burges, C. A Tutorial on Support Vector Machines for Pattern Recognition. Data Mining and Knowledge Discovery. 1998; 2: 121-167.
Cao, et al. Detecting dynamical changes in time series using the permutation entropy. Physical Review E. 2004; 70:046217-1-046217-7.
Carretero-Gonzalez, et al. Scaling and interleaving of subsystem Lyapunov exponents for spatio-temporal systems. Chaos. 1999; 9(2): 466-482.
Casdagli, et al. Characterizing nonlinearity in invasive EEG recordings from temporal lobe epilepsy. Physica D. 1996; 99 (2/3): 381-399.
Casdagli, et al. Nonlinear Analysis of Mesial Temporal Lobe Seizures Using a Surrogate Data Technique. Epilepsia. 1995; 36, suppl. 4, pp. 142.
Casdagli, et al. Non-linearity in invasive EEG recordings from patients with temporal lobe epilepsy. Electroencephalogr. Clin. Neurophysiol. 1997; 102(2): 98-105.
Cerf, et al. Criticality and synchrony of fluctuations in rhythmical brain activity: pretransitional effects in epileptic patients. Biol. Cybern. 2004; 90(4): 239-55.
Chaovalitwongse et al.; Reply to comments on “Performance of a seizure warning based on the dynamics of intracranial EEG”; Epilepsy Research, Elsevier Science Publishers, Amsterdam, NL; vol. 72; No. 1; pp. 82-84; Nov. 1, 2006.
Chaovalitwongse, et al. EEG Classification in Epilepsy. Annals. 2004; 2 (37): 1-31.
Chaovalitwongse, et al. Performance of a seizure warning algorithm based on the dynamics of intracranial EEG. Epilepsy Res. 2005; 64(3): 93-113.
Chavez, et al. Spatio-temporal dynamics prior to neocortical seizures: amplitude versphase couplings. IEEE Trans. Biomed. Eng. 2003; 50(5):571-83.
Chen et al.; Clinical utility of video-EEG monitoring; Pediatric Neurology; vol. 12; No. 3; pp. 220-224; 1995.
Crichton, Michael, “Terminal Man”, 1972, Ballantine Books, NY, NY, pp. 21-24, 32-33, 70-71, and 74-81.
D'Alessandro, et al. A multi-feature and multi-channel univariate selection process for seizure prediction. Clin. Neurophysiol. 2005; 116(3): 506-16.
D'Alessandro, et al. Epileptic seizure prediction using hybrid feature selection over multiple Intracranial EEG electrode contacts: a report of four patients. IEEE Trans. Biomed. Eng. 2003; 50(5): 603-15.
Drury, et al. Seizure prediction using scalp electroencephalogram. Exp. Neurol. 2003; 184 Suppl 1: S9-18.
Ebersole, J. S. Functional neuroimaging with EEG source models to localize epileptogenic foci noninvasively. Neurology, Available at http://www.uchospitals.edu/pdf/uch—001471.pdf. Accessed Feb. 28, 2006.
Ebersole, J. S. In search of seizure prediction: a critique. Clin. Neurophysiol. 2005; 116(3): 489-92.
Elbert et al. Chaos and Physiology: Deterministic Chaos in Excitable Cell Assemblies. Physiological Reviews. 1994; 74(1):1-47.
Elger, et al. Nonlinear EEG analysis and its potential role in epileptology. Epilepsia. 2000; 41 Suppl 3: S34-8.
Elger, et al. Seizure prediction by non-linear time series analysis of brain electrical activity. Eur. J. Neurosci. 1998; 10(2): 786-789.
Esteller, et al. A Comparison of Waveform Fractal Dimension Algorithms. IEEE Transactions on Circuits and Systems. 2001; vol. 48(2): 177-183.
Esteller, et al. Continuoenergy variation during the seizure cycle: towards an on-line accumulated energy. Clin. Neurophysiol. 2005; 116(3): 517-26.
Esteller, et al. Feature Parameter Optimization for Seizure Detection/prediction. Proceedings of the 23rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Istanbul, Turkey. Oct. 2001.
Faul, et al. An evaluation of automated neonatal seizure detection methods. Clin. Neurophysiol. 2005; 116(7): 1533-41.
Fein, et al. Common reference coherence data are confounded by power and phase effects. Electroencephalogr. Clin. Neurophysiol. 1988; 69:581-584.
Fell, et al. Linear inverse filtering improves spatial separation of nonlinear brain dynamics: a simulation study. J. Neurosci. Methods. 2000; 98(1): 49-56.
Firpi, et al. Epileptic seizure detection by means of genetically programmed artificial features. GECCO 2005: Proceedings of the 2005 conference on Genetic and evolutionary computation, vol. 1, pp. 461-466, Washington DC, USA, 2005. ACM Press.
Fisher et al. 1999. Reassessment: Vagnerve stimulation for epilepsy, A report of the therapeutics and technology assessment subcommitte of the American Academy of Neurology. Neurology.53: 666-669.
Franaszczuk et al.; An autoregressive method for the measurement of synchronization of interictal and ictal EEG signals; Biological Cybernetics, vol. 81; No. 1; pp. 3-9; 1999.
Gardner, A. B. A Novelty Detection Approach to Seizure Analysis from Intracranial EEG. Georgia Institute of Technology. Apr. 2004. A dissertation available at http://etd.gatech.edu/theses /available/etd-04122004-132404/unrestricted/gardner —andrew—b—200405 —phd.pdf. Accessed Feb. 28, 2006.
Geva, et al. Forecasting generalized epileptic seizures from the EEG signal by wavelet analysis and dynamic unsupervised fuzzy clustering. IEEE Trans. Biomed. Eng. 1998; 45(10): 1205-16.
Gigola, et al. Prediction of epileptic seizures using accumulated energy in a multiresolution framework. J. Neurosci. Methods. 2004; 138(1-2): 107-111.
Guyon, I. An introduction to variable and feature selection. Journal of Machine Learning Research. 2003; 3:1157-1182.
Guyon, I. Multivariate Non-Linear Feature Selection with Kernel Multiplicative Updates and Gram-Schmidt Relief. BISC FLINT-CIBI 2003 Workshop. Berkeley. 2003; p. 1-11.
Harrison, et al. Accumulated energy revised. Clin. Neurophysiol. 2005; 116 (3):527-31.
Harrison, et al. Correlation dimension and integral do not predict epileptic seizures. Chaos. 2005; 15(3): 33106-1-15.
Hearst M. Trends & Controversies: Support Vector Machines. IEEE Intelligent Systems. 1998; 13: 18-28.
Hively, et al. Channel-consistent forewarning of epileptic events from scalp EEG. IEEE Trans. Biomed. Eng. 2003; 50(5): 584-93.
Hively, et al. Detecting dynamical changes in nonlinear time series. Physics Letters A. 1999; 258: 103-114.
Hively, et al. Epileptic Seizure Forewarning by Nonlinear Techniques. ORNL/TM-2000/333 Oak Ridge National Laboratory. Nov. 2000. Available at http://computing.oml.gov/cse—home/staff/hively/NBICradaAnnualRpt FY00.pdf. Accessed Feb. 28, 2006.
Hjorth, B. Source derivation simplifies topographical EEG interpretation. Am. J. EEG Technol. 1980; 20: 121-132.
Hsu, et al. A practical guide to support vector classification. Technical report, Department of Computer Science and Information Technology, National Taiwan University, 2003. Available at http://www.csie.ntu.edu.tw/˜cjlin/papers/guide/guide.pdf. Accessed Feb. 28, 2006.
Huynh, J. A. Evaluation of Gene Selection Using Support Vector Machine Recursive Feature Elimination. Arizona State University. May 26, 2004. (28 pages).
Huynh, J. A. Evaluation of Gene Selection Using Support Vector Machine Recursive Feature Elimination. Presentation slides. (41 pages) (May 26, 2004).
Iasemidis, et al. Adaptive epileptic seizure prediction system. IEEE Trans. Biomed. Eng. 2003; 50(5):616-27.
Iasemidis, et al. Automated Seizure Prediction Paradigm. Epilepsia. 1998; vol. 39, pp. 56.
Iasemidis, et al, Chaos Theory and Epilepsy. The Neuroscientist. 1996; 2:118-126.
Iasemidis, et al. Comment on “Inability of Lyapunov exponents to predict epileptic seizures.” Physical Review Letters. 2005; 94(1):019801-1.
Iasemidis, et al. Detection of the Preictal Transition State in Scalp-Sphenoidal EEG Recordings. American Clinical Neurophysiology Society Annual Meeting, Sep. 1996. pp. C206.
Iasemidis, et al. Dynamical Interaction of the Epileptogenic Focwith Extrafocal Sites in Temporal Lobe Epilepsy (TLE). Ann. Neuro1.1997; 42, pp. 429. pp. M146.
Iasemidis, et al. Epileptogenic FocLocalization by Dynamical Analysis of Interictal Periods of EEG in Patients with Temporal Lobe Epilepsy. Epilepsia. 1997; 38, suppl. 8, pp. 213.
Iasemidis, et al. Localizing Preictal Temporal Lobe Spike Foci Using Phase Space Analysis. Electroencephalography and Clinical Neurophysiology. 1990; 75, pp. S63-64.
Iasemidis, et al. Long-term prospective on-line real-time seizure prediction. Clin. Neurophysiol. 2005; 116(3):532-44.
Iasemidis, et al. Long-Time-Scale Temporo-spatial Patterns of Entrainment of Preictal Electrocorticographic Data in Human Temporal Lobe Epilepsy. Epilepsia. 1990; 31(5):621.
Iasemidis, et al. Measurement and Quantification of Spatio-Temporal Dynamics of Human Epileptic Seizures. In: Nonlinear Signal Processing in Medicine, Ed. M. Akay, IEEE Press. 1999; pp. 1-27.
Iasemidis, et al. Modelling of ECoG in temporal lobe epilepsy. Biomed. Sci. Instrum. 1988; 24: 187-93.
Iasemidis, et al. Nonlinear Dynamics of EcoG Data in Temporal Lobe Epilepsy. Electroencephalography and Clinical Neurophysiology. 1998; 5, pp. 339.
Iasemidis, et al. Phase space topography and the Lyapunov exponent of electrocorticograms in partial seizures. Brain Topogr. 1990; 2(3): 187-201.
Iasemidis, et al. Preictal Entrainment of a Critical Cortical Mass is a Necessary Condition for Seizure Occurrence. Epilepsia. 1996; 37, suppl. 5. pp. 90.
Iasemidis, et al. Preictal-Postictal Versus Postictal Analysis for Epileptogenic Focus Localization. J. Clin. Neurophysiol. 1997; 14, pp. 144.
Iasemidis, et al. Quadratic binary programming and dynamic system approach to determine the predictability of epileptic seizures. Journal of Combinatorial Optimization. 2001; 5: 9-26.
Iasemidis, et al. Quantification of Hidden Time Dependencies in the EEG within the Framework of Non-Linear Dynamics. World Scientific. 1993; pp. 30-47.
Iasemidis, et al. Spatiotemporal dynamics of human epileptic seizures. World Scientific. 1996; pp. 26-30.
Iasemidis, et al. Spatiotemporal Evolution of Dynamical Measures Precedes Onset of Mesial Temporal Lobe Seizures. Epilepsia. 1994; 358, pp. 133.
Iasemidis, et al. Spatiotemporal Transition to Epileptic Seizures: A Nonlinear Dynamical Analysis of Scalp and Intracranial EEG Recordings. (In SILVA, F.L. Spatiotemporal Models in Biological and Artificial Systems. Ohmsha IOS Press. 1997; 37, pp. 81-88.).
Iasemidis, et al. The evolution with time of the spatial distribution of the largest Lyapunov exponent on the human epileptic cortex. World Scientific. 1991; pp. 49-82.
Iasemidis, et al. The Use of Dynamical Analysis of EEG Frequency Content in Seizure Prediction. American Electroencephalographic Society Annual Meeting, Oct. 1993.
Iasemidis, et al. Time Dependencies in Partial Epilepsy. 1993; 34, pp. 130-131.
Iasemidis, et al. Time dependencies in the occurrences of epileptic seizures. Epilepsy Res. 1994; 17(1): 81-94.
Iasemidis, L. D. Epileptic seizure prediction and control. IEEE Trans. Biomed. Eng. 2003; 50(5):549-58.
Jerger, et al. Early seizure detection. Journal of Clin. Neurophysiol. 2001; 18(3):259-68.
Jerger, et al. Multivariate linear discrimination of seizures. Clin. Neurophysiol. 2005; 116(3):545-51.
Jouny, et al. Characterization of epileptic seizure dynamics using Gabor atom density. Clin. Neurophysiol. 2003; 114(3):426-37.
Jouny, et al. Signal complexity and synchrony of epileptic seizures: is there an identifiable preictal period? Clin. Neurophysiol. 2005; 116(3):552-8.
Kapiris, et al. Similarities in precursory features in seismic shocks and epileptic seizures. Europhys. Lett. 2005; 69(4):657-663.
Katz, et al. Does interictal spiking change prior to seizures? Electroencephalogr. Clin. Neurophysiol. 1991; 79(2):153-6.
Kerem, et al. Forecasting epilepsy from the heart rate signal. Med. Biol. Eng. Comput. 2005; 43(2):230-9.
Khalilov, et al. Epileptogenic actions of GABA and fast oscillations in the developing hippocampus. Neuron. 2005; 48(5):787-96.
Korn, et al. Is there chaos in the brain? II. Experimental evidence and related models. C. R. Biol. 2003; 326(9):787-840.
Kraskov, A. Synchronization and Interdependence Measures and Their Application to the Electroencephalogram of Epilepsy Patients and Clustering of Data. Available at http://www.kfa-juelich.de/nic-series/volume24/nic-series-band24.pdf. Accessed Apr. 17, 2006 (106 pp).
Kreuz, et al. Measure profile surrogates: a method to validate the performance of epileptic seizure prediction algorithms. Phys. Rev. E. 2004; 69(6 Pt 1):061915-1-9.
Lachaux, et al. Measuring phase synchrony in brain signals. Hum. Brain Mapp. 1999; 8(4):194-208.
Lai, et al. Controlled test for predictive power of Lyapunov exponents: their inability to predict epileptic seizures. Chaos. 2004; 14(3):630-42.
Lai, et al. Inability of Lyapunov exponents to predict epileptic seizures. Phys. Rev. Lett. 2003; 91(6):068102-1-4.
Latka, et al. Wavelet analysis of epileptic spikes. Phys. Rev. E. 2003; 67(5 Pt 1):052902 (6 pages).
Le Van Quyen, et al. Anticipating epileptic seizures in real time by a non-linear analysis of similarity between EEG recordings. Neuroreport. 1999; 10(10):2149-55.
Le Van Quyen, et al. Author's second reply. The Lancet. 2003; 361:971.
Le Van Quyen, et al. Comparison of Hilbert transform and wavelet methods for the analysis of neuronal synchrony. J. Neurosci. Methods. 2001; 111(2):83-98.
Le Van Quyen, et al. Nonlinear analyses of interictal EEG map the brain interdependences in human focal epilepsy. Physica D. 1999; 127:250-266.
Le Van Quyen, et al. Preictal state identification by synchronization changes in long-term intracranial EEG recordings. Clin. Neurophysiol. 2005; 116(3):559-68.
Le Van Quyen, M. Anticipating epileptic seizures: from mathematics to clinical applications. C. R. Biol. 2005; 328(2):187-98.
Lehnertz, et al. Nonlinear EEG analysis in epilepsy: Its possible use for interictal focus localization, seizure anticipation, and prevention. J. Clin. Neurophysiol. 2001; 18(3):209-22.
Lehnertz, et al. Seizure prediction by nonlinear EEG analysis. IEEE Eng. Med. Biol. Mag. 2003; 22(1):57-63.
Lehnertz, et al. The First International Collaborative Workshop on Seizure Prediction: summary and data description. Clin. Neurophysiol. 2005; 116(3):493-505.
Lehnertz, K. Non-linear time series analysis of intracranial EEG recordings in patients with epilepsy—an overview. Int. J. Psychophysiol. 1999; 34(1):45-52.
Lemos, et al. The weighted average reference montage. Electroencephalogr. Clin. Neurophysiol. 1991; 79(5):361-70.
Li, et al. Fractal spectral analysis of pre-epileptic seizures in terms of criticality. J. Neural Eng. 2005; 2(2):11-16.
Li, et al. Linear and nonlinear measures and seizure anticipation in temporal lobe epilepsy. J. Comput. Neurosci. 2003; 15(3):335-45.
Li, et al. Non-linear, non-invasive method for seizure anticipation in focal epilepsy. Math. Biosci. 2003; 186(1):63-77.
Litt, et al. Prediction of epileptic seizures. Lancet Neurol. 2002; 1(1):22-30.
Litt, et al. Seizure prediction and the preseizure period. Curr. Opin. Neurol. 2002; 15(2):173-7.
Maiwald, et al. Comparison of three nonlinear seizure prediction methods by means of the seizure prediction characteristic. Physica D. 2004; 194:357-368.
Mangasarian, et al. Lagrangian Support Vector Machines. Journal of Machine Learning Research. 2001; 1:161-177.
Martinerie, et al. Epileptic seizures can be anticipated by non-linear analysis. Nat. Med. 1998; 4(10):1173-6.
McSharry, et al. Comparison of predictability of epileptic seizures by a linear and a nonlinear method. IEEE Trans. Biomed. Eng. 2003; 50(5):628-33.
McSharry, et al. Linear and non-linear methods for automatic seizure detection in scalp electro-encephalogram recordings. Med. Biol. Eng. Comput. 2002; 40(4):447-61.
McSharry, P. E. Detection of dynamical transitions in biomedical signals using nonlinear methods. Lecture Notes in Computer Science 2004; 3215:483-490.
Meng, et al. Gaussian mixture models of ECoG signal features for improved detection of epileptic seizures. Med. Eng. Phys. 2004; 26(5):379-93.
Mizuno-Matsumoto, et al. Wavelet-crosscorrelation analysis can help predict whether bursts of pulse stimulation will terminate after discharges. Clin. Neurophysiol. 2002; 113(1):33-42.
Mormann et al.; Seizure prediction: the long and winding road; Brain; vol. 130; No. 2; pp. 314-333; Sep. 28, 2006.
Mormann, et al. Automated detection of a preseizure state based on a decrease in synchronization in intracranial electroencephalogram recordings from epilepsy patients. Phys. Rev. E. 2003; 67(2 Pt 1):021912-1-10.
Mormann, et al. Epileptic seizures are preceded by a decrease in synchronization. Epilepsy Res. 2003; 53(3):173-85.
Mormann, et al. Mean phase coherence as a measure for phase synchronization and its application to the EEG of epilepsy patients. Physics D. 2000; 144:358-369.
Mormann, et al. On the predictability of epileptic seizures. Clin. Neurophysiol. 2005; 116(3):569-87.
Mormann, et al. Seizure anticipation: from algorithms to clinical practice. Curr. Opin. Neurol. 2006; 19(2):187-93.
Navarro, et al. Seizure anticipation in human neocortical partial epilepsy. Brain. 2002; 125:640-55.
Navarro, et al. Seizure anticipation: do mathematical measures correlate with video-EEG evaluation? Epilepsia. 2005; 46(3):385-96.
Niederhauser, et al. Detection of seizure precursors from depth-EEG using a sign periodogram transform. IEEE Trans. Biomed. Eng. 2003; 50(4):449-58.
Nigam, et al. A neural-network-based detection of epilepsy. Neurological Research. 2004; 26(1):55-60.
Osorio, et al. Automated seizure abatement in humans using electrical stimulation. Ann. Neurol. 2005; 57(2):258-68.
Osorio, et al. Performance reassessment of a real-time seizure-detection algorithm on long ECoG series. Epilepsia. 2002; 43(12):1522-35.
Osorio, et al. Real-time automated detection and quantitative analysis of seizures and short-term prediction of clinical onset. Epilepsia. 1998; 39(6):615-27.
Ossadtchi, et al. Hidden Markov modelling of spike propagation from interictal MEG data. Phys. Med. Biol. 2005; 50(14):3447-69.
Pflieger, et al. A noninvasive method for analysis of epileptogenic brain connectivity. Presented at the American Epilepsy Society 2004 Annual Meeting, New Orleans. Dec. 6, 2004. Epilepsia. 2004; 45 (Suppl. 7):70-71.
Pittman, V. Flexible Drug Dosing Produces Less Side-effects in People With Epilepsy. Dec. 29, 2005. Available at http://www.medicalnewstoday.com/medicalnews.php?newsid=35478. Accessed on Apr. 17, 2006.
Platt, et al. Large Margin DAGs for Multiclass Classification. S.A. Solla. T.K. Leen adn K. R. Muller (eds.). 2000; pp. 547-553.
Platt, J. C. Using Analytic QP and Sparseness to Speed Training of Support Vector Machines. Advances in Neural Information Processing Systems. 1999; 11:557-563.
Protopopescu, et al. Epileptic event forewarning from scalp EEG. J. Clin. Neurophysiol. 2001; 18(3):223-45.
Rahimi, et al. On the Effectiveness of Aluminum Foil Helmets: An Empirical Study. Available at http://people.csail.mit.edu/rahimi/helmet/. Accessed Mar. 2, 2006.
Rothman et al.; Local Cooling: a therapy for intractable neocortical epilepsy; Epilepsy Currents; vol. 3; No. 5; pp. 153-156; Sep./Oct. 2003.
Robinson, et al. Steady States and Global Dynamics of Electrical Activity in the Cerebral Cortex. Phys. Rev. E. 1998; (58):3557-3571.
Rudrauf, et al. Frequency flows and the time-frequency dynamics of multivariate phase synchronization in brain signals. Neurolmage. 2005. (19 pages.).
Saab, et al. A system to detect the onset of epileptic seizures in scalp EEG. Clin. Neurophysiol, 2005; 116:427-442.
Sackellares et al. Computer-Assisted Seizure Detection Based on Quantitative Dynamical Measures. American Electroencephalographic Society Annual Meeting, Sep. 1994.
Sackellares et al. Dynamical Studies of Human Hippocampin Limbic Epilepsy. Neurology. 1995; 45, Suppl. 4, pp. A 404.
Sackellares et al. Epileptic Seizures as Neural Resetting Mechanisms. Epilepsia. 1997; vol. 38, Sup. 3.
Sackellares et al. Measurement of Chaos to Localize Seizure Onset. Epilepsia. 1989; 30(5):663.
Sackellares et al. Relationship Between Hippocampal Atrophy and Dynamical Measures of EEG in Depth Electrode Recordings. American Electroencephalographic Society Annual Meeting, Sep. 1995. pp. A105.
Sackellares et al.; Predictability analysis for an automated seizure prediction algorithm; Journal of Clinical Neurophysiology; vol. 23; No. 6; pp. 509-520; Dec. 2006.
Sackellares, J. C. Epilepsy—when chaos fails. In: chaos in the brain? Eds. K. Lehnertz & C.E. Elger. World Scientific. 2000 (22 pages).
Salant, et al. Prediction of epileptic seizures from two-channel EEG. Med. Biol. Eng. Comput. 1998; 36(5):549-56.
Schelter et al.; Testing statistical significance of multivariate time series analysis techniques for epileptic seizure prediction; Chaos; vol. 16; pp. 013108-1-10; Jan. 2006.
Schelter, et al. Testing for directed influences among neural signals using partial directed coherence. J. Neurosci. Methods. 2006; 152(1-2):210-9.
Schindler, et al. EEG analysis with simulated neuronal cell models helps to detect pre-seizure changes. Clin. Neurophysiol. 2002; 113(4):604-14.
Schwartzkroin, P. Origins of the Epileptic State. Epilepsia. 1997; 38, supply. 8, pp. 853-858.
Sheridan, T. Humans and Automation. NY: John Wiley. 2002.
Shoeb et al. Patient-specific seizure detection. MIT Computer Science and Artificial Intelligence Laboratory. 2004; pp. 193-194.
Snyder et al; The statistics of a practical seizure warning system; Journal of Neural Engineering; vol. 5; pp. 392-401; 2008.
Staba, et al. Quantitative analysis of high-frequency oscillations (80-500 Hz) recorded in human epileptic hippocampand entorhinal cortex. J. Neurophysiol. 2002; 88(4):1743-52.
Stefanski, et al. Using chaos synchronization to estimate the largest Lyapunov exponent of nonsmooth systems. Discrete Dynamics in Nature and Society. 2000; 4:207-215.
Subasi, et al. Classification of EEG signals using neural network and logistic regression. Computer Methods Programs Biomed. 2005; 78(2):87-99.
Szoka et al. Procedure for preparation of liposomes with large internal aqueospace and high capture volume by reverse phase evaporation. 1978. Proc. Natl Acad. Sci. USA. 75: 4194-4198.
Tass, et al. Detection of n: m Phase Locking from Noisy Data: Application to Magnetoencephalography. Physical Review Letters. 1998; 81(15):3291-3294.
Terry, et al. An improved algorithm for the detection of dynamical interdependence in bivariate time-series. Biol. Cybern. 2003; 88(2):129-36.
Tetzlaff, et al. Cellular neural networks (CNN) with linear weight functions for a prediction of epileptic seizures. Int″l. J. of Neural Systems. 2003; 13(6):489-498.
Theiler, et al. Testing for non-linearity in time series: the method of surrogate data. Physica D. 1992; 58:77-94.
Tsakalis, K. S. Prediction and control of epileptic seizures: Coupled oscillator models. Arizona State University. (Slide: 53 pages) (No date).
Van Drongelen, et al. Seizure anticipation in pediatric epilepsy: use of Kolmogorov entropy. Pediatr. Neural. 2003; 29(3): 207-13.
Van Putten, M. Nearest neighbor phase synchronization as a measure to detect seizure activity from scalp EEG recordings. J. Clin. Neurophysiol. 2003; 20(5):320-5.
Venugopal, et al. A new approach towards predictability of epileptic seizures: KLT dimension. Biomed Sci. Instrum. 2003; 39:123-8.
Vonck, et al. Long-term amygdalohippocampal stimulation for refractory temporal lobe epilepsy. Ann. Neural, 2002; 52(5):556-65.
Vonck, et al. Long-term deep brain stimulation for refractory temporal lobe epilepsy. Epilepsia. 2005; 46(Suppl 5):98-9.
Vonck, et al. Neurostimulation for refractory epilepsy. Acta. Neurol. Belg. 2003; 103(4):213-7.
Weiss, P. Seizure prelude found by chaos calculation. Science News. 1998; 153(20):326.
Wells, R. B. Spatio-Temporal Binding and Dynamic Cortical Organization: Research Issues. Mar. 2005. Available at http://www.mrc.uidaho.edu/˜rwells/techdocs/Functional%20Column%20Research%20Issues.pdf. Accessed Mar. 2, 2006.
Widman, et al. Reduced signal complexity of intracellular recordings: a precursor for epileptiform activity? Brain Res. 1999; 836(1-2):156-63.
Winterhalder, et al. Sensitivity and specificity of coherence and phase synchronization analysis. (In Press) Phys. Lett. A. 2006.
Winterhalder, et al. The seizure prediction characteristic: a general framework to assess and compare seizure prediction methods. Epilepsy Behav. 2003; 4(3):318-25.
Wong et al.; A stochastic framework for evaluating seizure prediction algorithms using hiden markov models; Journal of Neurophysiology; vol. 97, No. 3; pp. 2525-2532; Oct. 4, 2006.
Yang et al.; Testing whether a prediction scheme is better than guess; Ch. 14 in Quantitative Neuroscience: Models, Algorithms, Diagnostics, and Therapeutic Applications; pp. 251-262; 2004.
Yang, et al. A supervised feature subset selection technique for multivariate time series. Available at http://infolab.usc.edu/DocsDemos/fsdm05.pdf. Accessed Mar. 2, 2006.
Yang, et al. CLe Ver: A feature subset selection technique for multivariate time series. T. B. Ho, D. Cheung, and H. Liu (Eds.): PAKDD. 2005; LNAI 3518: 516-522.
Yang, et al. Relation between Responsiveness to Neurotransmitters and Complexity of Epileptiform Activity in Rat Hippocampal CA1 Neurons. Epilepsia. 2002; 43(11):1330-1336.
Yatsenko, et al. Geometric Models, Fiber Bundles, and Biomedical Applications. Proceedings of Institute of Mathematics of NAS of Ukraine. 2004; 50 (Part 3):1518R1525.
Zaveri et al. Time-Frequency Analyses of Nonstationary Brain Signals. Electroencephalography and Clinical Neurophysiology. 1991; 79, pp. 28P-29P.
Zhang, et al. High-resolution EEG: cortical potential imaging of interictal spikes. Clin. Neurophysiol. 2003; 114(10):1963-73.
Bekas et al.; Low cost high performance uncertainty quantification; Conf. on High Performance Networking and Computing; Portland, Oregon; Article No. 8; (ISBN:978-1-60558-716-.
DiLorenzo, Daniel, U.S. Appl. No. 11/282,317 entitled “Closed-loop vagus nerve stimulation,” filed Nov. 17, 2005.
Related Publications (1)
Number Date Country
20100302270 A1 Dec 2010 US
Provisional Applications (1)
Number Date Country
61183449 Jun 2009 US