§7.1 假设检验的基本思想与概念
依赖于
被以下题目直接调用
正文部分
§7.1 假设检验的基本思想与概念
1. 假设
- 参数空间 Θ={θ} 的非空子集或有关参数 θ 的命题,称为统计假设,简称假设;
- 原假设,根据需要而设立的假设,常记为 H0:θ∈Θ0;
- 备择假设,在原假设被拒绝后而采用(接受)的假设,常记为 H1:θ∈Θ1。
注:要求 Θ0∩Θ1=∅,即原假设 H0 与备择假设 H1 不含公共参数。若假设 H0(或 H1)只含一个点,则称该假设是简单假设,否则称为复杂假设。
2. 检验 对原假设 H0:θ∈Θ0 作出是否拒绝 H0 的判断的法则称为检验法则,简称检验。检验有两个结果:
- “原假设不正确”,称为拒绝原假设,或称检验显著;
- “原假设正确”,称为接受(保留)原假设,或称检验不显著。
统计假设检验的着力点不在于说明原假设正确,而在于说明原假设不正确,因为用一个样本去说明原假设正确是根据不足的,而用一个样本去推翻原假设却是合理的,因此统计学在讨论假设检验时着力点在于建立检验的拒绝域。
3. 检验问题
- 由原假设 H0 和备择假设 H1 组成一个需要作判断的问题称为检验问题;
- 参数假设检验问题,两个假设都是有关参数的命题组成的检验问题;
- 非参数假设检验问题,两个假设都是有关分布的命题组成的检验问题。
常用的参数假设检验问题有如下三个,其中 θ0 是已知常数:
(1) H0:θ≤θ0(2) H0:θ≥θ0(3) H0:θ=θ0vsH1:θ>θ0;vsH1:θ<θ0;vsH1:θ=θ0.
其中 (1) 与 (2) 又称单侧检验问题,因为一个假设位于另一个假设的一侧;(3) 称为双侧检验问题,因为备择假设位于原假设的两侧。
4. 两类错误及其发生概率
- 原假设 H0 正确,但被拒绝,这种判断错误称为第一类错误,其发生概率称为犯第一类错误的概率,或称拒真概率,常记为 α;
- 原假设 H0 不真,但被接受,这种判断错误称为第二类错误,其发生概率称为犯第二类错误的概率,或称取伪概率,常记为 β。
5. 假设检验的基本步骤
(1) 建立假设,根据要求建立原假设 H0 和备择假设 H1.
(2) 选择检验统计量,给出拒绝域 W 的形式。
- 用于对原假设 H0 作出判断的统计量称为检验统计量;
- 使原假设被拒绝的样本观测值所在区域称为拒绝域,常用 W 表示;
- 一个拒绝域 W 唯一确定一个检验法则。
(3) 选择显著性水平 α (0<α<1).
(4) 给出拒绝域,由概率等式 P(W)=α 确定具体的拒绝域。
(5) 作出判断。
- 当样本 (x1,x2,…,xn)∈W,则拒绝 H0;
- 当样本 (x1,x2,…,xn)∈W,则接受 H0。
6. 势函数 设检验问题
H0:θ∈Θ0vsH1:θ∈Θ1
的拒绝域为 W,则样本观测值 x1,x2,…,xn 落在拒绝域 W 内的概率称为该检验的势函数,记为
g(θ)=Pθ((x1,x2,…,xn)∈W),θ∈Θ0∪Θ1.
由势函数 g(θ) 容易得到犯两类错误的概率
g(θ)={α(θ),1−β(θ),θ∈Θ0,θ∈Θ1.
7. 水平为 α 的检验 对检验问题
H0:θ∈Θ0vsH1:θ∈Θ1,
若一个检验犯第一类错误的概率 α(θ) 不超过事先给定的显著性水平 α,即
α(θ)≤α,θ∈Θ0,
则称该检验为显著性水平为 α 的显著性检验。
在实际使用中 α 不宜选得过小,α 过小会导致 β 过大,应在适当控制 α 中制约 β。最常用的选择是 α=0.05,有时也选用 α=0.10 或 α=0.01。
8. 检验的 p 值 在一个假设检验问题中,利用样本观察值能够做出拒绝原假设的最小显著性水平称为该检验的 p 值,引入检验的 p 值的好处是:
- 它比较客观,避免了事先确定显著性水平;
- 由检验的 p 值与人们心目中的显著性水平 α 进行比较:
- 若 α≥p,则在显著性水平 α 下拒绝 H0;
- 若 α<p,则在显著性水平 α 下接受 H0;
- 检验的 p 值的计算是复杂的,会涉及各种抽样分布。如今统计软件都有计算 p 值的功能,因此这对使用者反而方便,它不需要各用各种抽样分布的分位数表,而只需要观察计算机输出的 p 值多少就可以做出判断。
习题与解答 7.1
设 x1,x2,…,xn 是来自 N(μ,1) 的样本,考虑假设检验问题
H0:μ=2,H1:μ=3,
若检验由拒绝域
W={xˉ≥2.6}
确定。
(1) 当 n=20 时求检验犯两类错误的概率;
(2) 如果要使得检验犯第二类错误的概率 β≤0.01,n 最小应取多少?
(3) 证明:当 n→∞ 时,α→0,β→0。
解
(1) 由定义知,犯第一类错误的概率为
α=P(xˉ≥2.6∣H0)=P(1/20xˉ−2≥1/202.6−2)=1−Φ(2.68)=0.0037,
这是因为在 H0 成立下,xˉ∼N(2,1/20)。而犯第二类错误的概率为
β=P(xˉ<2.6∣H1)=P(1/20xˉ−3<1/202.6−3)=Φ(−1.79)=1−Φ(1.79)=0.0367,
这是因为在 H1 成立下,xˉ∼N(3,1/20)。
(2) 若使犯第二类错误的概率满足
β=P(xˉ<2.6∣H1)=P(1/nxˉ−3<1/n2.6−3)≤0.01,
即
1−Φ(0.4n)≤0.01,
或
Φ(0.4n)≥0.99.
查表得 0.4n≥2.33,由此给出
n≥33.93,
因而 n 最小应取 34,才能使检验犯第二类错误的概率 β≤0.01。
(3) 在样本量为 n 时,检验犯第一类错误的概率为
α=P(xˉ≥2.6∣H0)=P(1/nxˉ−2≥1/n2.6−2)=1−Φ(0.6n),
当 n→∞ 时,Φ(0.6n)→1,即 α→0。
检验犯第二类错误的概率为
β=P(xˉ<2.6∣H1)=P(1/nxˉ−3<1/n2.6−3)=Φ(−0.4n)=1−Φ(0.4n),
当 n→∞ 时,Φ(0.4n)→1,即 β→0。
注:从这个例子可以看出,要使得 α 与 β 都趋于 0,必须 n→∞ 才可实现,这一结论在一般场合仍成立,即要使得 α 与 β 同时很小,必须样本量 n 很大。由于样本量 n 很大在实际中常常是不可行的,故一般情况下人们不应要求 α 与 β 同时很小。
设 x1,x2,…,x10 是来自 0-1 总体 b(1,p) 的样本,考虑检验问题
H0:p=0.2vsH1:p=0.4,
取拒绝域为
W={xˉ≥0.5},
求该检验犯两类错误的概率。
解
x1,x2,…,x10∼b(1,p),则 10xˉ∼b(10,p),于是犯两类错误的概率分别为
α=P(xˉ≥0.5∣H0)=P(10xˉ≥5∣H0)=k=5∑10(k10)(51)k(54)10−k=0.0328,
β=P(xˉ<0.5∣H1)=P(10xˉ<5∣H1)=k=0∑4(k10)(52)k(53)10−k=0.6331.
讨论:这里 α=0.0328 已经很小了,但 β=0.6331 却很大。在样本量 n=10 固定下,若要使 α 更小则会导致 β 更大,为说明这一点,我们试着将拒绝域改变为
W={xˉ≥0.6},
则这时检验犯两类错误的概率分别为
α=P(xˉ≥0.6∣H0)=P(10xˉ≥6∣H0)=k=6∑10(k10)(51)k(54)10−k
=0.0328−(510)(51)5(54)5=0.0328−0.0264=0.0064,
β=P(xˉ<0.6∣H1)=P(10xˉ<6∣H1)=k=0∑5(k10)(52)k(53)10−k
=0.6331+(510)(52)5(53)5=0.6331+0.2007=0.8338.
这一现象在一般场合也是对的,即在样本量 n 固定下,减小 α 必导致增大 β,减小 β 也必导致增大 α。
设 x1,x2,…,x16 是来自正态总体 N(μ,4) 的样本,考虑检验问题
H0:μ=6vsH1:μ=6,
拒绝域取为
W={∣xˉ−6∣≥c},
试求 c 使得检验的显著性水平为 0.05,并求该检验在 μ=6.5 处犯第二类错误的概率。
解
在 H0 为真的条件下,
xˉ∼N(6,41),
因而由
P(∣xˉ−6∣≥c∣μ=6)=0.05,
得
P(0.5xˉ−6≥0.5c)=1−Φ(2c)=0.025,
也就是说
Φ(2c)=0.975,
2c=1.96,所以当 c=0.98 时,检验的显著性水平为 0.05。
该检验在 μ=6.5 处犯第二类错误的概率为
β=P(∣xˉ−6∣<0.98∣μ=6.5)
=P(−2×1.48<0.5xˉ−6.5<2×0.48)
=Φ(0.96)−Φ(−2.96)=Φ(0.96)+Φ(2.96)−1=0.83.
设总体为均匀分布 U(0,θ),x1,x2,…,xn 是样本,考虑检验问题
H0:θ≥3vsH1:θ<3,
拒绝域取为
W={x(n)≤2.5},
求检验犯第一类错误的最大值 α,若要使得该最大值 α 不超过 0.05,n 至少应取多大?
解
均匀分布 U(0,θ) 的最大次序统计量 x(n) 的密度函数为
fn(x)=⎩⎨⎧θnnxn−1,0,0<x<θ,其他,
因而检验犯第一类错误的概率为
α(θ)=P(x(n)≤2.5∣H0)=∫02.5nθnxn−1dx=(θ2.5)n.
它是 θ 的严格单调递减函数,故其最大值在 θ=3 处达到,即
α=α(3)=(32.5)n.
若要使得 α(3)≤0.05,则要求 nln(2.5/3)≤ln0.05,这给出 n≥16.43,即 n 至少为 17。
在假设检验问题中,若检验结果是接受原假设,则检验可能犯哪一类错误?若检验结果是拒绝原假设,则又有可能犯哪一类错误?
解
若检验结果是接受原假设,可能有两种情况:其一是原假设为真,此时检验是正确的,未犯错误;其二是原假设不真,此时检验结果就错了,这种错误是接受了不真的原假设,为第二类错误,故此时检验可能犯第二类错误。
若检验结果是拒绝原假设,也可能有两种情况:若原假设本身不真,检验是正确的;若原假设事实上是真的,则检验就犯了第一类错误,由此,在此种场合,检验可能会犯第一类错误。
设 x1,x2,…,x20 是来自 0-1 总体 b(1,p) 的样本,考虑检验问题
H0:p=0.2vsH1:p=0.2,
取拒绝域为
W={i=1∑20xi≥7 或 i=1∑20xi≤1},
- 求 p=0,0.1,0.2,…,0.9,1 时的势并由此画出势函数的图;
- 求在 p=0.05 时犯第二类错误的概率。
解
(1) 势函数的计算公式为
g(p)=Pp(i=1∑20xi≥7 或 i=1∑20xi≤1)=i=0∑1(i20)pi(1−p)20−i+i=7∑20(i20)pi(1−p)20−i.
则 p=0,0.1,0.2,…,0.9,1 时的势计算如下表:
pg(p)01.000.10.390.20.160.30.400.40.750.50.940.60.990.71.000.81.000.91.0011.00
可用软件计算,如 MATLAB 语句为 1 - binocdf(6,20,p) + binocdf(1,20,p);势函数图如图 7.1,它在 p=0.2 处达到最小。
\FigureSevenOne
(2) p=0.05 时,犯第二类错误的概率可采用如下 MATLAB 语句 binocdf(6,20,0.05) - binocdf(1,20,0.05) 计算给出 1−g(0.05),计算结果为 0.2641。
设一个单一观测的样本 x 取自密度函数为 p(x) 的总体,对 p(x) 考虑统计假设
H0:p0(x)=I[0,1](x)vsH1:p1(x)=2xI[0,1](x),
若其拒绝域的形式为 W={x∣x≥c},试确定一个 c,使得犯第一类、第二类错误的概率满足 α+2β 为最小,并求其最小值。
解
由 α=P0(x≥c)=1−c,β=1−P1(x≥c)=c2,可得
α+2β=1−c+2c2.
因此,当 c=41 时,α+2β 为最小,并且此时的最小值为 87。
设 x1,x2,…,x30 为取自泊松分布 P(λ) 的随机样本。
- 试给出单边假设检验问题
H0:λ≤0.1vsH1:λ>0.1
的显著性水平 α=0.05 的检验;
- 求此检验的势函数 g(λ) 在 λ=0.05,0.2,0.3,…,0.9 时的值,并据此画出 g(λ) 的图像。
解
(1) 选 ∑i=130xi 为检验统计量,其值愈大愈倾向于拒绝 H0,所以,该检验问题的拒绝域形式为
W={i=1∑30xi>c}.
注意到在 λ=0.1 时 ∑i=130xi∼P(3),从而第一类错误概率为
P0.1(i=1∑30xi>c)=1−i=0∑ci!3ie−3.
当 c=5 时,
P0.1(i=1∑30xi>5)=0.0839,
当 c=6 时,
P0.1(i=1∑30xi>6)=0.0335.
因此,该检验问题的拒绝域为
W={i=1∑30xi>6}.
(2) 势函数的计算公式为
g(λ)=Pλ(i=1∑30xi>6)=1−i=0∑6i!(30λ)ie−30λ,
则 λ=0.05,0.2,…,0.9 时的势计算如下表:
λg(λ)0.059.3×10−40.20.390.30.790.40.950.50.990.61.000.71.000.81.000.91.00
势函数图如图 7.2。
\FigureSevenTwo
设正态总体的方差 σ2 已知,均值 μ 只能取 μ0 或 μ1(μ1<μ0)两值之一,xˉ 为总体的容量为 n 的样本均值,考虑检验问题
H0:μ=μ0vsH1:μ=μ1,
若检验拒绝域取为
W={xˉ≤μ0+uασ/n},
则检验犯第二类错误的概率为
β=Pμ1(xˉ>μ0+uασ/n).
- 试验证:
β=Φ(u1−α+σ/nμ1−μ0),
从而在 α,β 给定时,有
n=(u1−α+u1−β)2(μ1−μ0)2σ2;
- 若 n 固定,当 α 减小时 β 怎样变化?当 β 减小时 α 怎样变化?
- 当 σ=0.12,μ1=μ0−0.02,并且要求 α≤0.05,β≤0.025 时,样本容量 n 至少应为多少?
解
(1) 由于 xˉ∼N(μ,σ2/n),故检验犯第二类错误的概率为
β=Pμ1(xˉ−μ1>μ0+σuα/n−μ1)=Pμ1(σ/nxˉ−μ1>uα+σ/nμ0−μ1)=1−Φ(uα+σ/nμ0−μ1)=Φ(−uα−σ/nμ0−μ1)=Φ(u1−α+σ/nμ1−μ0).
这给出
u1−α+σ/nμ1−μ0=uβ,
也即
u1−α+u1−β=σ/nμ1−μ0.
从而在 α,β 给定时,有
n=(u1−α+u1−β)2(μ1−μ0)2σ2.
(2) 若 n 固定,当 α 减小时,u1−α 就变大,由
u1−α+u1−β=σ/nμ1−μ0
为常量可知 u1−β 就变小,从而导致 β 增大。
同理可知:当 β 减小时 α 增大。
这说明,在样本量给定时,犯两类错误的概率一个变小另一个就会变大,不可能找到一个使得犯两类错误的概率都变小的检验方案。
(3) 由 α≤0.05,β≤0.025 查表可得 u1−α≥1.645,u1−β≥1.96,于是
n=(u1−α+u1−β)2(μ1−μ0)2σ2≥(μ1−μ0)2(1.645+1.96)2σ2.
将 σ=0.12,μ1=μ0−0.02 代入,有
n≥0.022(1.645+1.96)2⋅0.122=467.86,
即 n 至少应为 468。
假定 X 是连续随机变量,x 是 X 的一次观测值。关于总体密度函数 f(x) 有如下两个假设:
H0:f(x)={1/2,0,0≤x≤2,其他,vsH1:f(x)={x/2,0,0≤x≤2,其他.
检验的判断规则是:若 x≥2/3 则拒绝原假设 H0,试求检验犯两类错误的概率。
解
由所给条件,犯第一类错误的概率为
α=P(x≥32∣H0)=∫2/3221dx=32,
犯第二类错误的概率为
β=P(x<32∣H1)=∫02/32xdx=41×(32)2=91.
这个检验犯两类错误的概率都不小,不是一个好的检验,主要原因是样本量太小。
设 x1,x2,…,xn 是来自 U(0,θ) 的一个样本,对如下的检验问题
H0:θ≤21vsH1:θ>21
已给出拒绝域
W={x(n)≥c},
其中 x(n) 为样本的最大次序统计量。
- 求此检验的势函数;
- 若要求检验犯第一类错误概率不超过 0.05(即 α(θ)≤0.05),如何确定 c?
- 若在(2)的要求下进一步要求检验在 θ=43 处犯第二类错误的概率不超过 0.02(即 β(θ)≤0.02),n 至少要取多少?
- 如今 n=20,x(20)=0.48,对此检验问题作出判断。
解
(1) 此检验的势函数为
g(θ)=P(x(n)≥c)=1−P(x(n)<c)=1−P(x1<c,x2<c,…,xn<c)=⎩⎨⎧0,1−(θc)n,θ≤c,θ>c.
可见,在 θ>c 时,势函数 g(θ) 是 θ 的严格增函数。
(2) 在 H0 成立下,犯第一类错误的概率为 α(θ)=g(θ),故由题意知,应有
g(θ)=1−(θc)n≤0.05,θ≤21.
由于 g(θ) 是增函数,故 g(θ) 在 θ=21 处达到最大值,故只要使
g(21)=1−(2c)n=0.05
即可实现,由此解出
c=21(0.95)1/n.
譬如,在 n=5 时,c=0.4949;在 n=10 时,c=0.4974。
(3) 在备择假设 H1 成立下,犯第二类错误的概率为
β(θ)=1−g(θ)=(θc)n,θ>21.
现要求在 θ=43 处有 β(θ)≤0.02,即
(3/4c)n≤0.02.
若把(2)中的
c=21(0.95)1/n
代入,可得
n≥ln3−ln2ln95−ln2=9.52.
可见,若取 n=10 即可使 θ=43 处犯第二类错误的概率不超过 0.02。
(4) 若样本量 n=20,则其拒绝域为
W={x(n)≥c0},其中c0=21(0.95)1/20=0.4987.
如今 x(n)=0.48<c0,故应接受原假设 H0:θ≤21。
设 x1,x2,…,xn 是来自泊松分布 P(λ) 的一个样本。
- 利用泊松分布的充分统计量对检验问题
H0:λ=λ0vsH1:λ=λ1 (λ1>λ0)
在显著性水平 α 时给出其拒绝域;
- 证明(1)中的拒绝域也是检验问题
H0:λ≤λ0vsH1:λ>λ0
的显著性水平为 α 的显著性检验的拒绝域;
- 当样本量 n 较大时,由中心极限定理给出近似的拒绝域。
解
(1) 泊松分布 P(λ) 的充分统计量是 xˉ,它是 λ 的无偏估计。若原假设 H0:λ=λ0 成立,则 xˉ 不应该很大,因此,当 xˉ 较大时,就应该拒绝原假设 λ=λ0,所以此检验的拒绝域应有形式
W={xˉ≥c}={i=1∑nxi≥nc},
其中 c 应由给定的显著性水平 α 确定,即 c 由概率不等式
Pλ0(i=1∑nxi≥nc)≤α
或
Pλ0(i=1∑nxi<nc)≥1−α
确定。由于原假设成立下 xi∼P(λ0),故
i=1∑nxi∼P(nλ0).
若令泊松分布 P(nλ0) 的 1−α 分位数为 k1−α,则由 nc=k1−α 可得
c=nk1−α.
所以在给定 α 时,该检验的拒绝域为
W={xˉ≥nk1−α},
这里 k1−α 的寻求还不是一件易事。
(2) 若将上述拒绝域作为此处检验问题的拒绝域,我们只需要证明确该检验的势函数是单调增的即可说明它也是(2)的显著性水平为 α 的显著性检验。此处该检验的势函数为
g(λ)=Pλ(xˉ≥nk1−α)=1−Pλ(xˉ<nk1−α)=1−k=0∑mk!(nλ)ke−nλ,
其中 m 为如下整数:
m={k1−α−1,⌊k1−α⌋,k1−α 为整数,k1−α 为非整数.
考察 g(λ) 的单调性,为此求其导数
g′(λ)=−k=0∑mk!1[nk(nλ)k−1e−nλ−n(nλ)ke−nλ]=ne−nλ[k=0∑mk!(nλ)k−k=0∑m−1k!(nλ)k]=ne−nλm!(nλ)m>0.
所以势函数 g(λ) 是 λ 的严格增函数。由此可知,在原假设 H0:λ≤λ0 上,g(λ) 在 λ=λ0 处达到最大。
(3) 当样本量 n 较大时,由中心极限定理可得原假设成立时 xˉ 的渐近分布
λ0/nxˉ−λ0∼˙N(0,1).
对给定的显著性水平 α,有
P(λ0/nxˉ−λ0≥u1−α)=α,
P(xˉ≥λ0+u1−αnλ0)=α,
即拒绝域 W 中的临界值
c=λ0+u1−αnλ0.
譬如,λ0=1.5,n=10 和 α=0.05 时,有
c=1.5+1.645101.5=2.1371,
即当 n=10 时,若 xˉ≥2.1371,则应拒绝原假设 H0:λ≤1.5。
设 x1,x2,…,xn 是来自瑞利(Rayleigh)分布 Ra(θ) 的一个样本,瑞利分布的密度函数为
p(x)=⎩⎨⎧θ22xexp{−θ2x2},0,x>0,x≤0(θ>0).
- 求此分布的充分统计量;
- 利用充分统计量在给定显著性水平 α 下给出检验问题
H0:θ=1vsH1:θ>1
的拒绝域;
- 在样本量较大时,利用中心极限定理给出近似拒绝域。
解
(1) 样本的联合密度函数为
p(x1,x2,…,xn)=(i=1∏nθ22xi)exp{−θ21i=1∑nxi2}=exp{−θ21i=1∑nxi2−2nlnθ}i=1∏n2xi.
由因子分解定理知,θ 的充分统计量是
T=i=1∑nxi2.
(2) 注意到
E(x2)=∫0∞θ22x3exp{−θ2x2}dx=θ2∫0∞θ2x2exp{−θ2x2}d(θ2x2)=θ2∫0∞ye−ydy=θ2Γ(2)=θ2.
由此可见
n1i=1∑nxi2
是 θ2 的无偏估计。当
T=i=1∑nxi2
较大时,拒绝原假设 H0:θ=1 是合理的。故
H0:θ=1vsH1:θ>1
的拒绝域为
W={i=1∑nxi2≥c},
其中 c 由概率等式
P(i=1∑nxi2≥c)=α
确定。为了确定 c,需要充分统计量
T=i=1∑nxi2
的分布。
可以证明,当 X∼Ra(θ) 时,
X2∼Exp(θ21)=Ga(1,θ21),
由此可得
i=1∑nxi2∼Ga(n,θ21),
或者
θ22i=1∑nxi2∼Ga(n,21)=χ2(2n).
在原假设 H0:θ=1 成立下,有
2i=1∑nxi2∼χ2(2n),
利用 χ2 分布的分位数可确定临界值 c。由等式
P(i=1∑nxi2≥c)=α
可得
P(2i=1∑nxi2<2c)=1−α.
记 χ1−α2(2n) 是 χ2 分布的 1−α 分位数,可得
2c=χ1−α2(2n),或c=21χ1−α2(2n).
譬如,当 n=15,α=0.05 时,
χ0.952(30)=43.773,
所以
c=21.887,
即当检验统计量
i=1∑nxi2≥21.887
时,将拒绝原假设 H0:θ=1。
(3) 由
X2∼Exp(θ21)
可知
E(X2)=θ2,\Var(X2)=θ4,
从而有
E(i=1∑nxi2)=nθ2,\Var(i=1∑nxi2)=nθ4.
在原假设 H0:θ=1 成立下,有
E(i=1∑nxi2)=n,\Var(i=1∑nxi2)=n.
这里
i=1∑nxi2
可看作 n 个相互独立同分布随机变量之和,故由中心极限定理知
n∑i=1nxi2−n∼˙N(0,1),
从而有
P(i=1∑nxi2≥c)=P(n∑i=1nxi2−n≥nc−n)=1−Φ(nc−n).
故由等式
P(i=1∑nxi2≥c)=α
可得
Φ(nc−n)=1−α.
记 u1−α 为标准正态分布的 1−α 分位数,则有
nc−n=u1−α,
即
c=n+u1−αn.
若 n=15,α=0.05,查表得 u1−α=u0.95=1.645,从而
c=15+1.64515=21.371.
这个结果与(2)定出的精确值较为接近。
评论
支持 Markdown 和 LaTeX 数学公式。