Comparative study of four sigmoid models of pressurevolume curve in acute lung injury
 Guillermo M Albaiceta^{1, 2}Email author,
 Esteban Garcia^{1} and
 Francisco Taboada^{1, 3}
DOI: 10.1186/1475925X67
© Albaiceta et al; licensee BioMed Central Ltd. 2007
Received: 22 September 2006
Accepted: 14 February 2007
Published: 14 February 2007
Abstract
Background
The pressurevolume curve of the respiratory system is a tool to monitor and set mechanical ventilation in acute lung injury. Mathematical models of the static pressurevolume curve of the respiratory system have been proposed to overcome the inter and intraobserver variability derived from eyefitting. However, different models have not been compared.
Methods
The goodnessoffit and the values of derived parameters (upper asymptote, maximum compliance and points of maximum curvature) in four sigmoid models were compared, using pressurevolume data from 30 mechanically ventilated patients during the early phase of acute lung injury.
Results
All models showed an excellent goodnessoffit (R^{2} 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 withincase standard deviations of the pressures at the points of maximum curvature ranged from 2.33 to 6.08 cmH_{2}O.
Conclusion
There are substantial variabilities in relevant parameters obtained from the four different models of the static pressurevolume curve of the respiratory system.
Background
The static pressurevolume curve is one of the main measurements on lung mechanics [1]. 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 pressurevolume curves have been widely studied in acute lung injury and the acute respiratory distress syndrome [2]. They provide an assessment about the respiratory system mechanics, and their use has been proposed as a tool to set the mechanical ventilator [3].
Relevant data derived from the pressurevolume 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 [4]. 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 eyefitting has important interobserver variability [4, 5]. To overcome this limitation, pressurevolume data can be fitted to mathematical models and the parameters calculated in an objective manner. Although different equations have been proposed for modelling pressurevolume 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 pressurevolume 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.
Methods
Data from 30 patients (age 60 years [SD 14], APACHEII score 24 [SD 6] [10], lung injury score 3.04 [SD 0.29] [11], PaO_{2}/FiO_{2} 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 [14]. 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 volumecontrolled mode using an Evita 4 ventilator (Drager, Lubeck, Germany).
Pressurevolume curve maneuver
Pressure and flow were measured using the ventilator builtin heatedwire 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 pressurevolume curves were traced using changes in a Continuous Positive Airway Pressure (CPAP) [5], 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 cmH_{2}O. This pressure level was raised from 0 to 35 cmH_{2}O in 5 cmH_{2}O steps. After restoring basal ventilation for 3–5 minutes, three sighs more were applied and the ventilator switched again to CPAP mode at 0 cmH_{2}O, allowing the complete deflation of the lungs. After this, pressure was raised to 35 cmH_{2}O, and then decreased in 5 cmH_{2}O steps to 0 cmH_{2}O. 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 pressurevolume 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.
Volume  Pressure  

Model  Ref.  Equation  Lower asymptote  Upper asymptote  Inflection  LPMC  UPMC 
1  21  $V=a+\frac{b}{1+{e}^{(Pc)/d}}$  a  a + b  c  c  1.317d  c + 1.317d 
2  18  $V=\frac{b}{1+{e}^{(Pc)/d}}$  0  b  c  c  1.317d  c + 1.317d 
3  12  $V=\frac{{V}_{0}{V}_{0}{e}^{kP}}{1+{e}^{(Pc)/d}}$  0  V _{0}  *  *  * 
4  11  $V=a+\frac{b}{{\left[1+{e}^{(Pc)/d}\right]}^{s}}$  a  a + b  $cd\mathrm{ln}\frac{1}{s}$  $cd\mathrm{ln}\frac{(3s+1)\sqrt{5{s}^{2}+6s+1}}{2{s}^{2}}$  $cd\mathrm{ln}\frac{(3s+1)+\sqrt{5{s}^{2}+6s+1}}{2{s}^{2}}$ 
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 (HewlettPackard, Palo Alto, CA, USA).
Statistics
Data pairs of pressure and volume were fitted to each model using nonlinear regression with an iterative algorithm to minimize the sum of squared residuals (LevenbergMarquardt method). The convergence criteria for this sum were fixed at 10^{8}.
Data are presented as arithmetic mean and standard deviation. The R^{2} coeffients and the sums of squares (which represent the goodness of fit) were compared using an analysis of the variance [15], including the model and the curve limb as withincase factors. The values of the total lung capacity, maximum compliance and inflection points were compared using a oneway analysis of the variance. Posthoc test were done when appropriate using Bonferroni's correction [16]. The withincase standard deviation was used to estimate the magnitude of the differences between the models [17]. A p value lower than 0.05 was considered significant. All calculations were done using SPSS 11.0 software (SPSS inc, Chicago, IL, USA).
Results
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 posthoc tests). The withincase 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 withincase standard deviation was 117 ml.
There was a significant difference in the values of maximal inspiratory compliances in the four models (67 ml/cmH_{2}O [SD 18], 72 ml/cmH_{2}O [SD 20], 78 ml/cmH_{2}O [SD 24], 65 ml/cmH_{2}O [SD 18] for models 1 to 4 respectively, p = 0.019). Posthoc tests revealed a significant difference between model 4 and models 1 and 2 (p = 0.01 in both cases). The withincase standard deviation was 7 ml/cmH_{2}O. There were no differences in maximal expiratory compliances (59 ml/cmH_{2}O [SD 21], 59 ml/cmH_{2}O [SD 21], 57 ml/cmH_{2}O [SD 19], 58 ml/cmH_{2}O [SD 21] for models 1 to 4 respectively); p = 0.27). The withincase standard deviation was 4 ml/cmH_{2}O.
Points of maximum curvature
Withincase standard deviations of the values of pressure and volume at the points of maximum curvature.
Pressures (cmH_{2}O)  Volumes (ml)  

LPMC  UPMC  LPMC  UPMC  
Inspiratory limb  2.33  3.78  99  210 
Expiratory limb  6.08  1.82  184  370 
Discussion
Our results show that the goodness of fit of these four mathematical models of the static pressurevolume 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 pressurevolume curves
The diagnostic or prognostic utility of the static pressurevolume 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 [3], with the goals of improving oxygenation and oxygen transport [18] and minimizing ventilatorinduced lung injury [19]. 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 [20]. 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 [23]. Although there are experimental and clinical studies that challenge this interpretation of the pressurevolume curve [24], 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 [4] 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 [4]. 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 withincase standard deviation of 2.33 cmH_{2}O in the inspiratory LPMC implies that the difference between two models for the same subject is expected to be less than 6.45 cmH_{2}O (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 cmH_{2}O [25]. 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 pressurevolume curves in an automatized fashion. Previously published results suggest that the analysis of these curves should not be done by eye fitting [4], 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 goodnessoffit 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 [7] and patients [26] with injured lungs have confirmed this distribution. As suggested by Harris and coworkers [4], and demonstrated later [13], pressurevolume 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
 1.
It is a retrospective analysis of data. However, all the pressurevolume curves were performed using exactly the same methodology and the patients' characteristics in terms of course of the syndrome and severity were similar.
 2.
Limited pressure range. Only airway pressures up to 35 cmH_{2}O 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 [27].
 3.
Limited number of models: Other mathematical models were not considered. The selection was done because of the sigmoid nature of the equations.
Conclusion
Sigmoid equations allow fitting pressurevolume data acquired from patients with acute lung injury with an excellent goodnessoffit. 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.
Abbreviations
 CPAP:

Continuous positive airway pressure.
 LPMC:

Lower point of maximum curvature.
 UPMC:

Upper point of maximum curvature.
 PEEP:

Positive endexpiratory pressure.
Declarations
Acknowledgements
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 PI03/0833).
Authors’ Affiliations
References
 Harris RS: Pressurevolume curves of the respiratory system. Respir Care 2005, 50: 78–98. discussion 98–9Google Scholar
 Jonson B: Elastic pressurevolume curves in acute lung injury and acute respiratory distress syndrome. Intensive Care Med 2005, 31: 205–12. 10.1007/s0013400425179View 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 pressurevolume 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, LopezAguilar J, Taboada F, Blanch L: Application of continuous positive airway pressure to trace static pressurevolume 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 pressurevolume 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 pressurevolume 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 pressurevolume 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/0000324619851000000009View 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 pressurevolume curves to set endexpiratory pressure in acute lung injury. Intensive Care Med 2005, 31: 1370–8. 10.1007/s0013400527466View ArticleGoogle Scholar
 Albaiceta GM, Taboada F, Parra D, Luyando LH, Calvo J, Menendez R, Otero J: Tomographic study of the inflection points of the pressurevolume curve in acute lung injury. Am J Respir Crit Care Med 2004, 170: 1066–72. 10.1164/rccm.2003121644OCView ArticleGoogle Scholar
 Bernard GR, Artigas A, Brigham KL, Carlet J, Falke K, Hudson L, Lamy M, Legall JR, Morris A, Spragg R: The AmericanEuropean 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 endexpiratory 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 endexpiratory pressure during protective ventilation affects lung injury. Anesthesiology 2002, 97: 682–92. 10.1097/0000054220020900000023View ArticleGoogle Scholar
 Matamis D, Lemaire F, Harf A, BrunBuisson C, Ansquer JC, Atlan G: Total respiratory pressurevolume 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, LorenziFilho G, Kairalla RA, Deheinzelin D, Munoz C, Oliveira R, Takagaki TY, Carvalho CR: Effect of a protectiveventilation 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, PerezMendez L, AguirreJaime A: A high positive endexpiratory 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 pressurevolume curve in patients with acute respiratory distress syndrome. Curr Opin Crit Care 2002, 8: 32–8. 10.1097/0007519820020200000006View ArticleGoogle Scholar
 Brower RG, Lanken PN, MacIntyre N, Matthay MA, Morris A, Ancukiewicz M, Schoenfeld D, Thompson BT: Higher versus lower positive endexpiratory 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.200501048CPView ArticleGoogle Scholar
Copyright
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.