[論文レビュー] A Berry Esseen Theorem for the Lightbulb Process

Larry Goldstein, Haimeng Zhang|arXiv (Cornell University)|Jan 5, 2010

Random Matrices and Applications参考文献 11被引用数 7

ひとこと要約

本稿は、n個の電球が逐次切り替えられるマルコフ連鎖モデル「ライトバルブ過程」に対して、Berry-Esseen定理を確立し、正規近似誤差が $ O(n^{-1/2}) $ の速度で減少することを証明している。有界で単調なサイズバイアスカップリングを用いたシュタインの方法により、有限nに対して、標準化された点灯電球数と標準正規分布との間のコルモゴロフ-スミルノフ距離に対する明示的な上限が得られ、特に偶数nでは、過程のスペクトル分解によりより鋭い制御が可能である。

ABSTRACT

In the so called lightbulb process, on days $r=1,..., n$, out of $n$ lightbulbs, all initially off, exactly $r$ bulbs, selected uniformly and independent of the past, have their status changed from off to on, or vice versa. With $X$ the number of bulbs on at the terminal time $n$, an even integer, and $μ=n/2, σ^2=Var(X)$, we have $$ \sup_{z \in \mathbb{R}} |P(\frac{X-μ}σ \le z)-P(Z \le z)| \le \frac{n}{2σ^2} \barΔ_0 + 1.64 \frac{n}{σ^3}+ \frac{2}σ $$ where $Z$ is a standard normal random variable, and $$ \barΔ_0 = 1/2\sqrt{n}} + \frac{1}{2n} + 1/3 e^{-n/2} \qmq {for $n \ge 6$,} $$ yielding a bound of order $O(n^{-1/2})$ as $n o \infty$. A similar, though slightly larger bound holds for $n$ odd. The results are shown using a version of Stein's method for bounded, monotone size bias couplings. The argument for even $n$ depends on the construction of a variable $X^s$ on the same space as $X$ that has the $X$-size bias distribution, that is, that satisfies \beas E [X g(X)] =μE[g(X^s)] \quad for all bounded continuous $g$, \enas and for which there exists a $B \ge 0$, in this case B=2, such that $X \le X^s \le X+B$ almost surely. The argument for $n$ odd is similar to that for $n$ even, but one first couples $X$ closely to $V$, a symmetrized version of $X$, for which a size bias coupling of $V$ to $V^s$ can proceed as in the even case. In both the even and odd cases, the crucial calculation of the variance of a conditional expectation requires detailed information on the spectral decomposition of the lightbulb chain.

研究の動機と目的

• ライトバルブ過程（逐次的なランダムな切り替えを行うマルコフ連鎖モデル）における点灯電球数の有限標本に対するBerry-Esseen上限を確立すること。
• 終端時の点灯電球数の分布に関する漸近的正規性と近似精度の未解決問題を解消すること。
• 過程に内在する有界で単調な依存構造に特化したサイズバイアスカップリング技術の開発および適用すること。
• 偶数nおよび奇数nの両方に対して、コルモゴロフ-スミルノフ距離における明示的かつ非漸近的な誤差上限を提供すること。

提案手法

• 有界で単調なサイズバイアスカップリングを用いたシュタインの方法により、正規近似誤差を制限する。
• 終端カウント $ X $ のサイズバイアス版 $ X^s $ を構築し、$ X \leq X^s \leq X + 2 $ がほとんど確実に成り立ち、任意の有界連続関数 $ g $ に対して $ \mathbb{E}[Xg(X)] = \mu \mathbb{E}[g(X^s)] $ が成り立つようにする。
• 偶数 $ n $ の場合、$ X $ をそのサイズバイアス版に直接カップリングする。奇数 $ n $ の場合、まず $ X $ を対称化変数 $ V $ にカップリングし、その後 $ V $ をサイズバイアス化する。
• シュタインの方法の境界において重要な要素である条件付き期待値の分散を計算するために、ライトバルブ過程のスペクトル分解を用いる。
• 降べき階乗を用いて $ \lambda_{n,b,{\bf s}} $ の明示的表現を導出し、スイッチパターン $ {\bf s} $ の下での $ X_{{\bf s}} $ の平均と分散を特徴づける。
• スイッチの配置に関する組合せ的恒等式および条件付き期待値の計算を用い、シュタインの境界における重要な項の正当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 ライトバルブ過程における点灯電球数の正規分布への収束速度は何か？
• RQ2 シュタインの方法を用いて、この従属的かつi.i.d.でない和に対してBerry-Esseen上限を確立できるか？
• RQ3 ライトバルブ過程の構造、特に増加する数の電球が逐次的に切り替えられるという性質が、正規近似誤差にどのように影響するか？
• RQ4 基礎となるマルコフ連鎖のスペクトル分解が、サイズバイアスカップリング誤差の制御に果たす役割は何か？

主な発見

• 偶数 $ n $ の場合、標準化された終端カウント $ X $ と標準正規変数 $ Z $ の間のコルモゴロフ-スミルノフ距離は、$ \frac{n}{2\sigma^2}\overline{\Delta}_0 + 1.64\frac{n}{\sigma^3} + \frac{2}{\sigma} $ で抑えられ、$ n \geq 6 $ のとき $ \overline{\Delta}_0 = \frac{1}{2\sqrt{n}} + \frac{1}{2n} + \frac{1}{3}e^{-n/2} $ である。
• この上限は $ n \to \infty $ のとき $ O(n^{-1/2}) $ の速度で減少し、定量的な収束速度を伴う漸近的正規性が確認される。
• 奇数 $ n $ の場合、対称化カップリングを介して $ V $ に事前にカップリングすることで、類似した上限が得られ、定数はわずかに大きくなる。
• シュタインの方法フレームワークにおける条件付き期待値の分散は、ライトバルブ過程の遷移行列の詳細なスペクトル分解により計算されている。
• 解析により、正規近似誤差がスイッチパターン $ {\bf s} = (1,\ldots,n) $ の構造によって制御されることを確認し、平均と分散を支配する $ \lambda_{n,b,{\bf s}} $ の明示的表現が得られている。
• 証明では、スイッチ配置に関する条件付き重みに関する重要な組合せ的恒等式が検証されており、サイズバイアスカップリング項の正当性が保証されている。