BanditAlgorithm.least_squares_confidence_ellipsoid
Proved(Confidence ellipsoid; L&S Theorem 20.5) With the setup above and rewards , let
be the -regularized least-squares estimator (). For , with probability at least , simultaneously for all :
where .
Preamble
import Mathlib.Probability.Moments.SubGaussian import Mathlib.Probability.Martingale.Basic import Mathlib.Analysis.SpecialFunctions.Log.Basic import Mathlib.Data.Real.Sqrt import Definitions.Def_SelfNormalizedProcess open MeasureTheory ProbabilityTheory Matrix
Formal statement
theorem BanditAlgorithm.least_squares_confidence_ellipsoid
{Ω : Type} {mΩ : MeasurableSpace Ω} [StandardBorelSpace Ω]
{P : Measure Ω} [IsProbabilityMeasure P]
{d : ℕ} (ℱ : Filtration ℕ mΩ)
(A : ℕ → Ω → Fin d → ℝ) (η : ℕ → Ω → ℝ) (X : ℕ → Ω → ℝ) (θs : Fin d → ℝ)
(hA : ∀ t : ℕ, Measurable[ℱ t] (A (t + 1)))
(hη : ∀ t : ℕ, Measurable[ℱ (t + 1)] (η (t + 1)))
(hsg : ∀ t : ℕ, HasCondSubgaussianMGF (ℱ t) (ℱ.le t) (η (t + 1)) 1 P)
(hX : ∀ (t : ℕ) (ω : Ω), X (t + 1) ω = θs ⬝ᵥ A (t + 1) ω + η (t + 1) ω)
{lam : ℝ} (hlam : 0 < lam) {δ : ℝ} (hδ : δ ∈ Set.Ioo (0 : ℝ) 1) :
1 - δ ≤ P.real {ω | ∀ t : ℕ,
Real.sqrt ((regularizedLeastSquares d lam A X t ω - θs) ⬝ᵥ
regularizedDesignMatrix d lam A t ω *ᵥ
(regularizedLeastSquares d lam A X t ω - θs))
< Real.sqrt lam * Real.sqrt (θs ⬝ᵥ θs)
+ Real.sqrt (2 * Real.log (1 / δ)
+ Real.log ((regularizedDesignMatrix d lam A t ω).det / lam ^ d))} := by
sorry
Source
L&S Theorem 20.5, p.260