Which Spectrum? I. Kontoyiannis Athens U. of Econ & Business - - PowerPoint PPT Presentation

which spectrum
SMART_READER_LITE
LIVE PREVIEW

Which Spectrum? I. Kontoyiannis Athens U. of Econ & Business - - PowerPoint PPT Presentation

Which Spectrum? I. Kontoyiannis Athens U. of Econ & Business joint work with S.P. Meyn University of Illinois/Urbana-Champaign Athens Workshop on MCMC Convergence and Estimation 1 Outline Motivation In general: Geometric ergodicity


slide-1
SLIDE 1

Which Spectrum?

  • I. Kontoyiannis

Athens U. of Econ & Business joint work with S.P. Meyn University of Illinois/Urbana-Champaign Athens Workshop on MCMC Convergence and Estimation

1

slide-2
SLIDE 2

Outline

Motivation In general: Geometric ergodicity ⇔ spectral gap in LV

2

slide-3
SLIDE 3

Outline

Motivation In general: Geometric ergodicity ⇔ spectral gap in LV

Under reversibility: Geometric ergodicity ⇔ spectral gap in L2

3

slide-4
SLIDE 4

Outline

Motivation In general: Geometric ergodicity ⇔ spectral gap in LV

Under reversibility: Geometric ergodicity ⇔ spectral gap in L2 In the absence of reversibility Geometric ergodicity ⇐ spectral gap in L2 (explicit)

4

slide-5
SLIDE 5

Outline

Motivation In general: Geometric ergodicity ⇔ spectral gap in LV

Under reversibility: Geometric ergodicity ⇔ spectral gap in L2 In the absence of reversibility Geometric ergodicity ⇐ spectral gap in L2 (explicit) Geometric ergodicity ⇒ spectral gap in L2 (example)

5

slide-6
SLIDE 6

Outline

Motivation In general: Geometric ergodicity ⇔ spectral gap in LV

Under reversibility: Geometric ergodicity ⇔ spectral gap in L2 In the absence of reversibility Geometric ergodicity ⇐ spectral gap in L2 (explicit) Geometric ergodicity ⇒ spectral gap in L2 (example) Convergence rates Under reversibility: TV finite-n bound Without reversibility: Asymptotic V -norm bound

6

slide-7
SLIDE 7

The Setting

{Xn} Markov chain with general state space (Σ, S) X0 = x ∈ Σ initial state P(x, dy) transition kernel P(x, A) := P

x{X1 ∈ A} := Pr{Xn ∈ A|Xn−1 = x}

7

slide-8
SLIDE 8

The Setting

{Xn} Markov chain with general state space (Σ, S) X0 = x ∈ Σ initial state P(x, dy) transition kernel P(x, A) := P

x{X1 ∈ A} := Pr{Xn ∈ A|Xn−1 = x}

ψ-irreducibility and aperiodicity Assume that there exists σ-finite measure ψ on (Σ, S) such that P n(x, A) > 0 eventually for any x ∈ Σ and any A ∈ S with ψ(A) > 0

8

slide-9
SLIDE 9

The Setting

{Xn} Markov chain with general state space (Σ, S) X0 = x ∈ Σ initial state P(x, dy) transition kernel P(x, A) := P

x{X1 ∈ A} := Pr{Xn ∈ A|Xn−1 = x}

ψ-irreducibility and aperiodicity Assume that there exists σ-finite measure ψ on (Σ, S) such that P n(x, A) > 0 eventually for any x ∈ Σ and any A ∈ S with ψ(A) > 0 Recall Any kernel Q(x, dy) acts of functions F : Σ → R and measures µ on (Σ, S) as a linear operator: QF(x) =

  • Σ

Q(x, dy)F(y) µQ(A) =

  • Σ

µ(dx)Q(x, A)

9

slide-10
SLIDE 10

Geometric Ergodicity (GE) Equivalent Conditions ❀

There is an invariant measure π and functions ρ : Σ → (0, 1), C : Σ → [1, ∞): P n(x, ·) − πTV ≤ C(x)ρ(x)n n ≥ 0, π − a.s.

10

slide-11
SLIDE 11

Geometric Ergodicity (GE) Equivalent Conditions ❀

There is an invariant measure π and functions ρ : Σ → (0, 1), C : Σ → [1, ∞): P n(x, ·) − πTV ≤ C(x)ρ(x)n n ≥ 0, π − a.s.

There is an invariant measure π constants ρ ∈ (0, 1), B < ∞ and a π-a.s. finite V : Σ → [1, ∞]: P n(x, ·) − πV ≤ BV (x)ρn n ≥ 0, π − a.s. where FV := sup

x∈Σ

|F(x)| V (x) µV := sup

F : FV <∞

  • Fdµ
  • 11
slide-12
SLIDE 12

Geometric Ergodicity (GE) Equivalent Conditions ❀

There is an invariant measure π and functions ρ : Σ → (0, 1), C : Σ → [1, ∞): P n(x, ·) − πTV ≤ C(x)ρ(x)n n ≥ 0, π − a.s.

There is an invariant measure π constants ρ ∈ (0, 1), B < ∞ and a π-a.s. finite V : Σ → [1, ∞]: P n(x, ·) − πV ≤ BV (x)ρn n ≥ 0, π − a.s. where FV := sup

x∈Σ

|F(x)| V (x) µV := sup

F : FV <∞

  • Fdµ

Lyapunov condition (V4) There exist V : Σ → [1, ∞), δ > 0, b < ∞ and a “small” C ⊂ Σ: P V (x) ≤ (1 − δ)V (x) + bIC

12

slide-13
SLIDE 13

GE and the LV

∞ Spectrum Proposition 1: Geometric ergodicity ⇔ spectral gap in LV

[∼K-Meyn 2003] Suppose the chain {Xn} is ψ-irreducible and aperiodic. Then it is GE iff P admits a spectral gap in LV

∞ := {F : Σ → R s.t. FV < ∞}

13

slide-14
SLIDE 14

GE and the LV

∞ Spectrum Proposition 1: Geometric ergodicity ⇔ spectral gap in LV

[∼K-Meyn 2003] Suppose the chain {Xn} is ψ-irreducible and aperiodic. Then it is GE iff P admits a spectral gap in LV

∞ := {F : Σ → R s.t. FV < ∞}

Recall A set C ⊂ Σ is small if there exist n ≥ 1, ǫ > 0 and a probability measure ν

  • n (Σ, S) such that P n(x, A) ≥ ǫI

IC(x)ν(A) for all x ∈ Σ, A ∈ S The spectrum S(P) of P : LV

∞ → LV ∞ is the set of λ ∈ C s.t.

(I − λP)−1 : LV

∞ → LV ∞ does not exist

P : LV

∞ → LV ∞ admits a spectral gap if S(P) ∩ {z ∈ C : |z| ≥ 1 − ǫ}

contains only poles of finite multiplicity for some ǫ > 0

14

slide-15
SLIDE 15

Proof ideas (⇒)

Consider the potential operator Uz := [Iz − (P − IC ⊗ ν)]−1, z ∈ C Iterating the contraction provided (V4) gives a bound on | | |Uz| | |V for z ∼ 1 Use Uz to check that f0 ≡ 1 is an eigenfunction corresponding to λ0 = 1 Using an operator-inversion formula a la Nummelin [Iz − P]−1 = [Iz − (P − IC ⊗ ν)]−1

  • I +

1 1 − κIC ⊗ ν

  • show λ = 1 is maximal, isolated, and non-repeated

κ = ν[Iz − (P − IC ⊗ ν)]−1IC ✷

15

slide-16
SLIDE 16

GE, Reversibility and the L2(π) Spectrum

Proposition 2: Under reversibility: GE ⇔ spectral gap in L2 [Roberts-Rosenthal 1997] [Roberts-Tweedie 2001] [K-Meyn 2003] Suppose the chain {Xn} is reversible, ψ-irreducible and aperiodic. Then it is GE iff P admits a spectral gap in L2(π)

16

slide-17
SLIDE 17

GE, Reversibility and the L2(π) Spectrum

Proposition 2: Under reversibility: GE ⇔ spectral gap in L2 [Roberts-Rosenthal 1997] [Roberts-Tweedie 2001] [K-Meyn 2003] Suppose the chain {Xn} is reversible, ψ-irreducible and aperiodic. Then it is GE iff P admits a spectral gap in L2(π) Proof Analogous definitions, proof outline similar to Proposition 1 Big difference: In the Hilbert space setting, the spectral gap is simply 1 − sup νP2 ν2 : ν s.t. ν(Σ) = 1, ν2 = 0

  • where ν2 := dν/dπ2

17

slide-18
SLIDE 18

L2 Spectral Gap Always Implies GE

Theorem 1 Suppose the chain {Xn} is ψ-irreducible and aperiodic and that P admits a spectral gap in L2 Then the chain is geometrically ergodic [w.r.t. so some Lyapunov function V ]

18

slide-19
SLIDE 19

L2 Spectral Gap Always Implies GE

Theorem 1 Suppose the chain {Xn} is ψ-irreducible and aperiodic and that P admits a spectral gap in L2 Then for any h ∈ L2(π) there is a Vh ∈ L1(π) s.t.:

❀ (V4) holds w.r.t. Vh ❀ h ∈ LVh

19

slide-20
SLIDE 20

L2 Spectral Gap Always Implies GE

Theorem 1 Suppose the chain {Xn} is ψ-irreducible and aperiodic and that P admits a spectral gap in L2 Then for any h ∈ L2(π) there is a Vh ∈ L1(π) s.t.:

❀ (V4) holds w.r.t. Vh ❀ h ∈ LVh

Proof Prove “soft” GE Get explicit exponential bounds on explicit Kendall sets Let Vh(x) := Ex σC

  • n=0
  • 1 + |h(X(x))|
  • exp{1

2θn}

20

slide-21
SLIDE 21

L2 Spectral Gap ⇐ GE

Theorem 2 There is a (non-reversible) ψ-irreducible and aperiodic chain {Xn}

  • n a countable state space Σ, which is geometrically ergodic

but its transition kernel P does not admit a spectral gap in L2

21

slide-22
SLIDE 22

L2 Spectral Gap ⇐ GE

Theorem 2 There is a (non-reversible) ψ-irreducible and aperiodic chain {Xn}

  • n a countable state space Σ, which is geometrically ergodic

but its transition kernel P does not admit a spectral gap in L2 Proof Start with an example of H¨ aggstr¨

  • m or of Bradley:

GE chain {Xn} but CLT fails for some G ∈ L2 Spectral gap exists ⇒ autocorrelation function of {G(Xn)} decays exponentially ⇒ normalized partial sums of {G(Xn)} bdd in L2 ⇒ CLT ⇒ contradiction ✷

22

slide-23
SLIDE 23

Why do we care?

Theorem 3. [Roberts-Rosenthal 1997] Suppose the chain {Xn} is reversible, ψ-irreducible and aperiodic If P admits a spectral gap δ2 > 0 in L2 Then for any X0 ∼ µ: µP n − πTV ≤ µ − π2(1 − δ2)n

23

slide-24
SLIDE 24

Why do we care?

Theorem 3. [Roberts-Rosenthal 1997] Suppose the chain {Xn} is reversible, ψ-irreducible and aperiodic If P admits a spectral gap δ2 > 0 in L2 Then for any X0 ∼ µ: µP n − πTV ≤ µ − π2(1 − δ2)n Theorem 4. Suppose the chain {Xn} is ψ-irreducible and aperiodic If P admits a spectral gap δV > 0 in LV

Then for π-a.e. x: lim

n→∞

1 n log P n(x, ·) − πV = log(1 − δV ) In fact: lim

n→∞

1 n log

  • sup

x∈X, FV =1

|P nF(x) −

  • F dπ|

V (x)

  • = log(1 − δV )

24