CA2338843C - System and method for facilitating clinical decision making - Google Patents
System and method for facilitating clinical decision making Download PDFInfo
- Publication number
- CA2338843C CA2338843C CA002338843A CA2338843A CA2338843C CA 2338843 C CA2338843 C CA 2338843C CA 002338843 A CA002338843 A CA 002338843A CA 2338843 A CA2338843 A CA 2338843A CA 2338843 C CA2338843 C CA 2338843C
- Authority
- CA
- Canada
- Prior art keywords
- estimator
- estimators
- values
- physiological parameters
- physiological
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4821—Determining level or depth of anaesthesia
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/389—Electromyography [EMG]
Abstract
Disclosed is a system and method for providing information to the user of a medical monitoring or diagnostic device to aid in the clinical decision making process.
The preferred embodiment uses two estimators or predictors of the same physiological quantity, with each of the estimators being designed to detect specific states or artifacts in the estimated parameter and thus operating at a different point on its respective ROC curve; one chosen to provide high sensitivity, the other chosen to provide high specificity. The divergence between the estimators is indicated by the use of a shaded region between their respective time trends. The use of 2 estimators of the same parameters with different performance characteristics allows the system and method of the present invention to derive additional information about the underlying physiologic process over and above that which would be available from a single estimator. The system and method of the present invention can derive information from not only the instantaneous values of the estimators and the difference between them, but also from the time trend of the difference.
The preferred embodiment uses two estimators or predictors of the same physiological quantity, with each of the estimators being designed to detect specific states or artifacts in the estimated parameter and thus operating at a different point on its respective ROC curve; one chosen to provide high sensitivity, the other chosen to provide high specificity. The divergence between the estimators is indicated by the use of a shaded region between their respective time trends. The use of 2 estimators of the same parameters with different performance characteristics allows the system and method of the present invention to derive additional information about the underlying physiologic process over and above that which would be available from a single estimator. The system and method of the present invention can derive information from not only the instantaneous values of the estimators and the difference between them, but also from the time trend of the difference.
Description
System and Method for Facilitating Clinical Decision Making BACKGROUND OF THE INVENTION
This invention relates to physiological monitoring systems and methods, and more particularly to physiological monitoring systems and methods which provide improved measurements of physiological parameters.
A medical monitoring or diagnostic device generally includes the ability to perform a measurement to obtain an estimate of the true value of a physiological parameter. This estimate is then used by a clinician to derive information conceming a patient's physical state, for the purpose of determining or monitoring a course of action with regard to the patient's medical care. An implicit assumption in this process is that the estimate of the physiological parameter provided by the measurement is identical to the actual, true value.
In actuality, the estimate of a parameter provided by a measurement has an inherent error, or difference between the estimated (measured) value and the true value. A measurement system may be such that it produces a random error (with an average value of zero), or it may be biased, such that the average error of the estimates is not equal to zero.
Concepts well understood in medical science are that of the sensitivity and specificity of a test. A binary estimate is one that can have 1 of 2 values:
positive or negative, yes or no, present or absent, etc. For example, inq the case of a binary parameter whose value can be positive or negative, the probability that any one measurement or estimate yields a positive response or prediction when the true value actually is positive is the "sensitivity" of the measurement. The sensitivity is thus the
This invention relates to physiological monitoring systems and methods, and more particularly to physiological monitoring systems and methods which provide improved measurements of physiological parameters.
A medical monitoring or diagnostic device generally includes the ability to perform a measurement to obtain an estimate of the true value of a physiological parameter. This estimate is then used by a clinician to derive information conceming a patient's physical state, for the purpose of determining or monitoring a course of action with regard to the patient's medical care. An implicit assumption in this process is that the estimate of the physiological parameter provided by the measurement is identical to the actual, true value.
In actuality, the estimate of a parameter provided by a measurement has an inherent error, or difference between the estimated (measured) value and the true value. A measurement system may be such that it produces a random error (with an average value of zero), or it may be biased, such that the average error of the estimates is not equal to zero.
Concepts well understood in medical science are that of the sensitivity and specificity of a test. A binary estimate is one that can have 1 of 2 values:
positive or negative, yes or no, present or absent, etc. For example, inq the case of a binary parameter whose value can be positive or negative, the probability that any one measurement or estimate yields a positive response or prediction when the true value actually is positive is the "sensitivity" of the measurement. The sensitivity is thus the
-2-probability of detecting a true positive outcome. The probability that any one measurement or estimate yields a negative response or prediction when the true value actually is negative is the "specificity" of the measurement. The specificity is thus the probability of detecting a true negative outcome.
An ideal estimate or predictor would be one whose sensitivity and specificity are both 100%; that is, one in which there is zero probability of false detection, either false positive or false negative. Physically realizable measurements, however, rarely achieve this ideal. They instead typically have both sensitivity and specificity less than 100%.
It is well known in the art that the sensitivity and specificity of a measurement system can be biased by design either toward higher sensitivity at the expense of lower specificity, or toward higher specificity at the expense of lower sensitivity.
The possible operating points of a measurement system can be described by the Receiver Operating Characteristic (ROC) curve. A hypothetical ROC curve is shown in Figure 1. Actual ROC curves may be of different forms, but all exhibit a trade-off of sensitivity for specificity.
A binary estimate of a binary parameter translates directly to a predicted value for the parameter. It is also well understood in the art that predicted values can be obtained from continuous estimators. A threshold value may be applied to a continuous estimator of a binary parameter in order to provide a binary prediction.
For example, if the value of the estimate is less than the threshold, the parameter is predicted to have one of the binary states. If the value of the estimate is greater than or equal to the threshold, the parameter is predicted to have the other binary state.
The relationship between the estimator and the parameter may be modeled from experimental data as a probability of response, as shown in Figure 2. The value of the
An ideal estimate or predictor would be one whose sensitivity and specificity are both 100%; that is, one in which there is zero probability of false detection, either false positive or false negative. Physically realizable measurements, however, rarely achieve this ideal. They instead typically have both sensitivity and specificity less than 100%.
It is well known in the art that the sensitivity and specificity of a measurement system can be biased by design either toward higher sensitivity at the expense of lower specificity, or toward higher specificity at the expense of lower sensitivity.
The possible operating points of a measurement system can be described by the Receiver Operating Characteristic (ROC) curve. A hypothetical ROC curve is shown in Figure 1. Actual ROC curves may be of different forms, but all exhibit a trade-off of sensitivity for specificity.
A binary estimate of a binary parameter translates directly to a predicted value for the parameter. It is also well understood in the art that predicted values can be obtained from continuous estimators. A threshold value may be applied to a continuous estimator of a binary parameter in order to provide a binary prediction.
For example, if the value of the estimate is less than the threshold, the parameter is predicted to have one of the binary states. If the value of the estimate is greater than or equal to the threshold, the parameter is predicted to have the other binary state.
The relationship between the estimator and the parameter may be modeled from experimental data as a probability of response, as shown in Figure 2. The value of the
-3-estimator at which the probability of response is 50% is used as the threshold value in the preferred embodiment.
In a similar manner, multiple thresholds may be applied to a continuous estimator to predict the value of an ordinal parameter. For example, the Modified Observer's Assessment of Alertness/Sedation scale is commonly used to quantify an anesthetized surgical patient's level of sedation. A set of five ordinal thresholds may be applied to a single continuous estimator of the patient's sedation state to obtain a predicted sedation state.
MODIFIED OBSERVER'S ASSESSMENT OF ALERTNESS/SEDATION SCALE
RESPONSE SCORE LEVEL
Responds readily to name spoken in normal tone 5 Lethargic response to name spoken in normal tone 4 Responds only after name is called loudly and/or repeatedly 3 Responds only after mild prodding or shaking 2 Does not respond to mild prodding or shaking 1 Does not respond to test stimulus 0 The case of a continuous estimator used to predict an ordinal parameter may be extended to that of a continuous estimator used to predict the value of a continuous parameter. The relationship between the estimator and the parameter it predicts need not be linear, just as the threshold need not be evenly spaced in the case of an ordinal estimator. Any continuous parameter (e.g., a patient's state of sedation) may be considered to be a set of discrete ordinal states, each of which may be estimated by applying an appropriate set of threshold values to the continuous predictor.
While each state is distinct, they are all steps on a continuum of states from one end of the state range to the other. Thus the concept of a ROC curve may be applied to the
In a similar manner, multiple thresholds may be applied to a continuous estimator to predict the value of an ordinal parameter. For example, the Modified Observer's Assessment of Alertness/Sedation scale is commonly used to quantify an anesthetized surgical patient's level of sedation. A set of five ordinal thresholds may be applied to a single continuous estimator of the patient's sedation state to obtain a predicted sedation state.
MODIFIED OBSERVER'S ASSESSMENT OF ALERTNESS/SEDATION SCALE
RESPONSE SCORE LEVEL
Responds readily to name spoken in normal tone 5 Lethargic response to name spoken in normal tone 4 Responds only after name is called loudly and/or repeatedly 3 Responds only after mild prodding or shaking 2 Does not respond to mild prodding or shaking 1 Does not respond to test stimulus 0 The case of a continuous estimator used to predict an ordinal parameter may be extended to that of a continuous estimator used to predict the value of a continuous parameter. The relationship between the estimator and the parameter it predicts need not be linear, just as the threshold need not be evenly spaced in the case of an ordinal estimator. Any continuous parameter (e.g., a patient's state of sedation) may be considered to be a set of discrete ordinal states, each of which may be estimated by applying an appropriate set of threshold values to the continuous predictor.
While each state is distinct, they are all steps on a continuum of states from one end of the state range to the other. Thus the concept of a ROC curve may be applied to the
-4-overall continuum of states as well as any pair of states and the performance of the estimator may be tuned to vary its sensitivity and specificity with respect to the full range of the parameter. It is important to note that an explicit threshold need not be applied to a continuous estimator in order to use the ROC concept; the probability of response is implied through the correlation between the estimator and the underlying parameter.
If a physically realizable system does not have both 100% sensitivity and specificity, then an important design issue is the point on the ROC curve corresponding to the actual sensitivity and specificity of the measurement.
This is the operating point. In the case of a continuous estimator, varying the threshold value shifts the operating point along the ROC curve. The optimal location for the operating point on the ROC curve is preferably derived by means of a cost-benefit analysis.
Weighing the cost of a false detection (both false positive and false negative) against the benefit of a correct detection performs such an analysis. When the cost and benefit are not equivalent, it may be desirable to shift the operating point to bias the prediction toward increased sensitivity or specificity. For example, a system designed to screen patients for the presence of a fatal disease may be designed to have very high sensitivity at the expense of low specificity. The increased detection rate afforded by the elevated sensitivity may be well worth the extra testing required to rule out false positives. The design procedure, however, requires that a choice be made between high sensitivity and high specificity in order to specify an operating point.
US Patent No. 4,517,986 issued to Bilgutay teaches the use of a series of amplifiers with different sensitivities in order to correctly measure blood pressure from pressure pulse waves of different amplitudes. This patent, however, teaches only the use of a set of amplifiers with different transfer functions to amplify different parts of the pressure waveform. It does not teach selecting different points on the ROC
curve,
If a physically realizable system does not have both 100% sensitivity and specificity, then an important design issue is the point on the ROC curve corresponding to the actual sensitivity and specificity of the measurement.
This is the operating point. In the case of a continuous estimator, varying the threshold value shifts the operating point along the ROC curve. The optimal location for the operating point on the ROC curve is preferably derived by means of a cost-benefit analysis.
Weighing the cost of a false detection (both false positive and false negative) against the benefit of a correct detection performs such an analysis. When the cost and benefit are not equivalent, it may be desirable to shift the operating point to bias the prediction toward increased sensitivity or specificity. For example, a system designed to screen patients for the presence of a fatal disease may be designed to have very high sensitivity at the expense of low specificity. The increased detection rate afforded by the elevated sensitivity may be well worth the extra testing required to rule out false positives. The design procedure, however, requires that a choice be made between high sensitivity and high specificity in order to specify an operating point.
US Patent No. 4,517,986 issued to Bilgutay teaches the use of a series of amplifiers with different sensitivities in order to correctly measure blood pressure from pressure pulse waves of different amplitudes. This patent, however, teaches only the use of a set of amplifiers with different transfer functions to amplify different parts of the pressure waveform. It does not teach selecting different points on the ROC
curve,
-5-or using the divergence between simultaneous estimators as a measure of deriving information about blood pressure.
It is therefore a principal object of the present invention to provide a system and method for measuring physiological parameters that provide more information concerning the measured physiological parameter than known systems.
Another object of the present invention is to provide a system and method for measuring physiological parameters that uses multiple simultaneous estimators.
SUMMARY OF THE INVENTION
The system and method of the present invention provides information to the user of a medical monitoring or diagnostic device to aid in determining the true state of the measured parameter. The preferred embodiment uses two or more estimators or predictors of the same physiological process. Each of the estimators is designed to detect or predict either a specific state of the underlying parameter, or alternatively a type or class of artifact. Each estimator will therefore have a different sensitivity and specificity with respect to predicting a specific state or artifact, and thus will have a different operating point on the ROC curve corresponding to the prediction of that state or artifact. The estimators are designed such that they have similar values in the absence of the states or artifacts they are designed to predict. Furthermore, they are designed so that the degree of divergence between the estimators in the presence of specific states or artifacts will be similar across a population. The use of two or more estimators of the same parameter with different performance characteristics allows the user of the system and method of the present invention to derive information from not only the values of the estimators but also from the degree of their divergences. This enables the user to derive additional knowledge about the underlying measured parameter or state over and above that which would be available from a single = ' 69675-497
It is therefore a principal object of the present invention to provide a system and method for measuring physiological parameters that provide more information concerning the measured physiological parameter than known systems.
Another object of the present invention is to provide a system and method for measuring physiological parameters that uses multiple simultaneous estimators.
SUMMARY OF THE INVENTION
The system and method of the present invention provides information to the user of a medical monitoring or diagnostic device to aid in determining the true state of the measured parameter. The preferred embodiment uses two or more estimators or predictors of the same physiological process. Each of the estimators is designed to detect or predict either a specific state of the underlying parameter, or alternatively a type or class of artifact. Each estimator will therefore have a different sensitivity and specificity with respect to predicting a specific state or artifact, and thus will have a different operating point on the ROC curve corresponding to the prediction of that state or artifact. The estimators are designed such that they have similar values in the absence of the states or artifacts they are designed to predict. Furthermore, they are designed so that the degree of divergence between the estimators in the presence of specific states or artifacts will be similar across a population. The use of two or more estimators of the same parameter with different performance characteristics allows the user of the system and method of the present invention to derive information from not only the values of the estimators but also from the degree of their divergences. This enables the user to derive additional knowledge about the underlying measured parameter or state over and above that which would be available from a single = ' 69675-497
-6-estimator. The system and method of the present invention can derive information from not only the instantaneous values of the estimators and their divergences, but also from the time trend of the values of the estimators and their divergences.
According to one aspect of the present invention, there is provided a system for evaluating values of physiological parameters, said system comprising: a signal acquisition device for obtaining physiological signals from a human body; a processor for deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter and for computing relationships between said estimators; a display for displaying information concerning said relationships between said two or more estimators.
According to another aspect of the present invention, there is provided a system for evaluating values of physiological parameters, said system comprising: a signal acquisition device for obtaining physiological signals from a human body; a processor for deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter; a display for displaying said two or more estimators.
According to still another aspect of the present invention, there is provided a method for evaluating values of physiological parameters, said method comprising the steps of:
obtaining physiological signals from a human body; deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter and computing relationships between said estimators; displaying information concerning said relationships between said two or more estimators.
According to yet another aspect of the present invention, there is provided a method for evaluating values of physiological parameters, said method comprising the steps -6a-of: obtaining physiological signals from a human body; deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter; displaying said two or more estimators.
These and other objects and features of the present invention will be more fully understood from the following detailed description which should be read in light of the accompanying drawings.
BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 is a graph of a hypothetical Receiver Operating Characteristic (ROC) curve.
Figure 2 is a graph of a hypothetical Probability of Response curve.
Figure 3 is a block diagram of the functional components of a system of the present invention.
Figure 4 is a block diagram of the process for the determination of the first and second estimators utilized in the preferred embodiment of the system and method of the present invention.
Figure 5 is a graph of the ROC curves generated by the system and method of the present invention, in which each estimator (predictor) has a different ROC curve.
Figure 6 is a graph of the time trend of the two estimators (predictors) generated by the system and method of the present invention in which the shaded region indicates the divergence between the two estimators.
Figure 7 is a graph of the ROC curves generated by an alternate embodiment of the
According to one aspect of the present invention, there is provided a system for evaluating values of physiological parameters, said system comprising: a signal acquisition device for obtaining physiological signals from a human body; a processor for deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter and for computing relationships between said estimators; a display for displaying information concerning said relationships between said two or more estimators.
According to another aspect of the present invention, there is provided a system for evaluating values of physiological parameters, said system comprising: a signal acquisition device for obtaining physiological signals from a human body; a processor for deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter; a display for displaying said two or more estimators.
According to still another aspect of the present invention, there is provided a method for evaluating values of physiological parameters, said method comprising the steps of:
obtaining physiological signals from a human body; deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter and computing relationships between said estimators; displaying information concerning said relationships between said two or more estimators.
According to yet another aspect of the present invention, there is provided a method for evaluating values of physiological parameters, said method comprising the steps -6a-of: obtaining physiological signals from a human body; deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter; displaying said two or more estimators.
These and other objects and features of the present invention will be more fully understood from the following detailed description which should be read in light of the accompanying drawings.
BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 is a graph of a hypothetical Receiver Operating Characteristic (ROC) curve.
Figure 2 is a graph of a hypothetical Probability of Response curve.
Figure 3 is a block diagram of the functional components of a system of the present invention.
Figure 4 is a block diagram of the process for the determination of the first and second estimators utilized in the preferred embodiment of the system and method of the present invention.
Figure 5 is a graph of the ROC curves generated by the system and method of the present invention, in which each estimator (predictor) has a different ROC curve.
Figure 6 is a graph of the time trend of the two estimators (predictors) generated by the system and method of the present invention in which the shaded region indicates the divergence between the two estimators.
Figure 7 is a graph of the ROC curves generated by an alternate embodiment of the
-7-system and method of the present invention, in which each estimator (predictor) has an identical ROC curve.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT
It is well known in the state of the art that EEG data is predominantly contained in the frequency range of 0 - 40 Hz. Lighter states of anesthesia (i.e., those tending toward consciousness) are typically characterized by an increased amount of electrical power in the EEG frequency band of 20 - 40 Hz. Thus, an estimator intended to predict whether or not a patient was conscious might use the amount of power in this frequency band. An increased amount of power would predict an increased probability that the patient was conscious.
It is also well known that electromyographic (EMG) signals, the electrical activity arising from muscles, are predominantly contained in the frequency range of 30 - 300 Hz. EMG signals may thus appear in the same frequency range as EEG
signals. EMG activity tends to increase as an anesthetized patient approaches consciousness. This EMG activity tends to increase power in the 30-40 Hz frequency band in conjunction with that of the EEG, enhancing the ability of the power in the 20-40 Hz band to predict consciousness. Increasing EMG activity, however, may also occur in a deeply anesthetized .(deep) patient in conjunction with movement associated with reflexive motor activity. In this situation, the EMG activity tends to increase power in the 30-40 Hz frequency band, providing an effect on the estimator which is in opposition to that of the EEG. This increases the probability of a false prediction of consciousness and detracts from the overall performance of a diagnostic estimate of EEG activity based on the 20-40 Hz frequency band.
A preferred embodiment of this invention uses two (2) different estimators to predict the state of consciousness of a patient under anesthesia. Each of the estimators
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT
It is well known in the state of the art that EEG data is predominantly contained in the frequency range of 0 - 40 Hz. Lighter states of anesthesia (i.e., those tending toward consciousness) are typically characterized by an increased amount of electrical power in the EEG frequency band of 20 - 40 Hz. Thus, an estimator intended to predict whether or not a patient was conscious might use the amount of power in this frequency band. An increased amount of power would predict an increased probability that the patient was conscious.
It is also well known that electromyographic (EMG) signals, the electrical activity arising from muscles, are predominantly contained in the frequency range of 30 - 300 Hz. EMG signals may thus appear in the same frequency range as EEG
signals. EMG activity tends to increase as an anesthetized patient approaches consciousness. This EMG activity tends to increase power in the 30-40 Hz frequency band in conjunction with that of the EEG, enhancing the ability of the power in the 20-40 Hz band to predict consciousness. Increasing EMG activity, however, may also occur in a deeply anesthetized .(deep) patient in conjunction with movement associated with reflexive motor activity. In this situation, the EMG activity tends to increase power in the 30-40 Hz frequency band, providing an effect on the estimator which is in opposition to that of the EEG. This increases the probability of a false prediction of consciousness and detracts from the overall performance of a diagnostic estimate of EEG activity based on the 20-40 Hz frequency band.
A preferred embodiment of this invention uses two (2) different estimators to predict the state of consciousness of a patient under anesthesia. Each of the estimators
-8-is designed to have a value that is reflective of the probability that the patient is conscious. In addition, one of the estimators is designed to detect the presence of EMG activity. The estimators will have nearly identical values in the absence of EMG
activity, and will diverge in the presence of EMG activity. The estimators thus have different sensitivity and specificity with respect to the prediction of consciousness and the presence of EMG activity. The instantaneous value of each of the estimators is trended over time, and a shaded region is created on the trend display to graphically indicate the divergence between the two (2) estimators. The width of the shaded region is thus indicative of the degree to which EMG activity is influencing the estimation of the patient's state of consciousness, thereby imparting additional clinical information about the underlying state of consciousness over and above that of each of the estimators alone. It should be noted that this embodiment may be easily extended to N estimators, where N is a number greater than two. Each of these additional estimators would be designed to detect some specific consciousness or EEG
state or type of artifact, and to be similar to the others in the absence of the state or artifact it is designed to detect.
Referring to Figure 3, a system for generating multiple estimators and representations of relationships with such estimators is shown. Data acquisition device 12 can be any device suitable for obtaining physiological signals from a human body and converting them to a form for processing. Processor 14 processes such signals and computes the estimators as described below. In addition, processor 14 will calculate the relationships between such estimators. Display 16, which in a preferred embodiment is a CRT display, can be any video, hard copy or other similar display used to display computer generated information. In certain embodiments where the representations of the estimators and the relationship between such estimators is graphical in nature, the output device should be able to display graphical images.
The actual form of the estimators is derived using an empirical design process.
activity, and will diverge in the presence of EMG activity. The estimators thus have different sensitivity and specificity with respect to the prediction of consciousness and the presence of EMG activity. The instantaneous value of each of the estimators is trended over time, and a shaded region is created on the trend display to graphically indicate the divergence between the two (2) estimators. The width of the shaded region is thus indicative of the degree to which EMG activity is influencing the estimation of the patient's state of consciousness, thereby imparting additional clinical information about the underlying state of consciousness over and above that of each of the estimators alone. It should be noted that this embodiment may be easily extended to N estimators, where N is a number greater than two. Each of these additional estimators would be designed to detect some specific consciousness or EEG
state or type of artifact, and to be similar to the others in the absence of the state or artifact it is designed to detect.
Referring to Figure 3, a system for generating multiple estimators and representations of relationships with such estimators is shown. Data acquisition device 12 can be any device suitable for obtaining physiological signals from a human body and converting them to a form for processing. Processor 14 processes such signals and computes the estimators as described below. In addition, processor 14 will calculate the relationships between such estimators. Display 16, which in a preferred embodiment is a CRT display, can be any video, hard copy or other similar display used to display computer generated information. In certain embodiments where the representations of the estimators and the relationship between such estimators is graphical in nature, the output device should be able to display graphical images.
The actual form of the estimators is derived using an empirical design process.
9 PCT/US99/17688 A development database is first created which contains segments of data that serve as prototypes of the states or artifacts to be identified. The database must contain recordings from a large number of test subjects so that it is representative of a broad population. In the preferred embodiment described above, the development database contains segments of EEG signals recorded at various awake and anesthetized states and during intraoperative muscle reflex movements in a wide variety of subjects.
Estimators may be developed using various techniques that are common in the art, such as cluster or discriminant analysis. The estimators are designed such that they have similar values in the absence of the specific state or artifact they are designed to detect, and such that they have an approximately constant divergence across the entire population in the presence of the state or artifact they are designed to detect. It is important to recognize that other estimators may provide equivalent or even enhanced performance over those described in this prefen:ed embodiment. The present invention does not depend on the form of the estimators.
Figure 4 is a block diagram of the process for calculating the estimators in the first preferred embodiment. In step 20, the EEG signal is first sampled at 128 Hz and separated into two second epochs x, for processing.
xi where i=1,2,...N
The variable N is the number of samples in the epoch; in this case, N = 256.
The data is low pass filtered in step 22 using a fifth order Butterworth filter with a corner frequency (fc) at 50 Hz. This corner frequency was chosen to allow the entire EEG signal to pass through the filter, while excluding AC power line interference (60 Hz).
yi = b,xi + b,xi-, + b,xi-, + b4xi-3 + bsxi-4+ b6xi_s - azYi-, - a3Yi-2 - a4Yi-3 - aSYi-4 - a6Yi-s
Estimators may be developed using various techniques that are common in the art, such as cluster or discriminant analysis. The estimators are designed such that they have similar values in the absence of the specific state or artifact they are designed to detect, and such that they have an approximately constant divergence across the entire population in the presence of the state or artifact they are designed to detect. It is important to recognize that other estimators may provide equivalent or even enhanced performance over those described in this prefen:ed embodiment. The present invention does not depend on the form of the estimators.
Figure 4 is a block diagram of the process for calculating the estimators in the first preferred embodiment. In step 20, the EEG signal is first sampled at 128 Hz and separated into two second epochs x, for processing.
xi where i=1,2,...N
The variable N is the number of samples in the epoch; in this case, N = 256.
The data is low pass filtered in step 22 using a fifth order Butterworth filter with a corner frequency (fc) at 50 Hz. This corner frequency was chosen to allow the entire EEG signal to pass through the filter, while excluding AC power line interference (60 Hz).
yi = b,xi + b,xi-, + b,xi-, + b4xi-3 + bsxi-4+ b6xi_s - azYi-, - a3Yi-2 - a4Yi-3 - aSYi-4 - a6Yi-s
-10-where b, = 0.3198, b2 = 1.5990, b3 = 3.1980, b4 = 3.1980, b5 = 1.5990, b6 = 0.3198 a2=2.7872,a3=3.4091,a4=2.1972,a5=0.7377,a6=0.1023 In step 24, the filtered epoch y; is then multiplied by a Blackman window (w_blk) (J.N. Little, L. Shure, Matlab Signal Processing Toolbox User's Guide, The MathWorks, Inc. Natick, MA 1992, pg. 2-9) to improve spectral estimation, w_1ilk. = 0.42-0.5co 21t 1-1 + 0.08 4n =-1 where i=1,2,3,...N
y~'tlt = y w_blk The filtered, windowed data is transformed into the frequency domain in step 26 by use of the Fast Fourier Transform (FFT) technique N-, nxki Y_FILT = E y_~'ilt e N
k t-0 where In the preferred embodiment, n=256 and thus Y_FILT,, contains 129 complex values, corresponding to frequencies ranging from 0 to 64.0 Hz in 0.5 Hz increments.
y~'tlt = y w_blk The filtered, windowed data is transformed into the frequency domain in step 26 by use of the Fast Fourier Transform (FFT) technique N-, nxki Y_FILT = E y_~'ilt e N
k t-0 where In the preferred embodiment, n=256 and thus Y_FILT,, contains 129 complex values, corresponding to frequencies ranging from 0 to 64.0 Hz in 0.5 Hz increments.
-11-The power spectrum is computed as the squared magnitude of Y_FILT, in step 28:
P
k =I Y_FILT
Pk contains, 129 values, corresponding to power values at frequencies ranging from 0 to 64.0 Hz in 0.5 Hz increments.
The primary and secondary estimators are then calculated as ao EP
E =100 00 I ~ P~
E =100k-40 'E P
20 In steps 30 and 32, the estimators are calculated as a power in a specified frequency band relative to the power of the entire spectrum of interest, in percent. The estimator El is the percentage of the total power (2 - 40 Hz) that is contained in the 20 to 40 Hz frequency band. Estimator E2 is the percentage of the total power (2 - 30 Hz) that is contained in the 20 to 30 Hz frequency band. Estimator E, will be 25 sensitive to the presence of EMG activity while estimator E2 will not. In the absence of significant EMG activity, both estimators will have similar values since the 30 to 40 Hz frequency band will contain a negligible proportion of the total power. In this case, both will be reflective of the probability of consciousness as indicted by the power in the 20 to 30 Hz frequency band; higher power indicating an increased
P
k =I Y_FILT
Pk contains, 129 values, corresponding to power values at frequencies ranging from 0 to 64.0 Hz in 0.5 Hz increments.
The primary and secondary estimators are then calculated as ao EP
E =100 00 I ~ P~
E =100k-40 'E P
20 In steps 30 and 32, the estimators are calculated as a power in a specified frequency band relative to the power of the entire spectrum of interest, in percent. The estimator El is the percentage of the total power (2 - 40 Hz) that is contained in the 20 to 40 Hz frequency band. Estimator E2 is the percentage of the total power (2 - 30 Hz) that is contained in the 20 to 30 Hz frequency band. Estimator E, will be 25 sensitive to the presence of EMG activity while estimator E2 will not. In the absence of significant EMG activity, both estimators will have similar values since the 30 to 40 Hz frequency band will contain a negligible proportion of the total power. In this case, both will be reflective of the probability of consciousness as indicted by the power in the 20 to 30 Hz frequency band; higher power indicating an increased
-12-probability of consciousness, lower power indicating a lower probability of consciousness. In the presence of significant EMG activity, however, the value of E, will increase much more than that of E, and the values of the two estimators will diverge. Since increased EMG activity may be reflective of both increasing levels of consciousness as well as reflex movements in a deeply anesthetized (deeply unconscious) patient, E, will exhibit a higher sensitivity to an increased level of consciousness than E2. Conversely, E2 will exhibit a higher specificity to an increased level of consciousness than El. The effect of the design process is to produce two (or more) estimators of the same parameter that may have different ROC curves and different operating points (OP) on their respective curves. As shown in Figure 5, estimator Et operating at OP A provides predicts consciousness with a high sensitivity, while estimator E2 operating at OP B provides high specificity.
The divergence between the values of the estimators serves as an indicator of confidence in the prediction. When the estimators are in agreement (i.e., little or no divergence), the user can be confident that the specific states or artifacts that would be indicated by the divergence are not present and that the underlying parameter is being estimated with a high degree of accuracy. A divergence between the estimators indicates to the user the need for careful assessment of the situation. In a preferred embodiment, a divergence indicates the need for the user to utilize clinical judgement and clinical endpoints to determine which estimator more accurately reflects the state of the patient. If, for example, in the middle of surgery, E1 suddenly increases while E2 does not, an anesthesiologist must consider the level of sedative and analgesic agents administered to the patient. If these indicate that the patient should be deeply anesthetized and there is no concomitant increase in heart rate or blood pressure, then estimator E, has probably increased in response to a reflex movement. Since reflex movements are unrelated to the patient's anesthetic brain state, estimator E, will more accurately reflect the probability of consciousness. In contrast, if both estimators are increasing over time at the end of a case in response to a decreased dose or cessation
The divergence between the values of the estimators serves as an indicator of confidence in the prediction. When the estimators are in agreement (i.e., little or no divergence), the user can be confident that the specific states or artifacts that would be indicated by the divergence are not present and that the underlying parameter is being estimated with a high degree of accuracy. A divergence between the estimators indicates to the user the need for careful assessment of the situation. In a preferred embodiment, a divergence indicates the need for the user to utilize clinical judgement and clinical endpoints to determine which estimator more accurately reflects the state of the patient. If, for example, in the middle of surgery, E1 suddenly increases while E2 does not, an anesthesiologist must consider the level of sedative and analgesic agents administered to the patient. If these indicate that the patient should be deeply anesthetized and there is no concomitant increase in heart rate or blood pressure, then estimator E, has probably increased in response to a reflex movement. Since reflex movements are unrelated to the patient's anesthetic brain state, estimator E, will more accurately reflect the probability of consciousness. In contrast, if both estimators are increasing over time at the end of a case in response to a decreased dose or cessation
-13-of the administration of anesthetic agent and E, suddenly increases rapidly, the likely cause is increasing muscle tone due to returning consciousness. In this case, E, will be the more reliable predictor of the probability of consciousness. In all cases when the estimators diverge, it is a prompt to the clinician to consider "How did we get here?".
Figure 6 shows a sample trend of the two estimators and their divergence during a typical surgical case. At point A, the beginning of the surgical case, the patient is conscious. Both estimators are at high values, indicating a light anesthetic state and a high probability of consciousness. The high value of E2 (heavy line) 40 indicates a substantial contribution to the 20 - 40 Hz power from the EEG. The shaded divergence region 42 indicates that EMG activity contributes even more power, confirming the prediction of consciousness. At point B the two estimators are identical, both predicting a deep anesthetic state. The absence of a shaded divergence region indicates a high degree of confidence in the prediction of a low probability of consciousness.
At point C, estimator E, (thin line) 44 is indicating a light anesthetic state, while estimator E, (heavy line) 40 is indicating a deep anesthetic state. The wide divergence indicated by the shaded region 46 indicates a need to interpret this data in light of clinical judgement and other clinical endpoints. The sudden increase in E, occurred intraoperatively at a time when the patient had been administered doses of anesthetics generally considered adequate, and there was no sudden increase in surgical stimulation. For this reason, the increase in E, is most likely due to a burst of EMG
activity arising from an involuntary reflex movement, indicating that E, is the more accurate assessment of the patient's state of consciousness.
At point D, both estimators are trending upward, indicating a change in the patient to a lighter state of anesthesia. The absence of a divergence between the two
Figure 6 shows a sample trend of the two estimators and their divergence during a typical surgical case. At point A, the beginning of the surgical case, the patient is conscious. Both estimators are at high values, indicating a light anesthetic state and a high probability of consciousness. The high value of E2 (heavy line) 40 indicates a substantial contribution to the 20 - 40 Hz power from the EEG. The shaded divergence region 42 indicates that EMG activity contributes even more power, confirming the prediction of consciousness. At point B the two estimators are identical, both predicting a deep anesthetic state. The absence of a shaded divergence region indicates a high degree of confidence in the prediction of a low probability of consciousness.
At point C, estimator E, (thin line) 44 is indicating a light anesthetic state, while estimator E, (heavy line) 40 is indicating a deep anesthetic state. The wide divergence indicated by the shaded region 46 indicates a need to interpret this data in light of clinical judgement and other clinical endpoints. The sudden increase in E, occurred intraoperatively at a time when the patient had been administered doses of anesthetics generally considered adequate, and there was no sudden increase in surgical stimulation. For this reason, the increase in E, is most likely due to a burst of EMG
activity arising from an involuntary reflex movement, indicating that E, is the more accurate assessment of the patient's state of consciousness.
At point D, both estimators are trending upward, indicating a change in the patient to a lighter state of anesthesia. The absence of a divergence between the two
-14-estimators indicates that all the power in the 20-40 Hz region is coming from the EEG, and indicates a high level of confidence in the prediction.
At point E, the wide shaded divergence region 48 cautions the anesthesia provider to use clinical judgement to detenmine which estimator is more representative of the actual patient state. Both estimators have been trending up indicating an EEG
state tending toward increased probability of consciousness. Since this occurs at the end of the case, the divergence is most likely due to increased basal EMG
activity associated with returning consciousness and indicates that predictor E, more accurately reflects the state of the patient.
In this example, the differing sensitivity and specificity characteristics of the two estimators (predictors) which result from their differing operating points on their respective ROC curves are clearly demonstrated. The shaded divergence region directs the attention of the anesthetic provider to divergences between the patient state predicted by the different estimators and thus serves as an indicator of the level of confidence to be placed in the joint prediction and whether the provider needs to exercise extra clinical judgment in the interpretation.
The preceding embodiment utilizes two estimators (predictors) with different ROC curves to estimate the same parameter. An alternate embodiment utilizes two estimators (predictors) with identical ROC curves. The ROC curve for such an embodiment is shown in Figure 7. The operating points may be chosen so that one operating point (O.P. A) provides a high sensitivity at the expense of lower specificity, while the second operating point (O.P. B) provides high specificity at the expense of lower sensitivity. It should be noted that such an alternate embodiment, utilizing two estimators (predictors) with identical ROC curves, is functionally equivalent to the use of two different operating points on the ROC curve of a single estimator (predictor).
Such a system could be implemented using a single estimator by using two different
At point E, the wide shaded divergence region 48 cautions the anesthesia provider to use clinical judgement to detenmine which estimator is more representative of the actual patient state. Both estimators have been trending up indicating an EEG
state tending toward increased probability of consciousness. Since this occurs at the end of the case, the divergence is most likely due to increased basal EMG
activity associated with returning consciousness and indicates that predictor E, more accurately reflects the state of the patient.
In this example, the differing sensitivity and specificity characteristics of the two estimators (predictors) which result from their differing operating points on their respective ROC curves are clearly demonstrated. The shaded divergence region directs the attention of the anesthetic provider to divergences between the patient state predicted by the different estimators and thus serves as an indicator of the level of confidence to be placed in the joint prediction and whether the provider needs to exercise extra clinical judgment in the interpretation.
The preceding embodiment utilizes two estimators (predictors) with different ROC curves to estimate the same parameter. An alternate embodiment utilizes two estimators (predictors) with identical ROC curves. The ROC curve for such an embodiment is shown in Figure 7. The operating points may be chosen so that one operating point (O.P. A) provides a high sensitivity at the expense of lower specificity, while the second operating point (O.P. B) provides high specificity at the expense of lower sensitivity. It should be noted that such an alternate embodiment, utilizing two estimators (predictors) with identical ROC curves, is functionally equivalent to the use of two different operating points on the ROC curve of a single estimator (predictor).
Such a system could be implemented using a single estimator by using two different
-15-values of the estimator as thresholds defining a positive or negative prediction, corresponding to the different operating points.
While the preferred embodiment described above incorporates 2 estimators, the invention is easily extended to an alternate embodiment of N estimators, where N is a number greater than 2. In this case, additional estimators are incorporated which are designed to predict specific EEG or consciousness states or specific types or classes or signal artifact. For example, one estimator might be designed to detect the presence of low frequency EEG activity reflective of deep anesthetic states with a low probability of consciousness. A second might be designed to detect the presence of high frequency EEG reflective of lighter anesthetic states. A third might be designed to detect the presence of high frequency noise resulting from external electromagnetic interference. And a fourth might be designed to detect the presence of eye movement artifact. Each of these predictors would be designed to have similar values in the absence of the states or signals they are designed to detect. They would also be designed to have roughly constant divergences in response to similar signal situations across a broad population of subjects. Divergences may be indicated between each pair of estimators, or between selected estimators only.
Additionally, it should be recognized that many methods in addition to power spectral analysis may be used to compute estimators. These include, but- are not limited to, bispectral and higher order spectral analysis, time domain analysis methods, and stochastic and non-parametric modeling. This invention is not intended to be limited to any of the analysis methods described in the preferred or alternate embodiments.
The utility of estimating and displaying a range of predicted probabilities for a specific endpoint (e.g., consciousness) bounded by two different operating points is useful in the field of cardiology as well. For example, early detection of significant
While the preferred embodiment described above incorporates 2 estimators, the invention is easily extended to an alternate embodiment of N estimators, where N is a number greater than 2. In this case, additional estimators are incorporated which are designed to predict specific EEG or consciousness states or specific types or classes or signal artifact. For example, one estimator might be designed to detect the presence of low frequency EEG activity reflective of deep anesthetic states with a low probability of consciousness. A second might be designed to detect the presence of high frequency EEG reflective of lighter anesthetic states. A third might be designed to detect the presence of high frequency noise resulting from external electromagnetic interference. And a fourth might be designed to detect the presence of eye movement artifact. Each of these predictors would be designed to have similar values in the absence of the states or signals they are designed to detect. They would also be designed to have roughly constant divergences in response to similar signal situations across a broad population of subjects. Divergences may be indicated between each pair of estimators, or between selected estimators only.
Additionally, it should be recognized that many methods in addition to power spectral analysis may be used to compute estimators. These include, but- are not limited to, bispectral and higher order spectral analysis, time domain analysis methods, and stochastic and non-parametric modeling. This invention is not intended to be limited to any of the analysis methods described in the preferred or alternate embodiments.
The utility of estimating and displaying a range of predicted probabilities for a specific endpoint (e.g., consciousness) bounded by two different operating points is useful in the field of cardiology as well. For example, early detection of significant
-16-. coronary artery disease (CAD) in a patient is useful in minimizing the risk of a heart attack (i.e., early diagnosis enables early therapy or change in lifestyle that may inhibit the progression of heart disease). A cardiac index which predicts the presence of significant CAD can be constructed from parameters or features derived from the QRS
complex, the primary triangular component of the surface ECG waveform. For example, a set of features which may predict the presence of CAD include QRS
duration, average QRS bicoherence (600 < fl+f2 < 1000 Hz) and QRS power. Using a database of ECGs from healthy and diseased patients, a logistic model of a linear combination of these same features can be constructed which predicts the presence of significant CAD. An alternate embodiment of this invention, therefore, utilizes two models using these features, one model tuned for high sensitivity, the other tuned for high specificity. Subsequently, new patients can be evaluated for the presence of significant CAD by evaluating both predictive models. The result would be a lower and upper bound on the probability of CAD for the patient. In many instances, the range would be tight (e.g., for older men, or young women). In other cases, the range would be wide (e.g., for middle-aged women). The level of uncertainty would be reflected in the range of predicted probabilities. The estimation of the probability of CAD using each of the two models would produce a single value for each estimator, so these estimators (probabilities) would not be trended in time. However, the level of uncertainty indicated by the divergence between the estimators may be represented graphically by shading the region between them on a display device.
The notion of providing a range of predicted probabilities may be extended from using one estimator evaluated at two different operating points to using two different estimators (using different features, for example) in which one predictor would be highly sensitive while the other predictor would be highly specific.
For example, in the application of a set of predictors of consciousness, an accurate eye blink detector derived from the electrooculogram (EOG) would be a very specific detector for consciousness, but would be fairly insensitive. That is, the presence of
complex, the primary triangular component of the surface ECG waveform. For example, a set of features which may predict the presence of CAD include QRS
duration, average QRS bicoherence (600 < fl+f2 < 1000 Hz) and QRS power. Using a database of ECGs from healthy and diseased patients, a logistic model of a linear combination of these same features can be constructed which predicts the presence of significant CAD. An alternate embodiment of this invention, therefore, utilizes two models using these features, one model tuned for high sensitivity, the other tuned for high specificity. Subsequently, new patients can be evaluated for the presence of significant CAD by evaluating both predictive models. The result would be a lower and upper bound on the probability of CAD for the patient. In many instances, the range would be tight (e.g., for older men, or young women). In other cases, the range would be wide (e.g., for middle-aged women). The level of uncertainty would be reflected in the range of predicted probabilities. The estimation of the probability of CAD using each of the two models would produce a single value for each estimator, so these estimators (probabilities) would not be trended in time. However, the level of uncertainty indicated by the divergence between the estimators may be represented graphically by shading the region between them on a display device.
The notion of providing a range of predicted probabilities may be extended from using one estimator evaluated at two different operating points to using two different estimators (using different features, for example) in which one predictor would be highly sensitive while the other predictor would be highly specific.
For example, in the application of a set of predictors of consciousness, an accurate eye blink detector derived from the electrooculogram (EOG) would be a very specific detector for consciousness, but would be fairly insensitive. That is, the presence of
-17-eye blinks indicates consciousness; however, people can be conscious without creating eye blinks. In contrast, the previously mentioned EEG derived estimator E, is highly sensitive to conscious states, but is less specific. Separate estimates of the probability of consciousness derived from the EEG and EOG would provide a upper and lower bounds of the probability of consciousness in a manner similar to the two EEG
derived predictors described earlier.
While the foregoing invention has been described with reference to its preferred embodiments, various alterations and modifications will occur to those skilled in the art. All such alterations and modifications are intended to fall within the scope of the appended claims.
derived predictors described earlier.
While the foregoing invention has been described with reference to its preferred embodiments, various alterations and modifications will occur to those skilled in the art. All such alterations and modifications are intended to fall within the scope of the appended claims.
Claims (28)
1. A system for evaluating values of physiological parameters, said system comprising:
a signal acquisition device for obtaining physiological signals from a human body;
a processor for deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter and for computing relationships between said estimators;
a display l'or displaying information concerning said relationships between said two or more estimators.
a signal acquisition device for obtaining physiological signals from a human body;
a processor for deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter and for computing relationships between said estimators;
a display l'or displaying information concerning said relationships between said two or more estimators.
2. The system for evaluating values of physiological parameters of claim 1, wherein one estimator is derived using one physical or computational method and wherein at least one other estimator is derived using a different physical or computational method.
3. The system for evaluating values of physiological parameters of claim 1, wherein said relationships computed by said processor are instantaneous differences between the values of each estimator.
4. The system for evaluating values of physiological parameters of claim 1, in which said processor uses an instantaneous value of a relationship between estimators to determine the width of an identifiable region about, below, or above an instantaneous value of the first estimator displayed by said display.
The system for evaluating values of physiological parameters of claim 1, in which said processor uses an instantaneous value of a relationship between estimators to determine graphical attributes of Instantaneous values ol a first estimator displayed by said display
6. The system for evaluating values of physiological parameters of claim 1, in which said processor uses a measure independent of a first estimator of a parameter of interest to determine the width of an identifiable region about, below or above a first estimator displayed by said display.
7. The system for evaluating values of physiological parameters of claim 1, in which said processor uses a measure independent of a first estimator of a parameter of interest to determine graphical attributes of instantaneous values of a first estimator displayed by said display.
8. A system for evaluating values of physiological parameters, said system comprising:
a signal acquisition device for obtaining physiological signals from a human body;
a processor for deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter;
a display for displaying said two or more estimators.
a signal acquisition device for obtaining physiological signals from a human body;
a processor for deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter;
a display for displaying said two or more estimators.
9. The system for evaluating values of physiological parameters of claim 8.
wherein one estimator is derived using one physical or computational method and wherein at least one other estimator is derived using a different physical or computational method.
wherein one estimator is derived using one physical or computational method and wherein at least one other estimator is derived using a different physical or computational method.
The system for evaluating values ol physiological parameters of claim 8.
wherein said processor also computes relationships between said estimators and wherein said relationships are instantaneous differences between the values of each estimator.
wherein said processor also computes relationships between said estimators and wherein said relationships are instantaneous differences between the values of each estimator.
11. The system for evaluating values of physiological parameters of claim 8, in which said processor uses an instantaneous value of a relationship between estimators to determine the width of an identifiable region about, below, or above an instantaneous value of the first estimator displayed by said display.
12. The system for evaluating values of physiological parameters of claim 8, in which said processor uses an instantaneous value of a relationship between estimators to determine graphical attributes of instantaneous values of a first estimator displayed by said display.
13. The system for evaluating values of physiological parameters of claim 8, in which said processor uses a measure independent of a first estimator of a parameter of interest to determine the width of an identifiable region about, below or above a first estimator displayed by said display.
14. The system for evaluating values of physiological parameters of claim 8, in which said processor uses a measure independent of a first estimator of a parameter of interest to determine graphical attributes of instantaneous values of a first estimator displayed by said display.
15. A method for evaluating values of physiological parameters, said method comprising the steps of:
obtaining physiological signals from a human body;
deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter and computing relationships between said estimators;
displaying information concerning said relationships between said two or more estimators
obtaining physiological signals from a human body;
deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter and computing relationships between said estimators;
displaying information concerning said relationships between said two or more estimators
16. The method for evaluating values of physiological parameters of claim 15, wherein one estimator is derived using one physical or computational method and wherein at least one other estimator is derived using a different physical or computational method
17. The method for evaluating values of physiological parameters of claim 15, wherein said relationships are instantaneous differences between different estimators and instantaneous values of each estimator.
18. The method for evaluating values of physiological parameters of claim 15, in which an instantaneous value of a relationship between estimators is used to determine the width of an identifiable region about, below, or above an instantaneous value of the first estimator displayed by said display.
19. The method for evaluating values of physiological parameters of claim 15, in which an instantaneous value of a second estimator or a relationship between estimators is used to determine graphical attributes of an instantaneous values of a first estimator displayed by said display.
20. The method for evaluating values of physiological parameters of claim 15, in which a measure independent of a first estimator of a parameter of interest is used to determine the width of an identifiable region about, below or above a first estimator displayed by said display
21 The method for evaluating values of physiological parameters of claim 15, in which a measure independent of a first estimator of a parameter of interest is used to determine graphical attributes of instantaneous values of a first estimator displayed by said display.
22. A method for evaluating values of physiological parameters, said method comprising the steps of:
obtaining physiological signals from a human body;
deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter;
displaying said two or more estimators.
obtaining physiological signals from a human body;
deriving from said physiological signal values of two or more estimators of a pre-selected physiological parameter;
displaying said two or more estimators.
23. The method for evaluating values of physiological parameters of claim 22, wherein one estimator is derived using one physical or computational method and wherein at least one other estimator is derived using a different physical or computational method.
24. The method for evaluating values of physiological parameters of claim 22, further comprising the steps of computing relationships between said estimators and wherein said relationships are instantaneous differences between different estimators and instantaneous values of each estimator.
25. The method for evaluating values of physiological parameters of claim 22, in which an instantaneous value of a relationship between estimators is used to determine the width of an identifiable region about, below, or above an instantaneous value of the first estimator displayed by said display.
26. The method for evaluating values of physiological parameters of claim 22, in which an instantaneous value of a relationship between estimators is used to determine graphical attributes of an instantaneous values of a first estimator displayed by said display.
27. The method for evaluating values of physiological parameters of claim 22, in which a measure independent of a first estimator of a parameter of interest is used to determine the width of an identifiable region about, below or above a first estimator displayed by said display.
28. The method for evaluating values of physiological parameters of claim 22, in which a measure independent of a first estimator of a parameter of interest is used to determine graphical attributes of instantaneous values of a first estimator displayed by said display.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US09/135,931 | 1998-08-18 | ||
| US09/135,931 US6338713B1 (en) | 1998-08-18 | 1998-08-18 | System and method for facilitating clinical decision making |
| PCT/US1999/017688 WO2000011569A1 (en) | 1998-08-18 | 1999-08-04 | System and method for facilitating clinical decision making |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2338843A1 CA2338843A1 (en) | 2000-03-02 |
| CA2338843C true CA2338843C (en) | 2008-12-02 |
Family
ID=22470442
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002338843A Expired - Fee Related CA2338843C (en) | 1998-08-18 | 1999-08-04 | System and method for facilitating clinical decision making |
Country Status (8)
| Country | Link |
|---|---|
| US (1) | US6338713B1 (en) |
| EP (1) | EP1105812B1 (en) |
| JP (1) | JP4573149B2 (en) |
| AT (1) | ATE230863T1 (en) |
| AU (1) | AU773433B2 (en) |
| CA (1) | CA2338843C (en) |
| DE (1) | DE69904836T2 (en) |
| WO (1) | WO2000011569A1 (en) |
Families Citing this family (70)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19638738B4 (en) * | 1996-09-10 | 2006-10-05 | Bundesrepublik Deutschland, vertr. d. d. Bundesministerium für Wirtschaft und Technologie, dieses vertr. d. d. Präsidenten der Physikalisch-Technischen Bundesanstalt | Method for obtaining a diagnostic statement from signals and data of medical sensor systems |
| CA2404262C (en) * | 2000-03-29 | 2009-03-24 | University Of Virginia Patent Foundation | Method, system, and computer program product for the evaluation of glycemic control in diabetes from self-monitoring data |
| JP2004528935A (en) * | 2001-06-13 | 2004-09-24 | シーケーエム ダイアグノスティックス，インコーポレーテッド | Non-invasive detection method and detection device for tissue |
| US20030149597A1 (en) * | 2002-01-10 | 2003-08-07 | Zaleski John R. | System for supporting clinical decision-making |
| US20040122787A1 (en) * | 2002-12-18 | 2004-06-24 | Avinash Gopal B. | Enhanced computer-assisted medical data processing system and method |
| US20040122719A1 (en) * | 2002-12-18 | 2004-06-24 | Sabol John M. | Medical resource processing system and method utilizing multiple resource type data |
| US20040122705A1 (en) * | 2002-12-18 | 2004-06-24 | Sabol John M. | Multilevel integrated medical knowledge base system and method |
| US7187790B2 (en) * | 2002-12-18 | 2007-03-06 | Ge Medical Systems Global Technology Company, Llc | Data processing and feedback method and system |
| US20040122709A1 (en) * | 2002-12-18 | 2004-06-24 | Avinash Gopal B. | Medical procedure prioritization system and method utilizing integrated knowledge base |
| US20040122707A1 (en) * | 2002-12-18 | 2004-06-24 | Sabol John M. | Patient-driven medical data processing system and method |
| US20040122706A1 (en) * | 2002-12-18 | 2004-06-24 | Walker Matthew J. | Patient data acquisition system and method |
| US7490085B2 (en) * | 2002-12-18 | 2009-02-10 | Ge Medical Systems Global Technology Company, Llc | Computer-assisted data processing system and method incorporating automated learning |
| US20040122703A1 (en) * | 2002-12-19 | 2004-06-24 | Walker Matthew J. | Medical data operating model development system and method |
| US7367949B2 (en) * | 2003-07-07 | 2008-05-06 | Instrumentarium Corp. | Method and apparatus based on combination of physiological parameters for assessment of analgesia during anesthesia or sedation |
| JP4435092B2 (en) * | 2003-09-26 | 2010-03-17 | 東芝モバイルディスプレイ株式会社 | LCD panel |
| US20060085048A1 (en) * | 2004-10-20 | 2006-04-20 | Nervonix, Inc. | Algorithms for an active electrode, bioimpedance-based tissue discrimination system |
| WO2006044868A1 (en) * | 2004-10-20 | 2006-04-27 | Nervonix, Inc. | An active electrode, bio-impedance based, tissue discrimination system and methods and use |
| US7625340B2 (en) * | 2004-12-02 | 2009-12-01 | Instrumentarium Corporation | Identification of a dominant signal component in a biosignal |
| US20060129324A1 (en) * | 2004-12-15 | 2006-06-15 | Biogenesys, Inc. | Use of quantitative EEG (QEEG) alone and/or other imaging technology and/or in combination with genomics and/or proteomics and/or biochemical analysis and/or other diagnostic modalities, and CART and/or AI and/or statistical and/or other mathematical analysis methods for improved medical and other diagnosis, psychiatric and other disease treatment, and also for veracity verification and/or lie detection applications. |
| WO2006116529A2 (en) * | 2005-04-28 | 2006-11-02 | Katalytik, Inc. | System and method for managing healthcare work flow |
| JP5584413B2 (en) * | 2005-06-22 | 2014-09-03 | コーニンクレッカ フィリップス エヌ ヴェ | Patient monitoring system and monitoring method |
| JP5203215B2 (en) * | 2005-12-01 | 2013-06-05 | ネバ ヘルス エルエルシー | System and method for analyzing and evaluating depression and other mood disorders using electroencephalogram (EEG) measurements |
| JP5292104B2 (en) * | 2006-01-05 | 2013-09-18 | ユニバーシティ オブ バージニア パテント ファウンデーション | Computer-implemented method, system, and computer program for evaluating blood glucose variability in diabetes from self-monitoring data |
| US7676446B2 (en) * | 2006-01-11 | 2010-03-09 | Decision Command, Inc. | System and method for making decisions |
| US8442932B2 (en) * | 2006-01-11 | 2013-05-14 | Decision Command, Inc. | System and method for making decisions |
| US20070276197A1 (en) * | 2006-05-24 | 2007-11-29 | Lifescan, Inc. | Systems and methods for providing individualized disease management |
| US7751340B2 (en) * | 2006-11-03 | 2010-07-06 | Microsoft Corporation | Management of incoming information |
| US20080154513A1 (en) * | 2006-12-21 | 2008-06-26 | University Of Virginia Patent Foundation | Systems, Methods and Computer Program Codes for Recognition of Patterns of Hyperglycemia and Hypoglycemia, Increased Glucose Variability, and Ineffective Self-Monitoring in Diabetes |
| US20090012812A1 (en) * | 2007-03-06 | 2009-01-08 | Tracy Rausch | System and method for patient care |
| US20080235053A1 (en) * | 2007-03-20 | 2008-09-25 | Pinaki Ray | Communication medium for diabetes management |
| US20080234943A1 (en) * | 2007-03-20 | 2008-09-25 | Pinaki Ray | Computer program for diabetes management |
| US8758245B2 (en) * | 2007-03-20 | 2014-06-24 | Lifescan, Inc. | Systems and methods for pattern recognition in diabetes management |
| TWI324918B (en) * | 2007-04-04 | 2010-05-21 | Ind Tech Res Inst | Monitoring apparatus, system and method |
| US20090240127A1 (en) * | 2008-03-20 | 2009-09-24 | Lifescan, Inc. | Methods of determining pre or post meal time slots or intervals in diabetes management |
| US11382571B2 (en) | 2008-10-29 | 2022-07-12 | Flashback Technologies, Inc. | Noninvasive predictive and/or estimative blood pressure monitoring |
| US11395594B2 (en) | 2008-10-29 | 2022-07-26 | Flashback Technologies, Inc. | Noninvasive monitoring for fluid resuscitation |
| US11395634B2 (en) | 2008-10-29 | 2022-07-26 | Flashback Technologies, Inc. | Estimating physiological states based on changes in CRI |
| US20110172545A1 (en) * | 2008-10-29 | 2011-07-14 | Gregory Zlatko Grudic | Active Physical Perturbations to Enhance Intelligent Medical Monitoring |
| US11478190B2 (en) | 2008-10-29 | 2022-10-25 | Flashback Technologies, Inc. | Noninvasive hydration monitoring |
| US11857293B2 (en) | 2008-10-29 | 2024-01-02 | Flashback Technologies, Inc. | Rapid detection of bleeding before, during, and after fluid resuscitation |
| US8512260B2 (en) | 2008-10-29 | 2013-08-20 | The Regents Of The University Of Colorado, A Body Corporate | Statistical, noninvasive measurement of intracranial pressure |
| US11406269B2 (en) | 2008-10-29 | 2022-08-09 | Flashback Technologies, Inc. | Rapid detection of bleeding following injury |
| US8364254B2 (en) * | 2009-01-28 | 2013-01-29 | Brainscope Company, Inc. | Method and device for probabilistic objective assessment of brain function |
| JP2012532323A (en) * | 2009-06-30 | 2012-12-13 | ライフスキャン・インコーポレイテッド | Analyte testing method and system |
| WO2011025953A1 (en) * | 2009-08-28 | 2011-03-03 | Lexicor Medical Technology, Llc | Systems and methods to identify a subgroup of adhd at higher risk for complicating conditions |
| CA2762418A1 (en) | 2009-09-01 | 2011-03-10 | Exxonmobil Upstream Research Company | Method of using human physiological responses as inputs to hydrocarbon management decisions |
| WO2011038231A2 (en) * | 2009-09-25 | 2011-03-31 | Med-El Elektromedizinische Geraete Gmbh | Hearing implant fitting |
| DE102009053256A1 (en) * | 2009-11-06 | 2011-05-19 | Baars, Jan H., Dr. med. | Method for determining the analgesic level of a sedated or anaesthetized individual |
| US20130310422A1 (en) | 2010-09-01 | 2013-11-21 | The General Hospital Corporation | Reversal of general anesthesia by administration of methylphenidate, amphetamine, modafinil, amantadine, and/or caffeine |
| EP2621333B1 (en) | 2010-09-28 | 2015-07-29 | Masimo Corporation | Depth of consciousness monitor including oximeter |
| US9615781B2 (en) | 2011-02-03 | 2017-04-11 | Covidien Lp | Systems and methods for monitoring depth of consciousness |
| EP2734103B1 (en) | 2011-07-22 | 2020-12-23 | Flashback Technologies, Inc. | Hemodynamic reserve monitor and hemodialysis control |
| US20140081659A1 (en) | 2012-09-17 | 2014-03-20 | Depuy Orthopaedics, Inc. | Systems and methods for surgical and interventional planning, support, post-operative follow-up, and functional recovery tracking |
| US10383580B2 (en) | 2012-12-31 | 2019-08-20 | Abbott Diabetes Care Inc. | Analysis of glucose median, variability, and hypoglycemia risk for therapy guidance |
| US9351670B2 (en) | 2012-12-31 | 2016-05-31 | Abbott Diabetes Care Inc. | Glycemic risk determination based on variability of glucose levels |
| WO2014145335A1 (en) | 2013-03-15 | 2014-09-18 | Abbott Diabetes Care Inc. | System and method to manage diabetes based on glucose median, glucose variability, and hypoglycemic risk |
| ES2877551T3 (en) | 2013-04-24 | 2021-11-17 | Fresenius Kabi Deutschland Gmbh | Operating procedure of a control device to control an infusion device |
| WO2014176441A1 (en) * | 2013-04-24 | 2014-10-30 | The General Hospital Corporation | System and method for monitoring level of dexmedatomidine-induced sedation |
| US10314503B2 (en) | 2013-06-27 | 2019-06-11 | The General Hospital Corporation | Systems and methods for tracking non-stationary spectral structure and dynamics in physiological data |
| WO2014210527A1 (en) | 2013-06-28 | 2014-12-31 | The General Hospital Corporation | System and method to infer brain state during burst suppression |
| JP6586093B2 (en) | 2013-09-13 | 2019-10-02 | ザ ジェネラル ホスピタル コーポレイション | System for improved brain monitoring during general anesthesia and sedation |
| WO2016057553A1 (en) | 2014-10-07 | 2016-04-14 | Masimo Corporation | Modular physiological sensors |
| US10786168B2 (en) | 2016-11-29 | 2020-09-29 | The General Hospital Corporation | Systems and methods for analyzing electrophysiological data from patients undergoing medical treatments |
| US11723579B2 (en) | 2017-09-19 | 2023-08-15 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement |
| US11717686B2 (en) | 2017-12-04 | 2023-08-08 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to facilitate learning and performance |
| US11478603B2 (en) | 2017-12-31 | 2022-10-25 | Neuroenhancement Lab, LLC | Method and apparatus for neuroenhancement to enhance emotional response |
| US11364361B2 (en) | 2018-04-20 | 2022-06-21 | Neuroenhancement Lab, LLC | System and method for inducing sleep by transplanting mental states |
| EP3849410A4 (en) | 2018-09-14 | 2022-11-02 | Neuroenhancement Lab, LLC | System and method of improving sleep |
| US11918386B2 (en) | 2018-12-26 | 2024-03-05 | Flashback Technologies, Inc. | Device-based maneuver and activity state-based physiologic status monitoring |
| US11786694B2 (en) | 2019-05-24 | 2023-10-17 | NeuroLight, Inc. | Device, method, and app for facilitating sleep |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4517986A (en) * | 1984-03-05 | 1985-05-21 | Bilgutay Ilhan M | Four function vital sign monitor |
| JPH069546B2 (en) * | 1990-11-30 | 1994-02-09 | 株式会社コンピュータコンビニエンス | Diagnostic device using pulse wave and / or heartbeat collected from body surface |
| US5195531A (en) | 1991-03-01 | 1993-03-23 | Bennett Henry L | Anesthesia adequacy monitor and method |
| JP3369201B2 (en) * | 1991-10-02 | 2003-01-20 | マツダ株式会社 | Arousal maintenance device |
| US5458117A (en) | 1991-10-25 | 1995-10-17 | Aspect Medical Systems, Inc. | Cerebral biopotential analysis system and method |
| JPH07231880A (en) * | 1994-02-24 | 1995-09-05 | Sanyo Electric Co Ltd | Stress evaluation method and device therefor |
| US5590650A (en) * | 1994-11-16 | 1997-01-07 | Raven, Inc. | Non-invasive medical monitor system |
| JP3351160B2 (en) * | 1995-03-24 | 2002-11-25 | 松下電器産業株式会社 | Biological evaluation device |
| US5687733A (en) | 1995-10-26 | 1997-11-18 | Baxter International Inc. | System and method for estimating cardiac output |
| US5913308A (en) * | 1996-12-19 | 1999-06-22 | Hewlett-Packard Company | Apparatus and method for determining respiratory effort from muscle tremor information in ECG signals |
| US5924980A (en) * | 1998-03-11 | 1999-07-20 | Siemens Corporate Research, Inc. | Method and apparatus for adaptively reducing the level of noise in an acquired signal |
-
1998
- 1998-08-18 US US09/135,931 patent/US6338713B1/en not_active Expired - Lifetime
-
1999
- 1999-08-04 AU AU55474/99A patent/AU773433B2/en not_active Ceased
- 1999-08-04 AT AT99942007T patent/ATE230863T1/en active
- 1999-08-04 CA CA002338843A patent/CA2338843C/en not_active Expired - Fee Related
- 1999-08-04 WO PCT/US1999/017688 patent/WO2000011569A1/en active IP Right Grant
- 1999-08-04 DE DE69904836T patent/DE69904836T2/en not_active Expired - Lifetime
- 1999-08-04 EP EP99942007A patent/EP1105812B1/en not_active Expired - Lifetime
- 1999-08-04 JP JP2000566757A patent/JP4573149B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| DE69904836T2 (en) | 2003-11-13 |
| US6338713B1 (en) | 2002-01-15 |
| WO2000011569A1 (en) | 2000-03-02 |
| JP2002523163A (en) | 2002-07-30 |
| AU5547499A (en) | 2000-03-14 |
| AU773433B2 (en) | 2004-05-27 |
| EP1105812A1 (en) | 2001-06-13 |
| EP1105812B1 (en) | 2003-01-08 |
| CA2338843A1 (en) | 2000-03-02 |
| ATE230863T1 (en) | 2003-01-15 |
| DE69904836D1 (en) | 2003-02-13 |
| JP4573149B2 (en) | 2010-11-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2338843C (en) | System and method for facilitating clinical decision making | |
| Dolatabadi et al. | Automated diagnosis of coronary artery disease (CAD) patients using optimized SVM | |
| US11172835B2 (en) | Method and system for monitoring sleep | |
| AU2003206411B2 (en) | System and method of assessment of neurological conditions using EEG bispectrum | |
| US7309314B2 (en) | Method for predicting apnea-hypopnea index from overnight pulse oximetry readings | |
| US8155733B2 (en) | Lie detection method and system | |
| US20020143263A1 (en) | System and device for multi-scale analysis and representation of physiological data | |
| Lado et al. | Detecting sleep apnea by heart rate variability analysis: assessing the validity of databases and algorithms | |
| EP2203115A1 (en) | Electrocardiogram derived apnoea/hypopnea index | |
| Radhakrishna et al. | Nonlinear measures of heart rate time series: influence of posture and controlled breathing | |
| Hung | Detection of central sleep apnea based on a single-lead ECG | |
| WO2014091291A1 (en) | A device and method for determining the probability of response to pain and nociception of a subject t | |
| C. Dela Cruz et al. | Deriving heart rate and respiratory rate from ecg using wavelet transform | |
| Chase et al. | Quantifying agitation in sedated ICU patients using heart rate and blood pressure | |
| Griffel et al. | Path length entropy analysis of diastolic heart sounds | |
| Lin et al. | Real-time seismocardiogram feature extraction using adaptive gaussian mixture models | |
| Fischer et al. | Multi-and monofractal indices of short-term heart rate variability | |
| US20220142554A1 (en) | System and method for monitoring neural signals | |
| Ignácz et al. | NB-SQI: A novel non-binary signal quality index for continuous blood pressure waveforms | |
| US20160171170A1 (en) | Measuring respiration or other periodic physiological processes | |
| KR102554084B1 (en) | The assessment method of arterial stenosis by stethoscope sound analysis | |
| JP2003190109A (en) | Autonomic nervous system function evaluating method and system therefor | |
| Wesley | Approximate and Sample Entropy of Center of Pressure in Unperturbed Tandem Standing: Contribution of Embedding Dimension and Tolerance | |
| Yang | Recognition of sleep stages from sensor data | |
| EP2976016A1 (en) | Spectrum analysis of coronary artery turbulent blood flow |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| MKLA | Lapsed |
Effective date: 20190806 |