logo

多変量確率変数の分布収束 📂数理統計学

多変量確率変数の分布収束

定義1

pp次元ランダムベクトルX\mathbf{X}とランダムベクトルのシークエンス{Xn}\left\{ \mathbf{X}_{n} \right\}が次の条件を満たす時、nn \to \inftyで、Xn\mathbf{X}_{n}X\mathbf{X}に分布収束すると言い、XnDX\mathbf{X}_{n} \overset{D}{\to} \mathbf{X}と示される。 limnFXn(x)=FX(x),xCFX\lim_{n \to \infty} F_{\mathbf{X}_{n}} (x) = F_{\mathbf{X}} (x) \qquad, \forall x \in C_{F_{\mathbf{X}}}


  • FXF_{X}は確率変数XXの累積分布関数だ。
  • CFXC_{F_{\mathbf{X}}}は関数FXF_{\mathbf{X}}が連続である点の集合を表している。

多変量中心極限定理

{Xn}\left\{ \mathbf{X}_{n} \right\}が平均ベクトルμRp\mathbf{\mu} \in \mathbb{R}^{p}と共分散行列ΣRp×p\Sigma \in \mathbb{R}^{p \times p}を持つiidのランダムベクトルのシークエンスとする。ゼロベクトル0\mathbf{0}の近傍で、モーメント生成関数m(t)m ( \mathbf{t})の存在を仮定し、Yn\mathbf{Y}_{n}を以下のように定義する。 Yn:=1nk=1n(Xkμ)=n(Xμ)\mathbf{Y}_{n} := {{ 1 } \over { \sqrt{n} }} \sum_{k=1}^{n} \left( \mathbf{X}_{k} - \mathbf{\mu} \right) = \sqrt{n} \left( \overline{\mathbf{X}} - \mathbf{\mu} \right) すると、Yn\mathbf{Y}_{n}は多変量正規分布Np(0,Σ)N_{p} \left( \mathbb{0} , \mathbf{\Sigma} \right)に分布収束する。

証明

ゼロベクトル0\mathbf{0}の近傍のtRp\mathbf{t} \in \mathbb{R}^{p}に対して、YnY_{n}のモーメント生成関数は次の通り。Wk:=t(Xμ)W_{k} := \mathbf{t} ' \left( \mathbf{X} - \mathbf{\mu} \right)とする Mn(t)=E[exp{t1nk=1n(Xkμ)}]=E[exp{1nk=1nt(Xkμ)}]=E[exp{1nk=1nWk}] \begin{align*} M_{n} \left( \mathbf{t} \right) =& E \left[ \exp \left\{ \mathbf{t}’ { { 1 } \over { \sqrt{n} } } \sum_{k=1}^{n} \left( \mathbf{X}_{k} - \mathbf{\mu} \right) \right\} \right] \\ =& E \left[ \exp \left\{ { { 1 } \over { \sqrt{n} } } \sum_{k=1}^{n} \mathbf{t} ' \left( \mathbf{X}_{k} - \mathbf{\mu} \right) \right\} \right] \\ =& E \left[ \exp \left\{ { { 1 } \over { \sqrt{n} } } \sum_{k=1}^{n} W_{k} \right\} \right] \end{align*} ここでWkW_{k}たちはiidで平均が00、分散がVar(Wk)=tΣt\operatorname{Var} \left( W_{k} \right) = \mathbf{t} ' \mathbf{\Sigma} \mathbf{t}なので、一変量中心極限定理によると 1nk=1nWkDN(0,tΣt) { { 1 } \over { \sqrt{n} } } \sum_{k=1}^{n} W_{k} \overset{D}{\to} N \left( 0, \mathbf{t} ' \mathbf{\Sigma} \mathbf{t} \right) これは、nn \to \inftyの時、Mn(t)M_{n} \left( \mathbf{t} \right)Mn(t)etΣt/2 M_{n} \left( \mathbf{t} \right) \to e^{\mathbf{t} ' \mathbf{\Sigma} \mathbf{t} / 2} これは多変量正規分布Np(0,Σ)N_{p} \left( \mathbf{0}, \mathbf{\Sigma} \right)のモーメント生成関数である。


  1. Hogg et al. (2013). Introduction to Mathematical Statistcs(7th Edition): p322. ↩︎