Critical region of likelihood ratio test

This is problem # 5 from RSS's 2014 Graduate Diploma Module 2:

$$\mathbb{P}(X_j=k) = \begin{cases} (1-p)^3 & k=0\\ 3p(1-p) & k=1\\ p^3 & k=2\\ 0 & \text{otherwise} \end{cases}$$

$$Y_k = \sum_{i=1}^n I(X_i=k)\quad\text{for }\, k=0,1,2$$

It is straightforward to get an MLE estimate of $$p$$: $$\hat{p} = \frac{Y_1+3Y_2}{3Y_0+2Y_1+3Y_2}$$

Required to test: $$H_0: p=p_0 \quad \text{versus}\quad H_1: p\neq p_0$$.

I am stuck in the following part:

Using an asymptotic result, find the critical region of this test if its size is to be approximately 0.05

Attempt

Likelihood ratio(LR) test: $$P(LR < k) = \alpha = 0.05$$

\begin{align*} LR &= \frac{L(p|\mathbf{X})|_{p=p_0}}{L(p|\mathbf{X})|_{p=\hat{p}}}\\ &= \frac{(1-p_0)^{3Y_0}(3p_0(1-p_0))^{Y_1}(p_0^3)^{Y_2}}{(1-\hat{p})^{3Y_0}(3\hat{p}(1-\hat{p}))^{Y_1}(\hat{p}^3)^{Y_2}}\\ &= \left(\frac{p_0}{\hat{p}}\right)^{Y_1+3Y_2} \left(\frac{1-p_0}{1-\hat{p}}\right)^{Y_1+3Y_0} \end{align*}

EDIT 1:

Using Wilk's theorem: $$-2\log(LR) \stackrel{D}\longrightarrow \chi^2_1(0.95)$$ in distribution.

Thus critical region is

$$C = \left\{\mathbf{Y}: (Y_1+3Y_2) \log{\frac{p_0}{\hat{p}}} + (Y_1+3Y_0) \log\left(\frac{1-p_0}{1-\hat{p}}\right) \geq \frac{-\chi^2_1(0.95)}{2} \right\}$$

Is this further reducable since the final subpart of the question asks:

Show that in the case $$Y_0 = Y_2$$ this confidence interval can be written as $$p(1-p) \geq 0.25 \exp{\left(\frac{-1.92}{Y_1+3Y_2}\right)}$$

• @Slowloris, it is best not to edit someone's Q to add the [self-study] tag, but to ask them to add it themselves & read our policies. That way there is a greater chance they will be familiar w/ our policy. If they don't add the tag in a reasonable period of time, we can close the thread. There is example text you can use here. – gung - Reinstate Monica May 12 '16 at 17:31
• Have you come across Wilk's Theorem? – Greenparker May 12 '16 at 17:32
• Please add the [self-study] tag & read its wiki. – gung - Reinstate Monica May 12 '16 at 17:33
• @Greenparker Thanks! That takes me a step closer to finding the rejection region. See edit # 1 – user239438 May 12 '16 at 17:56
• Its a little confusing what confidence interval you are talking about, since $C$ is the rejection interval. Also, the "confidence interval" you have has $p$ instead of $p_0$, so again, it is unclear what you mean. In any case, if you plug in $Y_0 = Y_2$, you will notice that $\hat{p}$ turns out to be nice, and solving you get that very expression in terms of $p_0(1-p_0) \geq...$ – Greenparker May 12 '16 at 18:37

1 Answer

$\chi^2_1(0.95) = 3.841$

$C = \{\mathbf{Y}: (Y_1+3Y_2) \log{\frac{p_0}{\hat{p}}} + (Y_1+3Y_0) \log(\frac{1-p_0}{1-\hat{p}}) \geq \frac{-\chi^2_1(0.95)}{2} \}$

When $Y_0 =Y_2$, $\hat{p} = 1/2$

Thus, $\log(4p_0(1-p_0)) \geq \frac{-3.841}{2(Y_1+3Y_2)}$ $\implies p_0(1-p_0) \geq 0.25 \exp{(\frac{-1.92}{Y_1+3Y_2})}$