 Research
 Open Access
 Published:
Nonlinear parameters of specific resistance loops to characterise obstructive airways diseases
Respiratory Researchvolume 18, Article number: 9 (2017)
Abstract
Background
Specific resistance loops appear in different shapes influenced by different resistive properties of the airways, yet their descriptive ability is compressed to a single parameter  its slope. We aimed to develop new parameters reflecting the various shapes of the loop and to explore their potential in the characterisation of obstructive airways diseases.
Methods
Our study included 134 subjects: Healthy controls (N = 22), Asthma with nonobstructive lung function (N = 22) and COPD of all disease stages (N = 90). Different shapes were described by geometrical and secondorder transfer function parameters.
Results
Our parameters demonstrated no difference between asthma and healthy controls groups, but were significantly different (p < 0.0001) from the patients with COPD. Grouping mild COPD subjects by an open or notopen shape of the resistance loop revealed significant differences of loop parameters and classical lung function parameters. Multiple logistic regression indicated RV/TLC as the only predictor of loop opening with OR = 1.157, 95% CI (1.064–1.267), pvalue = 0.0006 and R^{2} = 0.35. Inducing airway narrowing in asthma gave equal shape measures as in COPD nonopeners, but with a decreased slope (p < 0.0001).
Conclusion
This study introduces new parameters calculated from the resistance loops which may correlate with different phenotypes of obstructive airways diseases.
Background
Obstructive lung diseases, predominantly asthma and COPD, are a group of respiratory diseases characterized by airflow limitation [1, 2]. The primary pathophysiologic impairment in these diseases is an increase of airways resistance that originates from mucosal and submucosal inflammation, bronchial constriction and airway collapse during expiration. The rise in airways resistance requires greater intrapleural pressure changes to provide sufficient pressure gradients for the initiation and maintenance of airflow. Airway resistance particularly increases during expiration when positive intrathoracic pressures are further compressing the intraluminal space of the airways [3].
From a diagnostic point of view spirometry is considered as the gold standard to diagnose COPD and asthma, as well as to assess the level of airway obstruction [4, 5]. Spirometry requires forced maximal manoeuvres but cannot quantify increased resistance of the airways at tidal breathing, which may be characteristic and specific for the underlying disease. Wholebody plethysmography, however, allows the computation of airways resistance by measuring alveolar pressure changes and corresponding airway flows during tidal breathing in a closed body box [6]. More specifically, plethysmography records the small changes in box pressure by compression and decompression of thoracic gas which correspond to small volume changes, also known as shift volume. Simultaneously, it records airflow at the mouth generated by these subtle pressure changes. The relationship between these two measurements is linearly expressed as specific airways conductance (sG_{AW}) or specific airways resistance (sR_{AW}) [7]. Hence, the diagnostic power of resistance loops is reduced to one single parameter, the slope of the linear relation, rather than to its differences in shape which are affected by central or peripheral airway obstruction, expiratory collapse or endexpiratory airway closing, and which may vary in function of the breathing manoeuver and the thoracic gas volume (Fig. 1) [6, 7]. As demonstrated in our previous work, the slope of these loops is lacking discriminatory ability. It only marginally contributes to the diagnosis and differentiation of COPD and asthma [8].
Different studies have identified some limitations of the current parameterization of specific resistance loops, mainly by its variability, the necessity for artefact corrections and the clear evidence that the slope is significantly determined by small changes in the breathing frequency [9–11]. Similar observations were done with the forced oscillation technique (FOT), demonstrating important within and betweenbreath variation for which should be corrected to identify real airway patency or obstruction [12–15]. Other researchers have described methodologies to optimize the slope estimation as a key parameter for the calculation of airways resistance with plethysmography [7, 16, 17]. However, to the best of our knowledge, no one has explored the potential of other parameters that may better reflect the curvilinear twodimensional shape of plethysmographic resistance loops.
We hypothesize that a more detailed mathematical modelling of the tidal breathing loops may yield new parameters which better represent the nonlinear dynamics of specific resistance curves. With an empirical modelbased approach our first objective was to develop geometrical parameters reflecting the shape of the loop. Secondly, we designed secondorder mathematical models linking shift volumes to airflows during breathing. Finally, we checked whether these newly developed parameters associated with other important lung function measures in different clinical phenotypes of obstructive airways diseases.
Methods
Study subjects
For development of different mathematical models, lung function data of 134 subjects were retrospectively collected from a databank of the outpatient clinic of the University Hospital of Leuven (Belgium) based on an established diagnosis of asthma (N = 22) or COPD (N = 90) and compared with a healthy control group (N = 22). All enrolled subjects were Caucasians between 19 and 84 years old who had performed complete pulmonary function testing (including postbronchodilator spirometry, wholebody plethysmography for lung volumes and airway resistance, and diffusing capacity). All patients provided informed consent for the use of lung function and clinical data and the protocol was approved by the ethics committee of the University Hospital of Leuven. Baseline characteristics of the subjects are presented in the Table 1: i) Healthy group is defined as a group of asymptomatic subjects without smoking history with all lung function parameters falling in the normal reference range; ii) Asthma group: subjects with or without smoking history, a previous clinical diagnosis of asthma based on symptoms and therapy response, with a nonobstructive lung function (post bronchodilator FEV_{1}/FVC ratio > 0.7) and positive reaction on inhaled methacholine (FEV_{1} drop of at least 20% at maximal concentration of 8 mg/ml of methacholine) [18]; and iii) COPD group: subjects with at least 10 packyears of smoking history and postbronchodilator FEV_{1}/FVC ratio below 0.7 [19]. In a subgroup of asthma subjects (N = 11), a methacholine challenge test was repeated with pre and post resistance and lung volume measurements.
Pulmonary function tests
All pulmonary function tests were performed according to the American Thoracic Society (ATS)/ European Respiratory Society (ERS) criteria [20] using standardized equipment (Masterscreen Jaeger, Carefusion, Germany). Spirometry data are postbronchodilator measures and are expressed, along with prebronchodilator plethysmography measurements of airway resistance and lung volume as percent predicted of normal reference values [21, 22]. Diffusing capacity (DL,_{CO}) was measured by the singlebreath carbon monoxide gas transfer method and expressed as percent predicted of reference values [23].
Geometrical modelling
To geometrically explain the opening of the resistance loops we developed a series of new parameters that may capture what is visually apparent to the observer. From each expiratory resistance loop, using the same electronic data and MATLAB software as in the transfer function development, the following shape descriptive parameters were derived:

1/
Area of the expiratory loop: it stands for total surface covered by the expiratory phase of the breathing manoeuver (Fig. 2, panel I.A). Area of loop (AOL) per se is effort depended, as deeper breathing will cause larger loop surfaces. Therefore, we normalized each loop to a range [0, 1] for both flows and shift volumes.

2/
Roundness of the expiratory loop: In essence, this parameter measures how closely the expiratory loop approaches to a circle (Fig. 2, panel I.B). It appears in an interval from 0 to 1, where 0 means complete closure of the loop (flat line) while 1 means a perfect circle. Smooth eclipses will have a low roundness and an open loop is never expected to reach a perfect circle. Roundness (Rnd) is defined based on (Eq. 1):
$$ Rnd=\frac{4\uppi \ast \mathrm{A}\mathrm{O}\mathrm{L}}{\mathrm{Perimete}{\mathrm{r}}^2} $$(1)where AOL represents the complete area of the normalised expiratory loop, and Perimeter is the perimeter of the same loop.

3/
Median point of the expiratory loop: it is a two dimensional measure, with the values coming from both, X axis (median volume shift = point X) and Y axis (median flow = point Y) in the expiratory phase (Fig. 2, panel I.C). Ordinarily, both axes have a negative value as somewhat longer expiration and larger opening of the loop pulls the point lower to the third quadrant.

4/
Asynchrony between volume shift and flow: The nonlinearity in the relationship between volume shift and flow lies in the asynchrony between this two factors. It corresponds to airflow drops despite increasing volume shifts (Fig. 2, panel I.D, it can be also depicted from Fig. 2, panel II.C). Asynchrony is a measure expressed as time difference of volume shift peak and peak of the flow.

5/
SG _{ 0.5 } : Similar to sG_{AW}, yet limited to the linear drop between righthandside inspiratory and expiratory flow rates of ±0.5 L/s (Fig. 2, panel I.E.) [6].
Transfer function modelling
Development of the databased input–output transfer function models was performed in an offline framework in MATLAB (8.3, The MathWorks, Natick, Massachusetts) using the System Identification Toolbox [24]. From the Masterlab system, measurements of tidal expiration were exported at a sampling rate of 100Hz. Shape differences of specific resistance loop were described via relationship of their two creating factors: flow and volume shift. Volume shift was used as a model input, while the result of that generated pressure, flow, was used as a model output (input–output relationship shown in Fig. 2, panel II).
To characterize the breathing process and explain the data in a parametrically efficient way and yet sustain simplicity in the sense of model parameters and model order, an iterative system identification procedure was used [25]. Based on this procedure we have chosen a second order discretetime transfer function (TF) model (Eq. 2) as most the appropriate.
where F stands for Flow and Vs for Volume shift. Coefficients a _{ 1 }, a _{ 2 } and coefficients b _{ 0 }, b _{ 1 } and b _{ 2 } are the TF denominator and numerator, respectively. Z is a discrete domain operator. From the defined second order model, we further extracted process descriptive parameters:

1/
the steady state gain (SSG), defined as a ratio of the steady state output and the input (Eq. 3), which in simplified sense represents volume shift multiplication to reach a certain flow level in steady state conditions.
$$ SSG=\frac{\varDelta F(z)}{\varDelta Vs(z)}=\frac{{\displaystyle {\sum}_{i=0}^2{b}_i}}{1+{\displaystyle {\sum}_{i=1}^2{a}_i}} $$(3) 
2/
using denominator coefficients, two dynamic components of tidal expiration were derived: namely Pole _{ 1 } and Pole _{ 2 } (i.e. the roots of the denominator) (Eq. 4).
$$ Pol{e}_{1,2} = \frac{a_1\pm \sqrt{{a_1}^24{a}_2}}{2} $$(4) 
3/
finally, the numerator coefficients (b _{ 0 } , b _{ 1 } , and b _{ 2 }), were used for group comparison.
Statistical analysis
JMP Pro version 12, (SAS Institute, Cary, USA) was used to perform statistical analysis. The ShapiroWilk test was used to inspect normality of the groups. To control differences between two groups (paired or unpaired, where appropriate) with parametric and nonparametric distribution Ttest and Mann–Whitney test were used, respectively. In the case of multiple group comparisons, ANOVA or KruskalWallis tests were performed, depending of group distribution. The logisticregression model was applied for binary variables analyses, where in addition stepwise selection was used to identify the subset of variables that had the strongest relation to outcome, using default criteria of significance at the 0.25 level to enter and 0.5 level to leave the model. The model consisted of lung function parameters: FVC, %pred, FEV_{1}, %pred, FEV_{1}/FVC, %pred, PEF, %pred. TLC, %pred, RV, %pred, FRC, %pred., DL,_{CO}, %pred, K_{CO}, %pred., RV/TLC. To determine normal value range, twosided prediction interval of 99% was applied.
Results
Models
The assessment of the geometrical parameters was possible with the data of all subjects. Due to instable results of TF model parameters (z > 1), one healthy and ten COPD subjects were excluded from further analysis [26]. In general, confirmation of the appropriate model selection was demonstrated with a high goodness of fit expressed as normalized root mean square error (NRMSE) of 90 (85–92)% (values are median and IQR) for complete dataset. A visualisation of the model performance with its accuracy is shown in Fig. 2, Panel III.
Table 2 depicts the median values of all resistive parameters in each group, revealing significant differences in the majority of parameters when comparing COPD with Asthma or Healthy controls, and no significant differences between nonobstructive asthmatics and healthy controls. As expected, all geometrical parameters reflecting the opening of the loop (Roundness, Area of loop and Asynchrony) are significantly different in COPD. Point X and b’s from TF confirm that with airways obstruction in COPD, an increased pressure (increased volume shift) has to be generated to maintain tidal breathing without significant changes in flow (Point Y). This is coherent with the measuring standard of wholebody plethysmography where the subjects are instructed for a tidal breathing at 1 Hz at FRC securing similar flows for everyone (disease independent).
Resistance loop in COPD
As shown in Table 2, most resistance loops of COPD opened to an imperfect circular shape with significant differences on roundness, normalized AOL and asynchrony. As multivariate analysis indicated that roundness from all resistive parameters was the best differentiator between nonobstructive subjects (asthma and healthy) and COPD, roundness was chosen as the key variable to describe opening of the loop (data not shown). When applying a 99% confidence interval to define the normal range of roundness in healthy and nonobstructive asthma subjects, a 0.18 cutoff came out as the upper limit of normality. When using this > 0.18 cutoff for the definition of loop opening in COPD, 67/90 patients were considered as “Openers” whereas 23 were defined as “Nonopeners”. When comparing common lung function parameters between the group of “Openers” and “Nonopeners” in COPD, it was found that 22 of 23 “Nonopeners” were present in the mild to moderate stages (I and II) of COPD and that those subjects were characterized by a significantly higher FEV_{1}, a lower RV and RV/TLC ratio compared to the “openers” of GOLD stage I and II (Table 3). Multiple logistic regression with the stepwise selection demonstrated that in COPD GOLD I and II, RV/TLC ratio was the only predictor for loop opening with an OR of 1.157 (95% CI of (1.064–1.267), p value of 0.0006 and model R^{2} = 0.35). When looking at all COPD stages, again RV/TLC associated best with loop opening with an OR of 1.168 (95% CI of (1.092–1.259), p value <0.0001 and model R^{2} = 0.61) and a small contribution of FEV_{1}/FVC ratio (OR (CL) 0.866 (0.746–0.949), p value 0.0004, adding to total R^{2} = 0.72). Data presented in Table 4.
Changes in resistance loops during bronchoconstriction
In a subgroup of nonobstructive asthma patients (N = 11), airway narrowing was induced by increasing levels of methacholine during a challenge test. In all patients FEV_{1} dropped below 20% of its initial value at methacholine concentrations <4 mg/ml, indicative for acute bronchoconstriction with subsequent rises in airways resistance. These changes were confirmed by a significant decrease of specific conductance (sG_{AW} and sG_{0.5}), a significant increase in median volume shift (point X) and b _{ 0 } and significant drop in b _{ 1 } and SSG (data shown in Table 5). Although the median roundness statistically increased in the postmethacholine challenge conditions (from 0.08 to 0.13), these small changes were considered not to be relevant as in 7/11 patients roundness stayed under the upper limit of “Nonopeners”. Additionally, normalized AOL and asynchrony did not change, indicating that differences in the slope and changes in input related variables (point X, b _{ 0 }, b _{ 1 }, SSG) were perfectly in phase with changes in output (flow).
Finally, when comparing resistive and lung function parameters of obstructed asthma (post methacholine) with mild to moderate COPD “Openers” and “Nonopeners”, it was clear that obstructive asthma patients presented with the lowest slopes and highest point X for only little increases of roundness, asynchrony and AOL. By contrast, COPD openers presented with the highest roundness, asynchrony and AOL. Consistent with previous findings, not FEV_{1} (p = 0.093) but RV/TLC ratio (p < 0.0001) was the main clinical lung function determinant that differentiated between obstructive asthma and COPD openers (see Table 6).
Discussion
In this study we describe specific resistance loops measured by bodyplethysmography by using geometrical analyses and secondorder transfer functions. The newly developed parameters reflect the curvilinear aspect and rotation of the resistance loops and typically associate with different lung function characteristics. Loops that are mathematically identified by an open appearance as measured by roundness typically occur in the majority of COPD subjects with hyperinflation, whereas rotation of the slope without opening is apparent during bronchoconstriction, as demonstrated in asthma.
Opening of the loops in COPD was significantly associated with RV/TLC ratio, much more than with FEV_{1} or FEV_{1}/FVC ratio. In numerous studies the RV/TLC ratio has been linked to air trapping and hyperinflation, together with increases of RV and FRC [27–29]. In fact, the RV/TLC ratio reflects the proportion of trapped lung volume that cannot be mobilized by maximal breathing. Increases in static RV/TLC ratio are inversely correlated with maximal inspiratory capacity, which often declines during exercise by dynamic hyperinflation and strongly associates with breathing discomfort and dyspnea. From a mechanistic point of view, we can only speculate on the reasons why RV/TLC ratio is the best predictor for loop opening and asynchrony between alveolar pressures and flows. One possible explanation may be found in the competition for space between lung areas that are emptied during expiration with areas that are trapped with air and progressively compress adjacent airways during expiration. As such the heterogeneity of airflow may contribute to a wide distribution of time constants for gas emptying and thus asynchrony, typically occurring with hyperinflation. Heterogeneity in the ventilation may also occur following bronchoconstriction, but with limited hyperinflation. As methacholine challenge in the asthma group did not induce important asynchrony, one may hypothesize that a more homogenous reduction in flow was obtained. A last explanation may be found in increased airway collapse by reduced airway tethering with the loss of alveolar tissue. We did not observe any significant relationship between roundness and DL_{CO} or K_{CO}, as indicative lung function markers of emphysema. Unfortunately, in the absence of CT measures we were not able to study the relationship with emphysema in more depth. From a clinical point of view the strong association between the roundness of resistance loops and hyperinflation is very attractive. The reduction of static hyperinflation is a main target of several COPD treatments including lung volume reduction surgery, endoscopic valve displacement and bronchodilators. Bronchodilator responses are usually evaluated on the expiratory volumes of spirometry, although different studies suggest that their impact on static lung volumes may be much larger [30–32]. Therefore, the monitoring of hyperinflation during tidal breathing may become a promising evaluation tool for treatment responses in specific phenotypes.
Our data clearly indicate that bronchoconstriction due to airway muscle contraction and mucosal oedema in asthma is resulting in significant changes in input related parameters such as SSG, point X, b _{ 0 } and b _{ 1 }. These changes correspond to an increased volume shift in phase with the flow, which obviously result in a reduction of the slope. Although most of our asthma patients did not show opening, 4 subjects had significant increases in roundness after the bronchial challenge, which is not surprising as hyperinflation can occur during bronchoconstriction and heterogeneity in airflow during exhalation may occur with increased bronchoconstriction [33]. When using sG_{AW} as best correlate for the slope, we accept that volume shift and airflow have a linear relationship, which is often not the case. SSG, which is highly correlated to sG_{AW}, accounts for pressure multiplication to reach the level of exhaled flow but under steady state condition. Both SSG and sG_{AW} ignore much of the system’s dynamics, but as they are computed differently they may still be influenced by other factors in the flow volume shift relationship. Looking at the slope per se, our data suggest to use sG_{0.5} instead of or with sG_{AW} (computed as sG_{mid}) due to its potential to discriminate obstructive asthma from mild COPD. As sG_{0.5} practically represents the slope of the right hand side of the loop, it will be little influenced by resistive mechanisms that open the loop.
Previous studies with FOT have demonstrated that resistance measures can vary within and between breath as they are flow and volume dependent. In particular, bronchomotor challenge can alter the ventilation which translates into changes of airway resistance and obscures the true changes in airway patency. Moreover, it has been shown that in patients with COPD the degree of flow limitation is varying between breaths [12, 34]. In the transition period of mild disease the detection of flow limitation may therefore require the monitoring of several breathing loops. For FOT different mathematical models have been proposed that correct for loop and flow dependency during tidal breathing [12, 15]. Our approach with plethysmography is adjusting flows for pressure changes during the whole breathing cycle and is taking the nonlinearity of these dynamics into account. Indirectly, it also provides correction for betweenbreath differences as a representative breathing cycle of overlapping loops was carefully selected. Taken together, our data imply that airways are not just simple tubes, but flexible structures that may enlarge or compress in function of the manoeuver.
The new parameters that we have developed are potentially important as they are able to identify different phenotypes within the spectrum of obstructive airways diseases. Hence, it is only with interventions and longitudinal followup that we will be able to identify their exact clinical relevance. One weakness of the study is the smaller number of healthy subjects, since a larger number would secure confidence in defining the normality range for each parameter. Similarly, comparisons of postbronchodilator with prebronchodilator measures, particularly in asthma may reveal more subtle differences between asthma and healthy controls. Another weakness may be the lack of data from CT scans, as that would allow the correlation of newly developed parameters with CT parameters of hyperinflation, airtrapping, bronchial inflammation and emphysema [35–37].
Conclusion
We developed new parameters which are reflecting the dynamic relationships between alveolar pressures and flows. Some of them are potentially important as they are linked to specific phenotypes within the spectrum of obstructive airways diseases. Mechanistic and prospective followup studies are now needed to determine the true validity of these parameters in respiratory medicine.
References
 1.
Decramer M, Janssens W, Miravitlles M. Chronic obstructive pulmonary disease. Lancet. 2012;379:1341–51.
 2.
Decramer M, Selroos O. Asthma and COPD: differences and similarities. With special reference to the usefulness of budesonide/formoterol in a single inhaler (Symbicort) in both diseases. Int J Clin Pract. 2005;59:385–98.
 3.
O’Leary JPTA, Capote LR. The physiologic basis of surgery. 4th ed. 2008.
 4.
Bateman ED, Hurd SS, Barnes PJ, Bousquet J, Drazen JM, FitzGerald M, Gibson P, Ohta K, O’Byrne P, Pedersen SE, et al. Global strategy for asthma management and prevention: GINA executive summary. Eur Respir J. 2008;31:143–78.
 5.
Rabe KF, Hurd S, Anzueto A, Barnes PJ, Buist SA, Calverley P, Fukuchi Y, Jenkins C, RodriguezRoisin R, van Weel C, et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary. Am J Respir Crit Care Med. 2007;176:532–55.
 6.
Goldman M, Smith H, Ulmer W. Wholebody plethysmography. Eur Res Monograph. 2005;31:15.
 7.
Criee C, Sorichter S, Smith H, Kardos P, Merget R, Heise D, Berdel D, Köhler D, Magnussen H, Marek W. Body plethysmography–its principles and clinical use. Respir Med. 2011;105:959–71.
 8.
Topalovic M, Derom E, Osadnik CR, Troosters T, Decramer M, Janssens W, Belgian Pulmonary Function Study I. Airways resistance and specific conductance for the diagnosis of obstructive airways diseases. Respir Res. 2015;16:88.
 9.
Peslin R, Duvivier C, Malvestio P, Benis AR. Correction of thermal artifacts in plethysmographic airway resistance measurements. J Appl Physiol (1985). 1996;80:2198–203.
 10.
Peslin R, Duvivier C, Malvestio P, Benis AR, Polu JM. Frequency dependence of specific airway resistance in a commercialized plethysmograph. Eur Respir J. 1996;9:1747–50.
 11.
Peslin R, Gallina C, Rotger M. Methodological factors in the variability of lung volume and specific airway resistance measured by body plethysmography. Bull Eur Physiopathol Respir. 1987;23:323–7.
 12.
Dellaca RL, Duffy N, Pompilio PP, Aliverti A, Koulouris NG, Pedotti A, Calverley PM. Expiratory flow limitation detected by forced oscillation and negative expiratory pressure. Eur Respir J. 2007;29:363–74.
 13.
Oostveen E, MacLeod D, Lorino H, Farre R, Hantos Z, Desager K, Marchal F. The forced oscillation technique in clinical practice: methodology, recommendations and future developments. Eur Respir J. 2003;22:1026–41.
 14.
Peslin R, Duvivier C, Gallina C. Total respiratory input and transfer impedances in humans. J Appl Physiol (1985). 1985;59:492–501.
 15.
Tomalak W, Peslin R, Duvivier C. Variations in airways impedance during respiratory cycle derived from combined measurements of input and transfer impedances. Eur Respir J. 1998;12:1436–41.
 16.
Islam M, Ulmer W. Diagnostic value of ‘closing volume’ in comparison to ‘airway resistance/lung volume plot’. Respiration. 1974;31:449–58.
 17.
Matthys H, Orth U. Comparative measurements of airway resistance. Respiration. 1975;32:121–34.
 18.
Crapo RO, Casaburi R, Coates AL, Enright PL, Hankinson JL, Irvin CG, MacIntyre NR, McKay RT, Wanger JS, Anderson SD, et al. Guidelines for methacholine and exercise challenge testing1999. This official statement of the American Thoracic Society was adopted by the ATS Board of Directors, July 1999. Am J Respir Crit Care Med. 2000;161:309–29.
 19.
Vestbo J, Hurd SS, Agusti AG, Jones PW, Vogelmeier C, Anzueto A, Barnes PJ, Fabbri LM, Martinez FJ, Nishimura M, et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary. Am J Respir Crit Care Med. 2013;187:347–65.
 20.
Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, Crapo R, Enright P, van der Grinten CP, Gustafsson P, et al. Standardisation of spirometry. Eur Respir J. 2005;26:319–38.
 21.
Quanjer PH. Standardized lung function testing. Report working party. Bull Eur Physiopathol Respir. 1983;19 Suppl 5:1–95.
 22.
Quanjer PH, Tammeling GJ, Cotes JE, Pedersen OF, Peslin R, Yernault JC. Lung volumes and forced ventilatory flows. Work Group on Standardization of Respiratory Function Tests. European Community for Coal and Steel. Official position of the European Respiratory Society. Rev Mal Respir. 1994;11 Suppl 3:5–40.
 23.
Miller A, Thornton JC, Warshaw R, Anderson H, Teirstein AS, Selikoff IJ. Single breath diffusing capacity in a representative sample of the population of Michigan, a large industrial state: predicted values, lower limits of normal, and frequencies of abnormality by smoking history 1–3. Am Rev Respir Dis. 1983;127:270–7.
 24.
Ljung L. System Identification Toolbox for Use with {MATLAB}. 2007.
 25.
Ljung L. Prediction error estimation methods. Circ Syst Signal Processing. 2002;21:11–21.
 26.
Schüssler HW. A stability theorem for discrete systems. Acoustics, Speech and Signal Processing, IEEE Transactions on. 1976;24:87–9.
 27.
Gibson G. Pulmonary hyperinflation a clinical overview. Eur Respir J. 1996;9:2640–9.
 28.
O’Donnell DE, Webb KA, Neder JA. Lung hyperinflation in COPD: applying physiology to clinical practice. COPD Research and Practice. 2015;1:1.
 29.
Rubinsztajn R, Przybyłowski T, MaskeyWarzęchowska M, PaplińskaGoryca M, Karwat K, NejmanGryz P, Chazan R. Correlation between hyperinflation defined as an elevated RV/TLC ratio and body composition and cytokine profile in patients with chronic obstructive pulmonary disease. Pneumologia i Alergologia Polska. 2015;83:120–5.
 30.
Ferguson GT. Why does the lung hyperinflate? Proc Am Thorac Soc. 2006;3:176–9.
 31.
Newton MF, O’Donnell DE, Forkert L. Response of lung volumes to inhaled salbutamol in a large population of patients with severe hyperinflation. CHEST. 2002;121:1042–50.
 32.
O’Donnell D, Flüge T, Gerken F, Hamilton A, Webb K, Aguilaniu B, Make B, Magnussen H. Effects of tiotropium on lung hyperinflation, dyspnoea and exercise tolerance in COPD. Eur Respir J. 2004;23:832–40.
 33.
Lougheed MD, Fisher T, O’Donnell DE. Dynamic hyperinflation during bronchoconstriction in asthma: implications for symptom perception. CHEST. 2006;130:1072–81.
 34.
Dellaca RL, Santus P, Aliverti A, Stevenson N, Centanni S, Macklem PT, Pedotti A, Calverley PM. Detection of expiratory flow limitation in COPD using the forced oscillation technique. Eur Respir J. 2004;23:232–40.
 35.
Baroni RH, FellerKopman D, Nishino M, Hatabu H, Loring SH, Ernst A, Boiselle PM. Tracheobronchomalacia: comparison between endexpiratory and dynamic expiratory CT for evaluation of central airway collapse 1. Radiology. 2005;235:635–41.
 36.
Wagnetz U, Roberts HC, Chung T, Patsios D, Chapman KR, Paul NS. Dynamic airway evaluation with volume CT: initial experience. Can Assoc Radiol J. 2010;61:90–7.
 37.
Wielpütz MO, Eberhardt R, Puderbach M, Weinheimer O, Kauczor HU, Heussel CP. Simultaneous assessment of airway instability and respiratory dynamics with lowdose 4DCT in chronic obstructive pulmonary disease: a technical note. Respiration. 2014;87:294–300.
Acknowledgements
We thank Mrs. Kristien De Bent, Mr. Willem De Wit and Mr. Rafi Zonanashvili for the provided technical support and help with data collection.
Funding
This work was supported by an Astra Zeneca Chair 2013–2015. WJ is supported the Flemish Research Foundation (FWO).
Availability of data and materials
Data will not be shared.
Authors’ contributions
MT performed the data analysis, contributed to the study design and wrote the manuscript. VE contributed in the data analysis and critically reviewed the manuscript. GC, TT and JA critically reviewed the manuscript. WJ takes responsibility for the content of the manuscript, contributed to the study design, assisted in the data analysis and critically reviewed the manuscript. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicable.
Ethics approval and consent to participate
All patients provided informed consent for the use of lung function and clinical data and the protocol was approved by the ethics committee of the University hospital of Leuven, Belgium.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Bodyplethysmography
 Airway resistance
 Pulmonary function tests
 Chronic obstructive pulmonary disease
 Asthma