Stochastic Optimal Control of an Epidemic Under Partial Information
Ibrahim Mbouandi Njiasse, Florent Ouabo Kamkumo, Ralf Wunderlich
TL;DR
The paper tackles a stochastic optimal control problem for an epidemic under partial information, modeling with an $SI^{\pm}R^{\pm}H$ compartmental structure and incorporating social distancing, testing, and vaccination. It converts the partially observed problem into a fully observed Markov decision process via diffusion approximation and an extended Kalman filter, then solves the resulting Bellman equation with two scalable numerical schemes: (i) state discretization with optimal quantization and linear interpolation, and (ii) quantization with value-function regression using educated basis functions. Numerical experiments show that the optimal policy typically combines lockdown, testing, and vaccination to curb infections while respecting hospital capacity, with comparable results between the two methods and clear insights into how undetected infections drive interventions. The work integrates filtering, dynamic programming, and quantization-based techniques to provide practical decision-support for epidemic management under uncertainty.**
Abstract
In this paper, we address a social planner's optimal control problem for a partially observable stochastic epidemic model. The control measures include social distancing, testing, and vaccination. Using a diffusion approximation for the state dynamics of the epidemic, we apply filtering arguments to transform the partially observable stochastic optimal control problem into an optimal control problem with complete information. This transformed problem is treated as a Markov decision process. The associated Bellman equation is solved numerically using optimal quantization methods for approximating the expectations involved to mitigate the curse of dimensionality. We implement two approaches, the first involves state discretization coupled with linear interpolation of the value function at non-grid points. The second utilizes a parametrization of the value function with educated ansatz functions. Extensive numerical experiments are presented to demonstrate the efficacy of both methods.
