Comparative study of four sigmoid models of pressure-volume curve in acute lung injury
© Albaiceta et al; licensee BioMed Central Ltd. 2007
Received: 22 September 2006
Accepted: 14 February 2007
Published: 14 February 2007
The pressure-volume curve of the respiratory system is a tool to monitor and set mechanical ventilation in acute lung injury. Mathematical models of the static pressure-volume curve of the respiratory system have been proposed to overcome the inter- and intra-observer variability derived from eye-fitting. However, different models have not been compared.
The goodness-of-fit and the values of derived parameters (upper asymptote, maximum compliance and points of maximum curvature) in four sigmoid models were compared, using pressure-volume data from 30 mechanically ventilated patients during the early phase of acute lung injury.
All models showed an excellent goodness-of-fit (R2 always above 0.92). There were significant differences between the models in the parameters derived from the inspiratory limb, but not in those derived from the expiratory limb of the curve. The within-case standard deviations of the pressures at the points of maximum curvature ranged from 2.33 to 6.08 cmH2O.
There are substantial variabilities in relevant parameters obtained from the four different models of the static pressure-volume curve of the respiratory system.
The static pressure-volume curve is one of the main measurements on lung mechanics . It is constructed by plotting the airway, transpulmonary or pleural pressure against lung volume in static conditions, this is, in the absence of gas flow. Values of pressure and volume define two limbs in the curve, depending on whether the data was acquired during stepwise increments or decrements in pressure and volume. These inspiratory and expiratory limbs are usually different due to the respiratory system hysteresis. The pressure-volume curves have been widely studied in acute lung injury and the acute respiratory distress syndrome . They provide an assessment about the respiratory system mechanics, and their use has been proposed as a tool to set the mechanical ventilator .
Relevant data derived from the pressure-volume curve include total lung capacity, the maximal compliance and the value of the points of maximum curvature (which have been incorrectly cited in medical literature as "inflection points") of both limbs of the curve. There are different methods to determine these parameters . Manual methods are based on the visual identification of three linear zones of the curve, defining the points as the intersections between these three lines, the maximal compliance as the slope of the medial section and total lung capacity as the maximum value reached. However, this eye-fitting has important interobserver variability [4, 5]. To overcome this limitation, pressure-volume data can be fitted to mathematical models and the parameters calculated in an objective manner. Although different equations have been proposed for modelling pressure-volume curves, the sigmoid models [6–9] have shown an excellent fitting in experimental and clinical fields. However, some models have not been tested with clinical data, and the different models have not been compared.
Hypothesis and objectives
The objective of this study was to assess the goodness of fit of four previously published sigmoid models of pressure-volume curve using data from patients with acute lung injury, and to compare the values of the calculated parameters. Our hypothesis was that the different models would yield similar results.
Data from 30 patients (age 60 years [SD 14], APACHE-II score 24 [SD 6] , lung injury score 3.04 [SD 0.29] , PaO2/FiO2 146 mmHg [SD 43]) were retrospectively analyzed. Some data of these patients were previously published [12, 13] Informed consent was obtained for each patient's next of kin.
All patients were studied during the first 72 hours after meeting acute lung injury criteria . Exclusion criteria were: Less than 18 years old, contraindications to hypercapnia or elevated intrathoracic pressures (i.e. brain injury), air leaks, more than 5 days of mechanical ventilation or known or suspected chronic respiratory diseases. They were kept under deep sedation and muscle paralysis during the procedure, and ventilated in volume-controlled mode using an Evita 4 ventilator (Drager, Lubeck, Germany).
Pressure-volume curve maneuver
Pressure and flow were measured using the ventilator built-in heated-wire pneumotachograph and pressure transducer, connected to the patient via an endotracheal tube. Volume was calculated by integration of flow. Data was acquired and stored in a computer via serial port, at a sampling rate of 125 Hz, using the VentView software (Drager, Lubeck, Germany). Before the maneuver, the system was checked for air leaks during a prolonged inspiratory pause.
All pressure-volume curves were traced using changes in a Continuous Positive Airway Pressure (CPAP) , as described in detail in previously published papers [12, 13]. After three sighs (tidal volume 10 ml/Kg predicted body weight) to standardize volume history, the ventilator was switched to the CPAP mode, with a pressure level of 0 cmH2O. This pressure level was raised from 0 to 35 cmH2O in 5 cmH2O steps. After restoring basal ventilation for 3–5 minutes, three sighs more were applied and the ventilator switched again to CPAP mode at 0 cmH2O, allowing the complete deflation of the lungs. After this, pressure was raised to 35 cmH2O, and then decreased in 5 cmH2O steps to 0 cmH2O. Each pressure change was done only after reaching static conditions (zero flow). Changes in lung volume were computed for each pressure change. Data pairs of pressure and volume during the pressure increments were used to trace the inspiratory limb, and those from the decremental pressure maneuver used to trace the expiratory limb of the static pressure-volume curve.
Mathematical models of the curve
Equations for the mathematical models studied, and the general terms of lower and upper asymptotes, inflection point and points of maximum curvature. *The second derivative of model 3 is a transcendental equation, so a general term cannot be calculated.
a + b
c - 1.317d
c + 1.317d
c - 1.317d
c + 1.317d
a + b
The second derivative of model 3 is a transcendental equation, so a general term for pressure values when V" equals zero cannot be found. The values of inflection point and points of maximum curvature were numerically calculated using an iterative algorithm incorporated in a HP48 scientific calculator (Hewlett-Packard, Palo Alto, CA, USA).
Data pairs of pressure and volume were fitted to each model using non-linear regression with an iterative algorithm to minimize the sum of squared residuals (Levenberg-Marquardt method). The convergence criteria for this sum were fixed at 10-8.
Data are presented as arithmetic mean and standard deviation. The R2 coeffients and the sums of squares (which represent the goodness of fit) were compared using an analysis of the variance , including the model and the curve limb as within-case factors. The values of the total lung capacity, maximum compliance and inflection points were compared using a one-way analysis of the variance. Post-hoc test were done when appropriate using Bonferroni's correction . The within-case standard deviation was used to estimate the magnitude of the differences between the models . A p value lower than 0.05 was considered significant. All calculations were done using SPSS 11.0 software (SPSS inc, Chicago, IL, USA).
Goodness of fit
Upper asymptote and maximal compliances
The upper asymptotes yielded by model 3 were outside the expected limits in most of the cases (with values higher than 60 L), so this model was discarded in this comparison. There were significant differences in the predicted upper asymptote from the inspiratory curves (1836 ml [SD 680], 1647 ml [SD 532] and 2100 ml [SD 797] for the models 1, 2 and 4 respectively, p < 0.001 in the ANOVA and all post-hoc tests). The within-case standard deviation was 319.7 ml. However, the predicted upper asymptotes derived from expiratory curves were similar in the three models (1523 ml [SD 589], 1517 ml [SD 486] and 1536 ml [SD 577] ml for models 1, 2 and 4 respectively, p = 0.54). The within-case standard deviation was 117 ml.
There was a significant difference in the values of maximal inspiratory compliances in the four models (67 ml/cmH2O [SD 18], 72 ml/cmH2O [SD 20], 78 ml/cmH2O [SD 24], 65 ml/cmH2O [SD 18] for models 1 to 4 respectively, p = 0.019). Post-hoc tests revealed a significant difference between model 4 and models 1 and 2 (p = 0.01 in both cases). The within-case standard deviation was 7 ml/cmH2O. There were no differences in maximal expiratory compliances (59 ml/cmH2O [SD 21], 59 ml/cmH2O [SD 21], 57 ml/cmH2O [SD 19], 58 ml/cmH2O [SD 21] for models 1 to 4 respectively); p = 0.27). The within-case standard deviation was 4 ml/cmH2O.
Points of maximum curvature
Within-case standard deviations of the values of pressure and volume at the points of maximum curvature.
Our results show that the goodness of fit of these four mathematical models of the static pressure-volume curve is remarkably high. However, in contrast to our hypothesis, different models yield different results when applied in mechanically ventilated patients with acute lung injury. These differences are limited to the parameters derived from the inspiratory limb of the curve and could be relevant for the interpretation of the data.
Utility of pressure-volume curves
The diagnostic or prognostic utility of the static pressure-volume cuve in acute lung injury has not been demonstrated. However, there is substantial interest in these curves as a tool to set mechanical ventilation in acute lung injury , with the goals of improving oxygenation and oxygen transport  and minimizing ventilator-induced lung injury . The lower point of maximum curvature on the inspiratory limb of the curve has been used as a marker of optimal PEEP, and the upper as a marker of alveolar overstretching . Three clinical trials have shown that a PEEP level selected according these points is related to a decrease in mortality [21, 22] and a lower inflammatory response casued by mechanical ventilation . Although there are experimental and clinical studies that challenge this interpretation of the pressure-volume curve , there is considerable interest in the use of lung mechanics to titrate therapy.
Some of the discrepancies between studies could be attributed to different approaches to determinate these points of maximum curvature. Harris et al  found substantial intra- and interobserver variability during the determination of the lower point of maximum curvature, especially with small data sets. The application of mathematical models could overcome these differences. Most studies have used the model described by Venegas and coworkers . This consists in a sigmoid equation, symmetric to the inflection point of the curve. Other models have been proposed in an attempt to improve the flexibility to fit asymmetrical data [8, 9], but some have not been tested in the clinical arena. Our results show that the fitting of data to different equations can lead to different values of the derived parameters. The magnitude of these differences is small in the majority of the results, but is substantial in the points of maximum curvature, which are the parameters more frequently used and studied. For instance, a within-case standard deviation of 2.33 cmH2O in the inspiratory LPMC implies that the difference between two models for the same subject is expected to be less than 6.45 cmH2O (2.77 times the standard deviation) for 95% of the pairs of observations. This is a considerable range of pressures if LPMC is used for PEEP setting. In a clinical trial in patients with acute lung injury, changes in PEEP were done using a table with steps of 2 cmH2O . This shows that the differences between the models are above the range of pressures used in clinical practice.
Currently, some ventilators include tools to trace pressure-volume curves in an automatized fashion. Previously published results suggest that the analysis of these curves should not be done by eye fitting , so a mathematical analysis software must be included. Our results suggest that data about the mathematical model should be specified and taken into account in order to compare the results.
Characteristics of the mathematical models
After the publication of a mathematical model by Venegas and coworkers (model 1 in this study), other models were published. Model 2 and 3 are based on the assumption that the lower asymptote is 0. This assumption is probably incorrect in some cases with air trapping, but it allows to reduce the number of parameters in the equation.
The other models were developed to overcome the symmetric nature of the original model. An asymmetric equation allows to improve the goodness of fit (in our study, the best result was obtained with model 4, although this model has 5 fitting parameters), but an increment in the number of parameters can lead to overfitting (an increase in the goodness-of-fit but a decrease in the significance and robustness of the parameters). In the range of pressures used in a clinical setting, symmetric models have shown an excellent fitting. However, in cases where data range reaches high values of pressure, an asymmetric model could be superior.
The symmetry of models 1 and 2 reflects a Gaussian distribution of airway opening and closing pressures. Tomographic studies in animals  and patients  with injured lungs have confirmed this distribution. As suggested by Harris and coworkers , and demonstrated later , pressure-volume curves should be considered as recruitment functions in acute lung injury, and the contribution of overstretching (which would be the prevalent phenomenon at lung volumes near the upper asymptote) to the shape of the curve is low. This could explain also that during deflation, as there is no ongoing recruitment, the four models yield similar results.
The main limitation of these models is that they are descriptive in nature, but without a mechanistic background. They come from the intention to fit the data and calculate parameters in a standardized fashion. However, they are not based in the behaviour of lung tissue during inflation and deflation. This descriptive nature could be one of the causes of the differences between the four models. Other could be related to the small number of data points used for fitting and the fact that what models yield as a concrete point is, from a biological point of view, a range of pressures where abrupt changes in lung mechanics occur. A mechanistic model, based in the mechanical properties of the injured lungs, should be developed to move beyond a simple description of the curve.
Limitations of the study
It is a retrospective analysis of data. However, all the pressure-volume curves were performed using exactly the same methodology and the patients' characteristics in terms of course of the syndrome and severity were similar.
Limited pressure range. Only airway pressures up to 35 cmH2O were studied. This pressure level is probably associated to an incomplete aeration in some cases. A wider pressure range could have improved the results of the upper asymptote, especially in model 3, although the influence in the points of maximum curvature should be low. The range of pressures studied was selected according to clinical recommendations of airway pressure limitation .
Limited number of models: Other mathematical models were not considered. The selection was done because of the sigmoid nature of the equations.
Sigmoid equations allow fitting pressure-volume data acquired from patients with acute lung injury with an excellent goodness-of-fit. However, there could be significant differences in the derived parameters from these models, specifically in the points of maximum curvature (which are used to set the PEEP level during mechanical ventilation). These differences should be taken into account when applying a model to a set of data or when comparing results from different studies.
Continuous positive airway pressure.
Lower point of maximum curvature.
Upper point of maximum curvature.
Positive end-expiratory pressure.
The authors thank Jose G. Venegas his comments about the study. Supported in part by Instituto de Salud Carlos III (Fondo de Investigación Sanitaria FIS PI-03/0833).
- Harris RS: Pressure-volume curves of the respiratory system. Respir Care 2005, 50: 78–98. discussion 98–9Google Scholar
- Jonson B: Elastic pressure-volume curves in acute lung injury and acute respiratory distress syndrome. Intensive Care Med 2005, 31: 205–12. 10.1007/s00134-004-2517-9View ArticleGoogle Scholar
- Gattinoni L, Eleonora C, Caironi P: Monitoring of pulmonary mechanics in acute respiratory distress syndrome to titrate therapy. Curr Opin Crit Care 2005, 11: 252–8. 10.1097/01.ccx.0000160773.43122.35View ArticleGoogle Scholar
- Harris RS, Hess DR, Venegas JG: An objective analysis of the pressure-volume curve in the acute respiratory distress syndrome. Am J Respir Crit Care Med 2000, 161: 432–9.View ArticleGoogle Scholar
- Albaiceta GM, Piacentini E, Villagra A, Lopez-Aguilar J, Taboada F, Blanch L: Application of continuous positive airway pressure to trace static pressure-volume curves of the respiratory system. Crit Care Med 2003, 31: 2514–9. 10.1097/01.CCM.0000090003.87219.AAView ArticleGoogle Scholar
- Venegas JG, Harris RS, Simon BA: A comprehensive equation for the pulmonary pressure-volume curve. J Appl Physiol 1998, 84: 389–95.Google Scholar
- Pelosi P, Goldner M, McKibben A, Adams A, Eccher G, Caironi P, Losappio S, Gattinoni L, Marini JJ: Recruitment and derecruitment during acute respiratory failure: an experimental study. Am J Respir Crit Care Med 2001, 164: 122–30.View ArticleGoogle Scholar
- Henzler D, Orfao S, Rossaint R, Kuhlen R: Modification of a sigmoidal equation for the pulmonary pressure-volume curve for asymmetric data. J Appl Physiol 2003, 95: 2183–4. author reply 2184View ArticleGoogle Scholar
- Heller H, Brandt S, Schuster KD: Development of an algorithm for improving the description of the pulmonary pressure-volume curve. J Appl Physiol 2002, 92: 1770. author reply 1770–1View ArticleGoogle Scholar
- Knaus WA, Draper EA, Wagner DP, Zimmerman JE: APACHE II: a severity of disease classification system. Crit Care Med 1985, 13: 818–29. 10.1097/00003246-198510000-00009View ArticleGoogle Scholar
- Murray JF, Matthay MA, Luce JM, Flick MR: An expanded definition of the adult respiratory distress syndrome. Am Rev Respir Dis 1988, 138: 720–3.View ArticleGoogle Scholar
- Albaiceta GM, Luyando LH, Parra D, Menendez R, Calvo J, Pedreira PR, Taboada F: Inspiratory vs. expiratory pressure-volume curves to set end-expiratory pressure in acute lung injury. Intensive Care Med 2005, 31: 1370–8. 10.1007/s00134-005-2746-6View ArticleGoogle Scholar
- Albaiceta GM, Taboada F, Parra D, Luyando LH, Calvo J, Menendez R, Otero J: Tomographic study of the inflection points of the pressure-volume curve in acute lung injury. Am J Respir Crit Care Med 2004, 170: 1066–72. 10.1164/rccm.200312-1644OCView ArticleGoogle Scholar
- Bernard GR, Artigas A, Brigham KL, Carlet J, Falke K, Hudson L, Lamy M, Legall JR, Morris A, Spragg R: The American-European Consensus Conference on ARDS. Definitions, mechanisms, relevant outcomes, and clinical trial coordination. Am J Respir Crit Care Med 1994, 149: 818–24.View ArticleGoogle Scholar
- Altman DG, Bland JM: Comparing several groups using analysis of variance. BMJ 1996, 312: 1472–3.View ArticleGoogle Scholar
- Bland JM, Altman DG: Multiple significance tests: the Bonferroni method. BMJ 1995, 310: 170.View ArticleGoogle Scholar
- Bland JM, Altman DG: Measurement error. BMJ 1996, 313: 744.View ArticleGoogle Scholar
- Suter PM, Fairley B, Isenberg MD: Optimum end-expiratory airway pressure in patients with acute pulmonary failure. N Engl J Med 1975, 292: 284–9.View ArticleGoogle Scholar
- Takeuchi M, Goddon S, Dolhnikoff M, Shimaoka M, Hess D, Amato MB, Kacmarek RM: Set positive end-expiratory pressure during protective ventilation affects lung injury. Anesthesiology 2002, 97: 682–92. 10.1097/00000542-200209000-00023View ArticleGoogle Scholar
- Matamis D, Lemaire F, Harf A, Brun-Buisson C, Ansquer JC, Atlan G: Total respiratory pressure-volume curves in the adult respiratory distress syndrome. Chest 1984, 86: 58–66.View ArticleGoogle Scholar
- Amato MB, Barbas CS, Medeiros DM, Magaldi RB, Schettino GP, Lorenzi-Filho G, Kairalla RA, Deheinzelin D, Munoz C, Oliveira R, Takagaki TY, Carvalho CR: Effect of a protective-ventilation strategy on mortality in the acute respiratory distress syndrome. N Engl J Med 1998, 338: 347–54. 10.1056/NEJM199802053380602View ArticleGoogle Scholar
- Villar J, Kacmarek RM, Perez-Mendez L, Aguirre-Jaime A: A high positive end-expiratory pressure, low tidal volume ventilatory strategy improves outcome in persistent acute respiratory distress syndrome: a randomized, controlled trial. Crit Care Med 2006, 34: 1311–8. 10.1097/01.CCM.0000215598.84885.01View ArticleGoogle Scholar
- Ranieri VM, Suter PM, Tortorella C, De Tullio R, Dayer JM, Brienza A, Bruno F, Slutsky AS: Effect of mechanical ventilation on inflammatory mediators in patients with acute respiratory distress syndrome: a randomized controlled trial. JAMA 1999, 282: 54–61. 10.1001/jama.282.1.54View ArticleGoogle Scholar
- Hickling KG: Reinterpreting the pressure-volume curve in patients with acute respiratory distress syndrome. Curr Opin Crit Care 2002, 8: 32–8. 10.1097/00075198-200202000-00006View ArticleGoogle Scholar
- Brower RG, Lanken PN, MacIntyre N, Matthay MA, Morris A, Ancukiewicz M, Schoenfeld D, Thompson BT: Higher versus lower positive end-expiratory pressures in patients with the acute respiratory distress syndrome. N Engl J Med 2004, 351: 327–36. 10.1056/NEJMoa032193View ArticleGoogle Scholar
- Crotti S, Mascheroni D, Caironi P, Pelosi P, Ronzoni G, Mondino M, Marini JJ, Gattinoni L: Recruitment and derecruitment during acute respiratory failure: a clinical study. Am J Respir Crit Care Med 2001, 164: 131–40.View ArticleGoogle Scholar
- Hager DN, Krishnan JA, Hayden DL, Brower RG: Tidal volume reduction in patients with acute lung injury when plateau pressures are not high. Am J Respir Crit Care Med 2005, 172: 1241–5. 10.1164/rccm.200501-048CPView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.