Dynamic decision-making under model uncertainty is central to many economic environments, yet existing bandit and reinforcement learning algorithms rely on the assumption of correct model specification. This paper studies the behavior and performance of one of the most commonly used Bayesian reinforcement learning algorithms, Thompson Sampling (TS), when the model class is misspecified. We first provide a complete dynamic classification of posterior evolution in a misspecified two-armed Gaussian bandit, identifying distinct regimes: correct model concentration, incorrect model concentration, and persistent belief mixing, characterized by the direction of statistical evidence and the model-action mapping. These regimes yield sharp predictions for limiting beliefs, action frequencies, and asymptotic regret. We then extend the analysis to a general finite model class and develop a unified stochastic stability framework that represents posterior evolution as a Markov process on the belief simplex. This approach characterizes two sufficient conditions to classify the ergodic and transient behaviors and provides inductive dimensional reductions of the posterior dynamics. Our results offer the first qualitative and geometric classification of TS under misspecification, bridging Bayesian learning with evolutionary dynamics, and also build the foundations of robust decision-making in structured bandits.
Figure 2: Comparison of Case 1 and Case 2 scenarios. In Case 1, both models agree arm A1 is optimal. In Case 2, models disagree: ν prefers A1 while γ prefers A2 . The three regions (I, II, III) for the true expected rewards are indicated on each action line. In Region I the true expected reward is such that model ν is closer to the true DGP in KL divergence ( Δi>0 ), in Region III model γ is closer ( Δi<0 ), and Region II is the neutral boundary ( Δi=0 ).
Case
Condition
Behavior
Limit Distribution
Agreement
Δ1>0
Submartingale
δ1
Agreement
Δ1=0
Martingale
Bernoulli(1/2)
Agreement
Δ1<0
Supermartingale
δ0
Self-Confirming
Δ1>0>Δ2
Two attractors
Bernoulli(p∗)
Uniform Dominance
Δ1,Δ2>0
Global attractor
δ1
Self-Defeating
Δ1<0<Δ2
Ergodic
Invariant measure on (0,1)
Table 1: Complete classification of Thompson Sampling behavior under misspecification in the two-arm case
Case
Prob. of
Prob. of
Average
Action 1
Action 2
Regret
Agreement ( ϕ(ν)=ϕ(γ)=1 )
1
0
0
Agreement ( ϕ(ν)=ϕ(γ)=2 )
0
1
g(1)−g(2)
Self-Confirming
Bernoulli (p∗)
Bernoulli (1−p∗)
{0 w.p. p∗g(1)−g(2) w.p. 1−p∗
Uniform Dominance
1
0
0
Self-Defeating
α∗
1−α∗
(1−α∗)[g(1)−g(2)]
Table 2: Limiting action frequencies and average regret of Thompson Sampling in the two-arm model under misspecification
Figure 3: Vector field visualizations. Left two columns show S-space; right two columns show simplex space. (a,b) Self-confirming equilibrium; (c,d) Self-defeating equilibrium. (e,f) Drift towards face interior; (g,h) Drift towards vertex.
Figure 4: Different equilibrium types in multi-model Thompson Sampling. In each panel, the arrows d(a1),d(a2),d(a3) represent the drift vectors associated with each model, which are the expected directions in which log-odds would move if that model were played with probability one. The shaded region is the convex hull conv{d(a1),d(a2),d(a3)} , representing all possible mean drift directions achievable by mixing over models. We list out 4 representative cases: (a) Interior equilibrium with self-confirming models. (b) Interior equilibrium with self-defeating models. (c) No equilibria: drift towards face interior. (d) No equilibria: drift towards vertex.
0∈/conv{d(aj)} for all subfaces, recursion ends at a single vertex em⋆
Point mass δem⋆
Uniform dominance
Face-Ergodic
Push to some face ΔI (either 0∈/conv or interior unstable); on ΔI : First/Second Ergodic Condition. Models j∈/I are uniformly dominated .
Unique invariant μI on ΔI
Self-defeating (on the face) + uniform dominance
Nested/Mixed Outcomes
General recursive reduction: irrespective of interior fixed point existence, repeated tests on faces either (i) yield ergodic faces or (ii) descend to vertices.
Table 4: Limiting action distributions and average expected regret across multi-model categories. Here pa(π)=∑jπj1{ϕ(θ(j))=a} , a⋆=argmaxag(a) , r(a)=g(a⋆)−g(a) is the expected per-period regret, V⋆ are terminal vertices, and F⋆ are ergodic terminal faces.
Figure 5: Agreement case simulation results.
Figure 6: Self-confirming case simulation results.
Figure 7: Uniform dominance case simulation results.