[論文レビュー] Small parts in the Bernoulli sieve

Alexander Gnedin, Alexander Iksanov+1 more|ArXiv.org|Apr 18, 2008

Bayesian Methods and Mixture Models参考文献 7被引用数 8

ひとこと要約

本稿は、棒分割過程に基づく確率的割り当てモデルであるベルヌーイシーブの占有可能数の漸近的挙動を研究し、小さな部分（シングルトン、ダブルトンなど）の連続的分布が、単位ポアソン点過程と自己相似再生過程を含む極限過程に収束することを示している。主な結果として、占有数 $ (K_{n,0}, K_{n,1}, \dots) $ が $ (K_0^*, K_1^*, \dots) $ に弱収束することを確立し、$ \mathbb{E}[K_r^*] = (\mu r)^{-1} $（$ r > 0 $）および $ \mathbb{E}[K_0^*] = \nu / \mu $ という明示的な極限期待値の式が得られる。ここで $ \nu = \mathbb{E}[-\log(1-W)] $ である。

ABSTRACT

Sampling from a random discrete distribution induced by a `stick-breaking' process is considered. Under a moment condition, it is shown that the asymptotics of the sequence of occupancy numbers, and of the small-parts counts (singletons, doubletons, etc) can be read off from a limiting model involving a unit Poisson point process and a self-similar renewal process on the halfline.

研究の動機と目的

• ベルヌーイシーブモデルにおける占有数の漸近的分布—特にシングルトンやダブルトンなどの小さな部分—を理解すること。
• 連続的占有数の分布が、単位ポアソン過程と自己相似再生過程を含む取り扱いやすいモデルに収束する極限的状態を確立すること。
• $ r = 0,1,2,\dots $ における極限期待値 $ \mathbb{E}[K_r^*] $ を特徴づけ、$ W $ のモーメント条件の下で、$ \nu = \mathbb{E}[-\log(1-W)] $ の役割を明確にすること。
• $ \nu = \infty $ の場合を扱い、$ K_0^* = \infty $ がほとんど確実に成り立つことを示し、未占有区間数の発散と関連付けること。
• ポアソン化と点過程の曖昧収束を用いて、$ (K_{n,0}, K_{n,1}, \dots) \to_d (K_0^*, K_1^*, \dots) $ の厳密な収束証明を提供すること。

提案手法

• 本稿では、固定標本の均等割り当てをポアソン過程（レート $ n $）に置き換えるポアソン化技術を用い、占有数の解析を簡素化している。
• 再生定理を確率ウォーク $ \{-\log P_j\} $ に適用することで、スケーリングされた棒分割過程 $ nB $ が $ \mathbb{R}_+ $ 上の自己相似点過程 $ \mathcal{B} $ に曖昧収束することが確立される。
• 極限占有数 $ Z^{(i)} $ は、$ \mathcal{B} $ の左から右への区間内に落ちるポアソン点の数として定義され、$ p(n:m) = \binom{n}{m} \mathbb{E}[W^{n-m}(1-W)^m] $ を含む積の公式から連続的分布が導かれる。
• 基礎となる点過程の収束を用いて、有限標本の占有数 $ (Z_n^{(1)}, Z_n^{(2)}, \dots) $ が極限列 $ (Z^{(1)}, Z^{(2)}, \dots) $ に収束することが証明される。
• 未占有区間数 $ K_{n,0} $ の期待値は、$ \mathcal{B} $ の強度測度 $ (\mu x)^{-1} dx $ を用いて分析され、$ \mathbb{E}[K_0^*] = \nu / \mu $ が得られる。
• ゼロに近い区間からの $ K_{n,[r]} $ の寄与を、元のモデルと $ \mathcal{B} $ のズームインしたバージョンと比較することで解体し、$ s $ が十分に大きいときには一様に小さいことが示される。
• コッヘン＝ストーンの補題の適用により、$ \nu = \infty $ と $ K_0^* = \infty $ がほとんど確実に同値であることが示される。

実験結果

リサーチクエスチョン

• RQ1 ベルヌーイシーブモデルにおいて、小さな部分（例：シングルトン、ダブルトン）の占有数 $ K_{n,r} $ はどのように漸近的に振る舞うか？
• RQ2 $ n \to \infty $ のとき、列 $ (K_{n,1}, K_{n,2}, \dots) $ の極限的連続的分布は何か？ そして、ポアソン過程と自己相似再生過程とどのように関係するか？
• RQ3 極限期待値 $ \mathbb{E}[K_r^*] $ の正確な形は何か？ また、$ W $ のモーメントとパrameter $ \mu = \mathbb{E}[-\log W] $ にどのように依存するか？
• RQ4 区間 $ [U_{n,1}, 1] $ 内の未占有区間数 $ K_{n,0} $ の分布および期待値の収束はどのように成り立つか？
• RQ5 $ \nu = \mathbb{E}[-\log(1-W)] $ の発散と、極限において $ K_0^* = \infty $ がほとんど確実に成り立つことの関係は何か？

主な発見

• 占有数 $ (Z_n^{(1)}, Z_n^{(2)}, \dots) $ の連続的分布は弱収束し、極限確率は $ p(n:m) = \binom{n}{m} \mathbb{E}[W^{n-m}(1-W)^m] $ を含む積の公式によって与えられる。
• 極限期待値としての $ r $-タプルの数は $ \mathbb{E}[K_r^*] = (\mu r)^{-1} $（$ r > 0 $）である。ここで $ \mu = \mathbb{E}[-\log W] $ である。
• $ r = 0 $ の場合、極限期待値は $ \mathbb{E}[K_0^*] = \nu / \mu $ であり、$ \nu = \mathbb{E}[-\log(1-W)] $ で、これは無限大である可能性がある。
• コッヘン＝ストーンの補題により、$ \nu = \infty $ と $ K_0^* = \infty $ がほとんど確実に同値であることが示された。
• $ (K_{n,0}, K_{n,1}, \dots) $ から $ (K_0^*, K_1^*, \dots) $ への収束は、元のモデルとポアソン化バージョンとの比較と、$ \mathcal{B} $ の自己相似性を用いたゼロに近い区間の寄与のバウンディングにより証明された。
• 区間 $ [s/n, 1] $ 内の $ \mathbb{E}[K_{n,[r]}] $ の寄与は、$ s $ が十分に大きいとき一様に小さい。これは $ \mathcal{B} $ の単調性と自己相似構造のおかげで、近似のタイトネスが保証される。