DOFY's Blog
DOFY's Blog

假设检验

问题的提法

H_0: \theta \in \Theta_0 \leftrightarrow H_1:\theta \in \Theta_1

H_0:原假设,H_1:备择假设。

根据数据\bm x=(x_1,\cdots,x_n),回答接受H_0或拒绝H_0

检验方法:给出一个否定域\mathcal W,若\bm x\in {\mathcal W}则拒绝。不能根据数据选择否定域

功效函数\beta_{\mathcal W}(\theta)\triangleq P_\theta(\bm x\in \mathcal W)

水平\alpha\sup_{\theta\in \Theta_0}\beta_{\mathcal W}(\theta)\le \alpha

一致最大功效否定域(UMP):

1.\mathcal W 水平为\alpha
2.\forall \tilde {\mathcal W}\tilde {\mathcal W}水平为\alpha,有

\beta_{\mathcal W}(\theta)\ge \beta_{\tilde {\mathcal W}}(\theta)\qquad \forall \theta\in \Theta_1

无偏否定域

\beta_{\mathcal W}(\theta_0)\le \alpha\le \beta_{\mathcal W}(\theta_1)\qquad \forall \theta_0\in \Theta_0,\theta_1\in \Theta_1

这是说否定域的功效在\Theta_1上总是高于\Theta_0

最优无偏否定域(UMPU):

  • $\mathcal W$水平为$\alpha$,且无偏。
  • $\forall \tilde {\mathcal W}$,$\tilde {\mathcal W}$水平为$\alpha$,且$\tilde{\mathcal W}$无偏,有

\beta_{\mathcal W}(\theta_1)\ge \beta_{\tilde {\mathcal W}}(\theta_1)\qquad \forall \theta_1\in \Theta_1

可以认为是无偏否定域中的最大一致功效。

似然比检验

解决的问题:

H_0: \theta = \theta_0 \leftrightarrow H_1:\theta = \theta_1

即原假设和备择假设都只有一个参数。

似然函数:L(\bm x,\theta)=\prod_{i=1}^n f(x_i,\theta)

似然比

\frac{L(\bm x,\theta_1)}{L(\bm x,\theta_0)}

于是就有一种基于似然比的否定域类:

{\mathcal W}_{\lambda_0}=\{\bm x:L(\bm x,\theta_1)>\lambda_0 L(\bm x,\theta_0)\}

\lambda_0的取值由给定的水平\alpha确定,即

\beta_{\mathcal W}(\theta_0)=\int_{\mathcal W}L(\bm x,\theta_0)\mathrm{d}\bm x=\alpha

(Neyman-Pearson引理) {\mathcal W}_{\lambda_0}是水平为\alpha的UMP否定域。

证明思路:

\begin{gathered}
\int_{{\mathcal W}\setminus {\tilde{\mathcal W}}} L(\bm x,\theta_1)\mathrm{d} \bm x> \lambda_0 \int_{{\mathcal W}\setminus {\tilde{\mathcal W}}} L(\bm x,\theta_1) \mathrm{d} \bm x\\
\int_{{\tilde{\mathcal W}}\setminus {\mathcal W}} L(\bm x,\theta_1)\mathrm{d} \bm x\le \lambda_0 \int_{{\tilde{\mathcal W}}\setminus {\mathcal W}} L(\bm x,\theta_1)\mathrm{d} \bm x\\
\int_{{\mathcal W}\setminus {\tilde{\mathcal W}}} L(\bm x,\theta_1)\mathrm{d} \bm x=\int_{{\tilde{\mathcal W}}\setminus {\mathcal W}} L(\bm x,\theta_1)\mathrm{d} \bm x
\end{gathered}

{\mathcal W}_{\lambda_0}是无偏否定域,从而也是UMPU否定域。

证明思路:

\begin{aligned}
\beta_{\mathcal W}(\theta_1)-\alpha &=\int_{\mathcal W} L(\bm x,\theta_1)\mathrm{d} \bm x -\alpha\\
&=(1-\alpha)\int_{\mathcal W}L(\bm x,\theta_1)\mathrm{d} \bm x-\alpha\int_{\mathcal W^c}L(\bm x,\theta_1)\mathrm{d} \bm x\\
&\ge \lambda_0(1-\alpha)\int_{\mathcal W}L(\bm x,\theta_0)\mathrm{d} \bm x-\lambda_0\alpha\int_{\mathcal W^c}L(\bm x,\theta_0)\mathrm{d} \bm x\\
&=0
\end{aligned}

单参数模型的检验

稍微复杂的问题:

H_0: \theta \in \Theta_0 \leftrightarrow H_1:\theta = \theta_1

H_0: \theta \in \Theta_0 \leftrightarrow H_1:\theta \in \Theta_1

\Theta_0,\Theta_1不在是单点集。

不难验证,若否定域\mathcal W对任意的\theta_0\in \Theta_0都是UMP否定域,且不依赖于\theta_1,那么它也是这个问题的UMP否定域

单参数指数分布族

单参数指数分布族

f(x,\theta)=S(\theta)h(x)\exp\{C(\theta)T(x)\},\qquad \theta \in \Theta

其中\Theta是一个区间,C(\theta)严格增。

对于否定域类

{\mathcal W}=\{\bm x:\sum_{i=1}^n T(x_i)>c\}

其中c为任一常数,则\beta_{\mathcal W}(\theta)\theta的不减函数。

证明思路:

\theta_1<\theta_2,对于简单假设检验问题:

H_0: \theta =\theta_1 \leftrightarrow H_1:\theta =\theta_2

考虑似然比

\frac{L(\bm x,\theta_1)}{L(\bm x,\theta_2)}=\left[\frac{S(\theta_2)}{S(\theta_1)} \right]^n\exp\left\{[C(\theta_2)-C(\theta_1)]\sum_{i=1}^n T(x_i)\right\}

所以\mathcal W是该问题的一个似然比否定域,那么\mathcal W是无偏否定域,从而

\beta_{\mathcal W}(\theta_2)\ge \alpha\ge \beta_{\mathcal W}(\theta_1)

形象理解:\theta越大,T(x)的密度越向右偏移。

单边假设检验问题

H_0:\theta\le \theta_0 \leftrightarrow \theta>\theta_0

否定域类型:

{\mathcal W}=\{\bm x:\sum_{i=1}^n T(x_i)>c\}

选取c使得\beta_{\mathcal W}(\theta_0)=\alpha,那么\mathcal W就是UMP否定域。

$\sigma^2$已知,检验$\mu$

H_0:\mu\le \mu_0 \leftrightarrow H_1:\mu >\mu_0

f(x,\mu)=S(\mu)h(x)e^{\frac{\mu}{\sigma^2}x}\implies T(x)=x

{\mathcal W}=\{\bm x:\sum_{i=1}^n x_i>c\}

需要验证在H_0的情况下,\beta_{\mathcal W}(\mu)\le \beta_{\mathcal W}(\mu_0)(即最值通常在边界取到,后续问题省略此步骤)。

\beta_{\mathcal W}(\mu_0)=P_{\mu_0}\left(\sqrt n\frac{\bar X-\mu_0}{\sigma}>\tilde c\right)=P(Z>\tilde c)=\alpha

查表确定\tilde c,代入数据计算即可。

$\mu$已知,检验$\sigma^2$

H_0:\sigma^2\le \sigma_0^2 \leftrightarrow H_1:\sigma^2 >\sigma_0^2

f(x,\sigma^2)=\frac{1}{\sqrt {2\pi\sigma^2}}\exp\{\frac{1}{2\sigma^2}(x-\mu)^2\}\implies T(x)=(x-\mu)^2

{\mathcal W}=\{\bm x:\sum_{i=1}^n (x_i-\mu)^2>c\}

\beta_{\mathcal W}(\sigma_0^2)=P_{\sigma_0^2}\left(\sum_{i=1}^n\left(\frac{X_i-\mu}{\sigma_0}\right)^2< \tilde c\right)=P(K_n<\tilde c)=\alpha

查表确定\tilde c,代入数据计算即可。

双边假设检验问题

H_0:\theta=\theta_0\leftrightarrow H_1:\theta\neq\theta_0

双面假设检验问题不存在UMP否定域(无法兼顾两边)。

但有可能存在UMPU否定域。

若当\theta=\theta_0时,\frac 1n\sum_{i=1}^nT(X_i)的分布关于某个r_0对称,则

{\mathcal W}=\left\{\bm x:\left|\sum_{i=1}^nT(x_i)-r_0\right|>c\right\}

是UMPU否定域

证明略。

$\sigma^2$已知,检验$\mu$

H_0:\mu=\mu_0\leftrightarrow H_1:\mu\neq\mu_0

T(x)=x,r_0=\mu_0

{\mathcal W}=\left\{\bm x:\left|\sum_{i=1}^nx_i-n\mu_0\right|>c\right\}

\beta_{\mathcal W}(\mu_0)=P_{\mu_0}(\sqrt{n}\frac{|\bar x-\mu_0|}{\sigma}>c)=P(|Z|>c)=\alpha

广义似然比检验

H_0:\theta\in \Theta_0\leftrightarrow H_1:\theta\in \Theta_1

广义似然比

\frac{L(\bm x,\hat\theta)}{L(\bm x,\hat \theta_0)}

其中L(\bm x,\hat\theta_0)=\max_{\theta\in \Theta_0}L(\bm x,\theta)L(\bm x,\hat\theta)=\max_{\theta\in \Theta}L(\bm x,\theta),注意后者是在全局上求的最值。

基于广义似然比的否定域

{\mathcal W}=\{\bm x:\lambda(\bm x)=\frac{L(\bm x,\hat\theta)}{L(\bm x,\hat \theta_0)}>c\}

$\sigma^2$已知,检验$\mu$,双边

H_0:\mu=\mu_0\leftrightarrow H_1:\mu\neq\mu_0

最大似然估计:

\hat \mu=\bar x,\hat \mu_0=\mu_0

广义似然比:

\lambda(\bm x)=\frac{L(\bm x,\hat \mu)}{L(\bm x,\hat \mu_0)}=\exp\left\{\frac{1}{2\sigma^2}\left(\sum_{i=1}^n(x_i-\mu_0)^2-\sum_{i=1}^n(x_i-\bar x)^2\right)\right\}

{\mathcal W}=\{\bm x:\lambda(\bm x)>c\}=\{\bm x:|\bar x-\mu_0|>c\}

\beta_{\mathcal W}(\mu_0)=P_{\mu_0}(\sqrt{n}\frac{|\bar x-\mu_0|}{\sigma}>c)=P(|Z|>c)=\alpha

$\sigma^2$已知,检验$\mu$,单边

H_0:\mu\le \mu_0 \leftrightarrow H_1:\mu >\mu_0

最大似然估计:

\hat \mu=\bar x,\hat \mu_0=\mu_0(\bar x>\mu_0)

广义似然比:

\lambda(\bm x)=\exp\left(\frac{1}{2}\left(\frac{\sqrt n(\bar x-\mu_0)}{\sigma}\right)^2\right)

{\mathcal W}=\{\bm x:\frac{\sqrt n(\bar x-\mu_0)}{\sigma}>c\}

\max_{\mu\le \mu_0}\beta_{\mathcal W}(\mu)=P_{\mu_0}(\frac{\sqrt n(\bar x-\mu_0)}{\sigma}>c)=P(Z>c)=\alpha

$\sigma^2$未知,检验$\mu$,双边

H_0:\mu=\mu_0\leftrightarrow H_1:\mu\neq\mu_0

计算可得

L(\hat \mu,\hat \sigma^2)=(2\pi\hat \sigma^2)^{-n/2}e^{-n/2},L(\hat \mu_0,\hat \sigma_0^2)=(2\pi\hat \sigma_0^2)^{-n/2}e^{-n/2}

orz写不完了5555

UPD: 果不其然,期末考试翻在了没写完的地方

发表回复

textsms
account_circle
email

DOFY's Blog

假设检验
概统背不完了。
扫描二维码继续阅读
2021-01-08