 Research
 Open Access
 Published:
Comparative study of four sigmoid models of pressurevolume curve in acute lung injury
BioMedical Engineering OnLine volume 6, Article number: 7 (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
Four mathematical models were studied [6–9]. All are sigmoid equations where volume and pressure are the dependent and independent variable respectively. Their equations and parameters are presented in Table 1.
The total lung capacity was calculated as the upper asymptote of the curve. The maximum compliance was calculated as the maximal slope of the pressurevolume curve (as compliance is the first derivative of the pressurevolume curve, maximal compliance is the value where the second derivative equals zero). The inflection point and points of maximum curvature were calculated as the pressure values where the second and third derivatives of the model equal zero respectively. Due to the sigmoid nature of the models, there are two points of maximum curvature in each limb (the socalled "lower" and "upper" inflection points [LPMC and UPMC]). The general values of the calculated parameters are presented in table 1. Figure 1 shows a pressurevolume curve from one representative patient, including the parameters obtained after fitting to the equations.
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
All sets of data were fitted to the models, R^{2} being always above 0.92. Model 4 provided the best fit (mean R^{2} 0.9985 [SD 0.0027]), followed by model 1 (mean R^{2} 0.9976 [SD 0.0035]), model 2 (mean R^{2} 0.9943 [SD 0.0061]) and model 3 (mean R^{2} 0.9849 [SD 0.0284]). These differences were significant (p = 0.011 in the ANOVA test, p < 0.01 in all posthoc tests). The differences between inspiratory and expiratory limb fitting (p = 0.023 in the ANOVA) and the interaction between model and limb (p = 0.005) were also significant. Posthoc tests revealed a significant difference between inspiratory and expiratory fitting only in model 3 (R^{2} 0.9982 [SD 0.0024] vs 0.9715 [SD 0.0497] for the inspiratory and expiratory limb respectively, p = 0.007). The same differences were found if the sums of squares instead of R^{2} values were taken into account (data not shown). Figure 2 shows data from one patient with the four fitted curves.
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
The models yielded different results for the maximum curvature points on the inspiratory limb of the pressurevolume curve, but not for those on the expiratory limb. In the same sense, lung volumes at these points were different only in the points from the inspiratory limb of the curve. All posthoc test were significant (p < 0.05) except the comparison between inspiratory UPMCs derived from models 1 and 4 (p = 0.30 and 0.53 for pressure and volume respectively). These results are presented in figure 3. The withincase standard deviations for these data are presented in table 2.
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
Our study is subjected to some limitations:

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.
References
Harris RS: Pressurevolume curves of the respiratory system. Respir Care 2005, 50: 78–98. discussion 98–9
Jonson B: Elastic pressurevolume curves in acute lung injury and acute respiratory distress syndrome. Intensive Care Med 2005, 31: 205–12. 10.1007/s0013400425179
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.35
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.
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.AA
Venegas JG, Harris RS, Simon BA: A comprehensive equation for the pulmonary pressurevolume curve. J Appl Physiol 1998, 84: 389–95.
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.
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 2184
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–1
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/0000324619851000000009
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.
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/s0013400527466
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.2003121644OC
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.
Altman DG, Bland JM: Comparing several groups using analysis of variance. BMJ 1996, 312: 1472–3.
Bland JM, Altman DG: Multiple significance tests: the Bonferroni method. BMJ 1995, 310: 170.
Bland JM, Altman DG: Measurement error. BMJ 1996, 313: 744.
Suter PM, Fairley B, Isenberg MD: Optimum endexpiratory airway pressure in patients with acute pulmonary failure. N Engl J Med 1975, 292: 284–9.
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/0000054220020900000023
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.
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/NEJM199802053380602
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.01
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.54
Hickling KG: Reinterpreting the pressurevolume curve in patients with acute respiratory distress syndrome. Curr Opin Crit Care 2002, 8: 32–8. 10.1097/0007519820020200000006
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/NEJMoa032193
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.
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.200501048CP
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).
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The author(s) declare that they have no competing interests.
Authors' contributions
GMA and FT designed the study. GMA acquired the data in patients and wrote the manuscript. GMA and EG made the mathematical and statistical analysis. All authors read and approved the manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
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.
About this article
Cite this article
Albaiceta, G.M., Garcia, E. & Taboada, F. Comparative study of four sigmoid models of pressurevolume curve in acute lung injury. BioMed Eng OnLine 6, 7 (2007). https://doi.org/10.1186/1475925X67
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1475925X67
Keywords
 Continuous Positive Airway Pressure
 Acute Lung Injury
 Total Lung Capacity
 Maximum Curvature
 Lung Mechanic