US 7,992,557 B2Grant
Gas supply device for sleep apnea
Issue Date:2011-08-09
•18 Claims
•13 Drawing Sheets
Abstract
An apparatus is provided for controlling air supplied under pressure to a patient suffering from sleep disorders such as apnoea. The apparatus measures and/or calculates the inspiratory air flow, volume and pressure to the airway. Controlled pressurized air is supplied to the patient's upper anatomical airway, wherein the apparatus measures the air flow and pressure to the airway. The apparatus determines whether to increase or decrease the pressure to the patient's airway based on a determination of one or more factors, such as the occurrence of hypopnoea, hyperventilation, obstructive or central apnoea, air leakage and acoustical vibrations. Occurrences of events representing a sleep problem may be stored by the apparatus and retrieved by a clinician at a later date.
Metadata
Assignee
- Covidien AG
Inventors
- Hossein Nadjafizadeh
- Pascal Nicolazzi
- Véronique Griller Lanoir
Application Information
Application Number:US 12/106,778
Filing Date:2008-04-21
Priority Date:1999-02-12
Art Unit:3771
Classifications
IPC:
F16K31/02A61M16/00A62B7/00
Patent Drawings (13 sheets)
Description
Cross Reference to Related Applications
[0001] This application is a continuation of U.S. application Ser. No. 10/968,403 filed Oct. 18, 2004 now U.S. Pat. No. 7,370,650, which is a continuation of U.S. application Ser. No. 09/913,237 filed Nov. 13, 2001, now U.S. Pat. No. 6,814,074, which is an application filed under 35 U.S.C. §371, and claims priority to WIPO International Application Serial No. PCT/FR00/00334 filed Feb. 10, 2000, which claims priority to the French Application Serial No. 99/06515 filed Feb. 12, 1999.
Background of the Invention
[0002] The invention concerns a method of controlling an apparatus for supplying air pressure to a patient suffering from sleep problems.
[0003] The invention also concerns an apparatus for supplying air pressure to a patient suffering from sleep problems.
[0004] These sleep problems are respiratory and tend to waken the patient inopportunely.
[0005] They are for example apnoeas, hypopnoeas, acoustic vibrations or snores, or limitation of the respiratory flow, due to a narrowing of the upper airways of the patient.
[0006] The document U.S. Pat. No. 5,458,137 describes a method and a device for controlling respiration in the case of sleep problems, which use multiple and variable pressure levels.
[0007] A pressure source supplies a breathable gas compressed at a relatively low pressure to the airways of the user.
[0008] Pressure sensors monitor the pressures and convert them into electrical signals.
[0009] The electrical signals are filtered and processed in order to extract specific characteristics such as the duration and energy levels.
[0010] If these characteristics exceed chosen duration and energy level thresholds beyond a minimum time period, the microprocessor indicates the presence of a sleep respiratory problem.
[0011] If a chosen number of these events appears during a chosen time period, the microprocessor adjusts the pressure supplied by the source.
[0012] The document U.S. Pat. No. 5,490,502 describes a method and an apparatus for optimizing the controlled positive pressure in order to minimize the air flow coming from a generator while ensuring that flow limitation in the airways of the patient does not take place.
[0013] Provision is made therein to detect flow limitation by analysing a respiratory flow wave.
[0014] As soon as the presence of a flow limitation has been analysed, the system determines an action to be performed for adjusting the controlled positive pressure.
[0015] The pressure is increased, reduced or maintained depending on whether flow limitation has been detected and according to the previous actions implemented by the system.
[0016] The documents U.S. Pat. No. 5,335,654, l EP-A-661 071 and EP-A-651 971 should also be cited.
Brief Summary of the Invention
[0017] The invention aims to improve the methods and devices of the state of the art, to automatically and continuously adapt the delivered pressure to the state of the patient and to anticipate and prevent the appearance of problems.
[0018] A first object of the invention is a method of controlling an apparatus for supplying air pressure to a patient suffering from sleep problems such as apnoeas.
[0019] A second object of the invention is an apparatus for supplying air pressure to a patient suffering from sleep problems such as apnoea, implementing the supply method.
[0020] The patient wears a mask by means of which air under pressure is supplied to his upper airways by the apparatus.
[0021] According to the invention, a control algorithm is provided using an output flow signal from the apparatus for detecting apnoea, hypopnoea, flow limitation events and leakages, and using the analysis of an item of pressure information for determining the presence of snoring, also referred to as acoustic vibrations.
[0022] The pressure supplied to the upper airways of the patient by the apparatus can be maintained constant, be increased or reduced according to the determination of the event which has been performed by the control algorithm.
[0023] Thus, if no respiration is detected by the control algorithm within a predetermined minimum time depending on a calculated mean respiration time, the presence of an apnoea is determined.
[0024] This predetermined minimum apnoea detection time is for example equal to a time constant, for example 10 seconds, added to a proportionality factor multiplied by the calculated mean respiration time, this factor being for example equal to 5/8.
[0025] For each apnoea, the output flow signal is amplified and filtered in order to determine the presence or absence of cardiac oscillations.
[0026] If cardiac oscillations were detected during the last elapsed time interval, for example equal to 5 seconds, then the apnoea is classified as being central and no control takes place in the algorithm.
[0027] If no cardiac oscillation was detected in this time interval, the apnoea is classified as being obstructive, and the pressure is increased by a predetermined value a first time and, during the same apnoea, twice more regularly, for example every 15 seconds.
[0028] The control algorithm compares peak-to-peak flow variations during the latest respiration of the patient with respect to a predetermined number of previous respirations, for example equal to 8.
[0029] After each respiration, a classification is performed into:
-
- [0030] normal respiration, if the last peak-to-peak flow value is within a given range with respect to the mean value over the previous 8 respirations, for example from 40% to 150% or 140% thereof;
- [0031] hypopnoeic respiration, if the last flow value is below this range;
- [0032] hyperpnoeic respiration, if the last flow value is above this range.
[0033] A hypopnoea determination is made if hypopnoeic respiration detection takes place during at least a given time, for example 10 seconds, and terminates after a given number of normal or hyperpnoeic respirations, for example equal to 2.
[0034] A hypopnoea determination causes a given increase in pressure, for example 1 cm H2O first, and then, during the same hypopnoea, an increase in pressure by another given value, regularly, for example 0.5 cm H2O every two hypopnoeic respirations.
[0035] The control algorithm analyses and compares, respiration by respiration, the waveform of the respiratory flow with a sinusoidal waveform of the same period and same gradient.
[0036] After the comparison based on two flow form criteria, each respiration is first classified as normal, intermediate or limited flow.
[0037] A final classification, based on the combination of the flow classification and the occurrence of snores, changes the classification of respirations from normal into intermediate, respectively from intermediate into limited flow respiration.
[0038] Processing is decided upon when a certain number, for example 2, of successive limited flow respirations or a certain number, for example 5, of successive intermediate respirations take place after for example two normal respirations.
[0039] This processing causes a given increase in pressure, repeated regularly a certain number of times, for example 0.3 cm H2O three times every two respirations.
[0040] For each respiration, the pressure signal is amplified and filtered in order to detect the presence or absence of acoustic vibrations or snoring.
[0041] A determination of a valid snore is made by the control algorithm if the detected acoustic vibration occurred at least for a certain time, for example 7% of the mean duration of the last three respirations, and with a period less than a factor proportional to this mean time, for example 120% thereof.
[0042] In the case of a valid snore, the algorithm increases the pressure by a given value, for example 1 cm H2O, if the last control due to a snore took place more than a given time previously, for example 1 minute.
[0043] A mean leakage is determined as being equal to the mean flow during respiration.
[0044] The control algorithm continuously compares the current leakage with a leakage limit, it being possible to regulate said limit from the pressure.
[0045] If the current leakage exceeds the limit, all pressure increase controls generated following event detections are disabled.
[0046] After detection of an apnoea or a snoring event or a hypopnoea control or a processing decision, the algorithm will reduce the pressure by a given value, for example 0.5 cm H2O, in a first step after a given time, for example 5 minutes, and regularly for the following reductions, for example every minute.
[0047] A given maintenance pressure, for example 8 cm H2O, is supplied by the apparatus if no respiration has been detected during a given time, for example two minutes, or if the pressure supplied has been greater than or equal to a given value for a given time, for example 17 cm H2O for 10 or 30 minutes.
[0048] One advantage of the method is an automatic adaptation of the detection criteria to the respiratory characteristics of the patient.
[0049] Thus, any modification of the respiratory rhythm is taken into account by the algorithm for performing the detection.
[0050] The fact of involving a mean value of respiratory cycle time over a certain number of previous respiratory cycles has the effect of variations in the cycle and respiratory amplitude being tracked regularly and better detection.
[0051] The invention will be better understood from a reading of the following description, given with reference to the figures.
Brief Description of the Several Views of the Drawings
[0052] FIG. 1 is a diagram showing the apparatus for supplying air pressure to the patient.
[0053] FIG. 2 depicts an algorithm for taking a decision with a view to a first pressure increase control.
[0054] FIG. 3 depicts an algorithm for indicating the appearance of problems.
[0055] FIG. 4 depicts a respiration designation algorithm.
[0056] FIG. 5 depicts an algorithm for detecting central and obstructive apnoea and for pressure control according to the result of these detections, as well as an algorithm for reducing pressure according to the previous appearance or not of events representing sleep problems.
[0057] FIG. 6 depicts an algorithm for designating cycles as normal ventilation, hyperventilation or hypoventilation.
[0058] FIG. 7 depicts a hypopnoeic respiration detection algorithm.
[0059] FIG. 8 depicts a hyperpnoeic respiration detection algorithm.
[0060] FIG. 9 depicts a normal respiration detection algorithm.
[0061] FIG. 10 depicts a high pressure detection algorithm.
[0062] FIG. 11 depicts a mask leakage detection algorithm.
[0063] FIG. 12 depicts an acoustic vibration detection algorithm.
[0064] FIG. 13 depicts an algorithm for reducing pressure in the event of acoustic vibration detection.
Detailed Description of the Invention
[0065] In FIG. 1 , the apparatus for supplying air pressure to a patient has a central processing and pressure control unit U, a controlled pressure supply module MPD, a mask MVA for the upper airways of the patient, and a tube CF for supplying air pressure from the module MPD to the mask MVA.
[0066] The air flow supplied to the patient and the air pressure prevailing in the mask MVA are measured by means of a supplied air flow sensor CDAF, connected to the central unit U, and by means of a sensor CPM of pressure in the mask MVA, connected to the central unit U.
[0067] It is determined from the measured variables whether or not events representing sleep problems appear.
[0068] The algorithms of the method according to the invention are implemented by software integrated in the central unit U.
[0069] In FIG. 2 , it is determined from the measured variables whether the current respiratory cycle of the patient corresponds to a predetermined valid respiratory cycle.
[0070] A problem appearance indicator BLN is set to a first problem appearance state ON, if the appearance of one or more of the events representing sleep problems is determined.
[0071] The indicator BLN is set to a second problem absence state OFF, if the appearance of events representing sleep problems is not determined.
[0072] A count is made of a first number CCAR of valid respiratory cycles determined since the last pressure control.
[0073] A count is made of a second number CCON of valid respiratory cycles determined since the last change of the indicator BLN to the first state ON.
[0074] A count is made of a third number RC of successive changes of the indicator BLN from the second state OFF to the first state ON.
[0075] When the indicator BLN is in the first state ON, a first given increase of supplied air pressure is controlled, by means of the control C1, when all the following are true:
-
- [0076] the current respiratory cycle has been determined as being valid;
- [0077] the first number CCAR is greater than a first predetermined integer number RP;
- [0078] the second number CCON corresponds to one or more other second predetermined integer numbers N;
- [0079] the third number RC is greater than or equal to a third predetermined integer number X.
[0080] When the indicator BLN changes from the second state OFF to the first state ON, the first given increase of supplied air pressure is controlled, by means of the control C1, when, solely all the following are true:
-
- [0081] the current respiratory cycle has been determined as being valid;
- [0082] the first number CCAR is greater than a first predetermined integer number RP;
- [0083] the third number RC is greater than or equal to a third predetermined integer number X.
[0084] In one embodiment, the second integer numbers N are between 1 and 300.
[0085] In another embodiment, the second integer numbers N are the first three multiples of a given integer N0.
[0086] In another embodiment, the second integer numbers N are respectively 2, 4 and 6, N0 being equal to 2.
[0087] In another embodiment, the first predetermined integer number RP is between 1 and 255.
[0088] In another embodiment, the first predetermined integer number RP is equal to 10.
[0089] In another embodiment, the third predetermined integer number X is between 1 and 100.
[0090] In another embodiment, the third predetermined integer number X is equal to 1.
[0091] In another embodiment, the first given pressure increase control C1 is less than +10 mbar.
[0092] In another embodiment, the first given pressure increase control C1 is substantially equal to +0.3 mbar.
[0093] The first and third numbers CCAR; RC of counted valid respiratory cycles and counted changes are reset to 0, after the second counted number CCON of valid cycles has reached the largest of the second predetermined integer numbers N.
[0094] The second counted number CCON is reset to 0 when the indicator BLN changes from the second state OFF to the first state ON.
[0095] The predetermined valid respiratory cycle corresponds to a maximum respiratory flow greater than a predetermined flow value such as 50 ml/s, an inspiratory volume greater than a predetermined volume value such as 0.05 litres and an absence of saturation at flow detection time.
[0096] In FIG. 3 , in order to give the state ON or OFF to the problem appearance indicator BLN,
-
- [0097] when the apparatus is started up, a state variable ER is initialized to a third processing absence state NIR and the indicator BLN is initialized to the second state OFF.
[0098] Then sequentially,
-
- [0099] from the measured variables, the respiratory cycles are designated as belonging to different categories such as limited flow cycle, intermediate cycle, normal cycle and invalid cycle, each corresponding respectively to weightings RSV0, REV0; RSV1, REV1; RSV2, REV2; 0,0;
- [0100] the weightings of the category of the currently designated cycle are assigned to first and second weighting accumulators SV; EV;
- [0101] if the designated cycle belongs to the invalid cycle category, the state variable ER is reset to the third state NIR and the indicator BLN is reset to the second state OFF and a first counter FLC is initialized to a predetermined value.
[0102] If the state of the state variable ER corresponds to the third state NIR:
-
- [0103] if the value of a first accumulator SV is less than a first comparative value, the counter FLC is reinitialized to its predetermined value;
- [0104] if the value of the first accumulator SV is substantially equal to its first comparative value, no action is taken and the next test is passed to;
- [0105] if the value of the first accumulator SV is greater than its first comparative value, the state variable ER is changed to a fourth processing possibility state PR and the indicator BLN is set to the second state OFF.
[0106] If the state of the state variable ER corresponds to the fourth state PR and
-
- [0107] if the value of the first accumulator SV is less than its first comparative value, the first counter FLC is reinitialized to its predetermined value, and the state variable ER and the indicator BLN are reset respectively to the third and second states NIR; OFF;
- [0108] if the value of the first accumulator SV is substantially equal to its first comparative value, no action is taken and the next test is passed to;
- [0109] if the value of the first accumulator SV is greater than its first comparative value, the first counter FLC is made to take its previous value with the value of the first accumulator SV added to it, and if then the value of the first counter FLC is greater than or equal to a predetermined high stop RMS:
- [0110] a second counter NC is reinitialized to a predetermined value;
- [0111] the state variable ER is changed to a fifth processing state IR; and
- [0112] the indicator BLN is changed to the first state ON.
[0113] If the state of the state variable ER corresponds to the fifth processing state IR:
-
- [0114] if the value of the second accumulator EV is greater than a second comparative value, the second counter NC is made to take its previous value with the value of the second accumulator EV added to it, and
- [0115] if then the value of the second counter NC is greater than or equal to a low stop RME, the state variable ER and the indicator BLN are reset respectively to the third and second states NIR; OFF and the first and second counters FLC; NC are reinitialized to their predetermined respective values;
- [0116] or otherwise, the indicator BLN is changed to its first state ON;
- [0117] if the value of the second accumulator EV is less than its second comparative value, the second counter NC is reinitialized to its predetermined respective value and the indicator BLN is changed to the first state ON;
- [0118] if the value of the second accumulator EV is substantially equal to its second comparative value, no action is taken.
[0119] In one embodiment, the weightings RSV2, REV2; RSV1, REV1; RSV0, REVO; 0,0 corresponding to the normal cycle, intermediate cycle, limited flow cycle and invalid cycle categories, are respectively substantially equal to −1; 1; 5 and 0 for the first accumulator SV and are respectively substantially equal to 1; −1; −1 and 0 for the second accumulator EV.
[0120] The first and second comparative values and the predetermined initialization values of the first and second counters FLC; NC are each substantially equal to 0.
[0121] The high and low stops RMS; RME are respectively substantially equal to 10 and 2.
[0122] In FIG. 4 , the measured respiratory cycles are designated.
[0123] The predetermined valid respiratory cycle corresponds to a maximum inspiratory flow greater than a predetermined flow value such as 50 ml/s, an inspiratory volume greater than a predetermined volume value such as 0.05 litres, an absence of saturation at flow detection time, a measured inspiratory time within a predetermined interval such as 0.5 seconds to 6 seconds and a measured respiratory cycle duration within another predetermined interval such as 1.5 seconds to 20 seconds.
[0124] If the measured respiratory cycle is determined as being valid, then
-
- [0125] a calculation is made of an equivalent sinusoidal curve meeting predetermined characteristics with respect to the inspiratory curve of the measured inspiratory cycle;
- [0126] a calculation is made of a surface criterion CS proportional to the ratio of the area delimited by the inspiratory curve to the area delimited by the equivalent sinusoidal curve, each being taken over the same time interval, within the inspiratory phase of the measured respiratory cycle;
- [0127] a calculation is made of a criterion of correlation CC between the inspiratory curve of the measured inspiratory cycle and the equivalent sinusoidal curve;
- [0128] if the calculated correlation criterion CC is greater than or equal to a first predetermined normal limit LN, and if the calculated surface criterion CS is greater than a second predetermined surface limit LS, the measured respiratory cycle is designated as normal and otherwise, it is designated as a limited flow cycle.
[0129] If the measured respiratory cycle was designated as a limited flow cycle,
-
- [0130] if the calculated surface criterion CS is greater than a third predetermined expert limit LE, the measured respiratory cycle is redesignated as normal,
- [0131] or otherwise,
- [0132] if the calculated surface criterion CS is greater than a fourth predetermined flow limit LD, the measured respiratory cycle is redesignated as intermediate,
- [0133] and in the contrary case, it is designated as a limited flow cycle.
[0134] The second surface limit LS, the fourth flow limit LD and the third expert limit LE are predetermined in an ascending order.
[0135] The predetermined characteristics of the equivalent sinusoidal curve comprise a half period substantially equal to the measured inspiratory time and a gradient at the origin substantially equal to that of the inspiratory curve when it reaches substantially one third of its maximum amplitude.
[0136] In one embodiment, the calculated surface criterion CS is substantially equal to one hundred times the ratio of the areas each taken from substantially one quarter to three quarters of the duration of the inspiratory phase of the measured respiratory cycle.
[0137] The calculated correlation criterion CC is substantially equal to the maximum of one hundred times the coefficients of correlation between the inspiratory curve and the equivalent sinusoidal curve taken respectively over the second half of the inspiratory phase and over the whole thereof.
[0138] The first, second, fourth and third limits LN; LS; LD; LE are respectively between 45 and 100; 0 and 100; 0 and 100; 0 and 100 and are for example substantially equal to 87; 40; 60 and 90 respectively.
[0139] In FIG. 5 , obstructive apnoeas and central apnoeas are detected.
[0140] The algorithm depicted in FIG. 5 is performed during each of a number (NINT) of predetermined consecutive time intervals TAC(j).
[0141] The predetermined consecutive time intervals TAC(j) are those within a predetermined apnoea detection period PDAC.
[0142] In this algorithm, there are detected, for example by hardware means such as analogue or digital filters, the oscillations of the measured flow curve, which are of frequencies within a frequency range P2.
[0143] Then it is detected whether the amplitude of the detected oscillations of the measured flow curve goes successively above and then below a first predetermined central apnoea threshold SAC or whether this amplitude remains less than the first central apnoea threshold SAC, as depicted schematically at the right of FIG. 5 by:
-
- [0144] the behaviour of an obstructive apnoea flow curve (curve constantly below the first threshold SAC);
- [0145] the behaviour of a central apnoea flow curve (curve going a number of times successively above and then below the first threshold SAC).
[0146] In the presence of at least one detection of a passage above and then below the first threshold SAC, a central apnoea detection CAC(D) is counted.
[0147] Then, at each apnoea detection period PDAC,
-
- [0148] the sum SIG is performed of the numbers CAC(i) of central apnoea detections counted, successively over the last (D+1) apnoea detection periods;
- [0149] a second predetermined increase of delivered air pressure is controlled C2 if the sum SIG of the numbers CAC(i) of counted detections is less than or equal to a second predetermined central apnoea designation threshold SQAC;
- [0150] a maintenance of delivered air pressure is controlled, if the sum SIG of the numbers CAC(i) of counted detections is greater than the second threshold SQAC.
[0151] In one embodiment, the second central apnoea designation threshold SQAC is between 0 and 50, and is for example substantially equal to 10.
[0152] The predetermined consecutive time intervals TAC(j) correspond to ten (NINT) consecutive time intervals each of substantially 100 ms, the apnoea detection period PDAC corresponding substantially to 1 second.
[0153] The second pressure increase control C2 is between 1 and 10 mbar and is for example substantially equal to +1 mbar.
[0154] The number (D+1) of apnoea detection periods PDAC, over which the sum of the counted central apnoea detection numbers CAC(i) is performed, is substantially equal to 5.
[0155] The second oscillation frequency range P2 is between substantially 2.5 and 47 Hz.
[0156] The counted central apnoea detection numbers CAC(i) are reset to 0 when the apparatus is started up.
[0157] FIG. 5 also depicts an algorithm for pressure reduction according to the previous appearance or not of events representing sleep problems.
[0158] According to this algorithm, depicted at the bottom of FIG. 5 , the measured pressure P is compared with a predetermined pressure value MPL.
[0159] After determination of the appearance of one or more events,
-
- [0160] if the measured pressure P is less than the predetermined value MPL, a third predetermined pressure reduction control C3 is performed;
- [0161] if the measured pressure P is greater than or equal to the predetermined value MPL, a fourth predetermined pressure reduction control C4 is performed;
- [0162] then, if no event appearance has been detected after one or more of the pressure reduction controls C3; C4, the fourth predetermined pressure reduction control C4 is performed.
[0163] The fourth pressure reduction control C4 is such that it causes a greater pressure reduction per unit of time than that caused by the third control C3.
[0164] In one embodiment, the fourth pressure reduction control C4 is substantially −0.5 mbar/1 minute and the third pressure reduction control C3 is substantially −0.5 mbar/5 minutes, the comparative pressure value MPL is between 4 and 19 mbar and is for example substantially equal to 17 mbar.
[0165] This algorithm for pressure reduction according to the appearance or not of events is implemented after the one for central and obstructive apnoea detection as depicted in FIG. 5 but is also implemented, in non-depicted embodiments, after the other algorithms such as:
-
- [0166] the one for processing decision taking, when the indicator BLN has changed from the first state ON to the second state OFF;
- [0167] the one for hypopnoeic respiration detection described below;
- [0168] the one for acoustic vibration detection described below.
[0169] In FIGS. 6 to 9 , the respiratory cycles are designated as hyperventilated, hypoventilated or normal ventilation cycles and pressure controls are generated according to the designations made.
[0170] At each measured respiratory cycle end, the mean amplitude AM over a fourth predetermined number Y4 of previous respiratory cycles is calculated.
[0171] As depicted in FIG. 7 , if the measured amplitude of the last respiratory cycle is less than the calculated mean amplitude AM multiplied by a first predetermined hypopnoea factor FHO, then the duration TC of the last measured respiratory cycle is added to a hypopnoea time counter CTHO,
-
- [0172] if the current value of the hypopnoea time counter CTHO is greater than or equal to a minimum hypopnoea time TMHO, a fifth predetermined pressure increase is controlled by means of a control C5;
- [0173] after the end of a fifth predetermined number Y5 of respiratory cycles following the fifth pressure increase control C5, a sixth predetermined pressure increase is controlled C6;
- [0174] after the end of a sixth predetermined number Y6 of respiratory cycles, greater than the fifth number Y5, following the fifth pressure increase control C5, a seventh pressure increase is controlled by means of a control C7.
[0175] The hypopnoea time counter CTHO is initialized to 0 when the apparatus is started up.
[0176] In one embodiment, the fourth given number Y4 of respiratory cycles for mean amplitude calculation is substantially equal to 8.
[0177] The first predetermined hypopnoea factor FHO is between 1 and 100% and is for example substantially equal to 40%.
[0178] The minimum hypopnoea time TMHO is between 1 second and 25 seconds and is for example substantially equal to 10 seconds.
[0179] The fifth and sixth predetermined numbers Y5; Y6 of respiratory cycles are substantially equal to respectively 2 and 4.
[0180] The fifth predetermined pressure increase C5 is between 0.1 mbar and 10 mbar and is for example substantially equal to +1 mbar.
[0181] The sixth and seventh predetermined pressure increases C6; C7 are each less than the fifth control C5 and are for example each substantially equal to half the fifth pressure increase C5.
[0182] As depicted in FIGS. 8 and 9 , if the measured amplitude of the last respiratory cycle is greater than or equal to the calculated mean amplitude AM multiplied by the first hypopnoea factor FHO, then the mean respiratory cycle time TCM over a seventh predetermined number Y7 of previous cycles is calculated.
[0183] If the measured duration TC of the last cycle is greater than an eighth predetermined number Y8 multiplied by the calculated mean respiratory cycle time TCM, the measured duration TC of the last cycle, multiplied by a second hypopnoea factor F2, is added to the hypopnoea time counter CTHO.
[0184] If the measured amplitude of the last measured respiratory cycle is greater than a third hyperventilation factor F3, greater than the first hypopnoea factor FHO, multiplied by the calculated mean amplitude AM, the last cycle is designated as hyperventilated, a hyperventilated cycle counter CCH is incremented by one unit, a normal ventilation cycle counter CCN is reset to 0 and
-
- [0185] if the value of the hyperventilated cycle counter CCH is greater than or equal to a ninth predetermined number Y9,
- [0186] if the duration of the last cycle TC is greater than or equal to the eighth number Y8 multiplied by the calculated mean cycle time TCM, the second factor F2 multiplied by the duration of the last respiratory cycle TC is added to the hypopnoea time counter CTHO;
- [0187] and otherwise, the hypopnoea time counter CTHO is reset to 0;
then a hypoventilated cycle counter CCHO is reset to 0 and the mean respiratory cycle amplitude AM over the predetermined number Y4 of previous respiratory cycles is calculated.
[0188] If the measured amplitude of the last measured respiratory cycle is less than or equal to the third factor F3 multiplied by the calculated mean amplitude AM, the last cycle is designated as a normal ventilation cycle, the hyperventilated cycle counter CCH is reset to 0 and the normal ventilation cycle counter CCN is incremented by one unit, and
-
- [0189] if the value of the normal ventilation cycle counter CCN is greater than or equal to a tenth predetermined number Y10,
- [0190] if the duration of the last cycle TC is greater than or equal to the eighth number Y8 multiplied by the calculated mean cycle time TCM, the second factor F2 multiplied by the duration of the last cycle TC is assigned to the hypopnoea time counter CTHO and the normal ventilation cycle counter CCN is reset to 0,
- [0191] and otherwise, the hypopnoea time counter CTHO is reset to 0;
then the hypoventilated cycle counter CCHO is reset to 0 and the mean amplitude of the respiratory cycle over the predetermined number Y4 of respiratory cycles is calculated.
[0192] In one embodiment, the second factor F2 is substantially equal to 5/8.
[0193] The third hyperventilation factor F3 is between 100% and 200% and is for example substantially equal to 140%.
[0194] The seventh, eighth, ninth and tenth predetermined numbers Y7; Y8; Y9; Y10 are respectively substantially equal to 3; 2; 2; and 2.
[0195] In FIG. 10 , it is detected whether the pressure is too high.
[0196] If the measured pressure P is less than a predetermined high pressure value PH, a high pressure time counter TPH is reset to 0.
[0197] If the value of the high pressure time counter TPH is greater than a maximum high pressure time TMPH and
-
- [0198] if the maximum regulated pressure value Pmaxi is less than a predetermined safety pressure value PSEC, the pressure P is controlled to this maximum regulated pressure value Pmaxi;
- [0199] if the minimum regulated pressure value Pmini is greater than a predetermined safety pressure value PSEC, the pressure P is controlled to this minimum regulated pressure value Pmini;
- [0200] if the previous two conditions are not fulfilled, the pressure P is controlled to the safety pressure value PSEC.
then - [0201] the high pressure time counter TPH is reset to 0.
[0202] In one embodiment, the high pressure value PH is between 10 mbar and 25 mbar and is for example substantially equal to 17 mbar.
[0203] The maximum high pressure time TMPH is between 1 and 100 minutes and is for example substantially equal to 10 minutes or 30 minutes.
[0204] The safety pressure value PSEC is substantially equal to 8 mbar.
[0205] In FIG. 11 , an air leakage is measured, substantially equal to the mean flow during respiration of the patient.
[0206] If the measured air leakage is greater than a predetermined leakage level NFM, the pressure increase controls are invalidated.
[0207] In one embodiment, NFM=A×Pfiltered+B.
[0208] According to this formula, the predetermined leakage level NFM is substantially equal to a leakage coefficient A multiplied by a filtered air pressure in the mask, added to an additive leakage coefficient B, the leakage coefficient A being between 0 and 10 litres/minute.mbar and being for example substantially equal to 2.5 litres/minute.mbar.
[0209] The additive leakage coefficient B is between 0 and 100 litres/min and is for example substantially equal to 50 litres/min.
[0210] In FIG. 12 , it is detected whether the measured pressure curve has oscillations, such as acoustic vibrations, within a frequency range P1.
[0211] This detection is performed for example by hardware means such as analogue or digital filters.
[0212] A measurement is made of the detected oscillation presence time RF1 between two successive absences of detected oscillations and the detected oscillation absence time RF0 between two successive presences of detected oscillations.
[0213] If the sum of the measured detected oscillation absence and presence times RF0; RF1 is within a prescribed time range BIP; BSP.
[0214] If the measured oscillation presence time RF1 is greater than or equal to a minimum oscillation time TMRH and if the value of a counter CTAR of elapsed time since the last but one time that the previous time conditions were fulfilled is greater than a prescribed waiting time TAR, an eighth predetermined pressure increase is controlled C8 and the elapsed time counter CTAR is reset to 0.
[0215] The algorithms for acoustic vibration detection and control in the case of acoustic vibrations are implemented at prescribed time intervals, notably regularly and for example every 100 ms.
[0216] At the start of the acoustic vibration detection algorithm depicted in FIG. 12 , if the value of the elapsed time counter CTAR is less than the prescribed waiting time TAR, this counter is incremented (INC CTAR) by the prescribed time interval mentioned above.
[0217] If the sum of the measured detected oscillation presence and absence times RF0; RF1 is below the prescribed time range BIP; BSP or if the measured detected oscillation presence time RF1 is less than the minimum oscillation time TMRH,
-
- [0218] the measured detected oscillation absence time RF0 is replaced by the sum of the measured detected oscillation absence and presence times RF0; RF1, and then
- [0219] the measured detected oscillation presence time RF1 is reset to 0.
[0220] If the sum of the measured detected oscillation absence and presence times RF0; RF1 is above the predetermined time range BIP; BSP or a predetermined maximum time TCMax, each of the measured detected oscillation absence and presence times RF0; RF1 is reset to 0.
[0221] If the two conditions mentioned above concerning the sum of the presence and absence times RF1, RF0 and the presence time RF1 are not fulfilled, each of the measured detected oscillation absence and presence times RF0; RF1 is reset to 0.
[0222] In one embodiment, the predetermined maximum time TCMax is substantially equal to twice the mean respiratory cycle time TCM over the last three measured cycles.
[0223] The prescribed time range BIP; BSP is substantially between 10% and 120% of the calculated mean cycle time TCM.
[0224] The minimum oscillation time TMRH is substantially equal to 7% of the calculated mean cycle time TCM.
[0225] The prescribed waiting time TAR is between 1 and 30 minutes and is for example substantially equal to 1 minute.
[0226] The eighth pressure increase control C8 is between 0.1 mbar and 10 mbar and is for example substantially equal to 1 mbar.
[0227] The oscillation detection frequency range P1 is between substantially 30 and 300 Hz.
[0228] The chronology of the detected events is stored and the stored chronology is read, for example after one night.
[0229] To that end, the central unit U of the apparatus has a memory, not depicted, capable of being written and read with the chronology of the detected events.
[0230] This chronology can be displayed, for example on a monitor, by reading the content of the memory, by means of a computer, not depicted.
Claims
The invention claimed is:
1. A method for managing gas delivery to a patient, comprising:
delivering gas to a patient;
determining a flow amplitude for each of a first plurality of breath cycles;
determining a mean of the determined flow amplitudes for the first plurality of breath cycles;
determining a flow amplitude for a particular breath cycle occurring subsequent to the first plurality of breath cycles;
identifying the particular breath cycle as a hypoventilated breath cycle if the flow amplitude for the particular breath cycle is less than the product of the determined mean flow amplitude for the first plurality of breath cycles and a first factor; and
determining whether to increase the pressure of gas delivered to the patient based at least on whether the particular breath cycle is identified as a hypoventilated breath cycle.
2. A method according to claim 1, further comprising identifying the particular breath cycle as a hyperventilated breath cycle if the flow amplitude for the particular breath cycle is greater than the product of the determined mean flow amplitude for the first plurality of breath cycles and a second factor greater than the first factor.
3. A method according to claim 1, wherein the first factor is substantially equal to 40%.
4. A method according to claim 1, comprising:
determining a duration of consecutive breath cycles identified as hypoventilated breath cycles; and
determining to increase the pressure of gas delivered to the patient if the duration of consecutive breath cycles identified as hypoventilated breath cycles exceeds a predetermined threshold duration.
5. A method according to claim 1, further comprising:
determining a mean breath cycle duration of a second plurality of breath cycles by the patient; and
determining whether the duration of the particular breath cycle is greater than the product of the mean breath cycle duration of the second plurality of breath cycles and a second factor.
6. A method according to claim 5, further comprising:
increasing a hypoventilated breath cycle time counter if the duration of the particular breath cycle is greater than the product of the mean breath cycle duration of the second plurality of breath cycles and the second factor; and
using the hypoventilated breath cycle time counter for determining whether to increase the pressure of gas delivered to the patient based on one or more hypoventilated breath cycles subsequent to the particular breath cycle.
7. A method according to claim 5, wherein the second plurality of breath cycles at least partially overlaps with the first plurality of breath cycles.
8. A method according to claim 1, further comprising:
identifying the particular breath cycle as a hyperventilated breath cycle if the flow amplitude for the particular breath cycle is greater than the product of the determined mean flow amplitude for the first plurality of breath cycles and a second factor greater than the first factor; and
identifying the particular breath cycle as a normal breath cycle if the flow amplitude for the particular breath cycle is neither (a) less than the product of the determined mean flow amplitude for the first plurality of breath cycles and the first factor, nor (b) greater than the product of the determined mean flow amplitude for the first plurality of breath cycles and the second factor.
9. A method according to claim 1, wherein each determined flow amplitude comprises a peak inspiratory flow amplitude.
10. A breathing assistance apparatus, comprising:
a pressure supply module configured to supply pressurized air toward a patient;
at least one sensor configured to monitor air supplied by the pressure supply module; and
a pressure control unit coupled to the pressure supply module, the pressure control unit configured to:
determine a flow amplitude for each of a first plurality of breath cycles;
determine a mean of the determined flow amplitudes for the first plurality of breath cycles;
determine a flow amplitude for a particular breath cycle occurring subsequent to the first plurality of breath cycles;
identify the particular breath cycle as a hypoventilated breath cycle if the flow amplitude for the particular breath cycle is less than the product of the determined mean flow amplitude for the first plurality of breath cycles and a first factor; and
determine whether to increase the pressure of gas delivered to the patient based at least on whether the particular breath cycle is identified as a hypoventilated breath cycle.
11. An apparatus according to claim 10, wherein the pressure control unit is further configured to identify the particular breath cycle as a hyperventilated breath cycle if the flow amplitude for the particular breath cycle is greater than the product of the determined mean flow amplitude for the first plurality of breath cycles and a second factor greater than the first factor.
12. An apparatus according to claim 10, wherein the pressure control unit is configured to:
determine a duration of consecutive breath cycles identified as hypoventilated breath cycles; and
determine to increase the pressure of gas delivered to the patient if the duration of consecutive breath cycles identified as hypoventilated breath cycles exceeds a predetermined threshold duration.
13. An apparatus according to claim 10, wherein the pressure control unit is further configured to:
determine a mean breath cycle duration of a second plurality of breath cycles by the patient; and
determine whether the duration of the particular breath cycle is greater than the product of the mean breath cycle duration of the second plurality of breath cycles and a second factor.
14. An apparatus according to claim 13, wherein the pressure control unit is further configured to:
increase a hypoventilated breath cycle time counter if the duration of the particular breath cycle is greater than the product of the mean breath cycle duration of the second plurality of breath cycles and the second factor; and
use the hypoventilated breath cycle time counter for determining whether to increase the pressure of gas delivered to the patient based on one or more hypoventilated breath cycles subsequent to the particular breath cycle.
15. An apparatus according to claim 13, wherein the second plurality of breath cycles at least partially overlaps with the first plurality of breath cycles.
16. An apparatus according to claim 10, wherein the pressure control unit is further configured to:
identify the particular breath cycle as a hyperventilated breath cycle if the flow amplitude for the particular breath cycle is greater than the product of the determined mean flow amplitude for the first plurality of breath cycles and a second factor greater than the first factor; and
identify the particular breath cycle as a normal breath cycle if the flow amplitude for the particular breath cycle is neither (a) less than the product of the determined mean flow amplitude for the first plurality of breath cycles and the first factor, nor (b) greater than the product of the determined mean flow amplitude for the first plurality of breath cycles and the second factor.
17. An apparatus according to claim 10, wherein each determined flow amplitude comprises a peak inspiratory flow amplitude.
18. Software for managing gas delivery to a patient, the software embodied in tangible computer-readable media and when executed by a processor, configured to:
determine a flow amplitude for each of a first plurality of breath cycles;
determine a mean of the determined flow amplitudes for the first plurality of breath cycles;
determine a flow amplitude for a particular breath cycle occurring subsequent to the first plurality of breath cycles;
identify the particular breath cycle as a hypoventilated breath cycle if the flow amplitude for the particular breath cycle is less than the product of the determined mean flow amplitude for the first plurality of breath cycles and a first factor; and
determine whether to increase the pressure of gas delivered to the patient based at least on whether the particular breath cycle is identified as a hypoventilated breath cycle.
Patent Citations (118)
| Patent | Date | Inventor | Cited By |
|---|---|---|---|
| US3675640(A) | 1972-07-01 | Gatts | Applicant |
| US3991304(A) | 1976-11-01 | Hillsman | Applicant |
| US4077404(A) | 1978-03-01 | Elam | Applicant |
| US4365636(A) | 1982-12-01 | Barker | Applicant |
| US4440177(A) | 1984-04-01 | Anderson et al. | Applicant |
| US4444201(A) | 1984-04-01 | Itoh | Applicant |
| US4448192(A) | 1984-05-01 | Stawitcke et al. | Applicant |
| US4602644(A) | 1986-07-01 | DiBenedetto et al. | Applicant |
| US4655213(A) | 1987-04-01 | Rapoport et al. | Applicant |
| US4686974(A) | 1987-08-01 | Sato et al. | Applicant |
| US4723543(A) | 1988-02-01 | Beran | Applicant |
| US4782832(A) | 1988-11-01 | Trimble et al. | Applicant |
| US4803997(A) | 1989-02-01 | Bowman | Applicant |
| US4811755(A) | 1989-03-01 | Bourdon et al. | Applicant |
| US4944310(A) | 1990-07-01 | Sullivan | Applicant |
| US4957107(A) | 1990-09-01 | Sipin | Applicant |
| US5065756(A) | 1991-11-01 | Rapoport | Applicant |
| US5107831(A) | 1992-04-01 | Halpern et al. | Applicant |
| US5134995(A) | 1992-08-01 | Gruenke et al. | Applicant |
| US5146918(A) | 1992-09-01 | Kallok et al. | Applicant |
| US5148802(A) | 1992-09-01 | Sanders et al. | Applicant |
| US5161541(A) | 1992-11-01 | Bowman et al. | Applicant |
| US5190048(A) | 1993-03-01 | Wilkinson | Applicant |
| US5199424(A) | 1993-04-01 | Sullivan et al. | Applicant |
| US5203343(A) | 1993-04-01 | Axe et al. | Applicant |
| US5231979(A) | 1993-08-01 | Rose et al. | Applicant |
| US5239995(A) | 1993-08-01 | Estes et al. | Applicant |
| US5245995(A) | 1993-09-01 | Sullivan et al. | Applicant |
| US5246995(A) | 1993-09-01 | Murakami et al. | Applicant |
| US5259373(A) | 1993-11-01 | Gruenke et al. | Applicant |
| US5269296(A) | 1993-12-01 | Landis | Applicant |
| US5275159(A) | 1994-01-01 | Griebel | Applicant |
| US5313937(A) | 1994-05-01 | Zdrojkowski | Applicant |
| US5335654(A) | 1994-08-01 | Rapoport | Applicant |
| US5353788(A) | 1994-10-01 | Miles | Applicant |
| US5365922(A) | 1994-11-01 | Raemer | Applicant |
| US5413111(A) | 1995-05-01 | Wilkinson | Applicant |
| US5433193(A) | 1995-07-01 | Sanders et al. | Applicant |
| US5438980(A) | 1995-08-01 | Phillips | Applicant |
| US5456264(A) | 1995-10-01 | Series et al. | Applicant |
| US5458137(A) | 1995-10-01 | Axe et al. | Applicant |
| US5490502(A) | 1996-02-01 | Rapoport et al. | Applicant |
| US5522382(A) | 1996-06-01 | Sullivan et al. | Applicant |
| US5535739(A) | 1996-07-01 | Rapoport et al. | Applicant |
| US5540733(A) | 1996-07-01 | Testerman et al. | Applicant |
| US5546933(A) | 1996-08-01 | Rapoport et al. | Applicant |
| US5549106(A) | 1996-08-01 | Gruenke et al. | Applicant |
| USRE35339(E) | 1996-10-01 | Rapoport | Applicant |
| US5598838(A) | 1997-02-01 | Servidio et al. | Applicant |
| US5630411(A) | 1997-05-01 | Holscher | Applicant |
| US5645053(A) | 1997-07-01 | Remmers et al. | Applicant |
| US5664562(A) | 1997-09-01 | Bourdon | Applicant |
| US5704345(A) | 1998-01-01 | Berthon-Jones | Applicant |
| US5732696(A) | 1998-03-01 | Rapoport et al. | Applicant |
| US5746201(A) | 1998-05-01 | Kidd | Applicant |
| US5746697(A) | 1998-05-01 | Swedlow et al. | Applicant |
| US5752509(A) | 1998-05-01 | Lachmann et al. | Examiner |
| US5794614(A) | 1998-08-01 | Gruenke et al. | Applicant |
| US5803066(A) | 1998-09-01 | Rapoport et al. | Applicant |
| US5823187(A) | 1998-10-01 | Estes et al. | Applicant |
| US5829441(A) | 1998-11-01 | Kidd et al. | Applicant |
| US5845636(A) | 1998-12-01 | Gruenke et al. | Applicant |
| US5921238(A) | 1999-07-01 | Bourdon | Applicant |
| US5924979(A) | 1999-07-01 | Swedlow et al. | Applicant |
| US5931160(A) | 1999-08-01 | Gilmore et al. | Applicant |
| US6032065(A) | 2000-02-01 | Brown | Applicant |
| US6135105(A) | 2000-10-01 | Lampotang et al. | Applicant |
| US6135106(A) | 2000-10-01 | Dirks et al. | Applicant |
| US6158432(A) | 2000-12-01 | Biondi et al. | Applicant |
| US6283119(B1) | 2001-09-01 | Bourdon | Applicant |
| US6299581(B1) | 2001-10-01 | Rapoport et al. | Applicant |
| US6341606(B1) | 2002-01-01 | Bordewick et al. | Applicant |
| US6347631(B1) | 2002-02-01 | Hansen et al. | Applicant |
| US6398739(B1) | 2002-06-01 | Sullivan et al. | Applicant |
| US6418928(B1) | 2002-07-01 | Bordewick et al. | Applicant |
| US6431172(B1) | 2002-08-01 | Bordewick | Applicant |
| US6463930(B2) | 2002-10-01 | Biondi et al. | Applicant |
| US6488634(B1) | 2002-12-01 | Rapoport et al. | Applicant |
| US6505623(B1) | 2003-01-01 | Hansen | Applicant |
| US6516802(B2) | 2003-02-01 | Hansen et al. | Applicant |
| US6530373(B1) | 2003-03-01 | Patron et al. | Applicant |
| US6546930(B1) | 2003-04-01 | Emerson et al. | Applicant |
| US6571795(B2) | 2003-06-01 | Bourdon | Applicant |
| US6584973(B1) | 2003-07-01 | Biondi et al. | Applicant |
| US6644310(B1) | 2003-11-01 | Delache et al. | Applicant |
| US6644316(B2) | 2003-11-01 | Bowman et al. | Applicant |
| US6668829(B2) | 2003-12-01 | Biondi et al. | Applicant |
| US6668830(B1) | 2003-12-01 | Hansen et al. | Applicant |
| US6718974(B1) | 2004-04-01 | Moberg | Applicant |
| US6761167(B1) | 2004-07-01 | Nadjafizadeh et al. | Applicant |
| US6761168(B1) | 2004-07-01 | Nadjafizadeh et al. | Applicant |
| US6793629(B2) | 2004-09-01 | Rapoport et al. | Applicant |
| US6814074(B1) | 2004-11-01 | Nadjafizadeh et al. | Applicant |
| US6866040(B1) | 2005-03-01 | Bourdon | Applicant |
| US6960854(B2) | 2005-11-01 | Nadjafizadeh et al. | Applicant |
| DE3306607(A1) | 1983-02-01 | Applicant | |
| DE4038871(A1) | 1992-06-01 | Applicant | |
| EP46570 | 1982-03-01 | Applicant | |
| EP504945 | 1992-09-01 | Applicant | |
| EP504725(A2) | 1992-09-01 | Applicant | |
| EP651971(A1) | 1995-05-01 | Applicant | |
| EP661071(A1) | 1995-07-01 | Applicant | |
| FR2789593(A1) | 2000-08-01 | Applicant | |
| GB2294400 | 1996-05-01 | Applicant | |
| WO88/10108 | 1988-12-01 | Applicant | |
| WO89/10768 | 1989-11-01 | Applicant | |
| WO89/11823 | 1989-12-01 | Applicant | |
| WO90/14121 | 1990-11-01 | Applicant | |
| WO91/06832 | 1991-05-01 | Applicant | |
| WO92/11054(A1) | 1992-07-01 | Applicant | |
| WO92/22244(A1) | 1992-12-01 | Applicant | |
| WO94/06499(A1) | 1994-03-01 | Applicant | |
| WO94/23780(A1) | 1994-10-01 | Applicant | |
| WO95/28874 | 1995-11-01 | Applicant | |
| WO96/08285 | 1996-03-01 | Applicant | |
| WO97/14462(A1) | 1997-04-01 | Applicant | |
| WO97/28838(A1) | 1997-08-01 | Applicant | |
| WO/47262 | 2000-02-01 | Applicant |
Non-Patent Literature (15)
- PCT International Search Report/PCT/FR00/00334, 4 pages, Apr. 20, 2000.Applicant
- Non-Final Office Action for U.S. Appl. No. 10/968,403, 8 pages, Aug. 31, 2007.Applicant
- Response to Non-Final Office Action dated Aug. 31, 2007 for U.S. Appl. No. 10/968,403, 16 pages, Nov. 30, 2007.Applicant
- Dupuis, Yvon G., Ventilators—Theory and Clinical Application, The C.V. Mosby Company, pp. 107-117, 1986.Applicant
- Schwartz, Alan R. et al., “Induction of upper airway occlusion in sleeping individuals with subatmospheric nasal pressure”, The American Physiological Society, pp. 535-542, 1988.Applicant
- Sullivan, Colin E., et al., “Reversal of Obstructive Sleep Apnoea by Continuous Positive Airway Pressure Applied through the Nares”, The Lancet, pp. 862-865, Apr. 18, 1981.Applicant
- Rapoport, David M., “Techniques for Administering Nasal CPAP”, Respiratory Management, pp. 17-21, Jul.-Aug. 1987.Applicant
- Garay, Stuart M., “Therapeutic Options for Obstructive Sleep Apnea”, Respiratory Management, pp. 11-15, Jul.-Aug. 1987.Applicant
- Southworth et al., “Digital Computation and Numerical Methods”, McGraw-Hill, Inc., pp. 6-10, 1965.Applicant
- “Sleep Apnea—Diagnosis and Treatment”, Topics in Primary Care Medicine, The Western Journal of Medicine, 145:249-250, Aug. 1986.Applicant
- Guilleminault, C. et al.; “Unattended CPAP Titration: Toward a Smart Machine”; 1992; pp. 342 [Sleep Research, vol. 21].Applicant
- Miles, L.E. et al.; “Development and Application of an Automatic Nasal CPAP Calibration Procedures for Use in the Unsupervised Home Environment”; pp. S118-S119 [Supplemental Issue: Sleep and Breathing, vol. 16.], Dec. 1993.Applicant
- Miles, L.E. et al.; “Different Roles for an Automatic Nasal CPAP Calibration Procedure and Smart Pap”, Sleep Research 21, p. 238, 1993.Applicant
- Miles, L.E. et al.; “Development and Application of an Automatic Nasal CPAP Calibration Procedure for Use in the Unsupervised Home Environment”; p. 352 [Sleep Research, vol. 21], 1992.Applicant
- Remmers, et al.; “Mechanics of the Pharynx in Patients With Obstructive Sleep Apnea”; pp. 261-271 [Sleep and Respiration], 1990.Applicant