2014-Kazemnejad

JRHS 2014; 14(2): 115-121

Copyright© Journal of Research in Health Sciences

Global Epidemic Trend of Tuberculosis during 1990-2010: Using Segmented Regression Model

Anoushiravan Kazemnejad (PhD)a, Shahram Arsang Jang (MSc)b*, Firouz Amani (PhD)c, Alireza Omidi (MSc)d

a Department of Biostatistics, Faculty of Medical Sciences, Tarbiat Modares University, Tehran, Iran

b Department of Epidemiology and Biostatistics, Faculty of Health, Qom University of Medical Sciences, Qom, Iran

c Department of Biostatistics, Faculty of Health, Ardabil University of Medical Sciences, Ardabil, Iran

d Department of Public Health, Faculty of Health, Qom University of Medical Sciences, Qom, Iran

* Correspondence: Shahram Arsang Jang(MSc), E-mail: shahramarsang@gmail.com

Received: 18 November 2013, Revised: 18 January 2014, Accepted: 05 March 2014, Available online: 12 March 2014

Abstract

Background: Tuberculosis (TB) is a pandemic disease. It is the second leading cause of death from infectious diseases after human immunodeficiency virus (HIV) in the world.The main objective of this paper was to determine and compare the epidemiology of TB incidence rate and its trend changes during 1990-2010 in six WHO regions regarding age, gender and income levels.

Methods: The Average Annual Percent Change (AAPC) and Annual Percent Change (APC) of TB incidence, mortality, treatment-successes, case detection rates, as well as change points of trend was estimated using segmented regression model. The number of change points was selected by the permutation procedure based on likelihood ratio test.

Results: Two change points for global TB incidence rate trend with AAPC5years equaling -1.4 % was estimated, the maximum AAPC5years of six regions was attributed to the American region (-3.5%). AACP of TB treatment-successes rate for Eastern Mediterranean (+2.2), the Americas (+1.6), south East Asia (+.8) and Global (+1.1) were significant (P<0.05). Moreover AACP5years of TB case detection rate for South East Asia (+7.5), Eastern Mediterranean (+4.9), Africa (+2.8) and the Americas (+1.7) were significant (P<0.05). Globally, all of income categories had descending trend of TB incidence and mortality rate, except the upper-middle income level that had ascending incidence trend (AAPC=+0.7%).

Conclusions: Globally, TB incidence and mortality rates have downturn trend and TB treatment successes and detection rates have upward trend, but their changes rate are insufficient to reach the goal of TB stop strategy. The economic levels have effect on trend, with no clear pattern, so it seems necessary that evaluation TB control programs based on characteristics of countries for reach TB control goals.

Keywords: Tuberculosis, Segmented regression, Change points, Annual Percent Change

Introduction

Tuberculosis (TB) is an infectious bacterial disease caused by Mycobacterium tuberculosis. It is the second leading cause of death worldwide amongst communicable diseases, striking 8.8 million people and killing 1.4 million in 20101.

Characterizing the trend of diseases and its changes synchronically can have an important role for evaluating the success of disease control strategies, health development indicators and health planning. Segmented linear regression is also known as breakpoint linear regression- piecewise linear regression- used for analysis and trend description of continuous data from mortality or incidence rate2-3. Understanding points of changes and number of changes are important. Permutation procedure, Bayesian Information Criteria (BIC), Akaike Information Criterion (AIC), and Generalized Cross Validation (GCV) can be used to select the number of change points4. Moreover, we can estimate the Annual Percent Changes (APC) and Average Annual Percent Changes (AAPC) using these methods. AAPC or APC provide summery statistics of trend. The estimate of change points in segmented regression is based on maximum likelihood, least square and Bayesian analysis. Segmented regression models fit with Lermans grid search (LGS) or Hudsons continuous fitting algorithm (HCA) approaches5-7

Several authors have examined segmented regression. Kim et al. compared cancer trend rate between two groups with segmented regression model and estimated the P-value using permutation test8. Arsang et al used the BIC and Joinpoint methods in order to determine the epidemiology of tuberculosis in Iran during 2001-20089; Nemes et al have used this method for detecting the relationship between abnormally expressed genes due to aberrant DNA copy numbers and subsequent altered gene expression profiles10.  Wijlaars et al have done the same to assess the effect of supervision of committee on safety of medicines on antidepressant prescription rates changes11; Bouadma et al has used the above-mentioned model to detect the impact of prevention program on ventilator associated pneumonia12. Muggeo et al. have used this model for fitting the trend of fertility data13.

The aims of this study were (a) estimating the number, location or time of change points for TB incidence, mortality, SP (Smear Positive) treatment-successes and case detection rate from 1990 to 2010  globally and six WHO regions (211 countries) including Eastern Mediterranean, Africa, The Americas, Europe, Western Pacific, South East Asia; (b) Pair wise comparison of parallelism and coincidence of trend between regions with the economic, age and gender categories as independent variables; (c) Providing the summery statistics of trend for each of segments for the last 5 and 10 years. The present research by determines TB trend changes can help to detect effective factors in TB control, improve current planning polices attributed to TB and prioritization. 

Methods

In order to determine the TB trend and the  number of change points, the segmented regression model was used with maximum three numbers of change points for n=20, in which n is the number of data points (1990 - 2010) as an independent variable. The maximum number of change points depends on the number of data points4. TB incidence rate per 105 populations, mortality incidence per 105 populations, percentage of SP treatment-successes and case detection rate of six WHO regions were considered response variables. Logarithmic transformation on response variable was performed. Parallelism of two combinations was tested in order to find out whether the two (trend of TB among regions were) regression mean functions were parallel. To do so, the years, six WHO regions, country income levels (low, lower-middle and upper-middle income countries), age groups and gender were inserted as independent variables. We used World Banks criterion that categorizes countries based on gross national income (GNI) per capita. The groups are: low income, $1,035 or less; lower middle income, $1,036 - $4,085; upper middle income, $4,086 or more.  Least squares methods for estimating regression parameters and Grid Search methods to determine the best fit for each individual model were used6.

We performed all the statistical analyses using Joinpoint software,V3.5.1 which is available at American Cancer Center14. Permutation test and BIC were also used for selecting the known number of change points4. Our data on global trends in TB came from WHO15 collected by WHO regional offices.

Segmented Regression

Segmented linear regression is one of the linear regression models used for segmenting nonlinear regression models into linear segments and points among segments called change points. For each segment there is a different fi(X) function, describing segmented regression curve for r segments as:

        

In which τ indicates the change point and fi(x,βi) regression function of ith segment. Regarding  (xij, yij)  pair of the jth observations for ith groups (i=1, 2; j=1,…,ni),e.g. yij can be TB incidence rate at the jth years, xij years, for the regional or group i, then regression response function mean for each segment or group i equals,

Where ki is the unknown number of change points, (i=1,.. , k) is the unknown change points, (i=1,2) and are the regression parameters and = for >0 8.

The Lermans grid search was used on, …,that uses  the least square estimation method, carrying out a grid search over trend to map minimum value of residuals sum of square function6. Parallelism test was used to find out whether two segmented mean regression functions are parallel, and to see whether two segmented mean regression functions are identical, coincident test was used accordingly. Permutation test was used for selection of models with i change points against j change points (i<j), the model with i change point was selected if

In which RSS denotes residual sum of squares, cn is a critical value obtained under the null model with i change points (cn equal the p-value), and is a non-decreasing penalty function for obtaining BIC4.

F statistic was used for testing H0: there is K0 change points against H1: there are KM change points, in which M is the maximum number of change points that is determined, a priori [by the researcher]. If H0 is rejected, then H0: k=1 against H1: k=M will be tested, otherwise H0: k=0 against H1: k=M-1 will be tested instead. And this procedure will be continued until hypothesis k=i stands against k=i+1, 0<i<M and the number of change point denoted by is estimated.

APC and AAPC

Annual Percent change is a way for describing the trend of incidence rate over time, with the assumption that incidence rates at year y with steady rate change into last y-1 year. To estimate the APC for a series of data, the following regression model is used:

Log (Ry)=b0+b1y where log (Ry) is the natural log of the rate in year y. the APC from year y to year

In order to present the summary of trend and determine the interval of years, AAPC can be used, providing single value of trend for the last 5 and 10 years, even if there are changes throughout the trend. This value is the average weight of APC segmented regression in which the weight equals the length of APC interval. To estimate the AAPC following equation is used:

That bis are the slope coefficients for each segment in the desired range of years, and the wis are the length of each segment runs in the range of years. In the Joinpoint software, the AAPC confidence interval is based on the normal distribution, and the APC confidence interval is based on a t distribution. If an AAPC lies entirely within a single Joinpoint segment, the AAPC is equal to the APC for that   segment 4,8-16.

Results

During 1995- 2010 there were 26,815,677 SP case notifications (64.3% males and 35.6% females). The maximum TB incidence rate was attributed to African region in 2004 with 303.36 per 105 populations. The maximum mortality rate was attributed to Western Pacific region in 1997 with 44.09 per 105 populations. TB incidence and mortality rate (exclude AIDS) globally for the last segment have declined, and two change points for incidence and one change point for mortality rate were estimated throughout the trend, so that the estimated mean segmented regression function for global TB incidence and mortality rate are:

   (Incidence rate)

     (Mortality rate)

TB incidence rate trend in all of the WHO regions for the last segment have declined (AAPC<0). Table 1 shows number and locale of change point(s). The maximum estimated AAPC of TB incidence rate for the last 5 years is attributed to the Americas region (AAPC=-3.5) and the minimum value is attributed to Eastern Mediterranean region (AAPC=-.8). The AAPC result of TB incidence and mortality incidence rate for the last 5 and 10 years are shown in Table 2.

While the global TB incidence rate during 1990-1996 has declined (APC=-.4, P=0.008), it has been constant during 1996-2004 (ACP=0.2, P=0.086) and has declined for the last segment (ACP=-1.4, P<0.001).

The parallelism hypotheses for TB incidence as well as mortality trend between WHO regions were rejected (P>0.05). The results of APC and estimated trends for TB mortality rate are shown in Figure 1.

The number and locale of change points of TB incidence, mortality, case detection, and SP treatment-success rate for WHO regions are shown in Table 1.

During 1999-2010 the maximum average value of treatment-success was attributed to South East Asia (89.07%) and the maximum average value of case detection rate was attributed to Europe (70.4%). The APC of TB SP treatment-success rates relating to the last segment were for Eastern Mediterranean (ACP=+2.2, P<0.001), the Americas (APC=+1.6, P=0.007), South East Asia (APC=0.8, P=0.016) and global (APC=1.1, P<0.001). The APC of TB case detection rates relating to the last segment were for Eastern Mediterranean (ACP=+4.9, P<0.001), Africa (ACP=+2.8, P<0.001), the Americas (ACP=+1.7, P=0.001) and South East Asia (ACP=+7.5, P=0.003). The TB case detection rate trend between the Americas and Europe was parallel (P=0.180) and coincident (P=0.171). The parallelism and coincident hypotheses between other WHO regions were rejected (P<0.05).

Trend of The TB incidence as well as mortality rate were analyzed separately based on country income categories, gender and age groups. AAPC results of TB incidence rate and mortality incidence rate for the last five years of low, lower- middle and upper middle income are shown in Figure 2 and 3.

  Pair-wise comparison demonstrates that AAPC5years differences of TB incidence rate between three segments (income level) for all regions were significant at 0.05 level, with the exclusion of the Americas (P>0.1), the low and lower-middle income countries globally and the low and lower-middle income countries of Europe. The global TB incidence has changed in 2003 from ascending to downturn trend for low (APC=-1.3; 95% CIAPC: -1.7, -1) and lower-middle income countries (APC=-1.1;95% CIAPC: -1.3, -1) but only its upward slopes is reduced for upper middle-income countries (APC=0.69; 95% CIAPC: 0.3, 1.1).

The global TB mortality trends of income categories were downturn for all segments and their slopes are increased after each change point.

We rejected the hypothesis that AAPC5years of TB incidence rate is different from zero for Eastern Mediterranean with low, Africa with upper and lower-middle, Europe with low and lower-middle and South East Asia with lower-middle income countries. The AAPCs5year of other regions and categories were different from zero (P<0.05) (Figure 2).  

During 1990-2010, the average proportion of TB notified in males in all regions for all age groups, except 0-14 age groups, was 2.018 times higher than females. With the passage of time (1990 to 2010) notification rate trend gap between female and male has been higher than that at previous times. Maximum TB notified during 1990-2010 was attributed to 25-44 age groups in male and 15-34 age groups in females. Globally, the smallest and the largest AAPC5 years differences between males and females were 1.2 (for age 0-14 years) and .4 (for age 35-45 years), respectively. We rejected the hypothesis that the AAPC of TB notified between two gender groups of all age groups for all regions are different, except some age groups of Eastern Mediterranean. In this region, AAPC differences of TB notified between males and females for 0-14 age groups (-6.2%; 95% CI: -9.4, -2.3), 15-24 age groups (-3.2%;95% CI: -4.4, -1.3) and 25-40 age groups (-3.5;95% CI:-5.1, -1.1) were significant, with AAPC for females being higher than males (P<0.001).

The output of BIC and permutation test for selection between ordinary or segmented regression and selection number of change points have produced similar results.

Figure 1: Trend and Annual Percent change (APC) of tuberculosis (TB) mortality incidence rate among WHO regions duration 1990-2010 (* means significant)

Figure 2: Average annual percent change (AAPC) of tuberculosis (TB) incidence rate for last 5 years. (NS: non-significant)

Figure 3: Average annual percent change (AAPC) of TB mortality rate for last 5 years. (NS: non-significant)

Table 1: Number and local of change points in tuberculosis (TB) incidence, mortality, case detection, and smear positive (SP) treatment-success rate trend

Table 2: The Average Annual Percent Changes (AAPC) for tuberculosis (TB) incidence, mortality case detection and smear positive (SP) treatment success rate for last 5 and 10 years

Discussion

The analysis of trend clearly indicates that compared to global trend, incidence reduction rate for all of the regions except Eastern Mediterranean (AAPC>-1.4), and mortality rate reduction for all of the regions except Europe and Africa were higher than those of global trend (AAPC>-4). In addition, mortality rate changes and decline were faster than TB incidence rate.

It can be said that the maximum rate of case detection for the last 5 and 10 years are attributed to South East Asia and the maximum rate of SP treatment-successes for the last 5 and 10 years to Eastern Mediterranean and Western Pacific regions. Since 2006 for the Americas and Western Pacific, 1996 for South East Asia and 2007 for globally, case detection has reached the goal of 85% TB SP treatment-success rate. AAPC of case detection rate during 2005-2010 for Europe and Western Pacific, also SP treatment successes rate for Africa, Europe and Western Pacific are not significant (fixed). Therefore, it seems necessary that attention to and re-evaluation of TB control programs be done for regions with fixed trend that have not reached the TB control goals (Africa, Europe and Eastern Mediterranean).

According to the findings; 1- TB incidence and mortality incidence rate of the WHO regions for the last segment have descent trend with different rates. 2- As regards the fact that more countries (a total of 180) implemented DOTS strategies from 1995 to 200317, and that change points occurred in 2003 (fix trend to downward) globally and in 2004 (decreasing incidence rate) for all lowincome countries, it can be concluded that TB control programs have been successful 18-19, but there is still a long way for the strategies to reach the goals of TB controls. If TB control programs continue with this trend and rate, it is difficult to reach the goal of stop TB and it is hoped that by 2050 only Americas and Western Pacific regions will reach the goal of incidence rate below one per 105 populations.

Like similar studies19, uniform effects of TB control programs would not be detected in this analysis, because social, biological and TB control program implementation variables are attributed to trend. TB incidence and mortality rate decline in most country were not remarkable; Only the American region has enjoyed the rate of decline (AAPC>4%), although the decline rate in some of them is the same but with regard to TB incidence rate for this region it is almost difficult to reach the goal of TB control at the same time. Therefore, it is vital, in the first step, to consider regions with high incidence and low change rate, and reinforce the DOTS strategy.

Upward trend in TB incidence rate and downward trend of mortality rate globally for uppermiddle income countries can be the result of improving the TB detection case that masked decline in TB incidence20. Moreover, upward trend in TB incidence and mortality rate for upper-middle Eastern Mediterranean can be as a detection case after long delays (regarding the fact that AAPC+ of mortality is greater than incidence). Unlike previous studies19, fixed trend of incidence and upward mortality trend for low income countries of Europe can be the result of decline or fixity of TB treatment successes, increase in extra pulmonary TB and multi- drug resistant 20.

As in other studies18, 21, there were low TB incidence in 0-14 age groups and high incidence in 15-44 age groups. It may be that exposure with TB risk factors are higher in younger people, compared to other age groups23. The ratio of TB cases notified in male were higher than females (except 0-14 age groups), which are expected to result from biological, behavioral and sociocultural components differences in male and female susceptibility to M. tuberculosis infection or the development of TB disease23-24.

Then incidence it seems that TB control strategy could deal with economic situation, so that regions with low income sometime have higher rate of reducing the TB incidence or mortality. This can be regarded as one of the successes in controlling TB , but not enough of TB stop strategies, because there are factors such as immigration, nutrition, early case detection, and access to health care service that impact the TB control programs. So to remove TB, attention must be paid to all forms of TB and the treatment processes should be improved (particularly early case detection).

Permutation test and BIC results demonstrate that when the relation between response and explanatory variable is nonlinear and there is/are change point(s) throughout trend, segmented regression provides unbiased estimation rather than nonlinear regression or Poisson models 4, 9. Segmented regression can be useful to detect the number and location of changes and also for comparing the trend of disease. Besides, this method can be used for dose-response study to detect the maximum effectiveness dose of drug and compare the parallelism and trend between drugs or groups. In this study we demonstrated TB trend for WHO regions, so for more details about trend, it is indispensable to determine the epidemiology of TB and effective factors such as co-morbidities, migration, drug resistance and severity of TB program implementation on trend changes, can be studied inside the regions.

Conclusions

As a result, the trend analyses provided useful measurement for comparing the successes the TB control strategies over time among the groups. These findings demonstrated that during 1990-2010, globally TB incidence and mortality rates have downturn trend. In addition, TB treatment successes and detection rates have upward trend, but there changes rate are insufficient to reach the goal of TB stop strategy. The economic levels have effect on AAPC of TB incidence and mortality rates, with no clear pattern, so it seems necessary that evaluation TB control programs based on characteristics of countries for reach TB control goals. The regions that have fixed or negative AAPC for TB incidence rate trend and positive AAPC of TB mortality rate trend, and also young age groups (15-44 years) need to more attention. As results, the America is in better situation than other regions.

Acknowledgments

The authors would like to thank the Department of Biostatistics, Faculty of Medical Sciences, Tarbiat Modares University, Tehran, Iran and all the study participants for their tremendous cooperation and support.

Conflict of interest statement

The authors have nothing to declare.

Funding

This study was funded by Tarbiat Modares University.

References

  1. World Health Organization. Tuberculosis key facts. WHO Web Site [cited March 2012]; Available from:  http://www.who.int/mediacentre/factsheets/fs104/en/index.html.
  2. Betts MG, Forbes GJ, Diamond AW.       Thresholds in Songbird Occurrence in Relation to Landscape Structure. Conserv Biol. 2007;21(4):1046-1058.
  3. Bang H, Mazumdar M, Spence D. Tutorial in biostatistics: analyzing associations between total plasma homocysteine and B vitamins using optimal categorization and segmented regression. Neuroepidemiology. 2006;27(4):188-200.
  4. Kim HJ, Yu B, Feuer EJ. Selecting the number of change points in segmented line regression. Stat Sin. 2009;19(2):597-609.
  5. Hyune JK, Binbing Y, Feuer EJ. Inference in segmented line regression: a simulation study. J Stat Comput Simul. 2008;78(11):1087-1103.
  6. Lerman PM. Fitting segmented regression models by grid search. J R Stat Soc Ser C Appl Stat. 1980;29(1):77-84.
  7. Feder PI. On asymptotic distribution theory in segmented regression problems identified case. Ann Stat. 1975;3(1):49-83.
  8. Kim HJ, Fay MP, Yu B, Barrett MJ, Feuer EJ. Comparability of segmented line regression models. Biometrics. 2004;60(4):1005-1014.
  9. Arsang SH, Kazemnejad A, Amani F. Applying segmented regression model to analysis the trend of tuberculosis incidence rate in Iran between 1964 -2008. Iranian Journal of Epidemiology. 2011;7(3):6-12.
  10. Nemes S, Parris TZ, Danielsson A, Kannius-Janson M, Jonasson JM, Steineck G, et al. Segmented regression, a versatile tool to analyze mRNA levels in relation to DNA copy number aberrations. Genes Chromosomes Cancer. 2012;51(1):77-82.
  11. Wijlaars LPMM, Nazareth I, Petersen I. Trends in Depression and Antidepressant Prescribing in Children and Adolescents: A Cohort Study in The Health Improvement Network (THIN). PLoS ONE. 2012;7(3):e33181.
  12. Bouadma L, Deslandes E, Lolom I, Corre1 BL, Mourvillier B, Regnier B, et al. Long-term impact of a multifaceted prevention program on ventilator-associated pneumonia in a medical intensive care unit. Clin Infect Dis. 2010;51(10):1115-1122.
  13. Muggeo VMR, Attanasio M, Porcu M. A segmented regression model for event history data: an application to the fertility patterns in Italy. J Appl Stat. 2009;36(9-10):973-988.
  14. Statistical Methodology and Applications Branch and Data Modeling Branch [computer program]. Version version 3.5: Surveillance Research Program National Cancer Institute; April 2011.
  15. Walls T, Shingadia D. Global epidemiology of paediatric tuberculosis. J Infect. 2004;48(1):13-22.
  16. Gallant AR, Fuller WA. Fitting segmented polynomial regression models whose join points have to be estimated. J Am Stat Assoc. 1973;68(341):144-147.
  17. Ukwaja K, Alobu I, Ifebunandu N, Osakwe C, Igwenyi C. From DOTS to the stop TB strategy: DOTS coverage and trend of tuberculosis notification in Ebonyi, southeastern Nigeria, 1998-2009. Pan Afr Med J. 2011;9:12.
  18. Guthmann J, Antoine D, Fonteneau L, Che D, Lévy-Bruhl D. Assessing BCG vaccination coverage and incidence of paediatric tuberculosis following two major changes in BCG vaccination policy in France. Euro Surveill. 2011;16:12.
  19. Dye C, Lönnroth K, Jaramillo E, Williams B, Raviglione M. Trends in tuberculosis incidence and their determinants in 134 countries. Bull World Health Organ. 2009;87(9):683-691.
  20. Dara M, Dadu A, Kremer K, Zaleskis R, Kluge HHP. Epidemiology of tuberculosis in WHO European Region and public health response. Eur Spine J. 2012:1-7.
  21. Kapata N, ChandaKapata P, OGrady J, Schwank S, Bates M, Mukonka V, et al. Trends of Zambias tuberculosis burden over the past two decades. Tropical Medicine & International Health. 2011;16(11):1404-1409.
  22. Wangdi K, Gurung MR. The epidemiology of tuberculosis in Phuentsholing General Hospital: a six-year retrospective study. BMC Res Notes. 2012;5(1):311.
  23. Neyrolles O, Quintana-Murci L. Sexual inequality in tuberculosis. PLoS Med. 2009;6(12):e1000199.
  24. Lönnroth K, Jaramillo E, Williams BG, Dye C, Raviglione M. Drivers of tuberculosis epidemics: The role of risk factors and social determinants. Soc Sci Med. 2009;68(12):2240-2246.


JRHS Office:

School of Public Health, Hamadan University of Medical Sciences, Shaheed Fahmideh Ave. Hamadan, Islamic Republic of Iran

Postal code: 6517838695, PO box: 65175-4171

Tel: +98 81 38380292, Fax: +98 81 38380509

E-mail: jrhs@umsha.ac.ir