Elsevier

Infection, Disease & Health

Volume 21, Issue 4, December 2016, Pages 184-191
Infection, Disease & Health

Research
Parameter estimation of tuberculosis transmission model using Ensemble Kalman filter across Indian states and union territories

https://doi.org/10.1016/j.idh.2016.11.001Get rights and content

Highlights

  • Deterministic TB model is used to model the TB transmission in India.

  • Parameters of model are estimated using Ensemble Kalman filter for the period 2006–2011.

  • Infection rate is found to be highest in Manipur.

Abstract

Background

Tuberculosis (TB) is one of the main causes of mortality on the globe. Besides the full implementation of Revised National Tuberculosis Control Programme (RNTCP), TB continues to be a major public health problem in India.

Methods

In the present study, parameters of a TB model are estimated using Ensemble Kalman filter (EnKf) approach. Infection rate and fraction of smear positive cases of TB are estimated in context of India.

Results and Conclusions

Results reveal that the infection rate is highest in Manipur and the ratio of smear positive cases is highest in Pondicherry. The infection rate of TB in Manipur is found to be 2.57 per quarter for the period 2006–2011.

Introduction

Tuberculosis (TB) is a well known infectious disease caused by bacterium M. tuberculosis which generally spreads through air. In 2011, 2.0–2.5 million new TB cases were estimated in India out of the global annual incidence of 9.4 million cases [1], [2]. A large number of factors contribute to the spread of TB; such as high prevalence of HIV/AIDS and diabetes, poor hygiene, crowding, illiteracy and lack of awareness make the TB situation critical in Indian context. All these factors directly contribute to high infection rates among the population. In 1997, the Government of India, with the help of World Bank, initiated RNTCP based on the internationally recommended Directly Observed Treatment Short-course (DOTS) strategy [2], [3]. RNTCP is the largest TB control program in terms of treatment of patients with full nationwide coverage.

Mathematical models and statistical techniques play a significant role in understanding the transmission dynamics of TB. In simple deterministic model of infectious disease, the number of susceptible persons who are infected by an infectious individual per unit of time is proportional to the total number of susceptible persons. This proportional coefficient is defined as infection rate. Estimation of parameters of mathematical model, for instance, infection rate contributes to better quantify the spread of disease.

Generally, inference of these parameters is a difficult task because of poor compatibility between observed data and models. Simulations and epidemiological data have been used to estimate the key parameters of deterministic models. Different techniques have been introduced and applied to estimate the parameters of TB models. Approximate Bayesian computation approach has been used to estimate TB transmission rate parameters for United States [4]. A synchronisation based method has been implemented to infer the parameters such as treatment rate, disease induced mortality rate and infection rate of a TB model. In particular, the infection rate in the study is estimated to be 2.04 for the quarterly data during 2003–2007 for Cameroon [5]. Liu et al. estimated the reactivation and infection rate of a TB model for China by assuming these rates as sinusoidal functions and infection rate is estimated to be 2.23 person per month for the period 2005–2009 [6]. A qualitative analysis of the TB model for Nigeria has been performed to analyse the effect of DOTS strategy [7]. Mandal et al. estimated parameters such as infection rate and treatment rate using annual prevalence and incidence data of TB [8]. In particular, the infection rate of TB have been estimated to be 11.03 per year for India [8]. Mishra et al. integrated quarantine compartment into TB model, which incorporates the multidrug-resistant TB patients. The model is further analysed and simulated in using TB data of Jharkhand, India [9].

In the present paper, we use Ensemble Kalman filter (EnKf) approach to estimate the parameters of a deterministic model of TB. Kalman filter has been extensively used to infer the parameters of models of various infectious diseases [10], [11], [12], [13], [14]. Parameters of an HIV/AIDS model has been estimated using Kalman filter approach [10]. An extension of Kalman filter has been implemented to analyse the spatio-temporal behaviour of measles outbreak using count data for the period 1960–1970 in Landon [12]. Influenza data of different cities within the United States has been analysed and demonstrated that ensemble filters were found to be more accurate than other filters in predicting the peaks of the influenza [13].

Section snippets

Dynamic model

In this paper, we use a variation of SIR (Susceptible-Infected-Recovered) model defined as SLIS (Susceptible-Latent-Infected-Susceptible). There are three exclusive groups of individuals; namely, susceptible, S, latently infected, L (infected with M. tuberculosis but not infectious), and actively infected with M. tuberculosis, I (infected and infectious). The model does not take into account genetic and demographic heterogeneity. The following are the governing differential equations for the

Results

The SLIS-EnKf framework explained in Section Methods and data is implemented in MATLAB and results are post processed and plotted in R-statistical programming language with the “spplot” package. In Fig. 1, reported smear positive cases of TB are plotted for Manipur state for the time period 2006 to 2011. It is observed that there is a declining trend in the number of cases. A seasonal pattern is observed in the number of cases over the years, with peak infection rate occurring in the mid of the

Discussion

We have implemented EnKf in conjunction with a deterministic model of TB. There are many alternate assimilation approaches that may be tested further [18]. The parameter estimation framework presented here captures seasonality well in the data which could not be expected from standard-likelihood methods. The technique is also computationally inexpensive when compared to Monte Carlo methods.

India's DOTS programme is the fastest expanding programme, placing more than 100,000 patients on treatment

Conclusion

To make RNTCP program more effective a shift in strategies is needed for states with higher transmission rates as compared with states with higher fraction of smear positive cases. For example, states with higher transmission rate require more focus on educational and awareness programs to mitigate the transmission of TB bacteria from infected to susceptible. Human behaviour plays an important role in the spread of infectious agent; transmission of the disease can be reduced by increasing the

Ethics

No ethical permission is required to conduct the present study.

Authorship contribution

SA, AB and PL designed the study. PN and VP carried out the analysis. All authors critically reviewed the manuscript.

Conflicts of interest

The authors of this paper have no conflicts of interest to declare.

Funding sources

The author(s) received no financial support for the research, authorship, and/or publication of this article.

References (24)

  • B.K. Mishra et al.

    Mathematical model on pulmonary and multidrug-resistant tuberculosis with vaccination

    J Egypt Math Soc

    (2013)
  • V. Dukic et al.

    Tracking epidemics with Google flu trends data and a state-space SEIR model

    J Am Stat Assoc

    (2012 Dec 1)
  • Cited by (9)

    • Analyzing the effects of observation function selection in ensemble Kalman filtering for epidemic models

      2021, Mathematical Biosciences
      Citation Excerpt :

      Compared to other filtering algorithms, ensemble Kalman filters avoid particle degeneracy issues associated with resampling and have shown to be computationally feasible for high-dimensional problems [4,5]. While these methods are commonly used in applications to weather prediction [6–8] and guidance, navigation, and control [9–11], ensemble Kalman-type filters have recently been utilized for parameter estimation and forecast prediction in a variety of epidemiological studies [12–19]. In the setting of epidemiology, while significant effort has been put towards developing useful forward models for a variety of infectious diseases, often less emphasis is placed on the choice of the observation function representing the data [20].

    • Short-term forecasts of the COVID-19 pandemic: a study case of Cameroon

      2020, Chaos, Solitons and Fractals
      Citation Excerpt :

      With the complexity of many phenomena, estimating of both states and parameters for dynamical systems are of high interest. EnKf has dealt with this interrogation [28,29]. In order to understand and control COVID-19 within a human community, several researchers have proposed different mathematical models for the dynamics of the new coronavirus [30–38].

    View all citing articles on Scopus
    View full text