Original Article

 

 

 

 

 

Reference Values for Lung Volumes in an Iranian Population: Introducing a New Equation Model

 

Mohammad Golshan MD*, Babak Amra MD*, Foroogh Soltani*, Robert O. Crapo MD**

 

 

Authors’ affiliations: *Bamdad Respiratory Research Institute, Isfahan University of Medical Sciences, Isfahan, Iran, **University of Utah and LDS Hospital, Pulmonary Division, LDS Hospital, 8th Avenue and C Street, Salt Lake City, UT, 84143 USA.

•Corresponding author and reprints:  Mohammad Golshan MD,

Bamdad Respiratory Research Institute, Bamdad Building, 11 Sheikh Mofid St., Isfahan 8164763831, Iran.

Fax:+98-311-662-3967, Tel: +98-311-663-7177-8

E-mail: golshan@med.mui.ac.ir

 

Accepted for publication: 30 December 2008



Background: Measurement of lung volumes, especially residual volume and total lung capacity are essential for assessment of restrictive lung disorders. Information regarding normative prediction values for lung volumes as measured by body plethysmography is scarce, and plethysmographic parameters are believed to be poorly reproducible. In this study, we   report a comprehensive set of predictive equations for static lung volumes from a general population sample of urban Iranians as measured by body plethysmography.

Methods: Standardized measurements were carried out on 1487 healthy nonsmoking volunteers (845 men and 642 women), aged six to 85 years, living in Isfahan, Iran. Nonlinear multiple regression analysis was used to calculate prediction equations based on subjects' ages and heights for the subdivisions of lung volumes [total lung capacity, functional residual capacity,  residual volume, and residual volume/total lung capacity (%)], separately for the two genders. The derived equations were used to calculate prediction values for the subjects. The two sets of predicted and measured values were compared by paired sample t-test.

Results: Prediction equations based on a new nonlinear model, (α1×age+α2×agen+β× height+c) which best fitted our data are presented. The measured and predicted values closely resemble and there is no significant difference between the two sets.

Since increments in total lung capacity, functional residual capacity, and residual volume disclose air trapping within the lungs, their upper limits of normal are as important as the lower limits. So, we have presented both for the equations.

Conclusion: Despite the usual beliefs we found rather reproducible prediction equations with high coefficient of determination (r2) and low standard error of estimate (SEE) in Iranian population.

 

 Archives of Iranian Medicine, Volume 12, Number 3, 2009: 256 – 261.

 

Keywords: Body plethysmography · Iran · lung volume(s) · prediction equation(s) · reference

values


 

 

Introduction
 

S

tatic lung volume (LV) measurements provide valuable information about the overall lung function and performance that can be fundamental in categorizing and staging of pulmonary disorders.1 Although forced vital  capacity  (FVC;  the  amount of air forcefully exhaled after a full inspiration) and static vital capacity (VC) and its components can be readily measured by simple spirometry,1 residual volume (RV; the volume of air remaining in the lungs after maximal exhalation), by definition, cannot be directly extracted from such a measurement. Measurement of RV allows functional residual capacity (FRC; the amount of air in the lungs at the end-tidal  position)  and  total  lung capacity (TLC; the amount of air in the chest after a maximum inspiration) to be derived when combined with the appropriate subdivisions of VC.2 Unfortunately, these more complex parameters which have greater physiologic and clinical implications than simpler variables are not simply accessible.1–3 Therefore, while simple spirometry will continue to be the first choice of pulmonary function study, in order to compensate for its failure in determining RV and TLC, body plethysmography is used as a complementary method. TLC and its subdivisions are usually measured by body plethysmography and/or gas dilution or nitrogen washout.2 

Reference values for pulmonary function tests are rather complex and variable: there are several potential sources of variability, some of which are related to individual characteristics such as sex, age, body size, ethnicity, level of regular physical activity, circadian rhythms, and the others are attributed to environmental factors including socio-economical status, exposures, altitude, smoking history, and technical issues such as posture, instrumentation, and the technique.4–7 Although it seems that much of this variability can be explained by the high degree of multi-colinearity among the mentioned factors, it should be noted that up to 20% of the total variability among populations cannot be explained at all.4,7 In an attempt to reduce this variability and improve accuracy, the use of reference values from a geographically-related population has been strongly recommended.8 Ideally, the population used for developing prediction equations should be asymptomatic general population.4 Unfortunately, to the best of our knowledge, there is no such source of reference values for static LVs obtained from a sample of the general population in Iran. Furthermore, it is usually stressed that values measured by body plethysmography are not exactly reproducible which could reflect technical difficulties with the method.9 

The purpose of this communication is to report a comprehensive set of predictive equations for static LVs from a general population sample of urban Iranians as measured by body plethysmography.

 

Materials and Methods

 

The Institutional Review Board for Medical Ethics at Isfahan Medical School approved the research protocol. During a 16-month period from January 2004 through April 2005, every 20th patient attending a general medical clinic in Isfahan, Iran was invited to bring her/his nonsmoker family members (parents, spouse, and children) for medical evaluation including a meticulous medical history, physical examination made by a general practitioner, and spirometry supplemented by body plethysmography. The exclusion criteria were: respiratory complaints (e.g., chronic or recurrent cough, dyspnea, sputum production, wheezy breathing, etc.), smoking one hundred cigarettes or more during lifetime, any history of hospitalization for pulmonary illness, physical findings suggesting cardiopulmonary disease, and any chest or skeletal deformity. Subjects were included in the study if they did not meet any of the exclusion criteria. Subject’s heights were measured, with bare feet, standing against a wall (buttocks, back, and head against the wall) with their heads erect in the Frankfort horizontal plane. A carpenter's square was placed against the wall and the subject’s head, and the subject was asked to step away from the wall, and the height was measured to the nearest centimeter from the floor to the bottom of the square with a metal ruler attached to the wall (Figure 1).

 

Figure 1. Height changes in correlation with age of the subjects in two sexes.

 

 

Age was asked (Figure 2), and in most cases, insurance cards, or identity documents were checked to confirm the accuracy of the stated age. Weight was measured with empty pockets and no shoes on.

 

 

Figure 2. Age distribution of the study population.

 

 

Barometric pressures, measured daily by Isfahan Airport, showed a range of 632 to 635 mmHg. Room temperature was monitored using a Brooklyn National Institute of Standards and Technology (NIST) Centigrade thermometer and kept between 21 and 25°C. Spirometry and body plethysmography results were automatically corrected for body temperature pressure saturation (BTPS) conditions by the device software. LVs and subdivision indices were measured, using a body plethysmography machine (ZAN500 Body Box) according to the published standards and criteria.7 The device was calibrated daily by its internal syringe. Spirograms were repeated until three acceptable tests were obtained. Studies were considered acceptable if the largest and second largest values for FVC and FEV1 were within 200 mL of each other. If the first maneuvers were not satisfactory, further maneuvers were used until the reproducibility criteria were satisfied or a maximum of eight maneuvers was reached. After obtaining spirometric measurements, the subject was allowed to rest for five minutes before proceeding to the measurement of lung subdivisions. Briefly, the subject was taught the panting technique while breathing through the mouthpiece of the body plethysmograph. The door was then closed and the measurements, including RV, FRC, and TLC, were made. The results were corrected for BTPS automatically by the Body Box software. All of the measurements were made by a single trained technician. The output report of each subject was scrutinized qualitatively by a pulmonologist (BA) before inclusion for statistical analysis.

 

Data analysis

Data were analyzed using the SPSS for Windows version 13 (SPSS Inc., Chicago, IL). Multiple regression analysis was applied to each of the LV indices and for each sex, separately. The relationship between LVs and anthropometric variables was examined first. Various regression models including quadratic and power functions, log-transformed, and linear relationships were compared.

The usefulness of various predictors of LVs was then compared. Age and standing height were considered as the best potential predictor variables. The dependent LV variables were first regressed individually against the independent variables. Stepwise multiple regression analyses were then used to determine which combination of parameters would fit the model best. Predictor variables were retained in the regression model only if their addition significantly improved the explained variance of the dependent variable. Wherever possible the most parsimonious model was chosen.

Using the derived equations, predicted parameter estimates were calculated for all subjects. Predictions from the new formula were compared with our measured values and also compared with predicted values extracted from previously published prediction equations7 the values predicted by Quanjer et al (Figure 3).3 Paired sample t-test was used to compare measured values versus predicted ones. To adjust for multiple comparisons, the level of statistical significance was set at 0.01.

 

Figure 3. Comparisons between measured versus predicted TLC in two sexes.

 

 

Results

 

Of the originally invited 1,953 individuals, 1,728 accepted to take part in the study (response rate=88.47%). The main reasons for refusal were: lack of time, negative attitude towards being a medical “research subject”, or being far from the clinic. Of the remaining, 198 were excluded by research criteria, and 146 could not provide satisfactory panting maneuvers. The remaining 1,386 body plethysmography recordings were analyzed to derive prediction equations.

Anthropometric characteristics for both sexes are presented in Table 1; height values ranged from 98 to 190 cm (169.54±11.80) in males, and from 104 to 179 cm (157.02±4.92) in females.

 

Table 1. Anthropometric characteristics for both sexes.

 

Females 

Males

Number of subjects

642

845

Age range

(Mean±SD)

7 – 73

(29.37±15.56)

6 – 85

(31.55±16.27)

Height range (Mean±SD)

119 – 179

(157.48±8.17)

114 – 190

(169.70±11.37)

 

 

In multiple regression approach, age and height were best fitted in all final models.

A nonlinear quadratic model was best fit for all parameters with a general form of:

 

α1 × age + α2 × agen + β × height + c

 

The derived equations for various lung function parameters are presented in Table 2. Ninety Five percent confidence intervals have been used to determine lower limits (LLN) and upper limits of normal (ULN).

 

 

Table 2. Derived  prediction equations based on the new model: ( α1 * age  + α2 *age **n + β * height + c).

Parameters

Model summaries for males

SEE for A1

SEE for B

r 2

TLC

LLN

ULN

41.228 *a - 40.872* a **1.00+0.34 * h - 2.547

41.228 *a - 40.872* a **1.00+0.030 * h - 2.547

41.228 *a - 40.872* a **1.00+0.037 * h - 2.547

0.01

.002

0.825

FRC

LLN

ULN

19.494*a -19.258 * a ** 1.002 +0.0202 * h –2.08

19.494*a -19.258 * a ** 1.002 +0.020 * h – 2.08

19.494*a -19.258 * a ** 1.002 +0.024 * h – 2.08

.002

.001

0.863

RV

LLN

ULN

95.09* age  -95.05*age**1.00+0.011×height- 0.521

95.09* age  -95.05*age **1.00+0.010*height- 0.521

95.09* age  -95.05*age **1.00+0.012*height- 0.521

.000

.001

.699

FRC/TLC

LLN

ULN

16.04×a- 16.04*a **1.00 + 0.001* h+ 0.415

16.04×a- 16.05*a **1.00 + 0.001* h+ 0.415

16.04 *a- 16.03*a **1.00 + 0.001* h+ 0.415

0.001

0.000

0.203

RV/TLC ratio

LLN

ULN

11.503 *a –11.512 *a ** 1.00+ 0.000 * h+ 0.335

11.503* a –11.512 * a ** 1.00 – 5.44 ×10-7* h+ 0.335

11.503 * a – 11.512 * a ** 1.00 + 0.000*h + 0.335

.000

.000

.240

 

Model summaries for females

 

 

r2

TLC

LLN

ULN

-12.416*a-12.07 *a **1.006 +0.049 * h – 5.388

-12.416* a-12.07 * a **1.006 +0.047 * h – 5.388

-12.416 * a-12.07 * a **1.006 +0.052 * h – 5.388

0.005

0.002

0.929

FRC

LLN

ULN

4.265* a-3.916 * a ** 1.01 + 0.016 * h – 2.567

4.265* a-3.916* a ** 1.01 + 0.013 * h – 2.567

4.265*a-3.916 * a ** 1.01 + 0.020 * h – 2.567

0.002

0.002

0.768

RV

LLN

ULN

-12.64 * a+12.821* a **0.997 + 0.005* h  - 0.77

-12.64 * a+12.821 * a **0.997+ 0.003* h  - 0.77

-12.64 * a+12.821 * a **0.997+ 0.007* h  - 0.77

.009

0.001

0.520

FRC/TLC ratio

LLN

ULN

2.387* a – 2.35* a **  1.003 – 0.002 * h+0.651

2.387* a – 2.35* a **  1.003 – 0.003* h+0.651

2.387* a – 2.35* a **  1.003 – 0.002* h+0.651

0.019

0.000

0.615

RV/TLC ratio

LLN

ULN

-41.08 * a+ 41.10 * a ** 1.00 -0.002 *h  +0.485

-41.08 * a+ 41.10 * a ** 1.00 -0.002 * h  +0.485

-41.08* a+ 41.10 * a ** 1.00 -0.002 * h  +0.485

0.002

.0002

0.210

a=age in years; h=height in centimeters; r2=coefficient of determination; LLN=lower limit of normal; ULN=upper limit of normal; **=to the power of; SEE= standard error of estimate.

 

As Table 3 shows our predicted values were significantly higher than published predictions.

 

Table 3.  Group comparisons.

 

Measured values (A)

Predicted by new equations (B)

P value*

( B vs. A)

Predicted by Stock’s

published equations (C)

P value*

( B vs. D)

Males 

 

 

 

 

 

TLC (L)

7.62±1.36**

7.64±1.18

0.93

6.15±1.10

0.000

FRC (L)

4.10±0.63

5.45±1.13

0.81

3.05±0.53

0.000

RV(L)

2.09±0.33

2.09±0.30

0.94

1.64±0.44

0.000

Females

 

 

 

 

 

TLC (L)

5.08±0.75

5.11±0.57

0.06

4.44±0.61

0.000

FRC (L)

2.95±0.58

2.94±0.39

0.15

2.36±0.37

0.000

RV(L)

1.42±0.32

1.41±0.21

0.26

2.83±1.56

0.000

*Paired sample t- test, **mean±SD.

 

Discussion
 

In many instances, the use of a widely recommended set of predictive equations for pulmonary function testes induces error.9 These results demonstrate the necessity of considering, as in most the biologic variables, reference values for LVs which are obtained from a racial, ethnic, and geographically-related population.4 Also, reference values for pulmonary function tests should derive from studies employing standardized procedures and equipment.4,6,7

Reported predicted values for static LVs generated with body plethysmography, using either volume-constant or volume-displacement plethysmography, vary markedly.10 Unfortunately, comparisons among studies are hindered because of some of the following factors: a) Equations for all the subdivisions of static LVs measured in the same group of reference subjects are not provided; b) Subjects with different smoking categories are not analyzed separately; c) The size of the reference population is small; d) A bias in the selection of the sample is present; e) Either the regression model is complex or age is not considered as a covariate.10 The present study provides a comprehensive description of the static LV values for a healthy, randomly selected sample of the Iranian  population.

Age has been consistently related to higher FRC and RV, independent of smoking.1,3 The aging process is associated with loss of lung elastance and increase in occlusion volume, both contributing to increased RV, FRC, and RV/TLC ratio. Since age-related changes in LVs are rather complex with an early increment, followed by a trivial flattening, and later decreasing slope of the curve, we used a new regression equation. We assumed that adding an extra age term helps to improve the prediction accuracy. From this one expects that the sign of α1 might be positive, whereas α2 will be negative, which means that, as the age increases, age contributes progressively less to the predicted value. But despite the fact that the derived equations work perfectly to predict the LV parameters, the signs for α1, α2, and n were not always the same as we had assumed. Actually various combinations of parameter signs reflect more complex correlation between the various independent variables with lung parameters.

In opposition to many previous studies,3,6,7,11–14 we obtained a rather high coefficient of determination (r2) and low standard error of estimate (SEE) in the LV prediction; however, still a substantial percentage of the LV variance could not be explained by the analyzed variables. These results are not consistent with the traditional notion that “LVs are far less reproducible than spirometric variables”, underlining the need for accurate and representative reference values. Additionally, the use of an adequate lower limit of normality is crucial.

Since increments in some parameters including RV and TLC are usually interpreted as findings in favor of air trapping,15 having ULN at hand, might be more important than their lower limits; in this study we are presenting the ULN to cope with this necessity. As much as we could find, this is the first time that ULN is introduced to the literature.

In summary, we have presented what we believe to be the first set of equations for prediction of static LVs in a random sample of an ethnically heterogeneous population from Iran and the Middle-East. Assuming that these values are at significant and nonparallel variance with a widely recommended set of equations from a North-American-European consensus,7 the use of equations obtained from foreign subjects with an "adjusting factor" is not advisable. Our results should ideally be applied to clinical and research contexts to evaluate the normalcy of static LV values in subjects aged three to 80 years with heights ranging 155 to 185 cm in males and 145 to 175 cm in females. The accuracy and validity of these equations, however, should be further confirmed in other samples of the Asian population with different ethnic and geographic backgrounds.

 

Acknowledgment

 

This work has been granted in part by Bamdad Respiratory Research Institute, Isfahan University of Medical Sciences, and ZAN Co. GMBH, Iranian Headquarter.

The authors cordially acknowledge the partial grants from the above institutions.

 

References

 

1           Ries AL. Measurements of lung volumes. Clin Chest Med. 1989; 10: 177 – 186.

2           Ruppel G. Lung volume tests. In: Ruppel G, ed. Manual of Pulmonary Function Testing. 6th ed. St. Louis: Mosby; 1994: 1 – 25.

3           Quanjer PH, Tammeling GJ, Cotes JE, Pedersen OF, Peslin R, Yernault JC. Lung volumes and forced ventilatory flows. Report Working Party Standardization of Lung Function Tests, European Community for Steel
and Coal. Official Statement of the European Respiratory Society.
Eur Respir J Suppl. 1993; 16: 5 – 40.

4           American Thoracic Society. Lung function testing. Selection of reference values and interpretative strategies. Am Rev Respir Dis. 1991; 144: 1202 – 1218.

5           Clausen JL. Prediction of normal values in pulmonary function testing. Clin Chest Med. 1989; 10: 135 – 143.

6           Quanjer PH, Dalhuijsen A, van Zomeren BC. Report working party: standardization of lung function tests. Bull Eur Soc Physiopathol Respir. 1983; 19 (suppl 5): S1 – S95.

7           Stocks J, Quanjer PH. Reference values for residual volume, functional residual capacity and total lung capacity. ATS Workshop on Lung Volume Measurements. Official Statement of The European Respiratory Society. Eur Respir J. 1995; 8: 492 – 506.

8           Golshan M, Nematbakhsh M, Amra B, Crapo RO. Spirometric reference values in a large Middle- Eastern population. Eur Respir J. 2003; 22: 529 – 534.

9           Wanger J, Clausen JL, Coates A, Pedersen OF, Brusasco V, Burgos F, et al. Standardization of the measurement of lung volumes. Eur Respir J. 2005; 26: 511 – 522.

10        Roca J, Burgos F, Barberà JA, Sunyer J, Rodriguez-Roisin R, Castellsagué J, et al. Prediction equations for plethysmographic lung volumes. Respir Med. 1998; 92: 454 – 460.

11        Donnelly PM, Yang TS, Peat JK, Woolcock AJ. What factors explain racial differences in lung volumes? Eur Respir J. 1991; 4: 829 – 838.

12        Boren HG, Kory RC, Syner JC. The Veteran's Administration-Army cooperative study of pulmonary function: the lung volume and its subdivisions in normal man. Am J Med. 1966; 41: 96 – 114.

13        Crapo RO, Morris AH, Clayton PD, Nixon CR. Lung volumes in healthy nonsmoking adults. Bull Eur Physiopathol Respir. 1982; 18: 419 – 425.

14        Ray CS, Sue DY, Bray G, Hansen JE, Wasserman K. Effects of obesity on respiratory function. Am Rev Respir Dis. 1983; 128: 501 – 506.

15        Pellegrino R, Viegi G, Brusasco V, Crapo RO, Burgos F, Casaburi R, et al. Interpretative strategies for lung function tests. Eur Respir J. 2005; 26: 948 – 968.

 


AIM Home | Table of Contents