14

Марковские цепи и MCMC

[26/0%]
Показать
LaTeX
Задача 14.1
?
(a)

Покажите, используя определение 14.1.1, что когда пространство состояний S\mathbb {S} счётно, для любого nn, при условии {Xn=an}\left\{ X_{n} = a_{n}\right\}, события {Xn+j=an+j,1≤j≤k}\left\{ X_{n+j} = a_{n+j}, 1 \leq j \leq k\right\} и {Xj=aj:0≤j≤n−1}\left\{ X_{j} = a_{j}: 0 \leq j \leq n-1\right\} независимы для всех выборов kk и {aj}j=0n+k\left\{ a_{j}\right\}_{j = 0}^{n+k}. Таким образом, при условии "настоящего" {Xn=an}\left\{ X_{n} = a_{n}\right\}, "прошлое" {Xj:j≤n−1}\left\{ X_{j}: j \leq n-1\right\} и "будущее" {Xj:j≥n+1}\left\{ X_{j}: j \geq n+1\right\} являются двумя семействами независимых случайных величин относительно условной вероятностной меры P(⋅∣Xn=an)P\left(\cdot \mid X_{n} = a_{n}\right), при условии P(Xn=an)>0P\left(X_{n} = a_{n}\right) > 0.

(b)

Докажите предложение 14.2.2, используя индукцию по nn (ср. главу 6).

Задача 14.2

В примере 14.1.1 (Лягушка в колодце) проверьте, что

?
(a)

если αi≡1−1ci,c>1,i≥1\alpha_{i} \equiv 1-\frac{1}{c i}, c > 1, i \geq 1, то состояние 1 нуль-возвратно,

(b)

если αi≡α,0<α<1\alpha_{i} \equiv \alpha , 0 < \alpha < 1, то состояние 1 положительно возвратно, и

(c)

если αi≡1−12i2\alpha_{i} \equiv 1-\frac{1}{2 i^{2}}, то состояние 1 невозвратно.

Задача 14.3

Рассмотрим простое симметричное случайное блуждание (ПССБ) в Z2\mathbb {Z}^{2}, где переходные вероятности p(i,j)(i′,j′)=14p_{(i, j)\left(i^{\prime }, j^{\prime }\right)} = \frac{1}{4} каждая, если (i′,j′)∈{(i+1,j),(i−1,j),(i,j+1),(i,j−1)}\left(i^{\prime }, j^{\prime }\right) \in \left\{ (i+1, j),(i-1, j),(i, j+1),(i, j-1)\right\}, и нулевые в противном случае. Проверьте, что при n=2kn = 2 k

p(0,0),(0,0)(2k)=142k(2kk)2∼1π1k p_{(0,0),(0,0)}^{(2 k)} = \frac{1}{4^{2 k}}\binom {2 k}{k}^{2} \sim \frac{1}{\pi } \frac{1}{k}

и заключите, что (0,0)(0,0) нуль-возвратно. Распространите это вычисление на ПССБ в Z3\mathbb {Z}^{3} и заключите, что (0,0,0)(0,0,0) невозвратно.

?
Задача 14.4

Покажите, что если ii поглощающее и j→ij \rightarrow i, то jj невозвратно, показав, что если j→ij \rightarrow i, то fji∗=P(Ti<Tj∣X0=j)>0f_{j i}^{*} = P\left(T_{i} < T_{j} \mid X_{0} = j\right) > 0 и 1−fjj≥fji∗1-f_{j j} \geq f_{j i}^{*}.

?
Задача 14.5
?
(a)

Пусть ii возвратно и i→ji \rightarrow j. Покажите, что jj возвратно, используя следствие 14.1.5. (Указание: покажите, что существуют n0n_{0} и m0m_{0} такие, что для всех nn, pjj(n0+n+m0)≥pji(n0)pii(n),pij(m0)p_{j j}^{\left(n_{0}+n+m_{0}\right)} \geq p_{j i}^{\left(n_{0}\right)} p_{i i}^{(n)}, p_{i j}^{\left(m_{0}\right)} с pji(n0)>0p_{j i}^{\left(n_{0}\right)} > 0 и pij(m0)>0p_{i j}^{\left(m_{0}\right)} > 0.)

(b)

Пусть ii и jj сообщаются. Покажите, что di=djd_{i} = d_{j}.

Задача 14.6

Покажите, что в конечной неприводимой цепи Маркова (S,\vectP)(\mathbb {S}, \vect {P}) все состояния положительно возвратны, показав

?
(a)

что для любых i,ji, j в S\mathbb {S} существует r,r≤Kr, r \leq K, такое, что pij(r)>0p_{i j}^{(r)} > 0, где KK — число состояний в S\mathbb {S},

(b)

для любого ii в S\mathbb {S} существует 0<α<10 < \alpha < 1 и c<∞c < \infty такие, что Pi(Ti>n)≤cαnP_{i}\left(T_{i} > n\right) \leq c \alpha^{n}.

Дайте альтернативное доказательство, показав, что если S\mathbb {S} конечно, то для любого начального распределения μ\mu меры пребывания

μn(⋅)≡1(n+1)∑j=0nPμ(Xj∈⋅) \mu _{n}^{(\cdot )} \equiv \frac{1}{(n+1)} \sum _{j = 0}^{n} P_{\mu }\left(X_{j} \in \cdot \right)

имеют подпоследовательность, сходящуюся к вероятностному распределению π\pi, стационарному для (S,\vectP)(\mathbb {S}, \vect {P}).

Задача 14.7

Докажите теорему 14.1.3, используя марковское свойство и индукцию.

?
Задача 14.8

Адаптируйте доказательство теоремы 14.1.9, чтобы показать, что для любых i,ji, j

1n∑j=1npij(k)→fijEjTj \frac{1}{n} \sum _{j = 1}^{n} p_{i j}^{(k)} \rightarrow \frac{f_{i j}}{E_{j} T_{j}}

если jj положительно возвратно, и 0 в противном случае. Заключите, что в случае конечного пространства состояний должно быть хотя бы одно положительно возвратное состояние.

?
Задача 14.9

Если j→ij \rightarrow i, то ζ1≡∑j=0Ti−1δXrj\zeta_{1} \equiv \sum_{j = 0}^{T_{i}-1} \delta_{X_{r} j} — число посещений jj до посещения ii — удовлетворяет Pi(ζ1>n)≤cαnP_{i}\left(\zeta_{1} > n\right) \leq c \alpha^{n} для некоторых 0<c<∞,0<α<10 < c < \infty , 0 < \alpha < 1 и всех n≥1n \geq 1.

?
Задача 14.10

Адаптируйте доказательство теоремы 14.1.9, чтобы установить следующие законы больших чисел. Пусть (S,\vectP)(\mathbb {S}, \vect {P}) неприводима и положительно возвратна со стационарным распределением π\pi.

?
(a)

Пусть h:S→Rh: \mathbb {S} \rightarrow \mathbb {R} такова, что ∑j∈S∣h(j)∣πj<∞\sum_{j \in \mathbb {S}}\left|h(j)\right| \pi_{j} < \infty. Тогда для любого начального распределения μ\mu,

1n+1∑j=0nh(Xj)→∑j∈Sh(j)πj п.н.  \frac{1}{n+1} \sum _{j = 0}^{n} h\left(X_{j}\right) \rightarrow \sum _{j \in \mathbb {S}} h(j) \pi _{j} \quad \text{ п.н. }

сначала проверив, что

Ei(∣∑j=0Ti−1h(Xj)∣)<∞ E_{i}\left(\left|\sum _{j = 0}^{T_{i}-1} h\left(X_{j}\right)\right|\right) < \infty
(b)

Пусть g:S×S→Rg: \mathbb {S} \times \mathbb {S} \rightarrow R такова, что ∑i,j∈S∣g(i,j)∣πipij<∞\sum_{i, j \in \mathbb {S}}\left|g(i, j)\right| \pi_{i} p_{i j} < \infty. Тогда для любого начального распределения μ\mu,

1n+1∑j=0ng(Xj,Xj+1)→∑i,j∈Sg(i,j)πipij п.н.  \frac{1}{n+1} \sum _{j = 0}^{n} g\left(X_{j}, X_{j+1}\right) \rightarrow \sum _{i, j \in \mathbb {S}} g(i, j) \pi _{i} p_{i j} \quad \text{ п.н. }
(c)

Зафиксируйте два непересекающихся подмножества AA и BB в S\mathbb {S}. Вычислите долгосрочную долю переходов из AA в BB.

(d)

Распространите (b), чтобы заключить, что хвостовая последовательность Zn≡{Xn+j:j≥0}Z_{n} \equiv \left\{ X_{n+j} : j \geq 0\right\} цепи Маркова {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0} сходится при n→∞n \rightarrow \infty в смысле конечномерных распределений к строго стационарной последовательности {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0}, являющейся цепью Маркова (S,\vectP)(\mathbb {S}, \vect {P}) с начальным распределением π\pi.

Задача 14.11

Пусть {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0} — неприводимая цепь Маркова, имеющая по крайней мере два состояния. Покажите, что п.н. траектории {Xn}\left\{ X_{n}\right\} не сходятся, т.е. п.н. lim⁡n→∞Xn\lim_{n \rightarrow \infty } X_{n} не существует.

?
Задача 14.12

Пусть {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0} — цепь Маркова с пространством состояний S\mathbb {S} и матрицей переходных вероятностей \vectP≡((pij))\vect {P} \equiv \left(\left(p_{i j}\right)\right). Вероятностное распределение π≡{πj:j∈S}\pi \equiv \left\{ \pi_{j}: j \in \mathbb {S}\right\} называется удовлетворяющим условию детального баланса или обратимости во времени относительно (S,\vectP)(\mathbb {S}, \vect {P}), если для всех i,j,πipij=πjpjii, j, \pi_{i} p_{i j} = \pi_{j} p_{j i}.

?
(a)

Покажите, что такая π\pi обязательно является стационарным распределением.

(b)

Для цепи рождения и гибели (пример 14.1.4) найдите условие в терминах интенсивностей рождения и гибели {αi,βi}i≥0\left\{ \alpha_{i}, \beta_{i}\right\}_{i \geq 0} для существования вероятностного распределения π\pi, удовлетворяющего условию детального баланса.

Задача 14.13

(Вероятности и времена поглощения). Пусть 0 — поглощающее состояние. Для любого i≠0i \neq 0 пусть θi=fi0≡Pi(T0<∞)\theta_{i} = f_{i 0} \equiv P_{i}\left(T_{0} < \infty \right) и ηi=EiT0\eta_{i} = E_{i} T_{0}. Покажите, используя марковское свойство, что для каждого i≠0i \neq 0,

θi=pi0+∑j≠0θjpijηi=1+∑j≠0ηjpij \begin{aligned} \theta _{i} & = p_{i 0}+\sum _{j \neq 0} \theta _{j} p_{i j} \\ \eta _{i} & = 1+\sum _{j \neq 0} \eta _{j} p_{i j} \end{aligned}

Примените это к задаче о разорении игрока с S={0,1,2,…,K}\mathbb {S} = \left\{ 0,1,2, \ldots , K\right\}, K<∞K < \infty и p00=1,pNN=1,pi,i+1=p,pi,i−1=1−p,0<p<1p_{00} = 1, p_{N N} = 1, p_{i, i+1} = p, p_{i, i-1} = 1-p, 0 < p < 1, 1≤i≤N−11 \leq i \leq N-1, и найдите вероятность и математическое ожидание времени ожидания разорения (поглощения в 0), начиная с начального капитала i,1≤i≤N−1i, 1 \leq i \leq N-1.

?
Задача 14.14

(Теория восстановления через цепи Маркова). Пусть {Xj}j≥1\left\{ X_{j}\right\}_{j \geq 1} — н.о.р. случайные величины, принимающие целые положительные значения. Пусть S0=0,Sn=∑j=1nXj,n≥1S_{0} = 0, S_{n} = \sum_{j = 1}^{n} X_{j}, n \geq 1, N(n)=kN(n) = k, если Sk≤n<Sk+1,k=0,1,2,…S_{k} \leq n < S_{k+1}, k = 0,1,2, \ldots — число восстановлений до момента nn, An=n−SN(n)A_{n} = n-S_{N(n)} — возраст текущего элемента в момент nn.

?
(a)

Покажите, что {An}n≥0\left\{ A_{n}\right\}_{n \geq 0} является цепью Маркова, и найдите её пространство состояний S\mathbb {S} и переходные вероятности.

(b)

Предполагая, что E[X1]<∞\mathbb {E}\left[X_{1} \right] < \infty, проверьте, что

πj=P(X1>j)E[X1]j=0,1,2,… \pi _{j} = \frac{P\left(X_{1} > j\right)}{\mathbb {E}\left[X_{1}\right] j} = 0,1,2, \ldots

является единственным стационарным распределением.

(c)

Предполагая, что X1X_{1} имеет апериодическое распределение и что выполняется теорема 14.1.18, покажите, что справедлива дискретная теорема восстановления.

Задача 14.15

Докажите предложение 14.2.1 для случая счётного пространства состояний.

?
Задача 14.16

Докажите предложение 14.2.2.

?
Задача 14.17

Установите утверждение (i) теоремы 14.2.3.

?
Задача 14.18

Покажите, что если P(⋅,⋅)P(\cdot , \cdot ) — переходная функция цепи Маркова {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0}, то для любого n≥0,Px(Xn∈A)=P(n)(x,A)n \geq 0, P_{x}\left(X_{n} \in A\right) = P^{(n)}(x, A), где P(n)(⋅,⋅)P^{(n)}(\cdot , \cdot ) определена итерацией

P(n+1)(x,A)=∫SP(n)(y,A)P(x,dy) P^{(n+1)}(x, A) = \int _{\mathbb {S}} P^{(n)}(y, A) P(x, d y)

с P(0)(x,A)=IA(x)P^{(0)}(x, A) = I_{A}(x).

?
Задача 14.19
?
(a)

Пусть {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0} — случайное блуждание, определённое схемой итерации Xn+1=Xn−+ηn+1X_{n+1} = X_{n}^{-}+\eta_{n+1}, где {ηn}n≥1\left\{ \eta_{n}\right\}_{n \geq 1} н.о.р. случайные величины, независимые от X0X_{0}. Предположим, что ν(⋅)=P(η1∈⋅)\nu (\cdot ) = P\left(\eta_{1} \in \cdot \right) имеет абсолютно непрерывную компоненту с плотностью, строго положительной п.в. на открытом интервале вокруг 0. Покажите, что {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0} неприводима по Харрису относительно меры Лебега на R\mathbb {R}. Покажите, что если, кроме того, E[η1]=0\mathbb {E}\left[\eta_{1}\right] = 0, то {Xn}\left\{ X_{n}\right\} также возвратна по Харрису.

(b)

Используйте теорему 14.2.11, чтобы установить второе утверждение в примере 14.2.10.

Задача 14.20

Покажите, что цепь времени ожидания (пример 14.2.6), определённая как Xn+1=max⁡{Xn+ηn+1,0}X_{n+1} = \max \left\{ X_{n}+\eta_{n+1}, 0\right\}, где {ηn}n≥1\left\{ \eta_{n}\right\}_{n \geq 1} н.о.р., неприводима с базовой мерой ϕ(⋅)≡δ0(⋅)\phi (\cdot ) \equiv \delta_{0}(\cdot ) — дельта-мерой в 0, при условии P(η1<0)>0P\left(\eta_{1} < 0\right) > 0. Покажите далее, что она ϕ\phi-возвратна, если Eη1<0E \eta_{1} < 0.

?
Задача 14.21

Докажите теорему 14.2.5 (i), используя лемму о CC-множестве.

?
Задача 14.22

Найдите h:[0,1]×[0,1]→[0,1]h:[0,1] \times [0,1] \rightarrow [0,1] такую, что h(x,y)h(x, y) разрывна по xx для почти всех yy в [0,1][0,1], и заключите, что функция P(x,A)=P(h(x,Y)∈A)P(x, A) = P(h(x, Y) \in A), где YY — равномерная на [0,1][0,1] случайная величина, не обязательно феллеровская.

?
Задача 14.23

Пусть (Ω,F,P)(\Omega , \mathcal{F}, P) — вероятностное пространство, а (S,S)(\mathbb {S}, \mathcal{S}) — измеримое пространство. Пусть h:S×Ω→Sh: \mathbb {S} \times \Omega \rightarrow \mathbb {S} совместно измерима. Покажите, что P(x,A)≡P(h(x,ω)∈A)P(x, A) \equiv P(h(x, \omega ) \in A) является переходной функцией.

?
Задача 14.24
?
(a)

Пусть {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0} — неприводимая цепь Маркова с пространством состояний S≡{0,1,2,…}\mathbb {S} \equiv \left\{ 0,1,2, \ldots \right\}. Предположим, что V:S→[0,∞)V: \mathbb {S} \rightarrow [0, \infty ) такова, что для некоторого K<∞,ExV(X1)≤V(x)K < \infty , E_{x} V\left(X_{1}\right) \leq V(x) для всех x>Kx > K и что lim⁡x→∞V(x)=∞\lim_{x \rightarrow \infty } V(x) = \infty. Покажите, что {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0} возвратна. (Указание: пусть {X~n}n≥0\left\{ \tilde{X}_{n}\right\}_{n \geq 0} — цепь Маркова с пространством состояний S≡{0,1,2,…}\mathbb {S} \equiv \left\{ 0,1,2, \ldots \right\} и переходными вероятностями, такими же, как у {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0}, за исключением того, что состояния {0,1,2,…,K}\left\{ 0,1,2, \ldots , K\right\} являются поглощающими. Проверьте, что {V(X~n)}n≥0\left\{ V\left(\tilde{X}_{n}\right)\right\}_{n \geq 0} — неотрицательный супермартингал, и, следовательно, что {X~n}n≥0\left\{ \tilde{X}_{n}\right\}_{n \geq 0} ограничена п.н. Теперь заключите, что должно существовать состояние xx, посещаемое бесконечно часто {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0}.)

(b)

Рассмотрим отражающееся неоднородное случайное блуждание на S≡{0,1,2,…}\mathbb {S} \equiv \left\{ 0,1,2, \ldots \right\} такое, что

pij={pi если j=i+11−pi если j=i−1 p_{i j} = \begin{cases} p_{i} & \text{ если } j = i+1 \\ 1-p_{i} & \text{ если } j = i-1 \end{cases}

с p0=1,0<pi≤qip_{0} = 1,0 < p_{i} \leq q_{i} для всех i≥k0i \geq k_{0} и некоторого 1≤k0<∞1 \leq k_{0} < \infty и 0<pi<10 < p_{i} < 1 для всех i≥1i \geq 1. Покажите, что {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0} неприводима и возвратна.

Задача 14.25

Пусть {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0} — неприводимая и возвратная цепь Маркова со счётным пространством состояний S\mathbb {S}. Пусть V:S→R+V: \mathbb {S} \rightarrow \mathbb {R}_{+}такова, что ExV(X1)≤V(x)E_{x} V\left(X_{1}\right) \leq V(x) для всех xx в S\mathbb {S}. Покажите, что V(⋅)V(\cdot ) постоянна на S\mathbb {S}.

?
Задача 14.26

Пусть {Cn}n≥1\left\{ C_{n}\right\}_{n \geq 1} — н.о.р. случайные величины со значениями в [0,4]. Пусть {Xn}n≥0\left\{ X_{n}\right\}_{n \geq 0} — цепь Маркова со значениями в [0,1], определённая схемой случайной итерации

Xn+1=Cn+1Xn(1−Xn),n≥0. X_{n+1} = C_{n+1} X_{n}\left(1-X_{n}\right), \quad n \geq 0.
?
(a)

Покажите, что если Elog⁡C1<0E \log C_{1} < 0, то Xn=O(λn)X_{n} = O\left(\lambda^{n}\right) п.н. для некоторого 0<λ<10 < \lambda < 1.

(b)

Покажите также, что если Elog⁡C1<0E \log C_{1} < 0 и 0<V(log⁡C1)<∞0 < V\left(\log C_{1}\right) < \infty, то существуют последовательности {an}n≥1\left\{ a_{n}\right\}_{n \geq 1} и {bn}n≥1\left\{ b_{n}\right\}_{n \geq 1} такие, что

log⁡Xn−anbn⟶dN(0,1) \frac{\log X_{n}-a_{n}}{b_{n}} \longrightarrow ^{d} N(0,1)