跳转至

单样本率非劣效检验

样本率用 \(\hat{p}\) 表示,总体率用 \(p\) 表示,非劣界值用 \(\delta\) 表示。

对于高优指标(\(\delta < 0\)),统计学假设如下:

\[ \begin{align} H_0 &: p - p_0 \leqslant \delta \\ H_1 &: p - p_0 \gt \delta \end{align} \]

对于低优指标(\(\delta > 0\)),统计学假设如下:

\[ \begin{align} H_0 &: p - p_0 \geqslant \delta \\ H_1 &: p - p_0 \lt \delta \end{align} \]

以下推导过程在边界条件 \(p - p_0 = \delta\) 下进行。

z-test using s(p0)

\(H_0\) 成立时,可构建 \(z\) 统计量:

\[ z = \frac{\hat{p}-p_0-\delta}{\sqrt{(p_0+\delta)(1-p_0-\delta)/n}} \sim N(0, 1) \]

\(H_1\) 成立时,可构建 \(z'\) 统计量:

\[ z' = \frac{\hat{p}-p_0-\delta}{\sqrt{(p_0+\delta)(1-p_0-\delta)/n}} \sim N\left(\frac{p-p_0-\delta}{\sqrt{(p_0+\delta)(1-p_0-\delta)/n}}, \frac{p(1-p)}{(p_0+\delta)(1-p_0-\delta)}\right) \]
\[ \begin{align} \text{Power} & = \operatorname{Pr}(z' > z_{1-\alpha}) \\ & = 1 - \Phi\left(\frac{z_{1-\alpha} - \frac{p-p_0-\delta}{\sqrt{(p_0+\delta)(1-p_0-\delta)/n}}}{\frac{\sqrt{p(1-p)}}{\sqrt{(p_0+\delta)(1-p_0-\delta)}}}\right) \\ & = 1 - \Phi\left(\frac{z_{1-\alpha}\sqrt{(p_0+\delta)(1-p_0-\delta)} - (p-p_0-\delta)\sqrt{n}}{\sqrt{p(1-p)}}\right) \end{align} \]
\[ \begin{align} \text{Power} & = \operatorname{Pr}(z' < z_{\alpha}) \\ & = \Phi\left(\frac{z_{\alpha} - \frac{p-p_0-\delta}{\sqrt{(p_0+\delta)(1-p_0-\delta)/n}}}{\frac{\sqrt{p(1-p)}}{\sqrt{(p_0+\delta)(1-p_0-\delta)}}}\right) \\ & = 1 - \Phi\left(\frac{z_{1-\alpha}\sqrt{(p_0+\delta)(1-p_0-\delta)} + (p-p_0-\delta)\sqrt{n}}{\sqrt{p(1-p)}}\right) \end{align} \]
样本量公式推导

根据标准正态分布分位数的定义:

\[ \frac{z_{1-\alpha}\sqrt{(p_0+\delta)(1-p_0-\delta)} \pm (p-p_0-\delta)\sqrt{n}}{\sqrt{p(1-p)}} = z_{\beta} \]

可解出:

\[ n = \frac{\left[z_{1-\alpha}\sqrt{(p_0+\delta)(1-p_0-\delta)} + z_{1-\beta}\sqrt{p(1-p)}\right]^2}{\left(p-p_0-\delta\right)^2} \]

z-test using s(p0) with continuity correction

z-test using s(p0) 的基础上加入校正项 \(c\)

\[ c = \begin{cases} - \frac{1}{2n} & , \text{if } p \gt p_0 + \delta \\ \frac{1}{2n} & , \text{if } p \lt p_0 + \delta \\ 0 & , \text{if } \left| p - p_0 - \delta \right| \leqslant \frac{1}{2n} \end{cases} \]

\(H_0\) 成立时,可构建 \(z\) 统计量:

\[ z = \frac{\hat{p}-p_0-\delta+c}{\sqrt{(p_0+\delta)(1-p_0-\delta)/n}} \sim N(0, 1) \]

\(H_1\) 成立时,可构建 \(z'\) 统计量:

\[ z' = \frac{\hat{p}-p_0-\delta+c}{\sqrt{(p_0+\delta)(1-p_0-\delta)/n}} \sim N\left(\frac{p-p_0-\delta+c}{\sqrt{(p_0+\delta)(1-p_0-\delta)/n}}, \frac{p(1-p)}{(p_0+\delta)(1-p_0-\delta)}\right) \]
\[ \begin{align} \text{Power} & = \operatorname{Pr}(z' > z_{1-\alpha}) \\ & = 1 - \Phi\left(\frac{z_{1-\alpha} - \frac{p-p_0-\delta+c}{\sqrt{(p_0+\delta)(1-p_0-\delta)/n}}}{\frac{\sqrt{p(1-p)}}{\sqrt{(p_0+\delta)(1-p_0-\delta)}}}\right) \\ & = 1 - \Phi\left(\frac{z_{1-\alpha}\sqrt{(p_0+\delta)(1-p_0-\delta)} - (p-p_0-\delta+c)\sqrt{n}}{\sqrt{p(1-p)}}\right) \end{align} \]
\[ \begin{align} \text{Power} & = \operatorname{Pr}(z' < z_{\alpha}) \\ & = \Phi\left(\frac{z_{\alpha} - \frac{p-p_0-\delta+c}{\sqrt{(p_0+\delta)(1-p_0-\delta)/n}}}{\frac{\sqrt{p(1-p)}}{\sqrt{(p_0+\delta)(1-p_0-\delta)}}}\right) \\ & = 1 - \Phi\left(\frac{z_{1-\alpha}\sqrt{(p_0+\delta)(1-p_0-\delta)} + (p-p_0-\delta+c)\sqrt{n}}{\sqrt{p(1-p)}}\right) \end{align} \]

z-test using s(phat)

\(H_0\) 成立时,可构建 \(z\) 统计量:

\[ z = \frac{\hat{p}-p_0-\delta}{\sqrt{p(1-p)/n}} \sim N(0, 1) \]

\(H_1\) 成立时,可构建 \(z'\) 统计量:

\[ z' = \frac{\hat{p}-p_0-\delta}{\sqrt{p(1-p)/n}} \sim N\left(\frac{p-p_0-\delta}{\sqrt{p(1-p)/n}}, 1\right) \]
\[ \text{Power} = \operatorname{Pr}(z' > z_{1-\alpha}) = 1 - \Phi\left(z_{1-\alpha} - \frac{p-p_0-\delta}{\sqrt{p(1-p)/n}}\right) \]
\[ \text{Power} = \operatorname{Pr}(z' < z_{\alpha}) = \Phi\left(z_{\alpha} - \frac{p-p_0-\delta}{\sqrt{p(1-p)/n}}\right) = 1 - \Phi\left(z_{1-\alpha} + \frac{p-p_0-\delta}{\sqrt{p(1-p)/n}}\right) \]
样本量公式推导

根据标准正态分布分位数的定义:

\[ z_{1-\alpha} \pm \frac{p-p_0-\delta}{\sqrt{p(1-p)/n}} = z_{\beta} \]

可解出:

\[ n = \frac{\left(z_{1-\alpha}+z_{1-\beta}\right)^2 p(1-p)}{\left(p-p_0-\delta\right)^2} \]

z-test using s(phat) with continuity correction

z-test using s(phat) 的基础上加入校正项 \(c\)

\[ c = \begin{cases} - \frac{1}{2n} & , \text{if } p \gt p_0 + \delta \\ \frac{1}{2n} & , \text{if } p \lt p_0 + \delta \\ 0 & , \text{if } \left| p - p_0 - \delta \right| \leqslant \frac{1}{2n} \end{cases} \]

\(H_0\) 成立时,构建 \(z\) 统计量:

\[ z = \frac{\hat{p}-p_0-\delta+c}{\sqrt{p(1-p)/n}} \sim N(0, 1) \]

\(H_1\) 成立时,构建 \(z'\) 统计量:

\[ z' = \frac{\hat{p}-p_0-\delta+c}{\sqrt{p(1-p)/n}} \sim N\left(\frac{p-p_0-\delta+c}{\sqrt{p(1-p)/n}}, 1\right) \]
\[ \text{Power} = \operatorname{Pr}(z' > z_{1-\alpha}) = 1 - \Phi\left(z_{1-\alpha} - \frac{p-p_0-\delta+c}{\sqrt{p(1-p)/n}}\right) \]
\[ \text{Power} = \operatorname{Pr}(z' < z_{\alpha}) = \Phi\left(z_{\alpha} - \frac{p-p_0-\delta+c}{\sqrt{p(1-p)/n}}\right) = 1 - \Phi\left(z_{1-\alpha} + \frac{p-p_0-\delta+c}{\sqrt{p(1-p)/n}}\right) \]