Risk and optimal policies in bandit experiments
Karun Adusumilli
TL;DR
The paper develops a diffusion-asymptotics framework for bandit experiments, proving that the minimal Bayes risk under $n^{-1/2}$-scaled rewards is characterized by a $2^{\textrm{nd}}$-order PDE and that this characterization extends from Gaussian to parametric and non-parametric reward models via a limit-of-experiments argument. It shows that, asymptotically, only two state variables per arm are needed, enabling practical computation of Bayes and minimax policies that systematically outperform Thompson sampling, with tuned MOSS offering near-minimax performance in several settings. The work provides existence/uniqueness results for the PDE, convergence of discrete schemes to the PDE solution, and constructive ways to implement near-optimal policies through piecewise-constant controls and Monte-Carlo or finite-difference computation. Empirically, the Bayes-optimal policies dominate standard bandit algorithms in both one- and two-armed illustrated examples, underscoring the practical impact of the PDE-based approach for adaptive experimentation and decision making. The framework also accommodates extensions such as discounting and best-arm identification, broadening its relevance to diverse sequential experimentation problems.
Abstract
We provide a decision theoretic analysis of bandit experiments under local asymptotics. Working within the framework of diffusion processes, we define suitable notions of asymptotic Bayes and minimax risk for these experiments. For normally distributed rewards, the minimal Bayes risk can be characterized as the solution to a second-order partial differential equation (PDE). Using a limit of experiments approach, we show that this PDE characterization also holds asymptotically under both parametric and non-parametric distributions of the rewards. The approach further describes the state variables it is asymptotically sufficient to restrict attention to, and thereby suggests a practical strategy for dimension reduction. The PDEs characterizing minimal Bayes risk can be solved efficiently using sparse matrix routines or Monte-Carlo methods. We derive the optimal Bayes and minimax policies from their numerical solutions. These optimal policies substantially dominate existing methods such as Thompson sampling; the risk of the latter is often twice as high.
