Thompson sampling: Precise arm-pull dynamics and adaptive inference

TL;DR

This paper investigates the detailed dynamics of arm-pull counts in Thompson sampling algorithms, revealing a dichotomy in stability and enabling new inference methods for both stable and unstable arms.

Contribution

It uncovers the precise asymptotic behavior of arm-pull counts in Thompson sampling, contrasting with UCB algorithms, and introduces novel approaches for inference in these regimes.

Findings

Stable arms have Gaussian limits for normalized means.

Unstable arms exhibit non-Gaussian limits, enabling inference beyond stability.

A unifying principle links arm stability to the interaction with statistical noise.

Abstract

Adaptive sampling schemes are well known to create complex dependence that may invalidate conventional inference methods. A recent line of work shows that this need not be the case for UCB-type algorithms in multi-armed bandits. A central emerging theme is a `stability' property with asymptotically deterministic arm-pull counts in these algorithms, making inference as easy as in the i.i.d. setting. In this paper, we study the precise arm-pull dynamics in another canonical class of Thompson-sampling type algorithms. We show that the phenomenology is qualitatively different: the arm-pull count is asymptotically deterministic if and only if the arm is suboptimal or is the unique optimal arm; otherwise it converges in distribution to the unique invariant law of an SDE. This dichotomy uncovers a unifying principle behind many existing (in)stability results: an arm is stable if and only if its interaction with statistical noise is asymptotically negligible. As an application, we show that normalized arm means obey the same dichotomy, with Gaussian limits for stable arms and a semi-universal, non-Gaussian limit for unstable arms. This not only enables the construction of confidence intervals for the unknown mean rewards despite non-normality, but also reveals the potential of developing tractable inference procedures beyond the stable regime. The proofs rely on two new approaches. For suboptimal arms, we develop an `inverse process' approach that characterizes the inverse of the arm-pull count process via a Stieltjes integral. For optimal arms, we adopt a reparametrization of the arm-pull and noise processes that reduces the singularity in the natural SDE to proving the uniqueness of the invariant law of another SDE. We prove the latter by a set of analytic tools, including the parabolic H\"ormander condition and the Stroock-Varadhan support theorem.