US 7,803,119 B2Grant
Respiration-based prognosis of heart disease
Issue Date:2010-09-28
•81 Claims
•8 Drawing Sheets
Abstract
A method for diagnosis includes receiving a signal associated with respiration of a patient during sleep, and processing the signal to detect occurrences of a pattern corresponding to multiple cycles of periodic breathing. A cumulative duration of the occurrences is computed, and a prognosis of the patient is determined based on the cumulative duration.
Metadata
Assignee
- Widemed Ltd.
Inventor
- Daniel Reisfeld
Application Information
Application Number:US 11/750,250
Filing Date:2007-05-17
Priority Date:2004-11-22
Art Unit:3735
Classifications
IPC:
A61B5/02A61B5/08
Patent Drawings (8 sheets)
Description
CROSS-REFERENCE TO RELATED APPLICATIONS
[0001] This application is a Continuation-In-Part of PCT patent application PCT/IL2005/001233, filed Nov. 22, 2005, which is a Continuation of Ser. No. 10/995,817, filed Nov. 22, 2004 and which claims the priority of U.S. Provisional Patent Application 60/843,107, filed Sep. 7, 2006, and is a Continuation-In-Part of PCT patent application PCT/IL2006/000148, filed Feb. 7, 2006, claiming the priority of U.S. Provisional Patent Application 60/651,295, filed Feb. 7, 2005. This application is related to three other U.S. patent applications, all filed on even date, which are entitled “Detection of Heart Failure Using a Photoplethysmograph,” “Respiration-Based Prognosis of Heart Disease,” and “Sleep Monitoring Using a Photoplethysmograph.” The disclosure of all of these related applications are incorporated herein by reference.
FIELD OF THE INVENTION
[0002] The present invention relates generally to physiological monitoring and diagnosis, and specifically to sleep recording and analysis.
BACKGROUND OF THE INVENTION
[0003] Human sleep is generally described as a succession of five recurring stages (plus waking, which is sometimes classified as a sixth stage). Sleep stages are typically monitored using a polysomnograph to collect physiological signals from the sleeping subject, including brain waves (EEG), eye movements (EOG), muscle activity (EMG), heartbeat (ECG), blood oxygen levels (SpO2) and respiration. The commonly-recognized stages include:
-
- [0004] Stage 1 sleep, or drowsiness. The eyes are closed during Stage 1 sleep, but if aroused from it, a person may feel as if he or she has not slept.
- [0005] Stage 2 is a period of light sleep, during which the body prepares to enter deep sleep.
- [0006] Stages 3 and 4 are deep sleep stages, with Stage 4 being more intense than Stage 3.
- [0007] Stage 5, REM (rapid eye movement) sleep, is distinguishable from non-REM (NREM) sleep by changes in physiological states, including its characteristic rapid eye movements.
Polysomnograms show brain wave patterns in REM to be similar to Stage 1 sleep. In normal sleep, heart rate and respiration speed up and become erratic, while the muscles may twitch. Intense dreaming occurs during REM sleep, but paralysis occurs simultaneously in the major voluntary muscle groups.
[0008] Sleep apneas commonly occur in conjunction with a variety of cardiorespiratory disorders. The relationship between sleep apnea and heart failure, for example, is surveyed by Bradley et al. in two articles entitled “Sleep Apnea and Heart Failure,” including “Part I: Obstructive Sleep Apnea,” Circulation 107, pages 1671-1678 (2003), and “Part II: Central Sleep Apnea,” Circulation 107, pages 1822-1826 (2003), which are incorporated herein by reference. The authors define “apnea” as a cessation of airflow for more than 10 sec. This term is distinguished from “hypopnea,” which is a reduction in but not complete cessation of airflow to less than 50% of normal, usually in association with a reduction in oxyhemoglobin saturation (commonly referred to as “desaturation”).
[0009] Sleep apneas and hypopneas are generally believed to fall into two categories: obstructive, due to collapse of the pharynx; and central, due to withdrawal of central respiratory drive to the muscles of respiration. Central sleep apnea (CSA) is commonly associated with Cheyne-Stokes respiration, which is a form of periodic breathing in which central apneas and hypopneas alternate with periods of hyperventilation, with a waxing-waning pattern of tidal volume. CSA is believed to arise as the result of heart failure, though obstructive sleep apnea (OSA) may also occur in heart failure patients.
[0010] Both OSA and CSA increase the strain on the cardiovascular system and thus worsen the prognosis of the heart failure patient. In some cases, both types of apneas may occur in the same patient, even at the same time (superposition). Classifying respiratory events as central or obstructive is considered to be a critical point, since treatment may differ according to the type of events, as pointed out by Pepin et al. in “Cheyne-Stokes Respiration with Central Sleep Apnea in Chronic Heart Failure: Proposals for a Diagnostic and Therapeutic Strategy,” Sleep Medicine Reviews 10, pages 33-47 (2006), which is incorporated herein by reference. Both CSA and OSA can be manifested in periodic breathing patterns.
[0011] Various methods have been proposed in the patent literature for automated apnea detection and diagnosis based on patient monitoring during sleep. For example, U.S. Patent Application Publication US 2004/0230105 A1 describes a method for analyzing respiratory signals using a Fuzzy Logic Decision Algorithm (FLDA). The method may be used to associate respiratory disorders with obstructive apnea, hypopnea, central apnea, or other conditions. As another example, U.S. Patent Application Publication US 2002/0002327 A1 and U.S. Pat. No. 6,839,581 describe methods for detecting Cheyne-Stokes respiration, which may be used on patients with heart failure. The methods involve performing spectral analysis of overnight oximetry recordings, from which a classification tree is generated. Another method, based on monitoring oxygen saturation and calculating the slope of desaturation events, is described in U.S. Pat. No. 6,760,608. Yet another method for classifying sleep apneas is described in U.S. Pat. No. 6,856,829. In this case, pulse waves from the body of a patient are detected, and the envelope of the pulse waves is created by connecting every peak of the pulse waves. The normalized amplitude and period of the envelope are used in determining whether the patient has OSA, CSA, or mixed sleep apnea syndrome. The disclosures of the patents and patent applications cited above are incorporated herein by reference.
[0012] It has been suggested that sleep monitoring can be used for assessing cardiorespiratory risk. For example, U.S. Pat. No. 5,902,250, whose disclosure is incorporated herein by reference, describes a home-based, wearable, self-contained system that determines sleep-state and respiratory pattern, and assesses cardiorespiratory risk. A respiratory disorder may be diagnosed from the frequency of eyelid movements and/or from ECG signals. Cardiac disorders (such as cardiac arrhythmia or myocardial ischemia) that are known to be linked to certain respiratory disorders also may be inferred upon detection of such respiratory disorders.
[0013] Photoplethysmograph devices, known commonly as pulse oximeters, provide instantaneous in vivo measurement of arterial oxygenation by determining the color of blood between a light source and a photodetector. To determine the blood oxygen saturation, light absorption measurement is carried out at two wavelengths in the red and infrared ranges. The difference between background absorption during diastole and peak absorption during systole at both wavelengths is used to compute the blood oxygen saturation.
[0014] Photoplethysmograph signals provide information not only on blood oxygenation, but also on other physiological signs. For example, U.S. Pat. No. 5,588,425 describes the use of a pulse oximeter in validating the heart rate and/or R—R intervals of an ECG, and for discriminating between sleep and wakefulness in a monitored subject. It also describes a method for distinguishing between valid pulse waveforms in the oximeter signal. U.S. Pat. No. 7,001,337 describes a method for obtaining physiological parameter information related to respiration rate, heart rate, heart rate variability, blood volume variability and/or the autonomic nervous system using photoplethysmography. U.S. Pat. No. 7,190,261 describes an arrhythmia alarm processor, which detects short-duration, intermittent oxygen desaturations of a patient using a pulse oximeter as a sign of irregular heartbeat. An alarm is triggered when the pattern of desaturations matches a reference pattern. The disclosures of the above-mentioned patents are incorporated herein by reference.
SUMMARY OF THE INVENTION
[0015] The photoplethysmograph signals that are output by a standard pulse oximeter can provide a wealth of information regarding the patient's vital signs and physiological condition. In embodiments of the present invention that are described hereinbelow, photoplethysmograph signals that are captured while the patient sleeps are analyzed in order to diagnose the patient's cardiorespiratory condition. In particular, the signals may be used to detect and assess the severity of conditions that are characteristic of heart failure (HF), such as premature ventricular contractions and Cheyne-Stokes breathing. The photoplethysmograph signals may also be used, even without monitoring other physiological parameters, to classify the sleep stages and “sleep quality” of the patient.
[0016] The power and versatility of the photoplethysmograph-based techniques that are described hereinbelow make it possible to monitor patients' oxygen saturation, heartbeat, respiration, sleep stages and autonomic nervous system during sleep using no more than a single pulse oximeter probe (which typically clips onto the patient's finger). As a result, the patient may be monitored comfortably and conveniently, at home or in a hospital bed, even without on-site assistance in setting up each night's monitoring.
[0017] In alternative embodiments, the principles of the present invention may be applied to analysis of respiration signals captured using monitors of other types.
[0018] There is therefore provided, in accordance with an embodiment of the present invention, a method for diagnosis, including:
[0019] receiving a signal associated with respiration of a patient during sleep;
[0020] processing the signal to detect occurrences of a pattern corresponding to multiple cycles of periodic breathing;
[0021] computing a cumulative duration of the occurrences; and
[0022] determining a prognosis of the patient based on the cumulative duration.
[0023] In disclosed embodiments, determining the prognosis includes determining the prognosis of heart failure (HF) in the patient. Determining the prognosis may include predicting a probability of survival of the patient or predicting a level of at least one marker, selected from a group of markers consisting of Brain Natriuretic Peptide (BNP) and N-terminal prohormone Brain Natriuretic Peptide (NT proBNP). Additionally or alternatively, determining the prognosis may include differentiating compensated from decompensated heart failure, predicting a hospitalization of the patient, or providing an indication of a need for a change in treatment of the patient.
[0024] In some embodiments, the signal includes a photoplethysmograph signal, and the method includes processing the photoplethysmograph signal so as to identify a cardiac arrhythmia of the patient. Typically, determining the prognosis includes finding the prognosis responsively both to the occurrences of the pattern and to the cardiac arrhythmia. Additionally or alternatively, the method may include processing the photoplethysmograph signal so as to measure a vasomodulation in the body or to measure a modulation of respiration by the patient.
[0025] In one embodiment, receiving the signal includes receiving the photoplethysmograph signal from a sensor on a finger of the patient, and processing the signal includes digitizing and filtering the signal so as to generate data using processing circuitry in a control unit that is fastened to a forearm of the patient, and storing the data in a memory in the control unit. Typically, processing the signal includes connecting the control unit to an external processor and uploading the data from the memory to the external processor for analysis, and storing the data includes receiving and recording the data in the control unit while the control unit is disconnected from the external processor.
[0026] Alternatively, receiving the signal includes receiving the photoplethysmograph signal from a sensor implanted in a body of the patient.
[0027] In a disclosed embodiment, treating the patient so as to reduce the cumulative duration.
[0028] Typically, determining the prognosis includes classifying the patient in a high-mortality group if the cumulative duration exceeds a predetermined threshold, wherein the threshold is at least 45 minutes.
[0029] In some embodiments, processing the signal includes identifying a periodic breathing event when an event duration of the multiple cycles is no less than a predetermined minimum duration or when a number of the multiple cycles is no less than a predetermined minimum number, such as three.
[0030] In disclosed embodiments, the signal is indicative of blood oxygen saturation, and processing the signal includes identifying a periodic breathing event when a representative level of the blood oxygen saturation during the cycles is less than the blood oxygen saturation before the cycles by at least a predetermined amount. Typically, the representative level of the blood oxygen saturation is one of a median saturation, a mean saturation, and a minimum saturation during the cycles, and the periodic breathing event is identified when the representative level of the blood oxygen saturation is less than the blood oxygen saturation before the cycles by at least 2%.
[0031] Typically, the periodic breathing includes Cheyne-Stokes breathing, and processing the signal includes distinguishing the Cheyne-Stokes breathing events from obstructive apnea events. In some embodiments, distinguishing the Cheyne-Stokes breathing events includes determining a slope of the signal, and distinguishing the Cheyne-Stokes breathing events responsively to the slope. Determining the slope may include finding a measure of the slope selected from a group of measures consisting of a median slope and a mean slope, and distinguishing the Cheyne-Stokes breathing events based on a maximal value of the measure of the slope in each of the cycles. In one embodiment, distinguishing the Cheyne-Stokes breathing events includes classifying a breathing event as Cheyne-Stokes only when the maximal value of the measure is no greater than 0.7%/sec.
[0032] In some embodiments, processing the signal includes extracting from the signal a measure of respiratory energy, and filtering the measure so as to eliminate signal components at frequencies equal to and greater than a respiratory frequency of the patient.
[0033] In an alternative embodiment, receiving the signal includes monitoring the respiration using a non-contact sensor.
[0034] There is also provided, in accordance with an embodiment of the present invention, apparatus for diagnosis, including:
[0035] a sensor, which is configured to be coupled to a body of a patient during sleep and to output a signal associated with respiration of the patient; and
[0036] a processor, which is coupled to process the signal to detect occurrences of a pattern corresponding to multiple cycles of periodic breathing, to compute a cumulative duration of the occurrences.
[0037] There is additionally provided, in accordance with an embodiment of the present invention, a computer software product, including a computer-readable medium in which program instructions are stored, which instructions, when read by a computer, cause the computer to receive a signal associated with respiration of a patient during sleep, and to process the signal to detect occurrences of a pattern corresponding to multiple cycles of periodic breathing, to compute and output a cumulative duration of the occurrences.
[0038] The present invention will be more fully understood from the following detailed description of the embodiments thereof, taken together with the drawings in which:
BRIEF DESCRIPTION OF THE DRAWINGS
[0039] FIG. 1 is a schematic, pictorial illustration of a system for sleep monitoring and diagnosis, in accordance with an embodiment of the present invention;
[0040] FIG. 2A is a schematic, pictorial illustration of apparatus for patient monitoring, in accordance with an embodiment of the present invention;
[0041] FIG. 2B is a block diagram that schematically shows functional elements of the apparatus of FIG. 2A , in accordance with an embodiment of the present invention;
[0042] FIGS. 3A and 3B are a flow chart that schematically illustrates a method for sleep monitoring and diagnosis, in accordance with an embodiment of the present invention;
[0043] FIG. 4 is a schematic plot of photoplethysmograph and ECG signals, illustrating detection of a cardiac arrhythmia in accordance with an embodiment of the present invention;
[0044] FIG. 5 is a Kaplan-Meier plot that schematically shows survival of heart failure patients as a function of Brain Natriuretic Peptide (BNP) levels;
[0045] FIG. 6 is a Kaplan-Meier plot that schematically shows survival of heart failure patients as a function of BNP levels, with a threshold level determined by severity of Cheyne-Stokes breathing, in accordance with an embodiment of the present invention;
[0046] FIG. 7 is a receiver operating characteristic (ROC) plot, which schematically shows the sensitivity and specificity of predicting heart failure prognosis in accordance with an embodiment of the present invention; and
[0047] FIG. 8 is a Kaplan-Meier plot that schematically shows survival of heart failure patients as a function of the severity of symptoms classified by a method in accordance with an embodiment of the present invention.
DETAILED DESCRIPTION OF EMBODIMENTS
System Overview
[0048] FIG. 1 is a schematic, pictorial illustration of a system 20 for sleep monitoring and diagnosis, in accordance with an embodiment of the present invention. In this embodiment, system 20 is used to monitor a patient 22 in a home, clinic or hospital ward environment, although the principles of the present invention may similarly be applied in dedicated sleep laboratories. System 20 receives and analyzes a photoplethysmograph signal from a suitable sensor, such as a pulse oximetry device 24. Device 24 provides a photoplethysmograph signal indicative of blood flow and a signal indicative of the level of oxygen saturation in the patient's blood. In the context of the present patent application and in the claims, the photoplethysmograph signal is thus considered to be a signal that is associated with blood oxygen saturation. Since the photoplethysmograph signal is modulated by both the heart rate and respiratory rate, it may also be used to provide a heart rate and respiratory rate signals. The sensor signals from device 24 are collected, digitized and processed by a console 28.
[0049] Optionally, system 20 may comprise sensors of other types (not shown), for collecting other physiological signals. For example, the system may receive an ECG signal, measured by skin electrodes, and a respiration signal measured by a respiration sensor. Additionally or alternatively, the techniques of monitoring and analysis that are described herein may be combined with EEG, EOG, leg motion sensors, and other sleep and/or cardiac monitoring modalities that are known in the art. As another example, console 28 may receive signals by telemetry from implantable cardiac devices, such as pacemakers and ICDs.
[0050] Console 28 may process and analyze the signals from pulse oximetry device 24 locally, using the methods described hereinbelow. In the present embodiment, however, console 28 is coupled to communicate over a network 30, such as a telephone network or the Internet, with a diagnostic processor 32. This configuration permits sleep studies to be performed simultaneously in multiple different locations. Processor 32 typically comprises a general-purpose computer processor (which may be embedded in a bedside or remote monitor) with suitable software for carrying out the functions described herein. This software may be downloaded to processor 32 in electronic form, or it may alternatively be provided on tangible media, such as optical, magnetic or non-volatile electronic memory. Processor 32 analyzes the signals conveyed by console 28 in order to analyze the physiological parameters, identify sleep stages, and extract prognostic information regarding patient 22, and to display the results of the analysis to an operator 34, such as a physician.
[0051] Alternatively, although the embodiments described herein relate mainly to methods and apparatus for monitoring and diagnosis during sleep, the principles of the present invention may also be applied, mutatis mutandis, to patients who are awake. In particular, these methods and apparatus may be used in monitoring patients who are reclining or otherwise at rest, even if they are not asleep.
[0052] Reference is now made to FIGS. 2A and 2B , which schematically illustrate apparatus 21 for patient monitoring, in accordance with an embodiment of the present invention. FIG. 2A is a pictorial illustration of the apparatus, while FIG. 2B is a block diagram showing functional components of the apparatus. Apparatus 21 is similar in functionality to elements of system 20, as shown in FIG. 1 , but apparatus 21 is particularly advantageous in that it can be worn comfortably by the patient during both sleep and waking hours and requires no wired connection to a console, except periodically (once a day, for example) for data upload and battery recharging.
[0053] As shown in FIG. 2A , pulse oximetry device 24 in apparatus 21 has the form of a ring, which fits comfortably over one of the fingers on a hand 23 of the patient (although other configurations of device 24 may alternatively be used in the apparatus). Device 24 is connected by a wire to a control unit 25, which may conveniently be fastened around the patient's wrist. Alternatively, the control unit may be fastened elsewhere on the patient's forearm, at any suitable location between the hand and the elbow, or elsewhere on the patient's body. The control unit may include a display 27, to present status information and/or readings of monitored parameters, such as heart rate, blood oxygen saturation and heart failure status. A connector 29 on the control unit is configured to connect to a console or docking station. In the illustrated embodiment, connector 29 comprises a receptacle for a cable with a standard plug, such as a USB cable. Alternatively, the connector may mate directly with a matching connector on a dedicated docking station.
[0054] As shown in FIG. 2B , a sensor 31 (typically comprising two light source/light detector subassemblies, as described below) in device 24 is connected via wire to signal processing circuitry 33 in control unit 25. The signal processing circuitry digitizes and filters the signals from sensor 31 and stores the results in a memory 35. (Alternatively or additionally, control unit 25 may transmit the results to a receiver using a suitable wireless communication protocol, such as Bluetooth® or ZigBee®). Optionally, circuitry 33 may also be configured to perform some of the additional processing functions that are shown in FIGS. 3A and 3B and described hereinbelow. The signal processing circuitry and peripheral components are powered by an internal power source, such as a battery 38, so that apparatus 21 can perform its data collection functions without wired connection to a console or to lines power.
[0055] Apparatus 21 may also comprise an actigraph 39, which is typically contained in control unit 25. The actigraph measures movement of the patient and typically comprises an accelerometer for this purpose. The measurements of patient movement are recorded together with the data from sensor 31 in memory 35 and may be used in subsequent analysis to determine the patient's state of sleep or arousal.
[0056] After apparatus 21 has recorded patient data in memory 35 for a sufficient period of time, the user (who may be the patient himself or herself) connects control unit 25 to the docking station or other console via connector 29. A controller 36 in the control unit is then able to communicate with the console or docking station via a suitable interface 37 (such as a USB interface in the example noted above). The controller reads out the data that are stored in memory 35 to a processor, such as processor 32, which analyzes the data, as described hereinbelow. In addition, interface 37 may comprise charging circuitry for recharging battery 38.
[0057] In the embodiments that are described below, pulse oximetry device 24 may be configured either as shown in FIG. 1 or as shown in FIG. 2A or 2B. Alternatively, substantially any suitable sort of photoplethysmograph may be used in these embodiments, including a photoplethysmographic sensor that is implanted in the body of the patient. An implantable oximeter that may be used for this purpose, for example, is described in U.S. Pat. No. 6,122,536, whose disclosure is incorporated herein by reference. Furthermore, the methods that are described hereinbelow may be used in conjunction with devices of other types that provide information on the breathing, oxygen saturation, and heart performances of the patient. In one embodiment, for instance, the methods described below are applied to the output of a non-contact respiratory monitor, such as the one described in U.S. Pat. No. 6,011,477, whose disclosure is incorporated herein by reference.
Diagnostic Method
[0058] FIGS. 3A and 3B are a flow chart that schematically illustrates a method for sleep monitoring and diagnosis, in accordance with an embodiment of the present invention. Pulse oximetry device 24 comprises two light source/light detector subassemblies 40 and 42. These subassemblies generate signals that are indicative of absorption and/or reflectance of light at two different wavelengths, typically one red and one infrared, as is known in the art. Each signal includes an AC component, which corresponds to the pulsatile change in the signal at the patient's heart rate, and a slow-changing DC component. Comparison of the AC components of the two signals gives a blood oxygen saturation signal 44. Alternatively, at least some of the methods described below can use the signals from only a single source/detector subassembly, or signals provided by other types of photoplethysmographic sensors.
[0059] The saturation signal is low-pass filtered to give a very-low-frequency (VLF) saturation signal 46. This filtering removes signal components at frequencies that are greater than or equal to the patient's respiratory frequency, so that the signal remaining reflects trends over multiple respiratory cycles. In some embodiments, the filtering is even more pronounced, and eliminates frequency components outside the Cheyne-Stokes cycle frequency, for example, components below 1/180 Hz or above 1/40 Hz.
[0060] Processor 32 analyzes shape characteristics of the VLF saturation signal in order to detect episodes of Cheyne-Stokes breathing (CSB). As noted earlier, this condition is characterized by a regular waxing and waning breathing pattern and occurs particularly among patients with heart failure and in patients who have experienced a stroke. CSB is present during sleep, and in more severe cases may also be observed during wakefulness. According to the American Academy of Sleep Medicine, Cheyne-Strokes breathing syndrome (CSBS) is characterized by the following criteria:
-
- [0061] 1. Presence of congestive heart failure or cerebral neurological disease.
- [0062] 2. Respiratory monitoring demonstrates:
- [0063] a. At least three consecutive cycles of a cyclical crescendo and decrescendo change in breathing amplitude. Cycle length is most commonly in the range of 60 seconds, although the length may vary.
- [0064] b. One or both of the following:
- [0065] i. Five or more central sleep apneas or hypopneas per hour of sleep.
- [0066] ii. The cyclic crescendo and decrescendo change in breathing amplitude has duration of at least 10 consecutive minutes.
[0067] The inventors have found the typical Cheyne-Stokes cycle length to be between 40 and 90 sec. The decrescendo phase is associated with decreased respiratory effort and rate (hypopnea/apnea); decrease in oxygen saturation; decrease in heart rate; and vasodilation, manifested in decreased blood pressure. The crescendo phase has the opposite effects: increase in respiratory effort and rate, i.e. hyperpnea; increase in heart rate; and vasoconstriction, leading to increased blood pressure. Sometimes the hyperpnea is accompanied by an arousal, which is manifested as a motion artifact in the photoplethysmograph signal. The changes in heart rate and vasomotion (dilation and constriction) depend on the severity of the heart failure, as discussed below.
[0068] The inventors have also found that decompensated heart failure patients nearly always present long sequences of periodic Cheyne-Stokes breathing episodes, with a cycle length between 55 and 180 seconds. In general, the longer the cycle length, the more severe is the state of the disease. Therefore, processor 32 uses the shape characteristics of the VLF saturation signal in measuring time characteristics 48 of the patient's Cheyne-Stokes episodes. Specifically, the processor detects desaturation episodes extending over multiple consecutive Cheyne-Stokes breathing cycles in order to identify the presence of CSBS.
[0069] In order to detect and measure the duration of multi-cycle Cheyne-Stokes episodes, processor 32 typically locates the local maxima and local minima of the VLF saturation signal. The processor may also compute the difference between the maximal and minimal saturation values (in the unsmoothed saturation signal 44), as well as the corresponding wavelengths. The processor extracts time sequences of cyclic breathing with similar desaturation values and similar wavelengths, falling in the range that is characteristic of Cheyne-Stokes cycles. (Typically only a certain percentage, such as 80%, of the desaturation and wavelength values are required to be close to the median values of the sequence, in order to avoid losing sequences due to intervening outliers. For example, a 50% deviation from the median value of 80% of the wavelength and desaturation values may be accepted for a sequence that is at least of a certain minimum duration, such as 5 min.) The processor chooses the longest segments that meet the above similarity criteria. Alternatively, a hysteresis procedure may be used to ensure robustness against outliers. The total Cheyne-Stokes time is then computed as the total duration of all the segments that are classified as Cheyne-Stokes breathing events.
[0070] In order to validate the automatic measurements of Cheyne-Stokes episodes described above, the inventors conducted a clinical trial, which included 91 full-night ambulatory polysomnography tests for patients with advanced heart failure. Cheyne-Stokes episodes were marked manually by an experienced scorer, and these manual results were compared to the results of the automatic process described above. The correlation between manual automatic scoring was 83%, which is as good as the typical correlation between different human scorers.
[0071] To ensure further that the sequences of cyclic breathing episodes are indeed associated with the severity of heart failure status, processor 32 evaluates the slope of the saturation signal (or of the DC component of the pulse oximeter signal) for each desaturation event. In central apnea, or when the heart failure state is grave, the slope of the exit from the cycle is moderate, i.e., it is similar to the typical (or specific) entry slope. Therefore, to identify a time sequence of cyclic breathing as Cheyne-Stokes, processor 32 requires that the sequence comprise mainly (typically at least 80%) events of moderate slope. (The above-mentioned PCT Patent Application PCT/IL2006/000148 defines formal criteria for assessing the symmetry of periodic breathing episodes, which may also be used in the present context for distinguishing Cheyne-Stokes events.) This requirement of moderate slope may be applied to the median slope value.
[0072] These observations with respect to the symmetry of the periodic breathing patterns apply both to the slowly-varying heart rate and saturation signals and to the envelopes of the other, rapidly-varying signals shown by the other traces. The term “envelope” in this context typically means a signal derived from the local minima and/or local maxima of another signal, with or without smoothing (by convolution or resampling, for example). “Envelopes” may also be derived by other mathematical operations known in the art, such as application of Hilbert transforms. The inventors have found that periodic breathing patterns associated with CSA generally tend to be more symmetrical than the patterns associated with OSA, presumably due to the different physiological mechanisms that are involved in the different types of apneas. Therefore, processor 32 may validate the prognostic value of the Cheyne-Stokes marker by considering only events with mild exit slope from desaturation events. The inventors found that computing the slope of the saturation curve by fitting a line (by the least-square method) to the curve over a nine-second epoch, and requiring that the slope of the line be less then 0.7 percent/second is a good implementation of this mild desaturation condition.
[0073] Processor 32 associates each segment with its segment duration and with its median desaturation value. The features of the Cheyne-Stokes segments are prognostic of patient outcome in cases of heart failure (and other illnesses). Long wavelength, in particular, is associated with bad prognosis. Thus, processor 32 typically detects signal components that have a period greater than a minimum period of at least 30 sec. In the marker validation experiments that are described herein, the inventors required the median cycle length to be above 55 seconds and the median desaturation value to be no less the 2% in order to classify a periodic breathing pattern as Cheyne-Stokes breathing.
[0074] On the other hand, time segments with steep exit saturation slope typically correspond to obstructive apnea/hypopnea events. Other features of obstructive apnea/hypopnea time segments include short wavelength, large vasomotion, and large heart rate modulations. These phenomena are generally associated with good prognosis, since they reflect the patient's ability to manifest enhanced sympathetic activity.
[0075] In addition to the saturation measurements and Cheyne-Stokes detection, processor 32 may also process an AC absorption or reflectance signal 52 that is output by device 24 in order to compute a heart rate 54, as is known in the art. Furthermore, the AC signal may be analyzed to detect a beat morphology 56. The processor identifies certain aberrations in this morphology as arrhythmias, such as premature ventricular contractions (PVCs) 58. It keeps a record of the occurrences of such arrhythmias, in a manner similar to a Holter monitor, but without requiring the use of ECG leads. The total number of abnormal heart beats and (specifically PVCs) that are accumulated in such a record, particularly during sleep, is indicative of bad prognosis. As the inventors have found that premature beats during sleep have the greatest prognostic value for advanced heart failure patients, the processor may be configured to count the number of premature beats only during sleep or during episodes of Cheyne-Stokes breathing.
[0076] FIG. 4 is a schematic plot of an AC photoplethysmograph signal 60, alongside a corresponding ECG signal 62, illustrating a method for detection of PVCs in accordance with an embodiment of the present invention. Signal 60 can be seen to comprise a series of regular waveforms, which are indicative of arterial blood flow. A PVC is manifested as an aberrant waveform 64 in signal 60, and likewise by an abnormal waveform 66 in signal 62. Processor 32 analyzes the shape, amplitude and timing of waveform 64 in the plethysmograph signal in order to determine that the aberrant wave represents PVC, even without the use of any sort of ECG monitoring.
[0077] In one embodiment, arrhythmias are identified in photoplethysmograph signal 60 based on the following features:
-
- [0078] 1. Local maxima and minima are extracted from the signal in segments of the signal whose length is less than the typical RR interval (i.e., the typical time difference between successive heart beats). For example, 0.3 seconds is an appropriate segment length for this purpose.
- [0079] 2. The width of each beat is defined, for example by measuring the time difference between successive locations of photoplethysmograph signal values whose energy is equal to the average (possibly a weighted average) of the local maximum and minimum.
- [0080] 3. Beats with short width typically correspond to PVCs, as shown in
FIG. 4 . The number of such beats is a measure of the severity of arrhythmia. - [0081] 4. An additional criterion for detecting an arrhythmia is that the time span of two beats, one of which has a short width, is roughly equal to the time span of two normal beats.
[0082] Although FIG. 4 and the above description relate specifically to PVCs, the principles of this embodiment may likewise be applied in detecting other types of premature heart beats, as well as various other types of heartbeat irregularities. Such irregularities are associated with reduced stroke volume, which in turn affect the amplitude, width and other features of the photoplethysmograph waveform.
[0083] Other types of aberrant waveforms in photoplethysmograph signal 60 may correspond to motion artifacts 80 (listed in FIG. 3A , but not shown in FIG. 4 ). Motion is characterized by local maxima well above normal beat range (for example, at least twice the normal value). The prevalence of motion artifacts can be used in detecting movement, which indicate whether the patient is in a sleep or waking state 82. Alternatively or additionally, a motion sensor may be used to detect arousals.
[0084] Referring further to FIGS. 3A and 3B , processor 32 may additionally extract other cardiorespiratory parameters from signal 52, either directly or indirectly. For example, the processor may apply very-low-frequency filtering to heart rate 54 in order to detect heart rate modulations 70. Additionally or alternatively, the envelope of signal 52 may be processed in order to detect characteristics of vasomodulation 72, i.e., arterial dilation and constriction.
[0085] Further additionally or alternatively, processor 32 may compute a respiration energy and/or rate characteristic 74 based on high-frequency components of signal 52. Respiratory sinus arrhythmia is a natural cycle of arrhythmia that occurs in healthy people through the influence of breathing on the flow of sympathetic and vagus impulses to the sinoatrial node in the heart. This effect may be used to calculate respiration from heart rate. Well-treated heart failure patients, however, are frequently under the control of cardiac pacemakers and often take beta-blockers and ACE inhibitors that suppress this phenomenon. High-frequency (10-30 cycles/min, i.e., 0.17-0.5 Hz) filtering of the photoplethysmograph signal enables the processor to determine respiration energy and/or rate characteristics in these cases, as well.
[0086] Very-low-frequency components of characteristic 74 are indicative of a respiration modulation 76. Processor 32 combines the various cardiac, respiratory and vasomodulation parameters described above in order to provide a general picture of cardiorespiratory effects 78, all on the basis of the photoplethysmograph signals.
[0087] Similar procedures to those described above can be applied to the detrended AC photoplethysmograph signal. One way to perform detrending is to replace the photoplethysmograph signal with its amplitude feature (maximum minus minimum signal). Other methods include subtracting a polynomial that approximates the signal, or using local maxima or local minima features. Following detrending, the processor applies very-low-frequency filtering followed by outlier rejection, and then computes the median vasomotion of each sequence.
[0088] The processor may perform similar analyses on heart rate and respiratory signals from other sources. Arousals can estimated from motion artifacts as described above or from other data if available (such as EEG alpha and beta frequencies, or scorer marking, or a motion sensor).
[0089] Information regarding sleep/wake state 82 is combined with Cheyne-Stokes time 48 to determine specific, cumulative Cheyne-Stokes time 84 during sleep. The total Cheyne-Stokes time and percentage of Cheyne-Stokes time during sleep have prognostic value: A large percentage of Cheyne-Stokes time is associated with mortality and high levels of brain natriuretic peptide (BNP), which are associated with severity of heart failure. Furthermore, information about sleep time can be used to ensure that low Cheyne-Stokes duration is not associated with little or no sleep. (The inventors have determined the prognostic value of total Cheyne-Stokes time only in patients who slept for at least a certain minimal duration, such as two hours.) The prognostic value of Cheyne-Stokes information derived in the above manner is illustrated in FIGS. 5-8 below.
[0090] This information regarding Cheyne-Stokes time 84 in turn is combined with the general picture of cardiorespiratory effects 78 in order to provide some or all of the following combined information 86 for each Cheyne-Stokes sequence during sleep:
-
- [0091] 1. duration
- [0092] 2. median wavelength
- [0093] 3. median desaturation
- [0094] 4. median vasomotion
- [0095] 5. median heart rate modulation
- [0096] 6. median respiratory modulation
- [0097] 7. number of PVCs and other premature beats.
- [0098] 8. arousal index: number of arousals
Alternatively, the above-mentioned median functions may be replaced by similar functions based on average values or average of values in the middle tertile, inter-quartile range, or any other appropriate segment. Each of the above parameters can also be computed separately for REM sleep and NREM sleep.
[0099] In an exemplary embodiment, the following criteria may be applied to the various processed outputs of oximetry device 24 in order to derive information 86 and measure the manifestations of Cheyne-Stokes breathing:
-
- [0100] 1. The saturation signal is filtered in the Cheyne-Stokes frequency range, typically 1/180-1/40 Hz.
- [0101] 2. A time segment is identified as a Cheyne-Stokes event (and the durations of such time segments are summed) if the segment contains a sequence of at least three cycles of desaturation for which:
- [0102] a. Median desaturation (compared to the previous saturation level) is at least 2%. Alternatively, another representative saturation level, such as the mean or minimum, may be used.
- [0103] b. Mean cycle length is long (55-180 sec). (Short cycle length is not associated with bad prognosis.)
- [0104] c. Cycle length fluctuation within each sequence may optionally be limited (to less then 10% fluctuation, for example).
- [0105] d. Moderate vasomotion, based on at least one of the following:
- [0106] i. Median of maximal desaturation slope in each cycle is less then a maximum slope limit, such as 0.7 percent/sec. (by least squares fit of a line to the desaturation curve). Alternatively, a measure of mean slope may be used.
- [0107] ii. Only moderate fluctuations (typically no greater than 10% of the normal range) occur in the VLF range of the detrended respiration signal.
- [0108] 3. Optionally, for a periodic breathing cycle to be identified as part of a Cheyne-Stokes event, respiration characteristic 74 may be required to reach a minimum indicating zero respiratory effort during the cycle. This minimum may be identified based on the VLF components of characteristic 74 in respiration modulation 76.
- [0109] 4. Outlier rejection procedures may be applied to the saturation and respiration values before classifying time segments. For example:
- [0110] a. As noted above, a certain fraction (typically up to 20%) of the desaturation and wavelength values may be far from the median values of the sequence, and extreme desaturations (for example, >50%) may be rejected as faulty readings.
- [0111] b. The mean cycle length can be calculated after discarding values that are far from the mean (for example, values in the top and bottom deciles.)
- [0112] c. Consistency may be enforced by permitting the relation between cycle length and desaturation to vary linearly within given bounds.
- [0113] 5. In addition to or instead of the above criteria, self-similarity measures can be used in identifying sequences of Cheyne-Stokes cycles. For example, a distinct peak in the 1/180-1/40 Hz range of the Fourier transform of the sequence of periodic breathing cycles or high autocorrelation of the cycles is an indicator of such self-similarity.
[0114] Similarly, when a respiration signal is obtained without a saturation signal, Cheyne-Stokes respiration segments can be found by applying the above criteria to the respiration signal (excluding the computations that relate to saturation values).
[0115] Referring again to FIG. 3B , sleep/wake state 82 may be combined with analysis of cardiorespiratory and Cheyne-Stokes effects in order to perform automatic sleep staging 88. All of the factors that are used in determining the sleep stage may be derived solely from the signals generated by oximetry device 24. Alternatively, other signals may be incorporated into the sleep staging calculation.
[0116] Sleep states are classified by processor 32 as light sleep, deep slow-wave sleep (SWS) 90 and REM 92. During REM sleep, the patient is partially paralyzed, so that there is no motion. Furthermore, due to the changes in autonomic control and the partial paralysis that characterize REM sleep, the Cheyne-Stokes wavelength tends to be longer and the desaturation deeper in REM sleep that in light sleep. On the other hand, there are no apnea episodes in deep sleep. Others factors characterizing deep sleep include regularity of respiratory cycle length and low vasomotion.
[0117] Processor 32 may use the distribution of sleep stages and of apnea events during sleep in computing a sleep quality index 94. Typically, high percentages of REM and SWS, as well as apnea-free (or nearly apnea-free) segments in non-SWS sleep, are indicative of good prognosis for heart failure patients. By contrast, low percentages of REM or SWS indicate a poor prognosis. Further aspects of sleep staging and sleep quality assessment are described in the above-mentioned PCT patent applications.
Clinical Application and Results
[0118] The methods described above for measuring and quantifying Cheyne-Stokes breathing and attendant heart failure prognosis may be used conveniently in performing frequent checks on patient status, both at home and in the hospital. Additionally or alternatively, occasional checks of this sort may be used for risk stratification and screening. As explained above, these methods may be implemented using measurements made solely by pulse oximetry device 24, or alternatively in conjunction with other sensors, as in a multi-monitor polysomnography system, or in an implantable device, or using other types of respiratory sensors.
[0119] The inventors compared Cheyne-Stokes time with heart failure status in 91 tests of advanced heart failure patients. Results of this study are presented below. The cumulative duration of Cheyne-Stokes breathing during a night's sleep was measured, wherein Cheyne-Stokes cycles were identified as described above (including the requirements of mild slope—up to 0.7 percent/sec, median desaturation of at least 2%, and median cycle length of 55 to 180 sec.) The status of the patients was determined by six-month survival and BNP levels, which are generally considered the best marker for heart failure status. For this purpose, a blood sample was drawn from each patient and tested for NT-proBNP on the night of the sleep study. Serum N-terminal prohormone Brain Natriuretic Peptide (NT-proBNP) was measured using the Elecsys® proBNP electro-chemiluminescence immunoassay, run on the Elecsys 1010 benchtop analyzer (Roche Diagnostics, Indianapolis, Ind.).
[0120] FIG. 5 is a Kaplan-Meier plot of patient survival according to the standard BNP kit values. According to accepted diagnostic standards, a state of decompensated heart failure is associated with a serum NT-proBNP level above 450, 900, or 1,800 pg/mL for patients whose age is less than 50, 50-75, or above 75, respectively. An upper trace 70 shows the rate of survival over time for the patients with low BNP (below the decompensation threshold), while an upper trace 72 shows the rate for patients with high BNP. The odd ratio for dying among patients with high (decompensated) BNP values was 3.5 times higher then for patients in the low-BNP group. This result, however, was not statistically significant in this trial (p=0.18).
[0121] FIG. 6 is a Kaplan-Meier plot of patient survival according to the automated Cheyne-Stokes marker described above, in accordance with an embodiment of the present invention. In this figure, an upper trace 74 shows the survival rate of the patients who had low cumulative duration of Cheyne-Stokes breathing episodes, while a lower trace 76 shows the survival rate for patients with high cumulative Cheyne-Stokes duration. The inventors have found that typically, a cumulative duration of Cheyne-Stokes breathing episodes in excess of 45 minutes during a night's sleep is indicative of poor prognosis. In the results shown in FIG. 6 , the Cheyne-Stokes cutoff (48 minutes) was selected to best predict BNP cutoff according to the standard guidelines described above.
[0122] Comparison of FIGS. 5 and 6 shows that the respiration-based measure of heart failure severity is superior to the standard BNP test in predicting mortality: The odd ratio for dying among patients with high Cheyne-Stokes duration is better then 5.2. (There was no mortality at all in the low Cheyne-Stokes group.) The results are statistically highly significant (p=0.017). For some higher cutoff points of Cheyne-Stokes time, the odd ratio improved even further. Furthermore, these results also demonstrate that Cheyne-Stokes duration is a strong predictor of BNP level and differentiates between compensated and decompensated levels of heart failure.
[0123] FIG. 7 is a receiver operating characteristic (ROC) plot, which schematically compares the sensitivity and specificity of predicting heart failure prognosis using BNP and duration of Cheyne-Stokes breathing episodes, as measured in accordance with an embodiment of the present invention. An upper trace 78 is the ROC curve for Cheyne-Stokes duration, while a lower trace 80 is the ROC curve for the BNP marker. To derive the results shown in the figure, both Cheyne-Stokes duration and BNP were tested against six-month mortality of the patients in the study. The plot shows that the Cheyne-Stokes marker gives greater sensitivity and specificity. The area under the curve (AUC) for BNP is 68% (p=0.082), while it is 75% (p=0.015) for the Cheyne-Stokes marker.
[0124] FIG. 8 is a Kaplan-Meier plot of six-month survival of the heart failure patients as a function of the severity of symptoms classified by the methods described above, in accordance with another embodiment of the present invention. In this case, patients were classified into two groups: one group with severe Cheyne-Stokes breathing coupled with cardiac arrhythmia, and the other with breathing and heart rhythm that showed mild or no symptoms of these kinds. For this purpose, patients who exhibited at least 200 premature beats in the course of the night's sleep were classified as suffering from cardiac arrhythmia. An upper trace 82 in the figure shows the survival rate for the group with mild or no symptoms, while a lower trace 84 shows the survival rate for patients in the severe/arrhythmic group. The plot demonstrates that the combined detection of Cheyne-Stokes breathing and cardiac arrhythmias is an even stronger predictor of survival than Cheyne-Stokes breathing by itself: The odd ratio in this case was 17.7 with p=0.000.
[0125] The results of FIGS. 5-8 show that photoplethysmographic monitoring during sleep may be used effectively for the prognosis of heart failure patients. As explained above, this sort of monitoring is simple to carry out in the patient's home or hospital bed and may be performed at regular intervals, at low cost and minimal discomfort to the patient. It provides physicians with an accurate prognostic indicator, which they can use in choosing the optimal treatment, such as determining whether a patient requires hospitalization, and adjusting treatment parameters (such as drug titration) based on changes in the patient's condition. For example, using the techniques described above, the physician may measure and quantify the patient's symptoms (Cheyne-Stokes duration and possibly arrhythmias) prior to initiating or making a change in treatment, and may repeat the measurement thereafter in order to assess the effectiveness of the treatment and possibly readjust treatment parameters.
[0126] As another alternative, the physician may fix a specific critical Cheyne-Stokes duration for individual patients, and then set a monitoring system to alarm whenever a specific duration is exceeded.
[0127] Although the embodiments described above relate mainly to signals captured by pulse oximetry device 24, the principles of the present invention may be applied to respiration signals captured by any other suitable type of sensor. Such sensors may be based, for example, on electrical measurements of thoracic and abdominal movement, using skin electrodes to make a plethysmographic measurement of the patient's respiratory effort, or a belt to sense changes in the body perimeter. Additionally or alternatively, air flow measurement, based on a pressure cannula, thermistor, or CO2 sensor, may be used for respiration sensing. In other embodiments of the present invention, a capnograph may be used in detecting sleep apneas, either in conjunction with or separately from the pulse oximeter used in the techniques described above.
[0128] It will thus be appreciated that the embodiments described above are cited by way of example, and that the present invention is not limited to what has been particularly shown and described hereinabove. Rather, the scope of the present invention includes both combinations and subcombinations of the various features described hereinabove, as well as variations and modifications thereof which would occur to persons skilled in the art upon reading the foregoing description and which are not disclosed in the prior art.
Claims
The invention claimed is:
1. A method for diagnosis, comprising:
receiving a signal associated with respiration of a patient during sleep;
operating a computer for processing the signal to distinguish Cheyne-Stokes breathing events from obstructive apnea events and to detect occurrences of a pattern corresponding to multiple cycles of periodic breathing based on the distinguished Cheyne-Stokes breathing events;
computing a cumulative duration of the occurrences;
determining a prognosis of at least a heart failure (HF) in the patient based on the cumulative duration.
2. The method according to claim 1, wherein determining the prognosis comprises predicting a probability of survival of the patient.
3. The method according to claim 1, wherein determining the prognosis comprises predicting a level of at least one marker, selected from a group of markers consisting of Brain Natriuretic Peptide (BNP) and N-terminal prohormone Brain Natriuretic Peptide (NT proBNP).
4. The method according to claim 1, wherein determining the prognosis comprises differentiating compensated from decompensated heart failure.
5. The method according to claim 1, wherein determining the prognosis comprises predicting a hospitalization of the patient.
6. The method according to claim 1, wherein determining the prognosis comprises providing an indication of a need for a change in treatment of the patient.
7. The method according to claim 1, wherein the signal comprises a photoplethysmograph signal.
8. The method according to claim 7, and comprising processing the photoplethysmograph signal so as to identify a cardiac arrhythmia of the patient.
9. The method according to claim 7, wherein determining the prognosis comprises finding the prognosis responsively both to the occurrences of the pattern and to the cardiac arrhythmia.
10. The method according to claim 7, and comprising processing the photoplethysmograph signal so as to measure a vasomodulation in the body.
11. The method according to claim 7, and comprising processing the photoplethysmograph signal so as to measure a modulation of respiration by the patient.
12. The method according to claim 7, wherein receiving the signal comprises receiving the photoplethysmograph signal from a sensor on a finger of the patient, and wherein processing the signal comprises digitizing and filtering the signal so as to generate data using processing circuitry in a control unit that is fastened to a forearm of the patient, and storing the data in a memory in the control unit.
13. The method according to claim 12, wherein processing the signal comprises connecting the control unit to an external processor and uploading the data from the memory to the external processor for analysis, and wherein storing the data comprises receiving and recording the data in the control unit while the control unit is disconnected from the external processor.
14. The method according to claim 7, wherein receiving the signal comprises receiving the photoplethysmograph signal from a sensor implanted in a body of the patient.
15. The method according to claim 1, and comprising treating the patient so as to reduce the cumulative duration.
16. The method according to claim 1, wherein determining the prognosis comprises classifying the patient in a high-mortality group if the cumulative duration exceeds a predetermined threshold.
17. The method according to claim 16, wherein the threshold is at least 45 minutes.
18. The method according to claim 1, wherein processing the signal comprises identifying a periodic breathing event when an event duration of the multiple cycles is no less than a predetermined minimum duration.
19. The method according to claim 1, wherein processing the signal comprises identifying a periodic breathing event when a number of the multiple cycles is no less than a predetermined minimum number.
20. The method according to claim 19, wherein the minimum number is three.
21. The method according to claim 1, wherein the signal is indicative of blood oxygen saturation.
22. The method according to claim 21, wherein processing the signal comprises identifying a periodic breathing event when a representative level of the blood oxygen saturation during the cycles is less than the blood oxygen saturation before the cycles by at least a predetermined amount.
23. The method according to claim 22, wherein the representative level of the blood oxygen saturation is one of a median saturation, a mean saturation, and a minimum saturation during the cycles.
24. The method according to claim 22, wherein the periodic breathing event is identified when the representative level of the blood oxygen saturation is less than the blood oxygen saturation before the cycles by at least 2%.
25. The method according to claim 1, wherein distinguishing the Cheyne-Stokes breathing events comprises determining a slope of the signal, and distinguishing the Cheyne-Stokes breathing events responsively to the slope.
26. The method according to claim 25, wherein determining the slope comprises finding a measure of the slope selected from a group of measures consisting of a median slope and a mean slope, and distinguishing the Cheyne-Stokes breathing events based on a maximal value of the measure of the slope in each of the cycles.
27. The method according to claim 26, wherein distinguishing the Cheyne-Stokes breathing events comprises classifying a breathing event as Cheyne-Stokes only when the maximal value of the measure is no greater than 0.7%/sec.
28. The method according to claim 1, wherein processing the signal comprises extracting from the signal a measure of respiratory energy, and filtering the measure so as to eliminate signal components at frequencies equal to and greater than a respiratory frequency of the patient.
29. The method according to claim 1, wherein receiving the signal comprises monitoring the respiration using a non-contact sensor.
30. Apparatus for diagnosis, comprising:
a sensor, which is configured to be coupled to a body of a patient during sleep and to output a signal associated with respiration of the patient; and
a processor, which is coupled to process the signal to distinguish Cheyne-Stokes breathing events from obstructive apnea events and to detect occurrences of a pattern corresponding to multiple cycles of periodic breathing based on the distinguished Cheyne-Stokes breathing events, to compute a cumulative duration of the occurrences, and to output an indication of a prognosis of at least a heart failure (HF) in the patient based on the cumulative duration.
31. The apparatus according to claim 30, wherein the prognosis comprises a prediction of a probability of survival of the patient.
32. The apparatus according to claim 30, wherein the prognosis comprises a prediction of a level of at least one marker, selected from a group of markers consisting of Brain Natriuretic Peptide (BNP) and N-terminal prohormone Brain Natriuretic Peptide (NT proBNP).
33. The apparatus according to claim 30, wherein the prognosis differentiates compensated from decompensated heart failure.
34. The apparatus according to claim 30, wherein the prognosis comprises a prediction of a hospitalization of the patient.
35. The apparatus according to claim 30, wherein the prognosis comprises an indication of a need for a change in treatment of the patient.
36. The apparatus according to claim 30, wherein the signal comprises a photoplethysmograph signal.
37. The apparatus according to claim 36, wherein the processor is configured to process the photoplethysmograph signal so as to identify a cardiac arrhythmia of the patient.
38. The apparatus according to claim 36, wherein the processor is configured to output an indication of a prognosis of the patient responsively both to the occurrences of the pattern and to the cardiac arrhythmia.
39. The apparatus according to claim 36, wherein the processor is configured to process the photoplethysmograph signal so as to measure a vasomodulation in the body.
40. The apparatus according to claim 36, wherein the processor is configured to process the photoplethysmograph signal so as to measure a modulation of respiration by the patient.
41. The apparatus according to claim 36, wherein the sensor is configured to be attached to a finger of the patient, and wherein the apparatus comprises a control unit, which contains the processor and a memory, and which is configured to be fastened to a forearm of the patient, and wherein the processor is configured to digitize and filter the signal so as to generate data and to store the data in the memory in the control unit.
42. The apparatus according to claim 41, wherein the control unit is configured to be connected to an external processor and to upload the data from the memory to the external processor for analysis, and wherein the processor is configured to receive and record the data in the memory while the control unit is disconnected from the external processor.
43. The apparatus according to claim 36, wherein the sensor is configured to be implanted in a body of the patient.
44. The apparatus according to claim 30, wherein the processor is configured to classify the patient in a high-mortality group if the cumulative duration exceeds a predetermined threshold.
45. The apparatus according to claim 44, wherein the threshold is at least 45 minutes.
46. The apparatus according to claim 30, wherein the processor is configured to identify a periodic breathing event when an event duration of the multiple cycles is no less than a predetermined minimum duration.
47. The apparatus according to claim 30, wherein the processor is configured to identify a periodic breathing event when a number of the multiple cycles is no less than a predetermined minimum number.
48. The apparatus according to claim 47, wherein the minimum number is three.
49. The apparatus according to claim 30, wherein the signal is indicative of blood oxygen saturation.
50. The apparatus according to claim 49, wherein the processor is configured to identify a periodic breathing event when a representative level of the blood oxygen saturation during the cycles is less than the blood oxygen saturation before the cycles by at least a predetermined amount.
51. The apparatus according to claim 50, wherein the representative level of the blood oxygen saturation is one of a median saturation, a mean saturation, and a minimum saturation during the cycles.
52. The apparatus according to claim 50, wherein the periodic breathing event is identified when the representative level of the blood oxygen saturation is less than the blood oxygen saturation before the cycles by at least 2%.
53. The apparatus according to claim 30, wherein the processor is configured to determine a slope of the signal, and to distinguish the Cheyne-Stokes breathing events responsively to the slope.
54. The apparatus according to claim 53, wherein the processor is configured to find a measure of the slope selected from a group of measures consisting of a median slope and a mean slope, and to distinguish the Cheyne-Stokes breathing events based on a maximal value of the measure of the slope in each of the cycles.
55. The apparatus according to claim 54, wherein the processor is configured to classify a breathing event as Cheyne-Stokes only when the maximal value of the measure is no greater than 0.7%/sec.
56. The apparatus according to claim 30, wherein the processor is configured to extract from the signal a measure of respiratory energy, and to filter the measure so as to eliminate signal components at frequencies equal to and greater than a respiratory frequency of the patient.
57. The apparatus according to claim 30, wherein receiving the sensor comprises a non-contact sensor.
58. A computer software product, comprising a computer-readable medium in which program instructions are stored, which instructions, when read by a computer, cause the computer to receive a signal associated with respiration of a patient during sleep, and to process the signal to distinguish Cheyne-Stokes breathing events from obstructive apnea events and to detect occurrences of a pattern corresponding to multiple cycles of periodic breathing based on the distinguished Cheyne-Stokes breathing events, to compute and output a cumulative duration of the occurrences, and to output an indication of a prognosis of at least a heart failure (HF) in the patient based on the cumulative duration.
59. The product according to claim 58, wherein the prognosis comprises a prediction of a probability of survival of the patient.
60. The product according to claim 58, wherein the prognosis comprises a prediction of a level of at least one marker, selected from a group of markers consisting of Brain Natriuretic Peptide (BNP) and N-terminal prohormone Brain Natriuretic Peptide (NT proBNP).
61. The product according to claim 58, wherein the prognosis differentiates compensated from decompensated heart failure.
62. The product according to claim 58, wherein the prognosis comprises a prediction of a hospitalization of the patient.
63. The product according to claim 58, wherein the prognosis comprises an indication of a need for a change in treatment of the patient.
64. The product according to claim 58, wherein the signal comprises a photoplethysmograph signal.
65. The product according to claim 64, wherein the instructions cause the computer to process the photoplethysmograph signal so as to identify a cardiac arrhythmia of the patient.
66. The product according to claim 64, wherein the instructions cause the computer to output an indication of a prognosis of the patient responsively both to the occurrences of the pattern and to the cardiac arrhythmia.
67. The product according to claim 64, wherein the instructions cause the computer to process the photoplethysmograph signal so as to measure a vasomodulation in the body.
68. The product according to claim 64, wherein the instructions cause the computer to process the photoplethysmograph signal so as to measure a modulation of respiration by the patient.
69. The product according to claim 58, wherein the instructions cause the computer to classify the patient in a high-mortality group if the cumulative duration exceeds a predetermined threshold.
70. The product according to claim 69, wherein the threshold is at least 45 minutes.
71. The product according to claim 58, wherein the instructions cause the computer to identify a periodic breathing event when an event duration of the multiple cycles is no less than a predetermined minimum duration.
72. The product according to claim 58, wherein the instructions cause the computer to identify a periodic breathing event when a number of the multiple cycles is no less than a predetermined minimum number.
73. The product according to claim 72, wherein the minimum number is three.
74. The product according to claim 58, wherein the signal is indicative of blood oxygen saturation.
75. The product according to claim 74, wherein the instructions cause the computer to identify a periodic breathing event when a representative level of the blood oxygen saturation during the cycles is less than the blood oxygen saturation before the cycles by at least a predetermined amount.
76. The product according to claim 74, wherein the representative level of the blood oxygen saturation is one of a median saturation, a mean saturation, and a minimum saturation during the cycles.
77. The product according to claim 74, wherein the periodic breathing event is identified when the representative level of the blood oxygen saturation is less than the blood oxygen saturation before the cycles by at least 2%.
78. The product according to claim 58, wherein the instructions cause the computer to determine a slope of the signal, and to distinguish the Cheyne-Stokes breathing events responsively to the slope.
79. The product according to claim 78, wherein the instructions cause the computer to find a measure of the slope selected from a group of measures consisting of a median slope and a mean slope, and to distinguish the Cheyne-Stokes breathing events based on a maximal value of the measure of the slope in each of the cycles.
80. The product according to claim 79, wherein the instructions cause the computer to classify a breathing event as Cheyne-Stokes only when the maximal value of the measure is no greater than 0.7%/sec.
81. The product according to claim 58, wherein the instructions cause the computer to extract from the signal a measure of respiratory energy, and to filter the measure so as to eliminate signal components at frequencies equal to and greater than a respiratory frequency of the patient.
Patent Citations (172)
| Patent | Date | Inventor | Cited By |
|---|---|---|---|
| US3835837(A) | 1974-09-01 | Peck | Applicant |
| US3863625(A) | 1975-02-01 | Viglione et al. | Applicant |
| US4258719(A) | 1981-03-01 | Lewyn | Applicant |
| US4545387(A) | 1985-10-01 | Balique | Applicant |
| US4589420(A) | 1986-05-01 | Adams et al. | Applicant |
| US4777962(A) | 1988-10-01 | Watson et al. | Applicant |
| US4955379(A) | 1990-09-01 | Hall | Applicant |
| US5101831(A) | 1992-04-01 | Koyama et al. | Applicant |
| US5187657(A) | 1993-02-01 | Forbes | Applicant |
| US5271411(A) | 1993-12-01 | Ripley et al. | Applicant |
| US5280791(A) | 1994-01-01 | Lavie | Applicant |
| US5385144(A) | 1995-01-01 | Yamanishi et al. | Applicant |
| US5513644(A) | 1996-05-01 | McClure et al. | Applicant |
| US5560369(A) | 1996-10-01 | McClure et al. | Applicant |
| US5588425(A) | 1996-12-01 | Sackner et al. | Applicant |
| US5605151(A) | 1997-02-01 | Lynn | Applicant |
| US5605158(A) | 1997-02-01 | Snell | Applicant |
| US5623933(A) | 1997-04-01 | Amano et al. | Applicant |
| US5685315(A) | 1997-11-01 | McClure et al. | Applicant |
| US5718242(A) | 1998-02-01 | McClure et al. | Applicant |
| US5720294(A) | 1998-02-01 | Skinner | Applicant |
| US5755229(A) | 1998-05-01 | Amano et al. | Applicant |
| US5769825(A) | 1998-06-01 | Lynn | Applicant |
| US5803066(A) | 1998-09-01 | Rapoport et al. | Applicant |
| US5819007(A) | 1998-10-01 | Elghazzawi | Applicant |
| US5846208(A) | 1998-12-01 | Pichlmayr et al. | Applicant |
| US5857978(A) | 1999-01-01 | Hively et al. | Applicant |
| US5865756(A) | 1999-02-01 | Peel | Applicant |
| US5891023(A) | 1999-04-01 | Lynn | Applicant |
| US5902250(A) | 1999-05-01 | Verrier et al. | Applicant |
| US5999846(A) | 1999-12-01 | Pardey et al. | Applicant |
| US6011477(A) | 2000-01-01 | Teodorescu | Applicant |
| US6032072(A) | 2000-02-01 | Greenwald et al. | Applicant |
| US6070098(A) | 2000-05-01 | Moore-Ede et al. | Applicant |
| US6095984(A) | 2000-08-01 | Amano et al. | Applicant |
| US6122536(A) | 2000-09-01 | Sun et al. | Applicant |
| US6142950(A) | 2000-11-01 | Allen et al. | Examiner |
| US6223064(B1) | 2001-04-01 | Lynn et al. | Applicant |
| US6304775(B1) | 2001-10-01 | Iasemidis et al. | Applicant |
| US6308094(B1) | 2001-10-01 | Shusterman et al. | Applicant |
| US6319205(B1) | 2001-11-01 | Goor et al. | Applicant |
| US6322515(B1) | 2001-11-01 | Goor et al. | Applicant |
| US6342039(B1) | 2002-01-01 | Lynn et al. | Examiner |
| US6361501(B1) | 2002-03-01 | Amano et al. | Applicant |
| US6375623(B1) | 2002-04-01 | Gavriely | Applicant |
| US6409675(B1) | 2002-06-01 | Turcott | Applicant |
| US6480733(B1) | 2002-11-01 | Turcott | Applicant |
| US6480734(B1) | 2002-11-01 | Zhang et al. | Applicant |
| US6511424(B1) | 2003-01-01 | Moore-Ede et al. | Applicant |
| US6519490(B1) | 2003-02-01 | Wiesel | Applicant |
| US6527729(B1) | 2003-03-01 | Turcott | Applicant |
| US6529752(B2) | 2003-03-01 | Krausman et al. | Applicant |
| US6549804(B1) | 2003-04-01 | Osorio et al. | Applicant |
| US6589188(B1) | 2003-07-01 | Street et al. | Applicant |
| US6595929(B2) | 2003-07-01 | Stivoric et al. | Applicant |
| US6600949(B1) | 2003-07-01 | Turcott | Applicant |
| US6608562(B1) | 2003-08-01 | Kimura et al. | Applicant |
| US6609016(B1) | 2003-08-01 | Lynn | Applicant |
| US6702752(B2) | 2004-03-01 | Dekker | Applicant |
| US6748252(B2) | 2004-06-01 | Lynn et al. | Examiner |
| US6752765(B1) | 2004-06-01 | Jensen et al. | Examiner |
| US6760608(B2) | 2004-07-01 | Lynn | Examiner |
| US6805673(B2) | 2004-10-01 | Dekker | Applicant |
| US6829501(B2) | 2004-12-01 | Nielsen et al. | Applicant |
| US6839581(B1) | 2005-01-01 | El Solh et al. | Applicant |
| US6856829(B2) | 2005-02-01 | Ohsaki et al. | Applicant |
| US6878121(B2) | 2005-04-01 | Krausman et al. | Applicant |
| US6881192(B1) | 2005-04-01 | Park | Examiner |
| US7001337(B2) | 2006-02-01 | Dekker | Applicant |
| US7020514(B1) | 2006-03-01 | Wiesel | Applicant |
| US7024233(B2) | 2006-04-01 | Al Ali et al. | Applicant |
| US7024234(B2) | 2006-04-01 | Margulies et al. | Applicant |
| US7039538(B2) | 2006-05-01 | Baker, Jr. | Applicant |
| US7070568(B1) | 2006-07-01 | Koh | Examiner |
| US7081095(B2) | 2006-07-01 | Lynn et al. | Applicant |
| US7082329(B2) | 2006-07-01 | Järverud | Applicant |
| US7160252(B2) | 2007-01-01 | Cho et al. | Applicant |
| US7177686(B1) | 2007-02-01 | Turcott | Applicant |
| US7190261(B2) | 2007-03-01 | Al Ali | Applicant |
| US7225013(B2) | 2007-05-01 | Geva et al. | Applicant |
| US7314451(B2) | 2008-01-01 | Halperin et al. | Examiner |
| US7324845(B2) | 2008-01-01 | Mietus et al. | Applicant |
| US7343198(B2) | 2008-03-01 | Behbehani et al. | Applicant |
| US7351206(B2) | 2008-04-01 | Suzuki et al. | Applicant |
| US7357775(B1) | 2008-04-01 | Koh | Examiner |
| US7386340(B2) | 2008-06-01 | Schlegel et al. | Applicant |
| US7396333(B2) | 2008-07-01 | Stahmann et al. | Applicant |
| US7398115(B2) | 2008-07-01 | Lynn | Applicant |
| US7413549(B1) | 2008-08-01 | Koh | Examiner |
| US7438686(B2) | 2008-10-01 | Cho et al. | Examiner |
| US7468040(B2) | 2008-12-01 | Hartley et al. | Applicant |
| US7479114(B2) | 2009-01-01 | Hartley et al. | Applicant |
| US7510531(B2) | 2009-03-01 | Lee et al. | Applicant |
| US7640055(B2) | 2009-12-01 | Geva et al. | Applicant |
| US7674230(B2) | 2010-03-01 | Reisfeld | Applicant |
| US2001/0018557(A1) | 2001-08-01 | Lynn et al. | Applicant |
| US2002/0002327(A1) | 2002-01-01 | Grant et al. | Applicant |
| US2002/0007124(A1) | 2002-01-01 | Woodward | Applicant |
| US2002/0095076(A1) | 2002-07-01 | Krausman et al. | Applicant |
| US2002/0138013(A1) | 2002-09-01 | Guerrero et al. | Applicant |
| US2002/0173707(A1) | 2002-11-01 | Lynn et al. | Applicant |
| US2002/0190863(A1) | 2002-12-01 | Lynn | Applicant |
| US2003/0000522(A1) | 2003-01-01 | Lynn et al. | Applicant |
| US2003/0004423(A1) | 2003-01-01 | Lavie et al. | Applicant |
| US2003/0004652(A1) | 2003-01-01 | Brunner et al. | Applicant |
| US2003/0055351(A1) | 2003-03-01 | Wiesel | Applicant |
| US2003/0120164(A1) | 2003-06-01 | Nielsen et al. | Applicant |
| US2003/0144597(A1) | 2003-07-01 | Bock | Applicant |
| US2003/0158466(A1) | 2003-08-01 | Lynn et al. | Applicant |
| US2004/0039292(A1) | 2004-02-01 | Schlegel et al. | Applicant |
| US2004/0059236(A1) | 2004-03-01 | Margulies et al. | Applicant |
| US2004/0068199(A1) | 2004-04-01 | Echauz et al. | Applicant |
| US2004/0073098(A1) | 2004-04-01 | Geva et al. | Applicant |
| US2004/0193068(A1) | 2004-09-01 | Burton et al. | Applicant |
| US2004/0230105(A1) | 2004-11-01 | Geva et al. | Applicant |
| US2005/0062609(A9) | 2005-03-01 | Lynn | Applicant |
| US2005/0076908(A1) | 2005-04-01 | Lee et al. | Applicant |
| US2005/0080349(A1) | 2005-04-01 | Okada et al. | Applicant |
| US2005/0148893(A1) | 2005-07-01 | Misczynski et al. | Applicant |
| US2005/0192508(A1) | 2005-09-01 | Lange et al. | Applicant |
| US2005/0240091(A1) | 2005-10-01 | Lynn | Applicant |
| US2006/0041201(A1) | 2006-02-01 | Behbehani et al. | Applicant |
| US2006/0111635(A1) | 2006-05-01 | Todros et al. | Applicant |
| US2006/0149144(A1) | 2006-07-01 | Lynn et al. | Applicant |
| US2006/0155206(A1) | 2006-07-01 | Lynn | Applicant |
| US2006/0155207(A1) | 2006-07-01 | Lynn et al. | Applicant |
| US2006/0161071(A1) | 2006-07-01 | Lynn et al. | Applicant |
| US2006/0189880(A1) | 2006-08-01 | Lynn et al. | Applicant |
| US2006/0195037(A1) | 2006-08-01 | Wiesel | Applicant |
| US2006/0195041(A1) | 2006-08-01 | Lynn et al. | Applicant |
| US2006/0235324(A1) | 2006-10-01 | Lynn | Applicant |
| US2006/0276695(A9) | 2006-12-01 | Lynn et al. | Applicant |
| US2006/0293602(A1) | 2006-12-01 | Clark | Applicant |
| US2007/0016095(A1) | 2007-01-01 | Low et al. | Applicant |
| US2007/0021979(A1) | 2007-01-01 | Cosentino et al. | Applicant |
| US2007/0093721(A1) | 2007-04-01 | Lynn et al. | Applicant |
| US2007/0118054(A1) | 2007-05-01 | Pinhas et al. | Applicant |
| US2007/0129647(A1) | 2007-06-01 | Lynn | Applicant |
| US2007/0149860(A1) | 2007-06-01 | Lynn et al. | Applicant |
| US2007/0149870(A1) | 2007-06-01 | Rosenthal | Applicant |
| US2007/0167843(A1) | 2007-07-01 | Cho et al. | Applicant |
| US2007/0191688(A1) | 2007-08-01 | Lynn | Applicant |
| US2007/0191697(A1) | 2007-08-01 | Lynn et al. | Applicant |
| US2007/0208269(A1) | 2007-09-01 | Mumford et al. | Examiner |
| US2007/0213620(A1) | 2007-09-01 | Reisfeld | Applicant |
| US2007/0213621(A1) | 2007-09-01 | Reisfeld et al. | Applicant |
| US2007/0213622(A1) | 2007-09-01 | Reisfeld | Applicant |
| US2007/0213624(A1) | 2007-09-01 | Reisfeld et al. | Applicant |
| US2007/0239057(A1) | 2007-10-01 | Pu et al. | Examiner |
| US2007/0244408(A1) | 2007-10-01 | Wingeier et al. | Applicant |
| US2008/0177195(A1) | 2008-07-01 | Armistead | Applicant |
| US2008/0200775(A1) | 2008-08-01 | Lynn | Applicant |
| US2008/0200819(A1) | 2008-08-01 | Lynn et al. | Applicant |
| US2008/0269583(A1) | 2008-10-01 | Reisfeld | Applicant |
| US2008/0287756(A1) | 2008-11-01 | Lynn | Applicant |
| US2008/0300500(A1) | 2008-12-01 | Reisfeld | Applicant |
| US2008/0312166(A1) | 2008-12-01 | Lynn et al. | Applicant |
| US2009/0281838(A1) | 2009-11-01 | Lynn et al. | Applicant |
| US2009/0281839(A1) | 2009-11-01 | Lynn et al. | Applicant |
| US2009/0292478(A1) | 2009-11-01 | Avinash et al. | Applicant |
| US2009/0299150(A1) | 2009-12-01 | Alberte, Jr. et al. | Applicant |
| US2009/0306482(A1) | 2009-12-01 | Davis et al. | Applicant |
| US2009/0312612(A1) | 2009-12-01 | Rantala | Applicant |
| WO96/16591 | 1996-06-01 | Applicant | |
| WO/10455 | 2000-03-01 | Applicant | |
| WO3/057025 | 2003-07-01 | Applicant | |
| WO2004/075746 | 2004-09-01 | Applicant | |
| WO2006/008743 | 2006-01-01 | Applicant | |
| WO2006/054306 | 2006-05-01 | Applicant | |
| WO2006/082589 | 2006-08-01 | Applicant | |
| WO2008/029399 | 2008-03-01 | Applicant | |
| WO2009/137682 | 2009-11-01 | Applicant |
Non-Patent Literature (98)
- Lempel et al., “A Universal Algorithm for Sequential Data Compression”, IEEE Transactions on Information Theory, IT-23:3 (1977), pp. 337-349.Applicant
- Akselrod et al., “Power Spectrum Analysis of Heart Rate Fluctuation: A Quantitative Probe of Beat-to-Beat Cardiovascular Control”, Science 213 (1981), pp. 220-222.Applicant
- Eamonn Keogh, et al., “An Online Algorithm for Segmenting Time Series”, pp. 289-296, 2001.Applicant
- Penzel, et al, “Computer Based Sleep Recording and Analysis”, Sleep Medicine Reviews 4:2 (2000), pp. 131-148.Applicant
- Tesler, et al, “Can One Detect Sleep Stage Transitions for On-Line Sleep Scoring by Monitoring the Heart Rate Variability?” Somnologie 8 (2004), pp. 33-41.Applicant
- A.K. Jain, et al, “Data Clustering: A Review”, ACM Computing Surveys, vol. 31, No. 3, Sep. 1999.Applicant
- http://www.sleepdisorderchannel.net, 2004.Applicant
- M.G. Terzano, et al, “Atlas, rules, and recording techniques for the scoring of cyclic alternating pattern (CAP) in human sleep”, Sleep Medicine 2 (2001), pp. 537-553.Applicant
- T.D. Bradley, et al., “Sleep apnea and heart failure”, Circulation 2003; 107:1671-1678.Applicant
- J.L. Peppin, et al., “Cheyne stokes respiration with central sleep apnoea in chronic heart failure: proposals for a diagnostic and therapeutic strategy”, Sleep medicine reviews (2006) 10, 33-47.Applicant
- U. Cora, et al., “Sleep and exertional periodic breathing in chronic heart failure”, Circulation 2006; 113:44-50.Applicant
- T. Penzel, et al., “Systematic comparison of different algorithms for apnoea detection based on electronic recordings”, Med. Biol. Eng. Comput., 2002, 40, 402-407.Applicant
- H.W. Duchna, et al., “Sleep-disordered breathing and cardio and cerebrovascular diseases: 2003 update of clinical significance and future perspectives”, Somnlogie 7: 101-121, 2003.Applicant
- M. Folke, et al., “Critical review of non-invasive respiratory monitoring in medical care”, Medical & Biol Eng & Comp 2003, vol. 41.Applicant
- U.S. Appl. No. 60/651,295.Applicant
- “Heart rate variability”, European Heart Journal 1996, 17, 354-381.Applicant
- J.J. Verbeek, et al., “Efficient greedy learning of Gaussian mixture models”, Neural Computation, 5(2): 469-485, Feb. 2003.Applicant
- L.R. Rabiner, “A tutorial on hidden markov models and selected applications in speech recognition”, Proc. of the IEEE, vol. 77, No. 2, Feb. 1989.Applicant
- P. Jourdain, et al., “Plasma brain natriuretic peptide-guided therapy to improve outcome in heart failure: the STARS-BNP multicenter study”, JACC vol. 49, No. 16, 2007.Applicant
- “Sleep related breathing disorders in adults: recommendations for syndrome definition and measurement techniques in clinical research”, SLEEP, vol. 22, No. 5, 1999.Applicant
- G.A. Wright, et al., “Natriuretic peptides as a prognostic marker and therapeutic target in heart failure”, Heart 2006; 92:149-151.Applicant
- W. Rosamond, et al., “Heart disease and stroke statistics—2007 update: A report from the American heart association statistics committee and stroke statistics subcommittee”, Circulation Feb. 6, 2007.Applicant
- C.A. Kushhida, et al., “Practice parameters for the indications for polysomnography and related procedures: an update for 2005”, SLEEP vol. 28, No. 4, 2005.Applicant
- G. Tusman, et al., “Effect of pulmonary perfusion on the slopes of single-breath test of CO2”, J Appl Physiol 99:650-655, 2005.Applicant
- International Preliminary Report on Patentability Dated Apr. 16, 2009 From the International Bureau of WIPO Re.: Application No. PCT/IL2007/001092.Applicant
- International Search Report Dated Jul. 22, 2008 From the International Searching Authority Re.: Application No. PCT/IL07/01092.Applicant
- Notice of Allowance Dated Feb. 2, 2007 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/678,773.Applicant
- Notice of Allowance Dated Sep. 3, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/677,176.Applicant
- Notice of Allowance Dated Nov. 9, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,221.Applicant
- Notice of Allowance Dated Apr. 15, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/995,817.Applicant
- Official Action Dated Oct. 1, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/677,176.Applicant
- Official Action Dated Apr. 3, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/677,176.Applicant
- Official Action Dated Sep. 3, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,173.Applicant
- Official Action Dated Feb. 5, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,173.Applicant
- Official Action Dated Feb. 5, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,250.Applicant
- Official Action Dated Jun. 6, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/677,176.Applicant
- Official Action Dated Oct. 8, 2009 From the US Patent and Trademark Office re.: U.S. Appl. No. 10/677,176.Applicant
- Official Action Dated Mar. 9, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,221.Applicant
- Official Action Dated Sep. 9, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,221.Applicant
- Official Action Dated Aug. 13, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,173.Applicant
- Official Action Dated May 14, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/995,817.Applicant
- Official Action Dated Oct. 14, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,250.Applicant
- Official Action Dated Aug. 15, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,250.Applicant
- Official Action Dated Jan. 21, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/995,817.Applicant
- Official Action Dated Jul. 22, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,221.Applicant
- Official Action Dated Jan. 23, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/995,817.Applicant
- Official Action Dated Jul. 31, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/995,817.Applicant
- Response Dated Jun. 1, 2009 to Official Action of Feb. 5, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,173.Applicant
- Response Dated Jun. 2, 2009 to Official Action of Feb. 5, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,250.Applicant
- Response Dated Jun. 3, 2009 to Official Action of Apr. 3, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/677,176.Applicant
- Response Dated Mar. 13, 2009 to Official Action of Jan. 21, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/995,817.Applicant
- Response Dated Oct. 17, 2008 to Official Action of Jul. 31, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/995,817.Applicant
- Response Dated Feb. 19, 2008 to Official Action of Jan. 23, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/995,817.Applicant
- Response Dated May 21, 2008 to Official Action of May 14, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/995,817.Applicant
- Response Dated Oct. 24, 2008 to Official Action of Aug. 13, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,173.Applicant
- Response Dated Oct. 24, 2008 to Official Action of Aug. 15, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,250.Applicant
- Response Dated Dec. 28, 2008 to Official Action of Oct. 1, 2008 From the US Patent and Trademark Office Re.: U.S. Appl. No. 10/677,176.Applicant
- Supplementary European Search Report and the European Search Opinion Dated Feb. 2, 2009 From the European Patent Office Re.: Application No. 05761327.5.Applicant
- Supplementary European Search Report and the European Search Opinion Dated Oct. 8, 2009 From the European Patent Office Re.: Application No. 06704643.3.Applicant
- Written Opinion Dated Jul. 22, 2008 From the International Searching Authority Re.: Application No. PCT/IL07/01092.Applicant
- Agarwal et al. “Automatic EEG Analysis During Long-Term Monitoring in the ICU”, Electroencephalography and Clinical Neurophysiology, XP002510061, 107(1): 44-58, Jul. 1998. p. 45, r-h col., Line 25-p. 50, l-h col., Line 5, p. 50, l-h col., Line 29-p. 51, r-h col., Line 24, p. 54, r-h col., Line 32-p. 55, r-h col., Line 21.Applicant
- Agarwal et al. “Computer-Assisted Sleep Staging Based on Segmentation and Clustering”, Proceedings of 23rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society, 2001 Conference Proceedings (EMBS), Istanbul, Turkey, Oct. 25-28, 2001, XP010594756, 2: 1695-1698, Oct. 25, 2001. p. 1695, r-h col., Line 27-p. 1697, l-h col., Line 36.Applicant
- Agarwal et al. “Computer-Assisted Sleep Staging”, IEEE Transactions on Biomedical Engineering, XP011007173, 48(12): 1412-1423, Dec. 1, 2001. p. 1413, l-h col., Line 26-p. 1418, r-h col., Line 14.Applicant
- Akselrod et al. “Power Spectrum Analysis of Heart Rate Fluctuation: A Quantitative Prode of Beat-to-Beat Cardiovascular Control”, Science, 213: 220-222, Jul. 10, 1981.Applicant
- Bodenstein et al. “Feature Extraction From the Electroencephalogram by Adaptive Segmentation”, Proceedings of the IEEE, 65(5): 642-652, May 1977.Applicant
- Gath et al. “Unsupervised Optimal Fuzzy Clustering”, IEEE Transactions on Pattern Analysis and Machine Intelligence, 11(7): 773-781, Jul. 1989.Applicant
- Geva et al. “Forecasting Generalized Epileptic Seizures From the EEG Signal by Wavelet Analysis ans Dynamic Unsupervised Fuzzy Clustering”, IEEE Transactions on Biomedical Engineering, 45(10): 1205, Oct. 1998. Abstract.Applicant
- Healthcommunities “Sleep Stages”, Sleepchannel.com, Retrieved From the Internet: http://www.sleepdisorderchannel.net/stages, 2004.Applicant
- Hutt et al. “Analysis and Modeling of Quasi-Stationary Multivariate Time Series and Their Application to Middle Latency Auditory Evoked Potentials”, Physica D, XP002510063, 177(1-4): 203-232, Mar. 15, 2003.Applicant
- Jain et al. “Data Clustering: A Review”, ACM Computing Surveys, 31(3): 264-323, Sep. 1999.Applicant
- Kaplan et al. “Macrostructural EEG Characterization Based on Nonparametric Change Point Segmentation: Application to Sleep Analysis”, Journal of Neurscience Methods, XP002510060, 106(1): 81-90, Mar. 30, 2001. p. 83, r-h col., Line 30-p. 86, l-h col., Line 2, p. 86, r-h col., Lines 18-54, p. 88, l-h col., Line 2-p. 90, l-h col., Line 38.Applicant
- Keogh et al. “An Online Algorithm for Segmenting Time Series”, Proceedings of the 2001 IEEE International Conference on Data-Mining, ICDM, p. 289-296, 2001.Applicant
- Klösch et al. “The SIESTA Project Polygraphic and Clinical Database. A New Approach to Studying Subjective and Objective Measurements of Human Sleep”, IEEE Engineering in Medicine and Biology Magazine: The Quarterly Magazine of the Engineering in Medicine & Biology Society, XP002510064, 20(3): 51-57, May 2001.Applicant
- Lempel et al. “A Universal Algorithm for Sequential Data Compression”, IEEE Transactions on Information Theory, 23(3): 337-349, May 1977.Applicant
- Malik et al. “Heart Rate Variability. Standards of Measurement, Physiological Interpretation, and Clinical Use”, European Heart Journal, 17: 354-381, 1996.Applicant
- Murphy “Hidden Markov Model (HMM) Toolbox for Matlab”, 1 P., 1998. Retrieved From the Internet: http://www.cs.ubc.ca/˜murphyk/Software/HMM/hmm.mtl.Applicant
- Paul et al. “Comparison of Quantitative EEG Characteristics of Quiet and Active Sleep in Newborns”, Sleep Medicine, XP002510062, 4(6): 543-552, Nov. 2003. p. 544, r-h col., Line 21-p. 548, l-h col., Line 38.Applicant
- Penzel et al. “Computer Based Sleep Recording and Analysis”, Sleep Medicine Review, 4(2): 131-148, 2000.Applicant
- Penzel et al. “Systematic Comparison of Different Algorithms for Apnoea Detection Based on Electrocardiogram Recordings”, Medical & Biological Engineering & Computing, 40: 402-407, 2002.Applicant
- Pépin et al. “Cheyne Stokes Respiration With Central Sleep Apnoea in Chronic Heart Failure: Proposals for A Diagnostic and Therapeutic Strategy”, Sleep Medicine Reviews, 10(1): 33-47, 2006.Applicant
- Roberts et al. “New Method of Automated Sleep Quantification”, Medical & Biological Engineering & Computing, 30: 509-517, Sep. 1992.Applicant
- Sackellares et al. “Epilepsy—When Chaos Fails in Chaos in the Brain”, Eds, 1999.Applicant
- Telser et al. “Can One Detect Sleep Stage Transitions for On-Line Sleep Scoring by Monitoring the Heat Rate Variability?”, Somnologie, 8: 33-41, 2004.Applicant
- Terzano et al. “Atlas, Rules, and Recording Techniques for the Scoring of Cyclic Alternating Pattern (CAP) in Human Sleep”, Sleep Medicine, 2: 537-553, 2001.Applicant
- Tusman et al. “Effect of Pulmonary Perfusion on the Slopes of Single-Breath Test of CO2”, Journal of Applied Physiology, 99(2): 650-655, 2005.Applicant
- Verbeek et al. “Efficient Greedy Learning of Gaussian Mixture Models”, Neural Computation, 15(2): 469-485, Feb. 2003.Applicant
- Xu-Sheng et al. “EEG Complexity As A Measure of Depth of Anesthesia for Patients”, IEEE Transactions on Biomedical Engineering, 48(12):1424-1433, Dec. 2001. Abstract.Applicant
- Ziv et al. “A Universal Algorithm for Sequential Data Compression”, IEEE Transactions on Information Theory, IT-23(3): 337-343, May 1977.Applicant
- Communication Pursuant to Article 94(3) EPC Dated Feb. 8, 2010 From the European Patent Office Re.: Application No. 05808274.4.Applicant
- Response Dated Jan. 4, 2010 to Official Action of Sep. 3, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,173.Applicant
- Supplementary European Search Report and the European Search Opinion Dated Nov. 3, 2009 From the European Patent Office Re.: Application No. 05808274.4.Applicant
- Covidien “Covidien Launches Alarm Management Systm for the Nellcor(™) OxiMax(™) N-600x(™) Pulse Oximeter”, Covidien Bussiness Wire, Apr. 16, 2009.Applicant
- Redmond et al. “Electrocardiogram-Based Automatic Sleep Staging in Sleep Disordered Breathing”, Computers in Cardiology, XP010698978, 30: 609-612, Sep. 21, 2003. Section 3, Feature Extraction, Section 5, Subject Specific Calssification, Section 6, Subject Independant Classification.Applicant
- Watanabe et al. “Estmation of Sleep Stages Based on Heart Rate Fluctuation and Body Movement”, SICE 2004 Annual Conference, IEEE, XP010824823, 3: 2153-2156, Aug. 4, 2004. Section 2.2 Measurement System, Section 4.1 Indices, Figs.1, 2.Applicant
- Response Dated Mar. 15, 2010 to Official Action of Oct. 14, 2009 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,222.Applicant
- Notice of Allowance Dated May 3, 2010 From the US Patent and Trademark Office Re.: U.S. Appl. No. 11/750,173.Applicant
- Notice of Allowance Dated Jun. 29, 2010 From the U.S. Patent and Trademark Office Re.: U.S. Appl. No. 11/750,222.Applicant
- Office Action Dated Jun. 7, 2010 From the Israel Patent Office Re.: Application No. 147502 and Its Translation Into English.Applicant