State estimation for nonlinear state-space transmission models of tuberculosis

Abstract

Given the high prevalence of tuberculosis (TB) and the mortality rate associated with the disease, numerous models, such as the Gammaitoni and Nucci (GN) model, were developed to model the risk of transmission. These models typically rely on a quanta generation rate as a measurement of infectivity. Since the quanta generation rate cannot be measured directly, the unique contribution of this work is to develop state estimators to estimate the quanta generation rate from available measurements. To estimate the quanta generation rate, the GN model is adapted into an augmented single-room GN model and a simplified two-room GN model. Both models are shown to be observable, i.e., it is theoretically possible to estimate the quanta generation rate given available measurements. Kalman filters are used to estimate the quanta generation rate. First, a continuous-time extended Kalman filter is used for both adapted models using a simulation and measurement sampling rate of 60 s. Accurate quanta generate rate estimates are achieved in both cases. A more realistic scenario is also considered with a measurement sampling rate of one day. For these estimates, a hybrid extended Kalman filter (HEKF) is used. Accurate quanta generation rate estimates are achieved for the more realistic scenario. Future work could potentially use the HEKFs, the adapted models, and real-time measurements in a control system feedback loop to reduce the transmission of TB in confined spaces such as hospitals.