Method and apparatus for predicting work of breathing

FIELD OF THE INVENTION

The present invention relates generally to the field of respiratory therapy and physiology, including ventilator and respiratory monitor technology, and, more particularly, to a method and apparatus for predicting a patient's physiologic work of breathing and imposed work of breathing.

BACKGROUND

Mechanical ventilatory support is widely accepted as an effective form of therapy and means for treating patients with respiratory failure. Ventilation is the process of delivering oxygen to and washing carbon dioxide from the alveoli in the lungs. When receiving ventilatory support, the patient becomes part of a complex interactive system that is expected to provide adequate ventilation and promote gas exchange to aid in the stabilization and recovery of the patient. Clinical treatment of a ventilated patient often calls for monitoring a patient's breathing to detect an interruption or an irregularity in the breathing pattern, for triggering a ventilator to initiate assisted breathing, and for interrupting the assisted breathing periodically to wean the patient off of the assisted breathing regime, thereby restoring the patient's ability to breathe independently.

In those instances in which a patient requires mechanical ventilation due to respiratory failure, a wide variety of mechanical ventilators are available. Most modern ventilators allow the clinician to select and use several modes of inhalation either individually or in combination via the ventilator setting controls that are common to the ventilators. These modes can be defined in three broad categories: spontaneous, assisted or controlled. During spontaneous ventilation without other modes of ventilation, the patient breathes at his own pace, but other interventions may affect other parameters of ventilation including the tidal volume and the baseline pressure, above ambient, within the system. In assisted ventilation, the patient initiates the inhalation by lowering the baseline pressure by varying degrees, and then the ventilator "assists" the patient by completing the breath by the application of positive pressure. During controlled ventilation, the patient is unable to breathe spontaneously or initiate a breath, and is therefore dependent on the ventilator for every breath. During spontaneous or assisted ventilation, the patient is required to "work" (to varying degrees) by using the respiratory muscles in order to breathe.

The total work of breathing (the work to initiate and sustain a breath) performed by a patient to inhale while intubated and attached to the ventilator may be divided into two major components: physiologic work of breathing (the work of breathing of the patient, "WOBp") and breathing apparatus (endotracheal tube and ventilator) imposed resistive work of breathing or imposed work of breathing ("WOBi"). The total work of breathing ("WOB") can be measured and quantified in Joules/L of ventilation. In the past, techniques have been devised to supply ventilatory therapy to patients for the purpose of improving patient's efforts to breathe by decreasing the total work of breathing to sustain the breath. Still other techniques have been developed that aid in the reduction of the patient's inspiratory work required to trigger a ventilator system "ON" to assist the patient's breathing. It is desirable to reduce the effort expended by the patient in each of these phases, because a high total work of breathing load can cause further damage to a weakened patient or be beyond the capacity or capability of small or disabled patients.

Furthermore, it is desirable to quantify the imposed work of breathing (WOBi) of a patient since this value is critical in the decision of when to extubate or remove ventilatory support from the patient. High work loads tend to create breathing patterns that are rapid and shallow (high frequency and low tidal volume). Rapid shallow breathing is typically a contra-indicator of extubation success. If this high work load, however, is caused mainly by a large imposed work of breathing such that most of the excess work is caused by the breathing apparatus, extubation success is much higher ("<nplcit id="ncit0001" npl-type="s"><text>Elevated imposed work of breathing masquerading as ventilator weaning intolerance." Chest. 1995 Oct.; 108(4): 1021-5</text></nplcit>).

The early generation of mechanical ventilators, prior to the mid-1960s, were designed to support alveolar ventilation and to provide supplemental oxygen for those patients who were unable to breathe due to neuromuscular impairment. Since that time, mechanical ventilators have become more sophisticated and complicated in response to increasing understanding of lung pathophysiology. In an effort to improve a patient's tolerance of mechanical ventilation, assisted or patient-triggered ventilation modes were developed. Partial positive pressure ventilation (PPV) support, in which mechanical support supplements spontaneous ventilation, became possible for adults outside the operating room when intermittent mandatory ventilation (IMV) became available in the 1970s. Varieties of "alternative" ventilation modes addressing the needs of severely impaired patients continue to be developed.

In recent years, microprocessors have been introduced into modem ventilators. Microprocessor ventilators are typically equipped with sensors that monitor breath-by-breath flow, pressure, volume, and derive mechanical respiratory parameters. Their ability to sense and transduce "accurately," combined with computer technology, makes the interaction between clinician, patient, and ventilator more sophisticated than ever. The prior art microprocessor controlled ventilators suffered from compromised accuracy due to the placement of the sensors required to transduce the data signals. Consequently, complicated algorithms were developed so that the ventilators could "approximate" what was actually occurring within the patient's lungs on a breath-by-breath basis. In effect, the computer controlled prior art ventilators were limited to the precise, and unyielding, nature of the mathematical algorithms that attempted to mimic cause-and-effect in the ventilator support provided to the patient.

The overall performance of the assisted ventilatory system is determined by both physiological and mechanical factors. The physiological determinants, which include the nature of the pulmonary disease, the ventilatory efforts of the patient, and many other physiological variables, changes with time and are difficult to diagnose. Moreover, the physician historically had relatively little control over these determinants. Mechanical input to the system, on the other hand, is to a large extent controlled and can be reasonably well characterized by examining the parameters of ventilator flow, volume, and/or pressure. Optimal ventilatory assistance requires both appropriately minimizing physiologic workloads to a tolerable level and decreasing imposed resistive workloads to zero. Doing both should ensure that the patient is neither overstressed nor oversupported. Insufficient ventilatory support places unnecessary demands upon the patient's already compromised respiratory system, thereby inducing or increasing respiratory muscle fatigue. Excessive ventilatory support places the patient at risk for pulmonary-barotrauma, respiratory muscle deconditioning, and other complications of mechanical ventilation.

In addition to total work of breathing (WOB), there are other measurements of patient effort, including power of breathing (POB), the rate at which total work of breathing is done, and the pressure time product (PTP), the integrated product of time multiplied by the decrease in pleural pressure during a breath. These methodologies are similar in their goal of measuring patient effort, but are calculated differently and provide different measures of the patient effort.

Although total work of breathing (and its alternatives) has been considered an important parameter for appropriately setting a ventilator, it has remained largely unused because of the difficulty in obtaining its value. Physiologic work of breathing is defined using a pleural pressure versus volume graph of a patient's breath. The pleura is a two-layered membrane that envelops the lung and contains lubricating fluid between its inner and outer layers. During breathing, the respiratory muscles either compress or expand the lungs by exerting forces on the pleura. The pressure in the pleural space therefore represents the respiratory effort. The patient's physiologic work of breathing is the area from the chest wall compliance line on the right to the pleural pressure versus volume loop on the left (see <figref idref="f0001">FIG. 1</figref>). Since the pleural pressure is very difficult to obtain and may be different at different positions in the pleural space, a typical surrogate for pleural pressure is esophageal pressure. The esophageal pressure is typically obtained by placing a balloon in the esophagus between the heart and the stomach.

Likewise, although imposed work of breathing has been considered an important parameter for appropriately setting a ventilator, it has also remained largely unused because of the difficulty in obtaining its value. Imposed work of breathing is the area below baseline pressure circumscribed within the tracheal pressure-tidal volume loop during spontaneous inhalation. Typically, this is done by using a catheter inserted into the tracheal tube or a lumen in the side of the tracheal tube that opens at the distal end of the tracheal tube. These devices are then attached to a pressure transducer to measure tracheal pressure. The greatest single difficulty with these devices is the harsh environment in which they exist and their propensity for becoming clogged. For this and other reasons, tracheal pressure is difficult to reliably measure and thus is not normally used.

<patcit id="pcit0001" dnum="US5316009A"><text>U.S. Patent 5,316,009</text></patcit> describes an apparatus for monitoring respiratory muscle activity based on measuring resistance and elastance of the lung and then calculating a value called Pmus from the standard equation of motion where Paw = Pmus + R*flow + V/C. It also discloses calculation of a PTP of pmus, which is not the standard PTP, and a "work" Wmus, but not necessarily real WOB. A problem with the method taught by the '009 patent is that Pmus is difficult to measure in a spontaneously breathing patient because the parameters R and C must be very accurately computed in order for Pmus to correlate with "work". Moreover, R and C in a spontaneously breathing patient with ventilator support are very difficult to obtain accurately.

Occlusion pressure at 0.1 seconds after breath initiation by the patient (P0.1) has also been proposed as an indicator of work of breathing. P0.1 can be based on esophageal pressure or airway pressure. An esophageal pressure P0.1 is invasive but correlates fairly well with work of breathing. An airway pressure P0.1 is non-invasive, but does not correlate nearly as well with work of breathing.

A number of other patents exist for respiratory systems including U. S. Patent Nos. <patcit id="pcit0002" dnum="US6439229B"><text>6,439,229</text></patcit>; <patcit id="pcit0003" dnum="US6390091B"><text>6,390,091</text></patcit>;<patcit id="pcit0004" dnum="US6257234B"><text> 6,257,234</text></patcit>; <patcit id="pcit0005" dnum="US6068602A"><text>6,068,602</text></patcit>; <patcit id="pcit0006" dnum="US6027498A"><text>6,027,498</text></patcit>,<patcit id="pcit0007" dnum="US6019732A"><text> 6,019,732</text></patcit>; <patcit id="pcit0008" dnum="US5941841A"><text>5,941,841</text></patcit>; <patcit id="pcit0009" dnum="US5887611A"><text>5,887,611</text></patcit>; <patcit id="pcit0010" dnum="US5876352A"><text>5,876,352</text></patcit>; <patcit id="pcit0011" dnum="US5807245A"><text>5,807,245</text></patcit>; and <patcit id="pcit0012" dnum="US5682881A"><text>5,682,881</text></patcit>. Also relevant are the article "<nplcit id="ncit0002" npl-type="s"><text>Breathing frequency and pattern are poor predictors of work of breathing in patients receiving pressure support ventilation" Chest. 1995 Jan.; 108(5): 1338-1344</text></nplcit> and the patent application <patcit id="pcit0013" dnum="WO0100265A1"><text>WO 01/00265 A1</text></patcit>.

Accordingly, there is a need in the art for a system and method to noninvasively and accurately predict physiologic work of breathing and imposed work of breathing in a patient. The present invention is designed to address this need.

SUMMARY OF THE INVENTION

The invention is defined in the claims. Broadly speaking, the present invention provides a method and apparatus for non-invasively predicting (estimating) physiologic work of breathing (the amount of effort expended by the patient to breathe) and imposed work of breathing (the amount of effort expended by the patient imposed by the breathing apparatus). This effort, typically invasively calculated as (imposed or physiologic) work of breathing (WOB), power of breathing (POB), or pressure time product (PTP), is useful in determining the most appropriate settings on a ventilator used to support the patient's breathing. Measuring patient effort allows for appropriate ventilatory support that avoids respiratory muscle fatigue and respiratory muscle deconditioning. Measuring imposed patient effort allows for more appropriate ventilatory support by allowing for the imposed effort to be driven to zero to simulate more natural breathing and also as an important extubation criteria.

In one aspect, the method comprises creating a mathematical model of the patient's inspiratory effort using predetermined parameters that are collected non-invasively, such as those collected with standard respiratory monitors. The respiratory monitors typically contain airway pressure and airway flow sensors that measure the flow going into and out of the patient, and often times a carbon dioxide sensor and pulse oximeter. From these time-waveforms, a variety of parameters are selectively derived that are used in characterizing different aspects of the patient's breathing and/or the patient's interaction with the ventilator. These parameters contain information that is extracted to accurately estimate the patient effort.

More specifically, the method comprises a method of estimating the actual patient effort parameter using a combination of multiple parameters derived from sensors that monitor the patient and/or ventilator. The patient effort parameter can be any parameter that represents the effort exerted by the patient to breathe, including but not limited to work of breathing, power of breathing, or pressure time product.

In this method, the parameters are preferably derived from the airway pressure, flow, and volume waveforms and the carbon dioxide and pulse oximeter waveforms normally collected by a respiratory monitor, including but not limited to tidal volume, breathing frequency, peak inspiratory pressure (PIP), inspiratory time, P0.1, trigger time, trigger depth, respiratory system resistance, respiratory compliance, end-tidal carbon dioxide, variations in the pulse oximeter plethysmogram, and the concavity/convexity of the pressure waveform.

This method includes using a linear combination of parameters or a nonlinear combination of parameters, including but not limited to a neural network, fuzzy logic, mixture of experts, or polynomial model. Moreover, multiple different models can be used to estimate the patient effort of different subsets of patients. These subsets can be determined by various means, including but not limited to patient condition (pathophysiology), patient physiologic parameters (lung resistance and compliance), or other parameters.

In an aspect of the invention, the method for estimating work of breathing in a patient comprises use of a neural network, wherein the neural network provides work of breathing information for the patient based upon input data, wherein the input data includes at least one of the following parameters: the airway pressure, flow, airway volume, carbon dioxide flow, and pulse oximeter plethysmogram waveforms normally collected by a respiratory monitor, including but not limited to tidal volume, breathing frequency, peak inspiratory pressure (PIP), inspiratory time, P01 (see reference), trigger time, trigger depth, respiratory system resistance, respiratory compliance, end-tidal carbon dioxide, variations in the pulse oximeter plethysmogram, and the concavity/convexity of the pressure waveform, wherein the work of breathing information is provided as an output variable.

In the above-noted method, the neural network is trained by clinical testing of a test population of patients to obtain teaching data, the teaching data which includes the above-noted input information. The teaching data is provided to the neural network, whereby the neural network is trained to provide an output variable corresponding to the work of breathing. Teaching data further includes esophageal pressure and/or tracheal pressure.

As a system for estimating work of breathing in a patient, the system comprises a neural network which first receives as input primary teaching data obtained from clinical testing of a test population of patients, whereby the neural network learns the teaching data and is trained to provide an output variable for work of breathing, such that when the neural network receives patient input data in the form of the above-noted parameters obtained from a patient, the neural network provides the output variable for estimating work of breathing for that patient.

The invention, as defined in the claims, can be implemented in numerous ways, including as a system (including a computer processing or database system), a method (including a computerized method of collecting and processing input data and a method for evaluating such data to provide an output(s)), an apparatus, a computer readable medium, a computer program product, or a data structure tangibly fixed in a computer readable memory. Several embodiments of the invention are discussed below.

As a system, an embodiment of the invention includes a processor unit having input and output devices. The processor unit operates to receive input parameters, process the input and provide an output corresponding to work of breathing. This output can be then used to control external devices, such as a ventilator. The processing of the data can be accomplished by various means such as neural networks, parallel distributed processing systems, neuromorphic systems, or the like.

As a method of predicting work of breathing, the method includes processing predetermined input variables (parameters), preferably through the use of a neural network.

As a computer readable media containing program instructions, an embodiment of the invention includes: computer readable code devices for receiving input variables, processing the input, and providing an output indicative of work of breathing. In a preferred embodiment, processing comprises utilizing a neural network. The method may further include controlling a ventilator in response to the output obtained.

The methods of the present invention may be implemented as a computer program product with a computer-readable medium having code thereon. The program product includes a program and a signal bearing media bearing the program.

As an apparatus, the present invention may include at least one processor, a memory coupled to the processor, and a program residing in the memory which implements the methods of the present invention.

Other aspects and advantages of the invention will become apparent from the following detailed description taken in conjunction with the accompanying drawings, illustrating, by way of example, the principles of the invention.

All patents, patent applications, provisional applications, and publications referred to or cited herein, or from which a claim for benefit of priority has been made, are referred to in their entirety to the extent they are not inconsistent with the explicit teachings of this specification.

BRIEF DESCRIPTION OF THE DRAWINGS

In order that the manner in which the above-recited and other advantages and objects of the invention are obtained, a more particular description of the invention briefly described above will be rendered by reference to specific embodiments thereof which are illustrated in the appended drawings. Understanding that these drawings depict only typical embodiments of the invention and are not therefore to be considered to be limiting of its scope, the invention will be described and explained with additional specificity and detail through the use of the accompanying drawings in which: <ul id="ul0001" list-style="none" compact="compact"><li><figref idref="f0001">FIG. 1</figref> is a graphical representation of pleural pressure showing elastic work, stored energy of chest wall and resistive work;</li><li><figref idref="f0002">FIG. 2</figref> depicts the method of one aspect of the invention for a patient on a ventilator;</li><li><figref idref="f0003">FIG. 3</figref> depicts a graph of the nonlinear prediction of Power of Breathing;</li><li><figref idref="f0004">FIG. 4</figref> depicts a neural network showing hidden layers; and</li><li><figref idref="f0005">FIG. 5</figref> depicts inputs and outputs of an adaptive system having backpropagation.</li></ul>

DESCRIPTION OF THE PREFERRED EMBODIMENT

Referring now to the drawings, the preferred embodiment of the present invention will be described.

In the embodiment depicted in <figref idref="f0002">FIG. 2</figref>, a patient <b>10</b> requiring respiratory support and connected to a ventilator <b>12</b> will have an airway flow and pressure sensor <b>14,</b> along with possibly a carbon dioxide detector attached at the y-piece of the standard ventilator circuit <b>16.</b> These sensors measure the flow, pressure, and partial pressure of carbon dioxide in the gases that pass to and from the patient. These raw signals <b>18</b> are preprocessed in a signal processor <b>20</b> using analog and digital signal processing to clean the signal, remove sensor biases and offsets, etc. These signals are then processed in a parameter extraction module <b>22</b> to calculate a variety of other parameters from the flow, pressure, and CO₂ data. For example, tidal volume is computed by integrating the flow into the patient over an inspiratory cycle; peak inspiratory pressure is calculated by determining the maximum pressure during a breath; P0.1 is calculated by measuring the change in airway pressure during the first tenth of a second of a breath; respiratory system resistance and compliance can be calculated by adaptively fitting a model, such as a linear model, to the airway pressure, flow, and volume signals; etc. In accordance with the invention, a neural network <b>24</b> is provided to model the parameters so that a ventilator may be controlled through controller <b>26.</b>

In order to create the model (training phase) patient information may be collected from one or more patients. The methodology of obtaining the test data (both input data and desired output data are used to create a model that will predict patient effort) is as follows:

An inflatable balloon is placed into the patient's esophagus and a catheter is placed in the patient's endotracheal tube. The balloon is partially inflated and creates a closed circuit that transmits esophageal pressure to a pressure sensor connected to the end of a catheter connected to the balloon. The pressure sensors read the tracheal pressure and the esophageal pressure, which is a surrogate for pleural pressure (<figref idref="f0001">FIG. 1</figref>). The esophageal pressure drop during normal breathing is plotted on a pressure-volume plot and a loop is created and integrated with a chest wall compliance line to calculate work of breathing. The esophageal pressure maybe first cleaned of cardiac interference with an adaptive noise cancellation technique that subtracts a scaled and shifted version of the ECG or plethysmogram from the esophageal pressure. The chest wall compliance can be calculated by paralyzing the patient and measuring the esophageal pressure rise in the circuit with a positive pressure breath. Many studies have found that 0.1 L/cm H20 is an appropriate average. This plot is called the Campbell diagram and is a method for calculating patient inspiratory or physiologic work of breathing from the measured esophageal pressure. The imposed work of breathing may be calculated by integrating the tracheal pressure drop versus volume waveform similar to the total work of breathing (with the exception of the chest wall compliance). For both total and imposed efforts, power of breathing is calculated as a per minute average of work of breathing. Pressure Time Product (PTP) is measured just like WOB except using a Pressure - Time axis integration instead of Pressure - Volume. The work of breathing, (both physiologic and imposed), the power of breathing, and the pressure time product may be used as output variable for the system. These are the three main "desired outputs" of the preferred system.

Some of the inputs are standard parameters that are calculated in a very straight forward manner, such as spontaneous breathing frequency (the number of breaths the patient takes per minute), tidal volume (the volume of air inspired per breath), etc. Other input parameters are more complex, such as respiratory system resistance, compliance, and airway pressure waveform characteristics. Resistance and compliance for example, are typically calculated using a least squares modeling of the pressure, volume, and flow waveforms (obtained by a standard respiratory monitor using airway pressure and flow sensors) using the equation: <maths id="math0001" num=""><math display="block"><mrow><mi mathvariant="italic">airway_pressure</mi><mo>=</mo><mi mathvariant="italic">flow</mi><mo>×</mo><mi mathvariant="italic">resistance</mi><mo>+</mo><mi mathvariant="italic">volume</mi><mo>×</mo><mi mathvariant="italic">compliance</mi><mo>+</mo><mi mathvariant="italic">positive_end_expiratory_pressure</mi></mrow></math><img file="EP1534131B1_D0001.tif" /></maths>

Using the flow, volume, and pressure obtained from the flow and pressure sensors, the only unknowns are compliance and resistance and these can be adaptively computed using techniques such as least squares optimization, as noted above.

The problem with these methods is that the equation is only valid if the patient is not exerting any effort and the ventilator is responsible for all work of breathing. It has been experimentally determined that better results may be obtained by using an end-inspiratory pause and analyzing the pressure drop with the cessation of flow (Δ<i>airway_pressure</i>/Δ<i>flow</i>) as resistance, and the <i>Tidal_volume</i>/(<i>airway_pressure-PEEP</i>) after flow has ceased for the compliance. This method, however, is also inaccurate since often the patient's effort will disrupt the inspiratory pause. A method of calculating resistance and compliance may be based on examining the initial pressure rise during the very first few hundredths of a second for each breath. This initial pressure rise is related to only the flow entering the system as volume has not yet accumulated in the lung, and thus reflects only the resistance. After resistance is accurately modeled, compliance can be obtained more accurately with the standard least squares technique mentioned above.

During data collection for the training data, the ventilator is adjusted, for example, the pressure support level may be adjusted and measured work of breathing and the other parameters are continuously monitored while changing the ventilator. Stable regions of the data are saved and parameters are averaged over 1 - 2 minutes to reduce the noise inherent in biological signals and the breath-to-breath variations. These averaged parameters are then used to create the model.

In an embodiment, the model, namely a neural network, is pretrained with clinical data and the input parameters can be collected non-invasively with a standard respiratory monitor. The neural network is trained to predict the physiologic and imposed WOB, POB, and PTP using the non-invasively acquired parameters described above (although invasive parameters may be added to the system, if desired.) Once a model having a desired degree of predictability has been achieved and verified, the esophageal pressure data (WOB, POB, PTP) is no longer needed and the network output, such as an actual breathing effort variable, may be used as an accurate predictor of patient effort.

Description of Neural Networks

Artificial neural networks loosely model the functioning of a biological neural network, such as the human brain. Accordingly, neural networks are typically implemented as computer simulations of a system of interconnected neurons. In particular, neural networks are hierarchical collections of interconnected processing elements (PEs). These elements are typically arranged in layers, where the input layer receives the input data, the hidden layers transform the data, and the output layer produces the desired output. Other embodiments of a neural network can also be used.

Each processing element in the neural network receives multiple input signals, or data values, that are processed to compute a single output. The inputs are received from the outputs of PEs in the previous layer or from the input data. The output value of a PE is calculated using a mathematical equation, known in the art as an activation function or a transfer function that specifies the relationship between input data values. As known in the art, the activation function may include a threshold, or a bias element. The outputs of elements at lower network levels are provided as inputs to elements at higher levels. The highest level element, or elements, produces a final system output, or outputs.

In the context of the present invention, the neural network is a computer simulation that is used to produce a noninvasive estimate of the quantified patient effort described previously. The neural network of the present invention may be constructed by specifying the number, arrangement, and connection of the processing elements which make up the network. A simple embodiment of a neural network consists of a fully connected network of processing elements. As shown in <figref idref="f0004">FIG. 4</figref>, the processing elements of the neural network are grouped into the following layers; an input layer <b>30</b> where the parameters collected and/or derived from the airway pressure and flow sensors are inputted to the network; a hidden layer or layers <b>32</b> of processing elements; and an output layer <b>34</b> where the resulting prediction of patient effort <b>36</b> is produced. The number of connections, and consequently the number of connection weights, is fixed by the number of elements in each layer <b>30, 32, 34.</b>

The most common training methodology for neural networks is based upon iterative improvement of the system parameters (normally called weights) by minimizing the mean squared difference between the desired output and the network output (mean squared error, MSE). The input is applied to the neural network, the neural network passes the data through its hierarchical structure, and an output is created. This network output is compared with the desired output corresponding to that input and an error is calculated. This error is then used to adjust the weights of the system so that the next time that particular input is applied to the system the network output will be closer to the desired output. There are many possible methodologies to adjust the weights, called the training algorithm. As shown in <figref idref="f0005">FIG. 5</figref>, the most common is called backpropagation that involves calculating each weight's responsibility for the error, and calculating a local gradient from this error in order to use a gradient descent learning rule for each weight.

Based on the foregoing specification, the invention may be implemented using computer programming or engineering techniques including computer software, firmware, hardware or any combination or subset thereof. Any such resulting program, having computer-readable code means, may be embodied or provided within one or more computer-readable media, thereby making a computer program product, i.e., an article of manufacture, according to the invention. The computer readable media may be, for instance, a fixed (hard) drive, diskette, optical disk, magnetic tape, semiconductor memory such as read-only memory (ROM), etc., or any transmitting/receiving medium such as the Internet or other communication network or link. The article of manufacture containing the computer code may be made and/or used by executing the code directly from one medium, by copying the code from one medium to another medium, or by transmitting the code over a network.

One skilled in the art of computer science will easily be able to combine the software created as described with appropriate general purpose or special purpose computer hardware to create a computer system or computer sub-system embodying the method of the invention. An apparatus for making, using or selling the invention may be one or more processing systems including, but not limited to, a central processing unit (CPU), memory, storage devices, communication links and devices, servers, I/O devices, or any sub-components of one or more processing systems, including software, firmware, hardware or any combination or subset thereof, which embody the invention. User input may be received from the keyboard, mouse, pen, voice, touch screen, or any other means by which a human can input data into a computer, including through other programs such as application programs.

Example 1:

With one particular subset of parameters (respiratory system resistance, compliance, tidal volume, and frequency), the present invention predicted power of breathing with a correlation coefficient well over 0.90. Similar results were obtained with other combinations of parameters. <figref idref="f0003">FIG. 3</figref> shows a plot of predicted versus actual, or measured, power of breathing. A perfect prediction would be indicated by all values having identical x- and y-axis values, and thus would fall on the diagonal black line 40 shown in the <figref idref="f0003">FIG. 3</figref>. The wider the spread from this line, the lower the correlation between predicted and actual power of breathing. In this particular prediction, data was obtained from 150 patients collected at Shands Hospital at the University of Florida. The data was collected with multiple ventilator settings for each patient, thus giving us approximately 500 data points in the study.

Example 2: Data on a Prediction (Table 1 includes the first 40-50 data points)

During data collection, respiratory measurements are made continuously and collected. These measurements are stored and later fed to the neural network for training. After the neural network model is trained successfully, the measurements from future are input to the "trained" model to produce a non-invasive prediction of actual power of breathing. This prediction of power of breathing is then used in a downstream fuzzy logic system that recommends changes to the ventilator.

In this example, the inputs/outputs to the system are: <ul id="ul0002" list-style="none"><li><i>Minute Ventilation -</i> the total volume delivered to the patient per minute (e.g. roughly tidal volume times breath frequency) in Liters.</li><li><i>Trigger Gradient -</i> the slope of the airway pressure drop due to the patient effort that triggers the breath (e.g. before the ventilator starts providing flow to the patient) in cm H20 per second.</li><li><i>Pressure Rise Time</i> - the time it takes for the flow to reach a maximum and then decay to 80% of the maximum during a PSV(pressure support ventilation) breath, as a proportion of total inspiratory time (no units).</li><li><i>Pmus using R Spike -</i> the muscle pressure (as described in the patent) using the resistance measurement from the initial time of pressure rise (as described in the patent) and using the least squares compliance (cm H20 *sec).</li><li><i>Actual POB -</i> the power of breathing (work per minute) measured (Joules/min).</li><li><i>Predicted POB -</i> the power of breathing predicted by the system (Joules/min).</li></ul>

These numbers were recorded continuously with data collection software. The values coming from the sensors are smoothed (time constant of roughly 2 minutes) and analyzed continuously as known in the art. When no errors, artifacts, or transients are detected, a set of parameters is labeled as "suitable for prediction", at one set of parameters per minute. These numbers are then collected and used to train the neural network. The neural network output is then used to optimally set the ventilator based on the predicted power of breathing and other relevant patient parameters (such as breathing frequency, tidal volume, and end-tidal CO2). <tables id="tabl0001" num="0001"><table frame="all"><title>Table 1</title><tgroup cols="6"><colspec colnum="1" colname="col1" colwidth="30mm" /><colspec colnum="2" colname="col2" colwidth="28mm" /><colspec colnum="3" colname="col3" colwidth="32mm" /><colspec colnum="4" colname="col4" colwidth="32mm" /><colspec colnum="5" colname="col5" colwidth="21mm" /><colspec colnum="6" colname="col6" colwidth="25mm" /><thead><row><entry align="center" valign="top">Minute Ventilation</entry><entry align="center" valign="top">Trigger Gradient</entry><entry align="center" valign="top">Pressure Rise Time</entry><entry align="center" valign="top">pmus using R Spike</entry><entry align="center" valign="top">Actual POB</entry><entry align="center" valign="top">Predicted POB</entry></row></thead><tbody><row><entry align="right">8.2</entry><entry align="right">14.2</entry><entry align="right">0.8</entry><entry align="right">63.4</entry><entry align="right">3.2</entry><entry align="right">4.1</entry></row><row><entry align="right">7.1</entry><entry align="right">6.4</entry><entry align="right">0.7</entry><entry align="right">40.7</entry><entry align="right">2.6</entry><entry align="right">2.7</entry></row><row><entry align="right">8.5</entry><entry align="right">15.3</entry><entry align="right">0.8</entry><entry align="right">72.9</entry><entry align="right">5.9</entry><entry align="right">6.2</entry></row><row><entry align="right">11.9</entry><entry align="right">9.5</entry><entry align="right">0.8</entry><entry align="right">81.2</entry><entry align="right">6.5</entry><entry align="right">6.2</entry></row><row><entry align="right">7.9</entry><entry align="right">8.0</entry><entry align="right">0.9</entry><entry align="right">37.7</entry><entry align="right">4.6</entry><entry align="right">5.8</entry></row><row><entry align="right">7.7</entry><entry align="right">12.4</entry><entry align="right">0.9</entry><entry align="right">64.3</entry><entry align="right">7.2</entry><entry align="right">6.5</entry></row><row><entry align="right">8.4</entry><entry align="right">13.4</entry><entry align="right">0.5</entry><entry align="right">51.5</entry><entry align="right">2.1</entry><entry align="right">2.1</entry></row><row><entry align="right">7.8</entry><entry align="right">8.1</entry><entry align="right">0.7</entry><entry align="right">79.4</entry><entry align="right">3.8</entry><entry align="right">3.4</entry></row><row><entry align="right">4.9</entry><entry align="right">15.2</entry><entry align="right">0.8</entry><entry align="right">41.5</entry><entry align="right">5.5</entry><entry align="right">5.0</entry></row><row><entry align="right">6.3</entry><entry align="right">26.3</entry><entry align="right">0.8</entry><entry align="right">41.4</entry><entry align="right">5.4</entry><entry align="right">5.7</entry></row><row><entry align="right">6.1</entry><entry align="right">2.3</entry><entry align="right">0.7</entry><entry align="right">37.3</entry><entry align="right">1.5</entry><entry align="right">2.0</entry></row><row><entry align="right">11.0</entry><entry align="right">15.4</entry><entry align="right">0.9</entry><entry align="right">83.5</entry><entry align="right">8.9</entry><entry align="right">9.0</entry></row><row><entry align="right">9.9</entry><entry align="right">7.5</entry><entry align="right">0.7</entry><entry align="right">74.2</entry><entry align="right">2.0</entry><entry align="right">3.2</entry></row><row><entry align="right">4.5</entry><entry align="right">8.1</entry><entry align="right">0.8</entry><entry align="right">47.3</entry><entry align="right">3.8</entry><entry align="right">3.2</entry></row><row><entry align="right">5.6</entry><entry align="right">3.4</entry><entry align="right">0.7</entry><entry align="right">21.8</entry><entry align="right">1.6</entry><entry align="right">2.3</entry></row><row><entry align="right">6.3</entry><entry align="right">12.0</entry><entry align="right">0.9</entry><entry align="right">77.0</entry><entry align="right">6.4</entry><entry align="right">5.8</entry></row><row><entry align="right">11.7</entry><entry align="right">11.2</entry><entry align="right">0.8</entry><entry align="right">66.4</entry><entry align="right">6.5</entry><entry align="right">7.0</entry></row><row><entry align="right">3.1</entry><entry align="right">12.9</entry><entry align="right">0.7</entry><entry align="right">31.3</entry><entry align="right">2.8</entry><entry align="right">2.4</entry></row><row><entry align="right">7.1</entry><entry align="right">2.0</entry><entry align="right">0.7</entry><entry align="right">15.0</entry><entry align="right">2.9</entry><entry align="right">2.3</entry></row><row><entry align="right">8.0</entry><entry align="right">7.4</entry><entry align="right">0.7</entry><entry align="right">79.5</entry><entry align="right">3.8</entry><entry align="right">3.4</entry></row><row><entry align="right">7.8</entry><entry align="right">4.4</entry><entry align="right">0.7</entry><entry align="right">22.1</entry><entry align="right">2.6</entry><entry align="right">2.5</entry></row><row><entry align="right">12.1</entry><entry align="right">26.3</entry><entry align="right">0.8</entry><entry align="right">179.2</entry><entry align="right">13.4</entry><entry align="right">11.4</entry></row><row><entry align="right">6.3</entry><entry align="right">6.4</entry><entry align="right">0.7</entry><entry align="right">21.4</entry><entry align="right">2.2</entry><entry align="right">2.7</entry></row><row><entry align="right">5.7</entry><entry align="right">18.3</entry><entry align="right">0.8</entry><entry align="right">31.1</entry><entry align="right">5.3</entry><entry align="right">5.4</entry></row><row><entry align="right">7.8</entry><entry align="right">10.9</entry><entry align="right">0.5</entry><entry align="right">73.9</entry><entry align="right">3.1</entry><entry align="right">2.5</entry></row><row><entry align="right">8.1</entry><entry align="right">6.1</entry><entry align="right">0.8</entry><entry align="right">56.5</entry><entry align="right">4.3</entry><entry align="right">3.7</entry></row><row><entry align="right">7.0</entry><entry align="right">16.6</entry><entry align="right">0.8</entry><entry align="right">69.5</entry><entry align="right">5.0</entry><entry align="right">5.6</entry></row><row><entry align="right">4.6</entry><entry align="right">13.5</entry><entry align="right">0.7</entry><entry align="right">26.9</entry><entry align="right">4.0</entry><entry align="right">3.0</entry></row><row><entry align="right">6.9</entry><entry align="right">10.9</entry><entry align="right">0.8</entry><entry align="right">29.4</entry><entry align="right">3.5</entry><entry align="right">5.0</entry></row><row><entry align="right">4.8</entry><entry align="right">14.0</entry><entry align="right">0.8</entry><entry align="right">39.7</entry><entry align="right">5.3</entry><entry align="right">4.8</entry></row><row><entry align="right">8.2</entry><entry align="right">33.5</entry><entry align="right">0.5</entry><entry align="right">95.5</entry><entry align="right">1.7</entry><entry align="right">2.8</entry></row><row><entry align="right">9.9</entry><entry align="right">11.4</entry><entry align="right">0.8</entry><entry align="right">35.6</entry><entry align="right">6.9</entry><entry align="right">5.9</entry></row><row><entry align="right">5.6</entry><entry align="right">18.5</entry><entry align="right">0.6</entry><entry align="right">30.2</entry><entry align="right">4.6</entry><entry align="right">3.7</entry></row><row><entry align="right">5.5</entry><entry align="right">12.1</entry><entry align="right">0.8</entry><entry align="right">49.5</entry><entry align="right">5.5</entry><entry align="right">4.9</entry></row><row><entry align="right">8.1</entry><entry align="right">12.4</entry><entry align="right">0.8</entry><entry align="right">82.0</entry><entry align="right">3.8</entry><entry align="right">4.3</entry></row><row><entry align="right">6.5</entry><entry align="right">16.0</entry><entry align="right">0.9</entry><entry align="right">46.0</entry><entry align="right">3.4</entry><entry align="right">5.4</entry></row><row><entry align="right">11.5</entry><entry align="right">16.7</entry><entry align="right">0.8</entry><entry align="right">45.9</entry><entry align="right">9.1</entry><entry align="right">9.7</entry></row><row><entry align="right">7.9</entry><entry align="right">9.3</entry><entry align="right">0.8</entry><entry align="right">59.0</entry><entry align="right">3.9</entry><entry align="right">4.7</entry></row><row><entry align="right">9.1</entry><entry align="right">7.6</entry><entry align="right">0.6</entry><entry align="right">47.4</entry><entry align="right">1.5</entry><entry align="right">2.5</entry></row><row><entry align="right">8.5</entry><entry align="right">4.0</entry><entry align="right">0.6</entry><entry align="right">32.0</entry><entry align="right">1.7</entry><entry align="right">2.2</entry></row><row><entry align="right">6.2</entry><entry align="right">13.6</entry><entry align="right">0.6</entry><entry align="right">33.1</entry><entry align="right">5.1</entry><entry align="right">3.3</entry></row><row><entry align="right">7.1</entry><entry align="right">21.6</entry><entry align="right">0.8</entry><entry align="right">77.6</entry><entry align="right">5.7</entry><entry align="right">5.2</entry></row><row><entry align="right">8.3</entry><entry align="right">13.7</entry><entry align="right">0.9</entry><entry align="right">87.6</entry><entry align="right">8.3</entry><entry align="right">6.8</entry></row><row><entry align="right">6.6</entry><entry align="right">25.2</entry><entry align="right">0.8</entry><entry align="right">38.8</entry><entry align="right">4.8</entry><entry align="right">5.7</entry></row><row><entry align="right">6.4</entry><entry align="right">19.2</entry><entry align="right">0.7</entry><entry align="right">36.5</entry><entry align="right">5.1</entry><entry align="right">5.1</entry></row><row><entry align="right">12.5</entry><entry align="right">12.8</entry><entry align="right">0.9</entry><entry align="right">93.4</entry><entry align="right">9.2</entry><entry align="right">9.9</entry></row><row><entry align="right">8.9</entry><entry align="right">14.7</entry><entry align="right">0.7</entry><entry align="right">71.2</entry><entry align="right">3.1</entry><entry align="right">3.9</entry></row></tbody></tgroup></table></tables>

It should be understood that the examples and embodiments described herein are for illustrative purposes only and that various modifications or changes in light thereof will be suggested to persons skilled in the art and are to be included within the scope of the claims.