두 개의 샘플 카이 제곱 테스트

이 질문은 Van der Vaart의 저서 Asymptotic Statistics, pg에서 발췌 한 것입니다. 253. # 3 :

및 이 매개 변수 및 갖는 독립 다항식 벡터 라고 가정합니다 . 귀무 가설에서 는 $\mathbf{X}_m$ $\mathbf{Y}_n$ $(m,a_1,\ldots,a_k)$ $(n,b_1,\ldots,b_k)$ $a_i=b_i$

\sum_{i = 1}^{k} \frac{(X_{m, i} - m {\hat{c}}_{i})^{2}}{m {\hat{c}}_{i}} + \sum_{i = 1}^{k} \frac{(Y_{n, i} - n {\hat{c}}_{i})^{2}}{n {\hat{c}}_{i}}

$\sum_{i=1}^k \dfrac{(X_{m,i} - m\hat{c}_i)^2}{m\hat{c}_i} + \sum_{i=1}^k \dfrac{(Y_{n,i} - n\hat{c}_i)^2}{n\hat{c}_i}$ 에는 분포가 있습니다. 여기서 입니다.

χ_{k - 1}^{2}

$\chi^2_{k-1}$

{\hat{c}}_{i} = (X_{m, i} + Y_{n, i}) / (m + n)

$\hat{c}_i = (X_{m,i} + Y_{n,i})/(m+n)$

시작하는 데 도움이 필요합니다. 여기서 전략은 무엇입니까? 나는 두 개의 summand를 다음과 같이 결합 할 수있었습니다.

\sum_{i = 1}^{k} \frac{(m Y_{n, i} - n X_{m, i})^{2}}{m n (m + n) {\hat{c}}_{i}}

$\sum_{i=1}^k \dfrac{(mY_{n,i} - nX_{m,i})^2}{mn(m+n)\hat{c}_i}$

그러나 과 의 가중치 조합이므로 CLT에서는 작동하지 않습니다 . 이것이 올바른 경로인지 확실하지 않습니다. 어떤 제안? $X_m$ $Y_n$

편집 : 경우 우리가 얻을 수 있기 때문에 매우 쉽습니다. $m=n$

\begin{aligned} \frac{m Y_{n} - n X_{m}}{\sqrt{m n (m + n)}} & = \frac{Y_{n} - X_{m}}{\sqrt{(m + n)}} \end{aligned}

$\begin{align*} \dfrac{mY_{n} - nX_{m}}{\sqrt{mn(m+n)}} &= \dfrac{Y_{n} - X_{m}}{\sqrt{(m+n)}} \end{align*}$

분자는 다항식 변수 의 차이의 합으로 볼 수 있으므로 CLT를 적용한 다음 같은 장의 정리 17.2로 마무리 할 수 있습니다. 그러나이 상황에서 다른 샘플 크기 로이 문제를 해결하는 방법을 알 수 없습니다. 어떤 도움? $(1,a_1,\ldots,a_k)$

van der Vaart의 Google 도서 17 장 링크

— 브데 오노 비치
소스

먼저 몇 가지 표기법. 하자가 및 와 관련된 범주 시퀀스 나타낸다 와 , 즉 . 이라고하자 . 이진화 여기서 는 크로네 커 델타입니다. 그래서 우리는 $\left\{X_t\right\}_{1,\ldots,m}$ $\left\{Y_t\right\}_{1,\ldots,n}$ $\mathbf{X}_m$ $\mathbf{Y}_n$ $\text{Pr}\left\{X_t = i\right\} = a_i, \text{Pr}\left\{Y_t = i\right\} = b_i$ $N=n+m$

\begin{aligned} X_{i}^{*} & = (X_{1, i}^{*}, \dots, X_{N, i}^{*}) = (δ_{i, X_{1}}, \dots, δ_{i, X_{n}}, 0, \dots, 0) \\ Y_{i}^{*} & = (Y_{1, i}^{*}, \dots, Y_{N, i}^{*}) = (0, \dots, 0, δ_{i, Y_{1}}, \dots, δ_{i, Y_{n}}) \end{aligned}

$\begin{align*} \mathbf{X}_{i}^* &= (X^*_{1,i},\ldots,X_{N,i}^*) = (\delta_{i,X_1},\ldots,\delta_{i,X_n},0,\ldots,0)\\ \mathbf{Y}_{i}^* &= (Y^*_{1,i},\ldots,Y_{N,i}^*)= (0,\ldots,0,\delta_{i,Y_1},\ldots,\delta_{i,Y_n})\\ \end{align*}$

δ_{i, j} \equiv 1_{i = j}

$\delta_{i,j}\equiv \mathbf{1}_{i=j}$

X_{m, i} = \sum_{t = 1}^{N} X_{t, i}^{*} = \sum_{t = 1}^{m} δ_{i, X_{t}} Y_{n, i} = \sum_{t = 1}^{N} Y_{t, i}^{*} = \sum_{t = 1}^{n} δ_{i, Y_{t}}

$X_{m,i} = \sum_{t=1}^{N} X_{t,i}^* = \sum_{t=1}^m \delta_{i,X_t} \qquad Y_{n,i} = \sum_{t=1}^{N} Y_{t,i}^* = \sum_{t=1}^n \delta_{i,Y_t}$

이제 우리는 증거를 시작합니다. 먼저 테스트 통계의 두 가지 요약을 결합합니다. 참고 따라서 테스트 통계를

\begin{aligned} X_{m, i} - m {\hat{c}}_{i} & = \frac{(n + m) X_{m, i} - m (X_{m, i} + Y_{n, i})}{n + m} \\ = \frac{n X_{m, i} - m Y_{n, i}}{n + m} \\ Y_{n, i} - n {\hat{c}}_{i} & = \frac{(n + m) Y_{n, i} - n (X_{m, i} + Y_{n, i})}{n + m} \\ = \frac{m Y_{n, i} - n X_{m, i}}{n + m} \end{aligned}

$\begin{align*} X_{m,i} - m\hat{c}_i &= \dfrac{(n+m)X_{m,i} - m(X_{m,i} + Y_{n,i})}{n+m}\\ &= \dfrac{nX_{m,i} - mY_{n,i}}{n+m}\\ Y_{n,i} - n\hat{c}_i &= \dfrac{(n+m)Y_{n,i} - n(X_{m,i} + Y_{n,i})}{n+m}\\ &= \dfrac{mY_{n,i} - nX_{m,i}}{n+m} \end{align*}$

\begin{aligned} S & = \sum_{i = 1}^{k} \frac{(X_{m, i} - m {\hat{c}}_{i})^{2}}{m {\hat{c}}_{i}} + \sum_{i = 1}^{k} \frac{(Y_{n, i} - n {\hat{c}}_{i})^{2}}{n {\hat{c}}_{i}} \\ = \sum_{i = 1}^{k} \frac{(n X_{m, i} - m Y_{n, i})^{2}}{(n + m)^{2} m {\hat{c}}_{i}} + \sum_{i = 1}^{k} \frac{(n X_{m, i} - m Y_{n, i})^{2}}{(n + m)^{2} n {\hat{c}}_{i}} \\ = \sum_{i = 1}^{k} \frac{(n X_{m, i} - m Y_{n, i})^{2}}{n m (n + m) {\hat{c}}_{i}} \end{aligned}

$\begin{align*} S &= \sum_{i=1}^k \dfrac{(X_{m,i} - m\hat{c}_i)^2}{m\hat{c}_i} + \sum_{i=1}^k \dfrac{(Y_{n,i} - n\hat{c}_i)^2}{n\hat{c}_i}\\ &= \sum_{i=1}^k \dfrac{(nX_{m,i} - mY_{n,i})^2}{(n+m)^2m\hat{c}_i} + \sum_{i=1}^k \dfrac{(nX_{m,i} - mY_{n,i})^2}{(n+m)^2n\hat{c}_i}\\ &= \sum_{i=1}^k \dfrac{(nX_{m,i} - mY_{n,i})^2}{nm(n+m)\hat{c}_i} \end{align*}$

다음 참고 와 다음 속성

n X_{m, i} - m Y_{n, i} = \sum_{t = 1}^{N} n X_{t, i}^{*} - m Y_{t, i}^{*} = Z_{i}

$nX_{m,i} - mY_{n,i} = \sum_{t=1}^N nX_{t,i}^* - mY_{t,i}^* = Z_{i}$

\begin{aligned} E [Z_{i}] & = n E [X_{m, i}] - m E [Y_{n, i}] \\ = n m a_{i} - n m a_{i} = 0 \\ Var [Z_{i}] & = Var [n X_{m, i} - m Y_{n, i}] \\ = n^{2} Var [X_{m, i}] - m^{2} Var [Y_{n, i}] Note X_{m, i} and Y_{n, i} are independent \\ = n^{2} m a_{i} (1 - a_{i}) + m^{2} n a_{i} (1 - a_{i}) \\ = n m (n + m) a_{i} (1 - a_{i}) \\ Cov [Z_{i}, Z_{j}] & = E [Z_{i} Z_{j}] - E [Z_{i}] E [Z_{j}] \\ = E [(n X_{m, i} - m Y_{n, i}) (n X_{m, j} - m Y_{n, j})] \\ = n^{2} (- m a_{i} a_{j} + m^{2} a_{i} a_{j}) - 2 n^{2} m^{2} a_{i} a_{j} + m^{2} (- n a_{i} a_{j} + n^{2} a_{i} a_{j}) \\ = - n m (n + m) a_{i} a_{j} \end{aligned}

$\begin{align*} \text{E}[Z_{i}] &= n\text{E}[X_{m,i}] - m\text{E}[Y_{n,i}]\\ &= nma_i - nma_i = 0\\ \text{Var}[Z_{i}] &= \text{Var}[nX_{m,i} - mY_{n,i}]\\ &= n^2\text{Var}[X_{m,i}] - m^2\text{Var}[Y_{n,i}] \qquad\text{Note $X_{m,i}$ and $Y_{n,i}$ are independent}\\ &= n^2ma_i(1-a_i) + m^2na_i(1-a_i)\\ &= nm(n+m)a_i(1-a_i)\\ \text{Cov}[Z_{i},Z_{j}] &= \text{E}[Z_{i}Z_{j}] - \text{E}[Z_{i}]\text{E}[Z_{j}]\\ &= \text{E}[(nX_{m,i} - mY_{n,i})(nX_{m,j} - mY_{n,j})]\\ &= n^2(-ma_ia_j + m^2a_ia_j) - 2n^2m^2a_ia_j + m^2(-na_ia_j+n^2a_ia_j)\\ &= -nm(n+m)a_ia_j \end{align*}$

따라서 다변량 CLT에 의해 여기서 , 의 번째 요소 입니다. 이후 Slutsky는 여기서 는 항등 행렬입니다.

\frac{1}{\sqrt{n m (n + m)}} Z = \frac{n X_{m} - m Y_{n}}{\sqrt{n m (n + m)}} \overset{D}{\to} N (0, Σ)

$\dfrac{1}{\sqrt{nm(n+m)}}\mathbf{Z} = \dfrac{n\mathbf{X}_m - m \mathbf{Y}_n}{\sqrt{nm(n+m)}}\overset{D}{\to} \text{N}(\mathbf{0},\Sigma)$

(i, j)

$(i,j)$

Σ

$\Sigma$

σ_{i j} = a_{i} (δ_{i j} - a_{j})

$\sigma_{ij} = a_i(\delta_{ij} - a_j)$

\hat{c} = ({\hat{c}}_{1}, \dots, {\hat{c}}_{k}) \overset{p}{\to} (a_{1}, \dots, a_{k}) = a

$\hat{\mathbf{c}} = (\hat{c}_1,\ldots,\hat{c}_k) \overset{p}{\to} (a_1,\ldots,a_k)=\mathbf{a}$

\frac{n X_{m} - m Y_{n}}{\sqrt{n m (n + m)} \hat{c}} \overset{D}{\to} N (0, I_{k} - \sqrt{a} {\sqrt{a}}^{'})

$\dfrac{n\mathbf{X}_m - m \mathbf{Y}_n}{\sqrt{nm(n+m)}\hat{\mathbf{c}}}\overset{D}{\to} \text{N}(\mathbf{0},\mathbf{I}_k - \sqrt{\mathbf{a}}\sqrt{\mathbf{a}}')$

I_{k}

$\mathbf{I}_k$

k \times k

$k\times k$

\sqrt{a} = (\sqrt{a_{1}}, \dots, \sqrt{a_{k}})

$\sqrt{\mathbf{a}} = (\sqrt{a_1},\ldots,\sqrt{a_k})$ . 이후 multiplicty 1 고유치 0 다중성의 고유 한 갖고 (연속 매핑 법칙에 의해 또는 참조 Lemma 17.1, van der Vaart의 정리 17.2)

I_{k} - \sqrt{a} {\sqrt{a}}^{'}

$\mathbf{I}_k - \sqrt{\mathbf{a}}\sqrt{\mathbf{a}}'$

k - 1

$k-1$

\sum_{i = 1}^{k} \frac{(n X_{m, i} - m Y_{n, i})^{2}}{n m (n + m) {\hat{c}}_{i}} \overset{D}{\to} χ_{k - 1}^{2}

$\sum_{i=1}^k \dfrac{(nX_{m,i} - mY_{n,i})^2}{nm(n+m)\hat{c}_i} \overset{D}{\to} \chi^2_{k-1}$

— 브데 오노 비치
소스