Abstract
Mathematical modeling is a powerful tool for understanding and predicting complex dynamical systems, ranging from gene regulatory networks to population-level dynamics. However, model predictions are highly sensitive to initial conditions, which are often unknown. In infectious disease models, for instance, the initial number of exposed individuals (E) at the time the model simulation starts is frequently unknown. This initial condition has often been estimated using an unrealistic, history-independent assumption for simplicity: the chance that an exposed individual becomes infectious is the same regardless of the timing of their exposure (i.e., exposure history). Here, we show that this history-independent method can yield serious bias in the estimation of the initial condition. To address this, we developed a history-dependent initial condition estimation method derived from a master equation expressing the time-varying likelihood of becoming infectious during a latent period. Our method consistently outperformed the history-independent method across various scenarios, including those with measurement errors and abrupt shifts in epidemics, for example, due to vaccination. In particular, our method reduced estimation error by 55% compared to the previous method in real-world COVID-19 data from Seoul, Republic of Korea, which includes likely infection dates, allowing us to obtain the true initial condition. This advancement of initial condition estimation enhances the precision of epidemic modeling, ultimately supporting more effective public health policies. We also provide a user-friendly package, Hist-D, to facilitate the use of this history-dependent initial condition estimation method.