Improving epidemic size prediction through stable reconstruction of disease parameters by reduced iteratively regularized Gauss–Newton algorithm

Abstract

AbstractClassical compartmental epidemic models of infectious diseases track the dynamic transition of individuals between different epidemiological states or risk groups. Reliable quantification of various transmission pathways in these models is paramount for optimal resource allocation and successful design of public health intervention programs. However, with limited epidemiological data available in the case of an emerging disease, simple phenomenological models based on a smaller number of parameters can play an important role in our quest to make forward projections of possible outbreak scenarios. In this paper, we employ the generalized Richards model for stable numerical estimation of the epidemic size (defined as the total number of infections throughout the epidemic) and its turning point using case incidence data of the early epidemic growth phase. The minimization is carried out by what we call the Reduced Iteratively Regularized Gauss–Newton (RIRGN) algorithm, a problem-oriented numerical scheme that takes full advantage of the specific structure of the non-linear operator at hand. The convergence analysis of the RIRGN method is suggested and numerical simulations are conducted with real case incidence data for the 2014–15 Ebola epidemic in West Africa. We show that the proposed RIRGN provides a stable algorithm for early estimation of turning points using simple phenomenological models with limited data.