Table of Contents
概率论与数理统计(Probability & Statistics I)
概率论与数理统计(Probability & Statistics II)
样本与抽样分布(Sample and Sampling Distribution)
数理统计(Statistics)
数理统计是以概率论为基础, 关于实验数据的收集、整理、分析与推断的一门科学与艺术。
总体(population):研究对象的全体;
个体(individual):总体中的每一个具体对象称为个体;
总体的容量(capacity):总体中包含的个体数;
有限总体(finite population):容量有限的总体;
无限总体(infinite population):容量无限的总体。通常将容量非常大的有限总体也按无限总体处理。
为了采用数理统计方法进行分析,首先要收集数据,数据收集方法一般有两种。
(1)通过调查、记录收集数据。(2)通过实验收集数据。
实际中人们通常只关注总体的某个(或几个)指标。
- 总体的某个指标 X X X, 对于不同的个体来说有不同的取值, 这些取值构成一个分布, 因此 X X X可以看成一个随机变量。
- 有时候直接将 X X X称为总体. 假设 X X X的分布函数为 F ( x ) F(x) F(x), 也称 X X X服从 F ( x ) F(x) F(x)分布 X ∼ F ( x ) X∼ F(x) X∼F(x)
如何推断总体分布的未知参数(或分布)?
需要从总体中抽取一部分个体, 根据这部分个体的数据,并利用概率论的知识等作出分析推断.被抽取的部分个体叫做总体的一个样本(sample)。
简单随机样本(simple random sample):满足以下两个条件的随机样本
(
X
1
,
X
2
,
…
,
X
n
)
(X_1,X_2,…,X_n)
(X1,X2,…,Xn) 称为容量是
n
n
n的简单随机样本,简称样本。
(1) 代表性:每个
X
i
X_i
Xi与
X
X
X同分布;
(2) 独立性:
(
X
1
,
X
2
,
…
,
X
n
)
(X_1,X_2,…,X_n)
(X1,X2,…,Xn)是相互独立的随机变量。
获得简单随机样本的抽样称为简单随机抽样。
- 对于有限总体,采用放回抽样.
- 但当总体容量很大的时候,放回抽样有时候很不方便, 因此在实际中当总体容量比较大时,通常将不放回抽样所得到的样本近似当作简单随机样本来处理.
- 对于无限总体, 一般采取不放回抽样.
对样本
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn 进行观测后,得到的数值:
x
1
,
x
2
,
…
,
x
n
x_1,x_2,…,x_n
x1,x2,…,xn 称为样本观测值(observed value)
观测前:
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn 是随机变量
观测后:
x
1
,
x
2
,
…
,
x
n
x_1,x_2,…,x_n
x1,x2,…,xn 是具体的数据
样本的联合分布:设总体
X
∼
F
(
x
)
X∼ F(x)
X∼F(x),则样本
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn 的联合分布函数为:
F
(
x
1
,
x
2
,
…
,
x
n
)
=
P
{
X
1
⩽
x
1
,
…
,
X
n
⩽
x
n
}
=
∏
i
=
1
n
F
(
x
i
)
F(x_1,x_2,…,x_n)=P\{X_1⩽ x_1,…,X_n⩽ x_n\}=\displaystyle\prod_{i=1}^{n}F(x_i)
F(x1,x2,…,xn)=P{X1⩽x1,…,Xn⩽xn}=i=1∏nF(xi)
若总体
X
X
X 的密度函数为
f
(
x
)
f(x)
f(x),则样本
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn 的联合密度函数为:
f
(
x
1
,
x
2
,
…
,
x
n
)
=
∏
i
=
1
n
f
(
x
i
)
f(x_1,x_2,…,x_n)=\displaystyle\prod_{i=1}^{n}f(x_i)
f(x1,x2,…,xn)=i=1∏nf(xi)
若总体
X
X
X 的分布律为
P
{
X
=
a
k
}
=
p
k
P\{X=a_k\}=p_k
P{X=ak}=pk,则样本
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn 的联合分布律为:
P
{
X
1
=
x
1
,
…
,
X
n
=
x
n
}
=
∏
i
=
1
n
P
{
X
=
x
i
}
P\{X_1=x_1,…,X_n=x_n\}=\displaystyle\prod_{i=1}^{n}P\{X=x_i\}
P{X1=x1,…,Xn=xn}=i=1∏nP{X=xi}
抽样分布(Sampling distribution)
从样本中提取有用的信息来研究总体的分布及各种特征数——构造统计量.
统计量(statistic):设
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn 是来自总体
X
∼
F
(
x
)
X∼ F(x)
X∼F(x) 的样本,
g
(
x
1
,
x
2
,
…
,
x
n
)
g(x_1,x_2,…,x_n)
g(x1,x2,…,xn) 是
n
n
n元实值连续函数,若函数
g
(
x
1
,
x
2
,
…
,
x
n
)
g(x_1,x_2,…,x_n)
g(x1,x2,…,xn) 不含未知参数, 则称随机变量
g
(
X
1
,
X
2
,
…
,
X
n
)
g(X_1,X_2,…,X_n)
g(X1,X2,…,Xn) 为统计量。
常用统计量
(1) 样本均值:
X
ˉ
=
1
n
∑
i
=
1
n
X
i
\bar X=\dfrac{1}{n}\displaystyle\sum_{i=1}^{n}X_i
Xˉ=n1i=1∑nXi
(2) 样本方差:
S
2
=
1
n
−
1
∑
i
=
1
n
(
X
i
−
X
ˉ
)
2
S^2=\dfrac{1}{n-1}\displaystyle\sum_{i=1}^{n}(X_i-\bar X)^2
S2=n−11i=1∑n(Xi−Xˉ)2
样本标准差:
S
=
S
2
S=\sqrt{S^2}
S=S2
(3) 样本 k 阶原点矩:
A
k
=
1
n
∑
i
=
1
n
X
i
k
A_k=\dfrac{1}{n}\displaystyle\sum_{i=1}^{n}X_i^k
Ak=n1i=1∑nXik
样本 k 阶中心矩:
B
k
=
1
n
∑
i
=
1
n
(
X
i
−
X
ˉ
)
k
B_k=\dfrac{1}{n}\displaystyle\sum_{i=1}^{n}(X_i-\bar X)^k
Bk=n1i=1∑n(Xi−Xˉ)k
(4) 将样本按由小到大的顺序排成
X
(
1
)
⩽
X
(
2
)
⩽
…
⩽
X
(
n
)
X_{(1)}⩽ X_{(2)}⩽ …⩽ X_{(n)}
X(1)⩽X(2)⩽…⩽X(n)
顺序统计量(order statistic):
X
(
1
)
,
X
(
2
)
,
…
,
X
(
n
)
X_{(1)}, X_{(2)}, …, X_{(n)}
X(1),X(2),…,X(n)
样本极小值(minimum):
X
(
1
)
=
min
{
X
(
1
)
,
X
(
2
)
,
…
,
X
(
n
)
}
X_{(1)}=\min \{X_{(1)}, X_{(2)}, …, X_{(n)}\}
X(1)=min{X(1),X(2),…,X(n)}
样本极大值(maximum):
X
(
n
)
=
max
{
X
(
1
)
,
X
(
2
)
,
…
,
X
(
n
)
}
X_{(n)}=\max \{X_{(1)}, X_{(2)}, …, X_{(n)}\}
X(n)=max{X(1),X(2),…,X(n)}
样本极差(range):
R
n
=
X
(
n
)
−
X
(
1
)
R_n=X_{(n)}-X_{(1)}
Rn=X(n)−X(1)
样本均值与样本方差的数字特征:设
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn 是来总体
X
X
X 的样本,且总体的均值与方差存在,记为
E
(
X
)
=
μ
,
D
(
X
)
=
σ
2
E(X)=μ,D(X)=σ^2
E(X)=μ,D(X)=σ2,则有
(1)
E
(
X
ˉ
)
=
E
(
1
n
∑
i
=
1
n
X
i
)
=
1
n
∑
i
=
1
n
E
(
X
i
)
=
μ
E(\bar X)=E(\dfrac{1}{n}\displaystyle\sum_{i=1}^{n}X_i)=\dfrac{1}{n}\displaystyle\sum_{i=1}^{n}E(X_i)=μ
E(Xˉ)=E(n1i=1∑nXi)=n1i=1∑nE(Xi)=μ
D
(
X
ˉ
)
=
D
(
1
n
∑
i
=
1
n
X
i
)
=
1
n
2
∑
i
=
1
n
D
(
X
i
)
=
σ
2
/
n
D(\bar X)=D(\dfrac{1}{n}\displaystyle\sum_{i=1}^{n}X_i)=\dfrac{1}{n^2}\displaystyle\sum_{i=1}^{n}D(X_i)=σ^2/n
D(Xˉ)=D(n1i=1∑nXi)=n21i=1∑nD(Xi)=σ2/n
(2)
E
(
S
2
)
=
σ
2
E(S^2)=σ^2
E(S2)=σ2
(3) 若
X
∼
N
(
μ
,
σ
2
)
X∼ N(μ,σ^2)
X∼N(μ,σ2),则
D
(
S
2
)
=
2
σ
4
n
−
1
D(S^2)=\dfrac{2σ^4}{n-1}
D(S2)=n−12σ4
抽样分布(sampling distribution):统计量的分布被称为抽样分布
当总体
X
X
X服从一般分布(如指数分布、均匀分布等),要得出统计量的分布是很困难的。当总体
X
X
X服从正态分布时,统计量
X
ˉ
,
S
2
\bar X,S^2
Xˉ,S2是可以计算的,本节主要介绍与标准正态总体相关的抽样分布。
χ
2
χ^2
χ2 分布(chi-square distribution):设
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn 相互独立,且都服从标准正态分布
N
(
0
,
1
)
N(0, 1)
N(0,1) ,则称随机变量
χ
2
=
∑
i
=
1
n
X
i
2
χ^2=\displaystyle\sum_{i=1}^{n}X_i^2
χ2=i=1∑nXi2 服从自由度为
n
n
n 的
χ
2
χ^2
χ2 分布,记为
χ
2
∼
χ
2
(
n
)
χ^2∼ χ^2(n)
χ2∼χ2(n),自由度指包含的独立变量的个数。
χ
2
χ^2
χ2 分布的概率密度为:
f
n
(
x
)
=
{
1
2
Γ
(
n
/
2
)
(
x
2
)
n
/
2
−
1
e
−
x
/
2
,
x
>
0
0
,
x
⩽
0
f_n(x)=\begin{cases} \dfrac{1}{2Γ(n/2)}(\dfrac{x}{2})^{n/2-1}e^{-x/2},\quad x>0 \\ 0,\quad x⩽0 \end{cases}
fn(x)=⎩⎨⎧2Γ(n/2)1(2x)n/2−1e−x/2,x>00,x⩽0 其中
Γ
(
x
)
=
∫
0
+
∞
t
x
−
1
e
−
t
d
t
Γ(x)=\displaystyle\int_{0}^{+∞}t^{x-1}e^{-t}\mathrm{d}t
Γ(x)=∫0+∞tx−1e−tdt
性质:
(1) 设
Y
∼
χ
2
(
n
)
Y∼ χ^2(n)
Y∼χ2(n),则
E
(
Y
)
=
n
,
D
(
Y
)
=
2
n
E(Y)=n,D(Y)=2n
E(Y)=n,D(Y)=2n
(2) 分布的可加性:设
Y
1
∼
χ
2
(
n
1
)
,
Y
2
∼
χ
2
(
n
2
)
Y_1∼ χ^2(n_1),Y_2∼ χ^2(n_2)
Y1∼χ2(n1),Y2∼χ2(n2),且
Y
1
,
Y
2
Y_1,Y_2
Y1,Y2相互独立,则
Y
1
+
Y
2
∼
χ
2
(
n
1
+
n
2
)
Y_1+Y_2∼ χ^2(n_1+n_2)
Y1+Y2∼χ2(n1+n2)
推广:设
Y
1
,
Y
2
,
⋯
,
Y
m
Y_1,Y_2,\cdots,Y_m
Y1,Y2,⋯,Ym相互独立,且
Y
i
∼
χ
2
(
n
i
)
Y_i∼ χ^2(n_i)
Yi∼χ2(ni),则
∑
i
=
1
m
Y
i
∼
χ
2
(
∑
i
=
1
m
n
i
)
\displaystyle\sum_{i=1}^{m}Y_i∼ χ^2(\sum_{i=1}^{m}n_i)
i=1∑mYi∼χ2(i=1∑mni)
t 分布(t-distribution):设
X
∼
N
(
0
,
1
)
,
Y
∼
χ
2
(
n
)
X∼ N(0,1),Y∼ χ^2(n)
X∼N(0,1),Y∼χ2(n),且 X与 Y 相互独立,则称随机变量
T
=
X
Y
/
n
T=\dfrac{X}{\sqrt{Y/n}}
T=Y/nX服从自由度为 n 的 t 分布,记为
T
∼
t
(
n
)
T∼ t(n)
T∼t(n)
t
(
n
)
t(n)
t(n)分布的概率密度为:
f
n
(
x
)
=
Γ
(
n
+
1
2
)
n
π
Γ
(
n
2
)
(
1
+
x
2
n
)
−
n
+
1
2
,
x
∈
R
f_n(x)=\dfrac{Γ(\frac{n+1}{2})}{\sqrt{nπ}Γ(\frac{n}{2})}(1+\dfrac{x^2}{n})^{-\frac{n+1}{2}},\quad x\in\R
fn(x)=nπΓ(2n)Γ(2n+1)(1+nx2)−2n+1,x∈R
特别地,
n
=
1
n=1
n=1 的 t 分布就是柯西分布
f
(
x
)
=
1
π
(
1
+
x
2
)
,
x
∈
R
f(x)=\dfrac{1}{π(1+x^2)},x\in\R
f(x)=π(1+x2)1,x∈R
lim
n
→
∞
f
n
(
x
)
=
1
2
π
e
−
x
2
/
2
,
x
∈
R
\lim\limits_{n\to∞}f_n(x)=\dfrac{1}{\sqrt{2π}}e^{-x^2/2},x\in\R
n→∞limfn(x)=2π1e−x2/2,x∈R
F 分布(F-distribution):设
X
∼
χ
2
(
n
1
)
,
Y
∼
χ
2
(
n
2
)
X∼ χ^2(n_1),Y∼ χ^2(n_2)
X∼χ2(n1),Y∼χ2(n2),且 X与 Y 相互独立,则称随机变量
F
=
X
/
n
1
Y
/
n
2
F=\dfrac{X/n_1}{Y/n_2}
F=Y/n2X/n1服从自由度为
(
n
1
,
n
2
)
(n_1,n_2)
(n1,n2) 的 F 分布,记为
F
∼
F
(
n
1
,
n
2
)
F∼ F(n_1,n_2)
F∼F(n1,n2),其中
n
1
n_1
n1称为第一自由度,
n
2
n_2
n2称为第二自由度。
F
(
n
1
,
n
2
)
F(n_1,n_2)
F(n1,n2)分布的概率密度为:
f
(
x
)
=
{
(
n
1
/
n
2
)
n
1
2
B
(
n
1
2
,
n
2
2
)
x
n
1
2
−
1
(
1
+
n
1
n
2
x
)
−
n
1
+
n
2
2
,
x
>
0
0
,
x
⩽
0
f(x)=\begin{cases} \dfrac{(n_1/n_2)^{\frac{n_1}{2}}}{B(\frac{n_1}{2},\frac{n_2}{2})}x^{\frac{n_1}{2}-1}(1+\frac{n_1}{n_2}x)^{-\frac{n_1+n_2}{2}},x>0 \\ 0,\quad x⩽ 0 \end{cases}
f(x)=⎩⎪⎨⎪⎧B(2n1,2n2)(n1/n2)2n1x2n1−1(1+n2n1x)−2n1+n2,x>00,x⩽0其中
B
(
a
,
b
)
=
∫
0
1
x
a
−
1
(
1
−
x
)
b
−
1
d
x
=
Γ
(
a
)
Γ
(
b
)
Γ
(
a
+
b
)
B(a,b)=\displaystyle\int_0^1x^{a-1}(1-x)^{b-1}\mathrm{d}x=\dfrac{Γ(a)Γ(b)}{Γ(a+b)}
B(a,b)=∫01xa−1(1−x)b−1dx=Γ(a+b)Γ(a)Γ(b)
重要性质:若
F
∼
F
(
n
1
,
n
2
)
F∼ F(n_1,n_2)
F∼F(n1,n2),则
1
F
∼
F
(
n
2
,
n
1
)
\dfrac{1}{F}∼ F(n_2,n_1)
F1∼F(n2,n1)
正态总体的抽样分布
设
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn 是来自正态总体
N
(
μ
,
σ
2
)
N(μ,σ^2)
N(μ,σ2)的样本,样本均值
X
ˉ
=
1
n
∑
i
=
1
n
X
i
\bar X=\dfrac{1}{n}\displaystyle\sum_{i=1}^{n}X_i
Xˉ=n1i=1∑nXi,样本方差
S
2
=
1
n
−
1
∑
i
=
1
n
(
X
i
−
X
ˉ
)
2
S^2=\dfrac{1}{n-1}\displaystyle\sum_{i=1}^{n}(X_i-\bar X)^2
S2=n−11i=1∑n(Xi−Xˉ)2,则
(1)
X
ˉ
∼
N
(
μ
,
σ
2
/
n
)
\bar X∼ N(μ,σ^2/n)
Xˉ∼N(μ,σ2/n)
(2)
(
n
−
1
)
S
2
σ
2
∼
χ
2
(
n
−
1
)
\dfrac{(n-1)S^2}{σ^2}∼χ^2(n-1)
σ2(n−1)S2∼χ2(n−1)
(3)
X
ˉ
\bar X
Xˉ与
S
2
S^2
S2相互独立
(4)
X
ˉ
−
μ
S
/
n
∼
t
(
n
−
1
)
\dfrac{\bar X-μ}{S/\sqrt{n}}∼ t(n-1)
S/nXˉ−μ∼t(n−1)
设
X
1
,
X
2
,
…
,
X
n
1
X_1,X_2,…,X_{n_1}
X1,X2,…,Xn1和
Y
1
,
Y
2
,
…
,
Y
n
2
Y_1,Y_2,…,Y_{n_2}
Y1,Y2,…,Yn2 分别来自正态总体
N
(
μ
1
,
σ
1
2
)
N(μ_1,σ_1^2)
N(μ1,σ12)和
N
(
μ
2
,
σ
2
2
)
N(μ_2,σ_2^2)
N(μ2,σ22)的样本,均值分别为
X
ˉ
,
Y
ˉ
\bar X,\bar Y
Xˉ,Yˉ,方差分别为
S
1
2
,
S
2
2
S_1^2,S_2^2
S12,S22,则
(1)
S
1
2
/
S
2
2
σ
1
2
/
σ
2
2
∼
F
(
n
1
−
1
,
n
2
−
1
)
\dfrac{S_1^2/S_2^2}{σ_1^2/σ_2^2}∼ F(n_1-1,n_2-1)
σ12/σ22S12/S22∼F(n1−1,n2−1)
(2)
(
X
ˉ
−
Y
ˉ
)
−
(
μ
1
−
μ
2
)
σ
1
2
n
1
+
σ
2
2
n
2
∼
N
(
0
,
1
)
\dfrac{(\bar X-\bar Y)-(μ_1-μ_2)}{\sqrt{\dfrac{σ_1^2}{n_1}+\dfrac{σ_2^2}{n_2}}}∼ N(0,1)
n1σ12+n2σ22(Xˉ−Yˉ)−(μ1−μ2)∼N(0,1)
(3) 当
σ
1
2
=
σ
2
2
=
σ
2
σ_1^2=σ_2^2=σ^2
σ12=σ22=σ2时
(
X
ˉ
−
Y
ˉ
)
−
(
μ
1
−
μ
2
)
S
w
1
n
1
+
1
n
2
∼
t
(
n
1
+
n
2
−
2
)
\dfrac{(\bar X-\bar Y)-(μ_1-μ_2)}{S_w\sqrt{\dfrac{1}{n_1}+\dfrac{1}{n_2}}}∼ t(n_1+n_2-2)
Swn11+n21(Xˉ−Yˉ)−(μ1−μ2)∼t(n1+n2−2)
其中
S
w
2
=
(
n
1
−
1
)
S
1
2
+
(
n
2
−
1
)
S
2
2
n
1
+
n
2
−
2
,
S
w
=
S
w
2
S_w^2=\dfrac{(n_1-1)S_1^2+(n_2-1)S_2^2}{n_1+n_2-2},S_w=\sqrt{S_w^2}
Sw2=n1+n2−2(n1−1)S12+(n2−1)S22,Sw=Sw2
参数估计(Parameter Estimation)
参数通常是刻画总体某些概率特征的数量。例如正态分布
N
(
μ
,
σ
2
)
N(μ,σ^2)
N(μ,σ2) 中的参数
μ
μ
μ 就是该分布的均值,参数 σ2 是该分布的方差。当该参数未知时,从总体中抽取一个样本,用某种方法对该未知参数进行估计,这就是参数估计。
参数估计的形式:先从该总体中抽样得到样本,然后构造样本函数,求出未知参数的估计值或取值范围,这就是点估计和区间估计。
假设总体 X ∼ F ( x ; θ 1 , θ 2 , … , θ m ) X ∼ F(x; θ_1, θ_2,…, θ_m) X∼F(x;θ1,θ2,…,θm),其中分布函数 F F F 的表达式已知,但参数 θ 1 , θ 2 , … , θ m θ_1, θ_2,…, θ_m θ1,θ2,…,θm 未知,若记 θ = ( θ 1 , θ 2 , … , θ m ) θ = (θ_1, θ_2,…, θ_m) θ=(θ1,θ2,…,θm),则总体分布可记为 X ∼ F ( x ; θ ) X∼ F(x; θ) X∼F(x;θ),参数 θ θ θ 的取值范围称为参数空间(parameter space),记为 Θ Θ Θ
点估计(Point estimation)
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn 是来自总体
X
∼
F
(
x
;
θ
1
,
θ
2
,
…
,
θ
m
)
X ∼ F(x; θ_1, θ_2,…, θ_m)
X∼F(x;θ1,θ2,…,θm)的一个样本,
θ
1
,
θ
2
,
…
,
θ
m
θ_1, θ_2,…, θ_m
θ1,θ2,…,θm 是未知参数,构造 m 个统计量:
随
机
变
量
{
θ
^
1
(
X
1
,
X
2
,
…
,
X
n
)
θ
^
2
(
X
1
,
X
2
,
…
,
X
n
)
⋮
θ
^
m
(
X
1
,
X
2
,
…
,
X
n
)
随机变量\begin{cases} \hat θ_1(X_1,X_2,…,X_n) \\ \hat θ_2(X_1,X_2,…,X_n) \\ \quad\vdots\\ \hat θ_m(X_1,X_2,…,X_n) \end{cases}
随机变量⎩⎪⎪⎪⎪⎨⎪⎪⎪⎪⎧θ^1(X1,X2,…,Xn)θ^2(X1,X2,…,Xn)⋮θ^m(X1,X2,…,Xn)
把样本观测值
x
1
,
x
2
,
…
,
x
n
x_1,x_2,…,x_n
x1,x2,…,xn 代入上统计量,就得到 m 个数值:
数
值
{
θ
^
1
(
x
1
,
x
2
,
…
,
x
n
)
θ
^
2
(
x
1
,
x
2
,
…
,
x
n
)
⋮
θ
^
m
(
x
1
,
x
2
,
…
,
x
n
)
数值\begin{cases} \hat θ_1(x_1,x_2,…,x_n) \\ \hat θ_2(x_1,x_2,…,x_n) \\ \quad\vdots\\ \hat θ_m(x_1,x_2,…,x_n) \end{cases}
数值⎩⎪⎪⎪⎪⎨⎪⎪⎪⎪⎧θ^1(x1,x2,…,xn)θ^2(x1,x2,…,xn)⋮θ^m(x1,x2,…,xn)
称
θ
^
k
(
X
1
,
X
2
,
…
,
X
n
)
\hat θ_k(X_1,X_2,…,X_n)
θ^k(X1,X2,…,Xn)为
θ
k
θ_k
θk的估计量(estimate),
θ
^
1
(
x
1
,
x
2
,
…
,
x
n
)
\hat θ_1(x_1,x_2,…,x_n)
θ^1(x1,x2,…,xn)为
θ
k
θ_k
θk的估计值(estimator)。
常用的估计方法:矩估计法,极大似然估计法,最小二乘估计法,贝叶斯方法
矩估计法(Moment Estimation)
矩估计思想:以样本矩估计总体矩,以样本矩的函数估计总体矩的函数。
由辛钦大数定律,有
A
k
=
1
n
∑
i
=
1
n
X
i
k
→
P
μ
k
(
θ
1
,
θ
2
,
…
,
θ
m
)
(
n
→
∞
)
A_k=\dfrac{1}{n}\displaystyle\sum_{i=1}^{n}X_i^k\xrightarrow{P}μ_k(θ_1, θ_2,…, θ_m)\quad (n\to∞)
Ak=n1i=1∑nXikPμk(θ1,θ2,…,θm)(n→∞)
因此当n较大时有
A
k
≈
μ
k
(
θ
1
,
θ
2
,
…
,
θ
m
)
A_k\approx μ_k(θ_1, θ_2,…, θ_m)
Ak≈μk(θ1,θ2,…,θm)假设总体的前m阶矩存在,令
A
k
=
μ
k
(
θ
1
,
θ
2
,
…
,
θ
m
)
,
k
=
1
,
2
,
⋯
,
m
A_k=μ_k(θ_1, θ_2,…, θ_m),k=1,2,\cdots,m
Ak=μk(θ1,θ2,…,θm),k=1,2,⋯,m得方程组,其解
θ
^
k
(
X
1
,
X
2
,
…
,
X
n
)
\hat θ_k(X_1,X_2,…,X_n)
θ^k(X1,X2,…,Xn) 称为
θ
k
θ_k
θk的矩估计量。
总体均值和方差的矩估计量表达式不因不同的总体分布而异,设总体
X
X
X的均值
μ
μ
μ及方差
σ
2
σ^2
σ2都存在
μ
^
=
X
ˉ
,
σ
^
2
=
1
n
∑
i
=
1
n
(
X
i
−
X
ˉ
)
2
≜
S
~
2
\hat μ=\bar X,\hatσ^2=\dfrac{1}{n}\displaystyle\sum_{i=1}^{n}(X_i-\bar X)^2 \triangleq \tilde S^2
μ^=Xˉ,σ^2=n1i=1∑n(Xi−Xˉ)2≜S~2
设总体
X
∼
U
(
a
,
b
)
⟹
a
^
=
X
ˉ
−
3
S
~
,
b
^
=
X
ˉ
+
3
S
~
X∼ U(a,b)\implies\hat a=\bar X-\sqrt{3}\tilde S,\hat b=\bar X+\sqrt{3}\tilde S
X∼U(a,b)⟹a^=Xˉ−3S~,b^=Xˉ+3S~
极大似然估计法(Maximum Likelihood Estimate,MLE)
Fisher的极大似然思想:
假设总体
X
X
X 是离散型随机变量,其分布律为
P
{
X
=
x
}
=
p
(
x
;
θ
)
,
θ
∈
Θ
P\{X=x\}=p(x;θ),θ\in Θ
P{X=x}=p(x;θ),θ∈Θ,θ未知,
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn为样本,其观察值为
x
1
,
x
2
,
…
,
x
n
x_1,x_2,…,x_n
x1,x2,…,xn,则事件
{
X
1
=
x
1
,
X
2
=
x
2
,
⋯
,
X
n
=
x
n
}
\{X_1=x_1,X_2=x_2,\cdots,X_n=x_n\}
{X1=x1,X2=x2,⋯,Xn=xn}发生的概率(联合分布律)为
L
(
θ
)
=
L
(
x
1
,
x
2
,
…
,
x
n
;
θ
)
=
∏
i
=
1
n
p
(
x
i
;
θ
)
,
θ
∈
Θ
L(θ)=L(x_1,x_2,…,x_n;θ)=\displaystyle\prod_{i=1}^{n} p(x_i;θ),θ\in Θ
L(θ)=L(x1,x2,…,xn;θ)=i=1∏np(xi;θ),θ∈Θ这一概率随θ的取值而改变,我们自然的认为使概率取得最大值的θ值较为合理,
L
(
θ
)
L(θ)
L(θ)称为样本的似然函数(likelihood function),方程
L
(
θ
^
)
=
max
θ
∈
Θ
L
(
θ
)
,
θ
^
=
θ
^
(
x
1
,
x
2
,
…
,
x
n
)
L(\hat θ)=\displaystyle\max_{θ\in Θ} L(θ),\hat θ=\hat θ(x_1,x_2,…,x_n)
L(θ^)=θ∈ΘmaxL(θ),θ^=θ^(x1,x2,…,xn)的解
θ
^
\hat θ
θ^称为极大似然估计值,统计量
θ
^
=
θ
^
(
X
1
,
X
2
,
…
,
X
n
)
\hat θ=\hat θ(X_1,X_2,…,X_n)
θ^=θ^(X1,X2,…,Xn)为θ 的极大似然估计量(MLE)
连续型随机变量
X
X
X的概率密度函数为
f
(
x
;
θ
)
,
θ
∈
Θ
f(x;θ),θ\in Θ
f(x;θ),θ∈Θ,可取得似然函数为
L
(
θ
)
=
L
(
x
1
,
x
2
,
…
,
x
n
;
θ
)
=
∏
i
=
1
n
f
(
x
i
;
θ
)
,
θ
∈
Θ
L(θ)=L(x_1,x_2,…,x_n;θ)=\displaystyle\prod_{i=1}^{n} f(x_i;θ),θ\in Θ
L(θ)=L(x1,x2,…,xn;θ)=i=1∏nf(xi;θ),θ∈Θ
在许多情况下,
p
(
x
;
θ
)
,
f
(
x
;
θ
)
p(x;θ),f(x;θ)
p(x;θ),f(x;θ)关于θ可微,这时
θ
^
\hat θ
θ^常可从方程
d
d
θ
L
(
θ
)
=
0
\dfrac{\mathrm{d}}{\mathrm{d}θ}L(θ)=0
dθdL(θ)=0解得,又因
L
(
θ
)
L(θ)
L(θ)和
ln
L
(
θ
)
\ln L(θ)
lnL(θ)在同一处取得极值,因此
θ
^
\hat θ
θ^也可从方程
d
d
θ
ln
L
(
θ
)
=
0
\dfrac{\mathrm{d}}{\mathrm{d}θ}\ln L(θ)=0
dθdlnL(θ)=0求得,这一方程被称作对数似然方程。
极大似然估计法同样适用于分布中含多个参数
θ
1
,
θ
2
,
…
,
θ
m
θ_1, θ_2,…, θ_m
θ1,θ2,…,θm的情况,这时似然函数L是这些未知参数的函数,分别令
∂
∂
θ
i
L
=
0
,
i
=
1
,
2
,
⋯
,
m
\dfrac{∂}{∂θ_i}L=0,i=1,2,\cdots,m
∂θi∂L=0,i=1,2,⋯,m或
∂
∂
θ
i
ln
L
=
0
,
i
=
1
,
2
,
⋯
,
m
\dfrac{∂}{∂θ_i}\ln L=0,i=1,2,\cdots,m
∂θi∂lnL=0,i=1,2,⋯,m解上述由m个方程组成的方程组即可,上式称为对数似然方程组。
设总体
X
∼
N
(
μ
,
σ
2
)
⟹
μ
^
=
X
ˉ
,
σ
^
2
=
1
n
∑
i
=
1
n
(
X
i
−
X
ˉ
)
2
≜
S
~
2
X∼ N(μ,σ^2)\implies\hat μ=\bar X,\hatσ^2=\dfrac{1}{n}\displaystyle\sum_{i=1}^{n}(X_i-\bar X)^2 \triangleq \tilde S^2
X∼N(μ,σ2)⟹μ^=Xˉ,σ^2=n1i=1∑n(Xi−Xˉ)2≜S~2
设总体
X
∼
U
(
a
,
b
)
⟹
a
^
=
min
{
X
1
,
X
2
,
⋯
,
X
n
}
=
X
(
1
)
,
b
^
=
max
{
X
1
,
X
2
,
⋯
,
X
n
}
=
X
(
n
)
X∼ U(a,b)\implies\hat a=\min\{X_1,X_2,\cdots,X_n\}=X_{(1)},\hat b=\max\{X_1,X_2,\cdots,X_n\}=X_{(n)}
X∼U(a,b)⟹a^=min{X1,X2,⋯,Xn}=X(1),b^=max{X1,X2,⋯,Xn}=X(n)
极大似然估计的不变性:设 θ ^ \hat θ θ^是 θ θ θ的极大似然估计,θ 的函数 u = u ( θ ) u=u(θ) u=u(θ)且有单值反函数 θ = θ ( u ) θ=θ(u) θ=θ(u),则 u ( θ ^ ) u(\hat θ) u(θ^)是 u u u 的极大似然估计。
估计量的评判标准(Criteria for estimators)
对同一个参数,不同方法得到的估计量可能不同。用什么标准来评价一个估计量的好坏?
(1) 无偏性(unbiased):设
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn为总体
X
∼
F
(
x
;
θ
)
(
θ
∈
Θ
)
X∼ F(x;θ) \quad (θ\inΘ)
X∼F(x;θ)(θ∈Θ) 的样本,若参数
θ
θ
θ 的估计量
θ
^
=
θ
^
(
X
1
,
X
2
,
…
,
X
n
)
\hat θ=\hat θ(X_1,X_2,…,X_n)
θ^=θ^(X1,X2,…,Xn)满足
E
(
θ
^
)
=
θ
E(\hat θ)=θ
E(θ^)=θ,则
θ
^
\hat θ
θ^是
θ
θ
θ的一个无偏估计量。
若
E
(
θ
^
)
≠
θ
E(\hat θ)\neq θ
E(θ^)=θ,则
∣
E
(
θ
^
)
−
θ
∣
|E(\hat θ)-θ|
∣E(θ^)−θ∣称为估计量
θ
θ
θ 的偏差。
若
lim
n
→
∞
E
(
θ
^
)
=
θ
\lim\limits_{n\to∞}E(\hat θ)= θ
n→∞limE(θ^)=θ,则称
θ
^
\hat θ
θ^是
θ
θ
θ的渐近无偏估计量
估计量 θ ^ ( X 1 , X 2 , … , X n ) \hat θ(X_1,X_2,…,X_n) θ^(X1,X2,…,Xn)是随机变量,而估计值会有波动性,无偏性的统计意义是指在大量重复试验下,保证了没有系统误差。
设
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn是总体
X
X
X 的样本,总体的k阶矩为
μ
k
μ_k
μk,期望和方差为
μ
,
σ
2
μ,σ^2
μ,σ2
E
(
A
k
)
=
μ
k
E(A_k)=μ_k
E(Ak)=μk,样本 k阶矩都是总体k阶矩的无偏估计
特别的,
E
(
X
ˉ
)
=
μ
,
E
(
S
2
)
=
σ
2
E(\bar X)=μ,E(S^2)=σ^2
E(Xˉ)=μ,E(S2)=σ2
而
E
(
S
~
2
)
=
n
−
1
n
E
(
S
2
)
=
n
−
1
n
σ
2
E(\tilde S^2)=\dfrac{n-1}{n}E(S^2)=\dfrac{n-1}{n}σ^2
E(S~2)=nn−1E(S2)=nn−1σ2不是总体方差的无偏估计
lim
n
→
∞
E
(
S
~
2
)
=
σ
2
\lim\limits_{n\to∞}E(\tilde S^2)=σ^2
n→∞limE(S~2)=σ2,
S
~
2
\tilde S^2
S~2是方差的渐进无偏估计。
纠偏方法:若 θ ^ \hat θ θ^是有偏估计, E ( θ ^ ) = a θ + b ( a ≠ 0 , θ ∈ Θ ) E(\hat θ)=aθ+b\quad (a\neq 0,θ\inΘ) E(θ^)=aθ+b(a=0,θ∈Θ),则 θ − b a \dfrac{θ-b}{a} aθ−b是θ 的无偏估计,上面就是用了此方法。
(2) 有效性(effective):设 θ ^ 1 , θ ^ 2 \hat θ_1,\hat θ_2 θ^1,θ^2 都是 θ 的无偏估计量,若有 ∀ θ ∈ Θ , D ( θ ^ 1 ) ⩽ D ( θ ^ 2 ) ∀ θ\inΘ,D(\hat θ_1) ⩽ D(\hat θ_2) ∀θ∈Θ,D(θ^1)⩽D(θ^2),且不等号至少对某一 θ ∈ Θ θ\inΘ θ∈Θ 成立,则 θ ^ 1 \hat θ_1 θ^1比 θ ^ 1 \hat θ_1 θ^1 有效。
均方误差准则:定义
E
(
θ
^
−
θ
)
2
E(\hat θ-θ)^2
E(θ^−θ)2是估计量
θ
^
\hat θ
θ^的均方误差(mean square error;MSE),记作
M
s
e
(
θ
^
)
Mse(\hat θ)
Mse(θ^),若
θ
^
\hat θ
θ^是θ 的无偏估计量,则
M
s
e
(
θ
^
)
=
D
(
θ
^
)
Mse(\hat θ)=D(\hat θ)
Mse(θ^)=D(θ^)。
设设
θ
^
1
,
θ
^
2
\hat θ_1,\hat θ_2
θ^1,θ^2 都是 θ的点估计,若有
∀
θ
∈
Θ
,
M
s
e
(
θ
^
1
)
⩽
M
s
e
(
θ
^
2
)
∀ θ\inΘ,Mse(\hat θ_1) ⩽ Mse(\hat θ_2)
∀θ∈Θ,Mse(θ^1)⩽Mse(θ^2),且不等号至少对某一
θ
∈
Θ
θ\inΘ
θ∈Θ 成立,则在均方误差条件下
θ
^
1
\hat θ_1
θ^1优于
θ
^
1
\hat θ_1
θ^1 。
在实际应用中,均方误差准则比无偏性准则更重要
(3) 相合性(consistence):设
θ
^
(
X
1
,
X
2
,
…
,
X
n
)
\hat θ(X_1,X_2,…,X_n)
θ^(X1,X2,…,Xn)是参数θ的估计量
若
∀
θ
∈
Θ
,
∀
ϵ
>
0
,
lim
n
→
∞
P
{
∣
θ
^
−
θ
∣
<
ϵ
}
=
1
∀ θ\inΘ,∀ϵ>0,\lim\limits_{n\to∞}P\{|\hat θ-θ|<ϵ\}=1
∀θ∈Θ,∀ϵ>0,n→∞limP{∣θ^−θ∣<ϵ}=1
即
θ
^
n
→
P
θ
(
n
→
∞
)
\hat θ_n\xrightarrow{P}θ\quad (n\to∞)
θ^nPθ(n→∞),称
θ
^
\hat θ
θ^为θ的相合估计量。
相合性的相关结论:
(1) 样本 k 阶矩
A
k
A_k
Ak是总体 k 阶矩
μ
k
μ_k
μk 的相合估计量;-----辛钦大数定律证明
(2)
S
2
,
S
~
2
S^2,\tilde S^2
S2,S~2是总体方差的相合估计量。-----切比雪夫不等式证明
(3) 矩估计量一般是相合估计量。
应用
(1) 设
X
∼
N
(
μ
,
σ
2
)
⟹
M
s
e
(
S
2
)
=
D
(
S
2
)
=
2
σ
4
n
−
1
,
M
s
e
(
S
~
2
)
=
2
n
−
1
n
2
σ
4
X∼ N(μ,σ^2)\implies Mse(S^2)=D(S^2)=\dfrac{2σ^4}{n-1},Mse(\tilde S^2)=\dfrac{2n-1}{n^2}σ^4
X∼N(μ,σ2)⟹Mse(S2)=D(S2)=n−12σ4,Mse(S~2)=n22n−1σ4
当
n
>
1
n>1
n>1时,
M
s
e
(
S
2
)
>
M
s
e
(
S
~
2
)
Mse(S^2)>Mse(\tilde S^2)
Mse(S2)>Mse(S~2)因此在均方误差条件下
S
~
2
\tilde S^2
S~2优于
S
2
S^2
S2
(2) 设
X
1
,
X
2
,
…
,
X
n
X_1,X_2,…,X_n
X1,X2,…,Xn是总体指数分布
X
∼
E
x
p
(
θ
)
X∼ Exp(θ)
X∼Exp(θ)的样本
θ 的矩估计量
X
ˉ
\bar X
Xˉ和极大似然估计量
n
Z
,
Z
=
min
{
X
1
,
X
2
,
…
,
X
n
}
nZ,Z=\min\{X_1,X_2,…,X_n\}
nZ,Z=min{X1,X2,…,Xn}均是无偏估计;
但
D
(
X
ˉ
)
=
θ
2
/
n
,
D
(
n
Z
)
=
θ
2
D(\bar X)=θ^2/n,D(nZ)=θ^2
D(Xˉ)=θ2/n,D(nZ)=θ2,当
n
>
1
n>1
n>1时,
D
(
X
ˉ
)
<
D
(
n
Z
)
D(\bar X)<D(nZ)
D(Xˉ)<D(nZ),所以
X
ˉ
\bar X
Xˉ比
n
Z
nZ
nZ有效;
区间估计(Interval estimation)
对于一个未知量,根据具体样本观测值,点估计提供一个明确的数值。还希望根据所给的样本确定一个随机区间, 使其包含参数真值的概率达到指定的要求。
置信区间(Confidence Intervals):设总体
X
∼
F
(
x
;
θ
)
,
θ
∈
Θ
X∼ F(x;θ),θ\inΘ
X∼F(x;θ),θ∈Θ,若对
给定的
α
(
0
<
α
<
1
)
α(0<α<1)
α(0<α<1),由来自
X
X
X的样本确定的两个统计量
θ
^
L
=
θ
^
L
(
X
1
,
X
2
,
…
,
X
n
)
,
θ
^
U
=
θ
^
U
(
X
1
,
X
2
,
…
,
X
n
)
(
θ
^
L
<
θ
^
U
)
\hat θ_L=\hat θ_L(X_1,X_2,…,X_n),\hat θ_U=\hat θ_U(X_1,X_2,…,X_n)(\hatθ_L<\hatθ_U)
θ^L=θ^L(X1,X2,…,Xn),θ^U=θ^U(X1,X2,…,Xn)(θ^L<θ^U),对于
∀
θ
∈
Θ
∀ θ\inΘ
∀θ∈Θ,满足
P
{
θ
^
L
<
θ
<
θ
^
U
}
⩾
1
−
α
P\{\hatθ_L<θ<\hatθ_U\}⩾ 1-α
P{θ^L<θ<θ^U}⩾1−α则称随机区间
(
θ
^
L
,
θ
^
U
)
(\hatθ_L,\hatθ_U)
(θ^L,θ^U)为θ的置信度为
1
−
α
1-α
1−α 的置信区间(confidence interval)。
θ
^
L
,
θ
^
U
\hatθ_L,\hatθ_U
θ^L,θ^U分别称为置信下限和置信上限,
1
−
α
1-α
1−α称为置信度或置信水平(confidence level)。
说明:
(1) 参数θ虽然未知,但是确定的值。置信区间
(
θ
^
L
,
θ
^
U
)
(\hatθ_L,\hatθ_U)
(θ^L,θ^U)是随机的,依赖于样本,样本不同,算出的区间也不同。对于有些样本观察值,区间覆盖θ,但对于另一些样本观察值,区间则不能覆盖。反复抽样多次(各次样本容量相同),按伯努利大数定律,在这些区间中,包含真值的比例约为
1
−
α
1-α
1−α
(2) 相同的置信水平也可以得到不同的区间估计。称置信区间
(
θ
^
L
,
θ
^
U
)
(\hatθ_L,\hatθ_U)
(θ^L,θ^U)的平均长度
E
(
θ
^
U
−
θ
^
L
)
E(\hatθ_U-\hatθ_L)
E(θ^U−θ^L)为区间的精确度(accuracy)(长度最短,精度最高),精确度的一半为误差限(error limit)。在给定的样本容量下,置信水平和精确度是相互制约的。
(3) Neyman原则:在置信水平达到
1
−
α
1-α
1−α的置信区间中,选精确度尽可能高的置信区间。
求解置信区间的一般方法
(1) 构造样本的一个函数
G
=
G
(
X
1
,
X
2
,
…
,
X
n
;
θ
)
G=G(X_1,X_2,…,X_n;θ)
G=G(X1,X2,…,Xn;θ)它含有待估参数θ,不含其它未知参数,其分布已知,且分布不依赖于待估参数(常由θ的点估计出发考虑),具有这种性质的函数称为枢轴量(pivot)。
(2) 对于给定的置信水平
1
−
α
1-α
1−α,求出两个分位点
a
,
b
a,b
a,b,使得
P
{
a
<
G
<
b
}
⩾
1
−
α
P\{a<G<b\}⩾ 1-α
P{a<G<b}⩾1−α
(3) 由
a
<
G
<
b
a<G<b
a<G<b得到与之等价的θ不等式
θ
^
L
<
θ
<
θ
^
U
\hatθ_L<θ<\hatθ_U
θ^L<θ<θ^U,那么
(
θ
^
L
,
θ
^
U
)
(\hatθ_L,\hatθ_U)
(θ^L,θ^U)就是置信度为
1
−
α
1-α
1−α 的置信区间。
分位点(percentile point):设连续型随机变量
X
X
X的概率密度函数为
f
(
x
)
f(x)
f(x),对给定的
α
(
0
<
α
<
1
)
α (0<α<1)
α(0<α<1),称满足条件
P
{
X
>
x
α
}
=
∫
x
α
∞
f
(
x
)
d
x
=
α
P\{X>x_α\}=\displaystyle\int_{x_α}^{∞}f(x)\mathrm{d}x=α
P{X>xα}=∫xα∞f(x)dx=α 的点
x
α
x_α
xα为 α 分位点。
标准正态分布
N
(
0
,
1
)
N(0,1)
N(0,1)的α 分位点记作
z
α
,
z
α
=
−
z
1
−
α
z_α,z_α=-z_{1-α}
zα,zα=−z1−α
χ
2
(
n
)
χ^2(n)
χ2(n)分布的α 分位点记作
χ
α
2
(
n
)
χ_α^2(n)
χα2(n)
t
(
n
)
t(n)
t(n)分布的α 分位点记作
t
α
(
n
)
,
t
α
(
n
)
=
−
t
1
−
α
(
n
)
t_α(n),t_α(n)=-t_{1-α}(n)
tα(n),tα(n)=−t1−α(n)
F
(
n
1
,
n
2
)
F(n_1,n_2)
F(n1,n2)分布的α 分位点记作
F
α
(
n
1
,
n
2
)
F_α(n_1,n_2)
Fα(n1,n2)
单侧置信区间(one-sided confidence interval)
如果
P
{
θ
^
L
<
θ
}
⩾
1
−
α
P\{\hatθ_L<θ\}⩾ 1-α
P{θ^L<θ}⩾1−α则称随机变量
θ
^
L
\hatθ_L
θ^L为θ的置信度为
1
−
α
1-α
1−α 的单侧置信下限(one-side confidence lower limit)。
如果
P
{
θ
<
θ
^
U
}
⩾
1
−
α
P\{θ<\hatθ_U\}⩾ 1-α
P{θ<θ^U}⩾1−α则称随机变量
θ
^
U
\hatθ_U
θ^U为θ的置信度为
1
−
α
1-α
1−α 的单侧置信上限(two-side confidence lower limit)。
正态总体参数的区间估计示例
设总体
X
∼
N
(
μ
,
σ
2
)
,
X
1
,
X
2
,
…
,
X
n
X∼ N(μ,σ^2),X_1,X_2,…,X_n
X∼N(μ,σ2),X1,X2,…,Xn为样本,
X
ˉ
,
S
2
\bar X,S^2
Xˉ,S2分别为样本均值和样本方差,置信区间为
1
−
α
1-α
1−α
(1)
σ
2
σ^2
σ2已知时,
X
ˉ
\bar X
Xˉ是 μ 的极大似然估计,取枢轴量
G
=
X
ˉ
−
μ
σ
/
n
∼
N
(
0
,
1
)
G=\dfrac{\bar X-μ}{σ/\sqrt{n}}∼ N(0,1)
G=σ/nXˉ−μ∼N(0,1)
设常数
a
,
b
a,b
a,b满足:
P
{
a
<
G
<
b
}
⩾
1
−
α
P\{a<G<b\}⩾ 1-α
P{a<G<b}⩾1−α
等价于
P
{
X
ˉ
−
σ
n
b
<
μ
<
X
ˉ
−
σ
n
a
}
⩾
1
−
α
P\{\bar X-\dfrac{σ}{\sqrt{n}}b<μ<\bar X-\dfrac{σ}{\sqrt{n}}a\}⩾ 1-α
P{Xˉ−nσb<μ<Xˉ−nσa}⩾1−α
此时区间的长度为
(
b
−
a
)
σ
/
n
(b-a)σ/\sqrt{n}
(b−a)σ/n
由正态分布的对称性知,
a
=
−
b
=
−
z
α
/
2
a=-b=-z_{α/2}
a=−b=−zα/2时,区间的长度达到最短
L
=
z
α
/
2
σ
/
n
L=z_{α/2}σ/\sqrt{n}
L=zα/2σ/n
n固定,置信水平提高,即
1
−
α
1-α
1−α增大,则
z
α
/
2
z_{α/2}
zα/2增大,所以
L
L
L 变大,精确度降低;反之亦然.
所得 μ 双侧置信区间为
(
X
ˉ
−
σ
n
z
α
/
2
,
X
ˉ
+
σ
n
z
α
/
2
)
(\bar X-\dfrac{σ}{\sqrt{n}}z_{α/2},\bar X+\dfrac{σ}{\sqrt{n}}z_{α/2})
(Xˉ−nσzα/2,Xˉ+nσzα/2)
单侧置信下限为
X
ˉ
−
σ
n
z
α
\bar X-\dfrac{σ}{\sqrt{n}}z_{α}
Xˉ−nσzα
单侧置信上限为
X
ˉ
+
σ
n
z
α
\bar X+\dfrac{σ}{\sqrt{n}}z_{α}
Xˉ+nσzα
(2)
σ
2
σ^2
σ2已知时,以
S
2
S^2
S2估计
σ
2
σ^2
σ2,的枢轴量
G
=
X
ˉ
−
μ
S
/
n
∼
t
(
n
−
1
)
G=\dfrac{\bar X-μ}{S/\sqrt{n}}∼ t(n-1)
G=S/nXˉ−μ∼t(n−1)
所得 μ 双侧置信区间为
(
X
ˉ
−
S
n
t
α
/
2
(
n
−
1
)
,
X
ˉ
+
S
n
t
α
/
2
(
n
−
1
)
)
(\bar X-\dfrac{S}{\sqrt{n}}t_{α/2}(n-1),\bar X+\dfrac{S}{\sqrt{n}}t_{α/2}(n-1))
(Xˉ−nStα/2(n−1),Xˉ+nStα/2(n−1))
单侧置信下限为
X
ˉ
−
σ
n
t
α
(
n
−
1
)
\bar X-\dfrac{σ}{\sqrt{n}}t_{α}(n-1)
Xˉ−nσtα(n−1)
单侧置信上限为
X
ˉ
+
σ
n
t
α
(
n
−
1
)
\bar X+\dfrac{σ}{\sqrt{n}}t_{α}(n-1)
Xˉ+nσtα(n−1)
(3) 其他总体均值的区间估计
当n足够大时(一般 n>30),由中心极限定理知,近似
X
ˉ
−
μ
σ
/
n
∼
N
(
0
,
1
)
\dfrac{\bar X-μ}{σ/\sqrt{n}}∼ N(0,1)
σ/nXˉ−μ∼N(0,1)
正态总体的均值、方差的置信区间和单侧置信限(置信水平为 1 − α 1-α 1−α)
一个正态总体
待估参数 | 其他参数 | 枢轴量G的分布 | 置信区间 | 单侧置信上/下限 |
---|---|---|---|---|
μ μ μ | σ 2 σ^2 σ2已知 | X ˉ − μ σ / n ∼ N ( 0 , 1 ) \dfrac{\bar X-μ}{σ/\sqrt{n}}∼ N(0,1) σ/nXˉ−μ∼N(0,1) | ( X ˉ ± σ n z α / 2 ) (\bar X±\dfrac{σ}{\sqrt{n}}z_{α/2}) (Xˉ±nσzα/2) | X ˉ ± σ n z α / 2 \bar X±\dfrac{σ}{\sqrt{n}}z_{α/2} Xˉ±nσzα/2 |
μ μ μ | σ 2 σ^2 σ2未知 | X ˉ − μ S / n ∼ t ( n − 1 ) \dfrac{\bar X-μ}{S/\sqrt{n}}∼ t(n-1) S/nXˉ−μ∼t(n−1) | ( X ˉ ± S n t α / 2 ( n − 1 ) ) \left(\bar X± \dfrac{S}{\sqrt{n}}t_{α/2}(n-1)\right) (Xˉ±nStα/2(n−1)) | X ˉ ± σ n t α ( n − 1 ) \bar X±\dfrac{σ}{\sqrt{n}}t_{α}(n-1) Xˉ±nσtα(n−1) |
σ 2 σ^2 σ2 | μ μ μ未知 | ( n − 1 ) S 2 σ 2 ∼ χ 2 ( n − 1 ) \dfrac{(n-1)S^2}{σ^2}∼χ^2(n-1) σ2(n−1)S2∼χ2(n−1) | ( ( n − 1 ) S 2 χ α / 2 2 ( n − 1 ) , ( n − 1 ) S 2 χ 1 − α / 2 2 ( n − 1 ) ) \left(\dfrac{(n-1)S^2}{χ^2_{α/2}(n-1)},\dfrac{(n-1)S^2}{χ^2_{1-α/2}(n-1)}\right) (χα/22(n−1)(n−1)S2,χ1−α/22(n−1)(n−1)S2) |
(
n
−
1
)
S
2
χ
α
/
2
2
(
n
−
1
)
\dfrac{(n-1)S^2}{χ^2_{α/2}(n-1)}
χα/22(n−1)(n−1)S2 ( n − 1 ) S 2 χ 1 − α / 2 2 ( n − 1 ) \dfrac{(n-1)S^2}{χ^2_{1-α/2}(n-1)} χ1−α/22(n−1)(n−1)S2 |
两个正态总体
待估参数 | 其他参数 | 枢轴量G的分布 | 置信区间 | 单侧置信上/下限 |
---|---|---|---|---|
μ 1 − μ 2 μ_1-μ_2 μ1−μ2 | σ 1 2 , σ 2 2 σ^2_1,σ^2_2 σ12,σ22 已知 | ( X ˉ − Y ˉ ) − ( μ 1 − μ 2 ) σ 1 2 n 1 + σ 2 2 n 2 ∼ N ( 0 , 1 ) \dfrac{(\bar X-\bar Y)-(μ_1-μ_2)}{\sqrt{\dfrac{σ_1^2}{n_1}+\dfrac{σ_2^2}{n_2}}}∼ N(0,1) n1σ12+n2σ22(Xˉ−Yˉ)−(μ1−μ2)∼N(0,1) | ( ( X ˉ − Y ˉ ) ± z α / 2 σ 1 2 n 1 + σ 2 2 n 2 ) \left((\bar X-\bar Y)± z_{α/2}\sqrt{\dfrac{σ_1^2}{n_1}+\dfrac{σ_2^2}{n_2}}\right) ((Xˉ−Yˉ)±zα/2n1σ12+n2σ22) | ( X ˉ − Y ˉ ) ± z α / 2 σ 1 2 n 1 + σ 2 2 n 2 (\bar X-\bar Y)± z_{α/2}\sqrt{\dfrac{σ_1^2}{n_1}+\dfrac{σ_2^2}{n_2}} (Xˉ−Yˉ)±zα/2n1σ12+n2σ22 |
μ 1 − μ 2 μ_1-μ_2 μ1−μ2 | σ 1 2 = σ 2 2 = σ 2 σ^2_1=σ^2_2=σ^2 σ12=σ22=σ2 未知 |
(
X
ˉ
−
Y
ˉ
)
−
(
μ
1
−
μ
2
)
S
w
1
n
1
+
1
n
2
∼
t
(
n
1
+
n
2
−
2
)
\dfrac{(\bar X-\bar Y)-(μ_1-μ_2)}{S_w\sqrt{\dfrac{1}{n_1}+\dfrac{1}{n_2}}}∼ t(n_1+n_2-2)
Swn11+n21(Xˉ−Yˉ)−(μ1−μ2)∼t(n1+n2−2) S w 2 = ( n 1 − 1 ) S 1 2 + ( n 2 − 1 ) S 2 2 n 1 + n 2 − 2 S_w^2=\dfrac{(n_1-1)S_1^2+(n_2-1)S_2^2}{n_1+n_2-2} Sw2=n1+n2−2(n1−1)S12+(n2−1)S22 | ( ( X ˉ − Y ˉ ) ± t α / 2 t ( n 1 + n 2 − 2 ) S w 1 n 1 + 1 n 2 ) \left((\bar X-\bar Y)± t_{α/2}t(n_1+n_2-2)S_w\sqrt{\dfrac{1}{n_1}+\dfrac{1}{n_2}}\right) ((Xˉ−Yˉ)±tα/2t(n1+n2−2)Swn11+n21) | ( X ˉ − Y ˉ ) ± t α / 2 t ( n 1 + n 2 − 2 ) S w 1 n 1 + 1 n 2 (\bar X-\bar Y)± t_{α/2}t(n_1+n_2-2)S_w\sqrt{\dfrac{1}{n_1}+\dfrac{1}{n_2}} (Xˉ−Yˉ)±tα/2t(n1+n2−2)Swn11+n21 |
σ 1 2 / σ 1 2 σ^2_1/σ^2_1 σ12/σ12 | μ μ μ未知 | S 1 2 / S 2 2 σ 1 2 / σ 2 2 ∼ F ( n 1 − 1 , n 2 − 1 ) \dfrac{S_1^2/S_2^2}{σ_1^2/σ_2^2}∼ F(n_1-1,n_2-1) σ12/σ22S12/S22∼F(n1−1,n2−1) | ( S 1 2 / S 2 2 F α / 2 ( n 1 − 1 , n 2 − 1 ) , S 1 2 / S 2 2 F 1 − α / 2 ( n 1 − 1 , n 2 − 1 ) ) \left(\dfrac{S_1^2/S_2^2}{F_{α/2}(n_1-1,n_2-1)},\dfrac{S_1^2/S_2^2}{F_{1-α/2}(n_1-1,n_2-1)}\right) (Fα/2(n1−1,n2−1)S12/S22,F1−α/2(n1−1,n2−1)S12/S22) | S 1 2 / S 2 2 F α ( n 1 − 1 , n 2 − 1 ) , S 1 2 / S 2 2 F 1 − α ( n 1 − 1 , n 2 − 1 ) \dfrac{S_1^2/S_2^2}{F_{α}(n_1-1,n_2-1)},\dfrac{S_1^2/S_2^2}{F_{1-α}(n_1-1,n_2-1)} Fα(n1−1,n2−1)S12/S22,F1−α(n1−1,n2−1)S12/S22 |
假设检验(Hypothesis Testing)
假设检验的基本思想(The basic idea of hypothesis testing)
假设检验是抽样推断中的一项重要内容。它是根据原资料作出一个总体指标是否等于某一个数值,某一随机变量是否服从某种概率分布的假设,然后利用样本资料采用一定的统计方法计算出有关检验的统计量,依据一定的概率原则,以较小的风险来判断估计数值与总体数值(或者估计分布与实际分布)是否存在显著差异,是否应当接受原假设选择的一种检验方法。
某工厂生产袋装葡萄糖,每袋糖的净重是个随机变量,服从正态分布,机器正常时,每袋糖净重
μ
0
μ_0
μ0 ,标准差为
σ
σ
σ 。为了检验机器是否正常,随机抽取 n 袋,称得净重:
x
1
,
x
2
,
⋯
,
x
n
x_1,x_2,\cdots,x_n
x1,x2,⋯,xn
由于长期实践表明方差比较稳定,设每袋糖的净重
X
∼
N
(
μ
,
σ
2
)
X∼ N(μ,σ^2)
X∼N(μ,σ2),接下来我们将根据样本值来判断
μ
=
μ
0
μ=μ_0
μ=μ0 或
μ
≠
μ
0
μ\neqμ_0
μ=μ0。
(1) 为此我们提出两个完全对立的假设:
H
0
:
μ
=
μ
0
,
H
1
:
μ
≠
μ
0
H_0: μ=μ_0,\quad H_1: μ\neqμ_0
H0:μ=μ0,H1:μ=μ0
观察到样本均值
X
ˉ
\bar X
Xˉ是
μ
μ
μ 无偏估计,若假设
H
0
H_0
H0为真,则偏差
∣
X
ˉ
−
μ
0
∣
|\bar X-μ_0|
∣Xˉ−μ0∣一般不应太大,考虑到
Z
=
X
ˉ
−
μ
0
σ
/
n
∼
N
(
0
,
1
)
Z=\dfrac{\bar X-μ_0}{σ/\sqrt{n}}∼ N(0,1)
Z=σ/nXˉ−μ0∼N(0,1),基于以上想法,我们选一正数
k
k
k,使得观察值满足
∣
X
ˉ
−
μ
0
σ
/
n
∣
⩾
k
\left|\dfrac{\bar X-μ_0}{σ/\sqrt{n}}\right|⩾ k
∣∣∣∣σ/nXˉ−μ0∣∣∣∣⩾k时,就拒绝原假设
H
0
H_0
H0,反之接受原假设。
(2) 然而,由于样本的随机性,当实际上
H
0
H_0
H0为真时,仍有可能拒绝原假设(这种可能性是无法消除的),犯这种错误的概率记为
P
μ
∈
H
0
{
拒
绝
H
0
}
=
P
{
t
y
p
e
I
e
r
r
o
r
}
P_{μ\in H_0}\{拒绝H_0\}=P\{\rm type\ I\ error\}
Pμ∈H0{拒绝H0}=P{type I error}。
为了确定常数
k
k
k,我们希望犯这类错误的最大概率为
α
∈
(
0
,
1
)
α\in(0,1)
α∈(0,1),令
P
{
t
y
p
e
I
e
r
r
o
r
}
=
P
{
∣
X
ˉ
−
μ
0
σ
/
n
∣
⩾
k
}
=
α
P\{\rm type\ I\ error\}=P\{\left|\dfrac{\bar X-μ_0}{σ/\sqrt{n}}\right|⩾ k\}=α
P{type I error}=P{∣∣∣∣σ/nXˉ−μ0∣∣∣∣⩾k}=α,由于
H
0
H_0
H0为真时,
Z
=
X
ˉ
−
μ
0
σ
/
n
∼
N
(
0
,
1
)
Z=\dfrac{\bar X-μ_0}{σ/\sqrt{n}}∼ N(0,1)
Z=σ/nXˉ−μ0∼N(0,1)
由正态分布分位点的定义知,
k
=
z
α
/
2
k=z_{α/2}
k=zα/2,当
∣
z
∣
⩾
z
α
/
2
|z|⩾ z_{α/2}
∣z∣⩾zα/2时,拒绝原假设。
综述
上例统计量
Z
=
X
ˉ
−
μ
0
σ
/
n
Z=\dfrac{\bar X-μ_0}{σ/\sqrt{n}}
Z=σ/nXˉ−μ0 称为检验统计量(test statistic),数
α
α
α 称为显著水平(significance level),一般
α
α
α很小
(
α
=
0.01
,
α
=
0.05
,
α
=
0.1
)
(α=0.01,α=0.05,α=0.1)
(α=0.01,α=0.05,α=0.1)。当检验统计量取某个区域 C 中的值,我们拒绝原假设,则称区域 C 为拒绝域(rejection region),他的补集称为接受域,拒绝域的边界点称为临界点(critical point)。
前面的检验问题叙述成:在显著性水平
α
α
α 下,检验
H
0
:
μ
=
μ
0
,
H
1
:
μ
≠
μ
0
H_0: μ=μ_0,\quad H_1: μ\neqμ_0
H0:μ=μ0,H1:μ=μ0
H
0
H_0
H0称为原假设(零假设 null hypothesis) ,
H
1
H_1
H1称为备择假设(对立假设 alternative hypothesis) 。
原假设与备择假设是不对称的!决定谁是原假设,依赖于立场、惯例、方便性。
(1) 保护原假设。如果错误地拒绝假设A比错误地拒绝假设B带来更严重的后果——A选作原假设。例如:2013年禽流感期间,一旦出现高烧一般先假定为禽流感患者
(2) 原假设为维持现状。为解释某些现象或效果的存在性,原假设常取为“无效果”、“无改进”、“无差异”等,拒绝原假设表示有较强的理由支持备择假设。
(3) 原假设取简单假设。只有一个参数(或分布)的假设称为简单假设。如果只有一个假设是简单假设,将其取为原假设。
参数假设的形式:设θ是反映总体指标某方面特征的量, 一般参数θ的假设有三种情形:
H
0
:
θ
=
θ
0
,
H
1
:
θ
<
θ
0
H_0: θ=θ_0,\quad H_1: θ<θ_0
H0:θ=θ0,H1:θ<θ0 (左边检验)
H
0
:
θ
=
θ
0
,
H
1
:
θ
>
θ
0
H_0: θ=θ_0,\quad H_1: θ>θ_0
H0:θ=θ0,H1:θ>θ0 (右边检验)
H
0
:
θ
=
θ
0
,
H
1
:
θ
≠
θ
0
H_0: θ=θ_0,\quad H_1: θ\neq θ_0
H0:θ=θ0,H1:θ=θ0 (双边检验)
任一检验规则在应用时,都有可能发生错误的判断——两类错误。
为此在确定检验法则时,应尽量避免犯两类错误的概率。而不考虑犯第Ⅱ类错误的检验称为显著性检验(significance test)。犯两类错误的概率记为
P
μ
∈
H
0
{
拒
绝
H
0
}
=
P
{
t
y
p
e
I
e
r
r
o
r
}
P_{μ\in H_0}\{拒绝H_0\}=P\{\rm type\ I\ error\}
Pμ∈H0{拒绝H0}=P{type I error}
P
μ
∈
H
1
{
接
受
H
0
}
=
P
{
t
y
p
e
I
I
e
r
r
o
r
}
P_{μ\in H_1}\{接受H_0\}=P\{\rm type\ II\ error\}
Pμ∈H1{接受H0}=P{type II error}
\quad | 原假设为真 | 原假设不真 |
---|---|---|
根据样本拒绝原假设 | type Ⅰ error | 正确 |
根据样本接受原假设 | 正确 | type Ⅱ error |
设总体服从
X
∼
N
(
μ
,
1
)
X∼ N(μ,1)
X∼N(μ,1),则
X
ˉ
∼
N
(
μ
,
1
/
n
)
\bar X∼ N(μ,1/n)
Xˉ∼N(μ,1/n),
H
0
:
μ
=
μ
0
,
H
1
:
μ
=
μ
1
>
μ
0
H_0: μ=μ_0,\quad H_1: μ=μ_1>μ_0
H0:μ=μ0,H1:μ=μ1>μ0,拒绝域
X
ˉ
⩾
C
\bar X⩾ C
Xˉ⩾C
正态总体下的假设检验(Hypothesis testing in a normal population)
正态总体均值和方差的检验法(显著性水平为 α α α)
假设检验 | 原假设 H 0 H_0 H0 | 检验统计量 | 备择假设 H 1 H_1 H1 | 拒绝域 |
---|---|---|---|---|
Z 检验 ( σ 2 σ^2 σ2已知) | μ = μ 0 μ ⩽ μ 0 μ ⩾ μ 0 μ=μ_0 \\ μ⩽ μ_0 \\ μ⩾μ_0 μ=μ0μ⩽μ0μ⩾μ0 | Z = X ˉ − μ 0 σ / n Z=\dfrac{\bar X-μ_0}{σ/\sqrt{n}} Z=σ/nXˉ−μ0 | μ ≠ μ 0 μ > μ 0 μ < μ 0 μ\neq μ_0 \\ μ> μ_0 \\ μ<μ_0 μ=μ0μ>μ0μ<μ0 | ∣ z ∣ ⩾ z α / 2 z ⩾ z α z ⩽ − z α \vert z\vert ⩾ z_{α/2} \\ z ⩾ z_{α} \\ z ⩽ -z_{α} ∣z∣⩾zα/2z⩾zαz⩽−zα |
t 检验 ( σ 2 σ^2 σ2未知) | μ = μ 0 μ ⩽ μ 0 μ ⩾ μ 0 μ=μ_0 \\ μ⩽ μ_0 \\ μ⩾μ_0 μ=μ0μ⩽μ0μ⩾μ0 | t = X ˉ − μ 0 S / n t=\dfrac{\bar X-μ_0}{S/\sqrt{n}} t=S/nXˉ−μ0 | μ ≠ μ 0 μ > μ 0 μ < μ 0 μ\neq μ_0 \\ μ> μ_0 \\ μ<μ_0 μ=μ0μ>μ0μ<μ0 | ∣ t ∣ ⩾ t α / 2 ( n − 1 ) t ⩾ t α ( n − 1 ) t ⩽ − t α ( n − 1 ) \vert t\vert ⩾ t_{α/2}(n-1) \\ t ⩾ t_{α}(n-1) \\ t ⩽ -t_{α}(n-1) ∣t∣⩾tα/2(n−1)t⩾tα(n−1)t⩽−tα(n−1) |
Z 检验 ( σ 1 2 , σ 2 2 σ_1^2,σ_2^2 σ12,σ22已知) | μ 1 − μ 2 = δ μ 1 − μ 2 ⩽ δ μ 1 − μ 2 ⩾ δ μ_1-μ_2=δ \\ μ_1-μ_2⩽ δ \\ μ_1-μ_2⩾δ μ1−μ2=δμ1−μ2⩽δμ1−μ2⩾δ | Z = ( X ˉ − Y ˉ ) − δ σ 1 2 n 1 + σ 2 2 n 2 Z=\dfrac{(\bar X-\bar Y)-δ}{\sqrt{\dfrac{σ_1^2}{n_1}+\dfrac{σ_2^2}{n_2}}} Z=n1σ12+n2σ22(Xˉ−Yˉ)−δ | μ 1 − μ 2 ≠ δ μ 1 − μ 2 > δ μ 1 − μ 2 < δ μ_1-μ_2\neqδ \\ μ_1-μ_2> δ \\ μ_1-μ_2<δ μ1−μ2=δμ1−μ2>δμ1−μ2<δ | ∣ z ∣ ⩾ z α / 2 z ⩾ z α z ⩽ − z α \vert z\vert ⩾ z_{α/2} \\ z ⩾ z_{α} \\ z ⩽ -z_{α} ∣z∣⩾zα/2z⩾zαz⩽−zα |
t 检验 ( σ 1 2 = σ 2 2 = σ 2 σ_1^2=σ_2^2=σ^2 σ12=σ22=σ2未知) | μ 1 − μ 2 = δ μ 1 − μ 2 ⩽ δ μ 1 − μ 2 ⩾ δ μ_1-μ_2=δ \\ μ_1-μ_2⩽ δ \\ μ_1-μ_2⩾δ μ1−μ2=δμ1−μ2⩽δμ1−μ2⩾δ |
t
=
(
X
ˉ
−
Y
ˉ
)
−
δ
S
w
1
n
1
+
1
n
2
t=\dfrac{(\bar X-\bar Y)-δ}{S_w\sqrt{\dfrac{1}{n_1}+\dfrac{1}{n_2}}}
t=Swn11+n21(Xˉ−Yˉ)−δ S w 2 = ( n 1 − 1 ) S 1 2 + ( n 2 − 1 ) S 2 2 n 1 + n 2 − 2 S_w^2=\dfrac{(n_1-1)S_1^2+(n_2-1)S_2^2}{n_1+n_2-2} Sw2=n1+n2−2(n1−1)S12+(n2−1)S22 | μ 1 − μ 2 ≠ δ μ 1 − μ 2 > δ μ 1 − μ 2 < δ μ_1-μ_2\neqδ \\ μ_1-μ_2> δ \\ μ_1-μ_2<δ μ1−μ2=δμ1−μ2>δμ1−μ2<δ | ∣ t ∣ ⩾ t α / 2 ( n 1 + n 2 − 2 ) t ⩾ t α ( n 1 + n 2 − 2 ) t ⩽ − t α ( n 1 + n 2 − 2 ) \vert t\vert ⩾ t_{α/2}(n_1+n_2-2) \\ t ⩾ t_{α}(n_1+n_2-2) \\ t ⩽ -t_{α}(n_1+n_2-2) ∣t∣⩾tα/2(n1+n2−2)t⩾tα(n1+n2−2)t⩽−tα(n1+n2−2) |
χ
2
χ^2
χ2检验 ( μ μ μ未知) | σ 2 = σ 0 2 σ 2 ⩽ σ 0 2 σ 2 ⩾ σ 0 2 σ^2=σ_0^2 \\ σ^2⩽ σ^2_0 \\ σ^2⩾σ^2_0 σ2=σ02σ2⩽σ02σ2⩾σ02 | χ 2 = ( n − 1 ) S 2 σ 2 χ^2=\dfrac{(n-1)S^2}{σ^2} χ2=σ2(n−1)S2 | σ 2 ≠ σ 0 2 σ 2 > σ 0 2 σ 2 < σ 0 2 σ^2\neq σ^2_0 \\ σ^2> σ^2_0 \\ σ^2<σ^2_0 σ2=σ02σ2>σ02σ2<σ02 | χ 2 ⩾ χ α / 2 2 ( n − 1 ) 或 χ 2 ⩽ χ 1 − α / 2 2 ( n − 1 ) χ 2 ⩾ χ α 2 ( n − 1 ) χ 2 ⩽ − χ α 2 ( n − 1 ) χ^2 ⩾ χ^2_{α/2}(n-1) 或\\ χ^2 ⩽ χ^2_{1-α/2}(n-1) \\ χ^2 ⩾ χ^2_{α}(n-1) \\ χ^2 ⩽ -χ^2_{α}(n-1) χ2⩾χα/22(n−1)或χ2⩽χ1−α/22(n−1)χ2⩾χα2(n−1)χ2⩽−χα2(n−1) |
F检验 ( μ 1 , μ 2 μ_1,μ_2 μ1,μ2未知) | σ 2 = σ 0 2 σ 2 ⩽ σ 0 2 σ 2 ⩾ σ 0 2 σ^2=σ_0^2 \\ σ^2⩽ σ^2_0 \\ σ^2⩾σ^2_0 σ2=σ02σ2⩽σ02σ2⩾σ02 | F = S 1 2 S 2 2 F=\dfrac{S_1^2}{S_2^2} F=S22S12 | σ 2 ≠ σ 0 2 σ 2 > σ 0 2 σ 2 < σ 0 2 σ^2\neq σ^2_0 \\ σ^2> σ^2_0 \\ σ^2<σ^2_0 σ2=σ02σ2>σ02σ2<σ02 | F ⩾ F α / 2 ( n 1 − 1 , n 2 − 1 ) 或 F ⩽ F 1 − α / 2 ( n 1 − 1 , n 2 − 1 ) F ⩾ F α ( n 1 − 1 , n 2 − 1 ) F ⩽ − F α ( n 1 − 1 , n 2 − 1 ) F ⩾ F_{α/2}(n_1-1,n_2-1) 或\\ F ⩽ F_{1-α/2}(n_1-1,n_2-1) \\ F ⩾ F_{α}(n_1-1,n_2-1) \\ F ⩽ -F_{α}(n_1-1,n_2-1) F⩾Fα/2(n1−1,n2−1)或F⩽F1−α/2(n1−1,n2−1)F⩾Fα(n1−1,n2−1)F⩽−Fα(n1−1,n2−1) |
t 检验 (成对数据D=X-Y) | μ D = 0 μ D ⩽ 0 μ D ⩾ 0 μ_D=0 \\ μ_D⩽ 0 \\ μ_D⩾ 0 μD=0μD⩽0μD⩾0 | t = D ˉ − 0 S D / n t=\dfrac{\bar D-0}{S_D/\sqrt{n}} t=SD/nDˉ−0 | μ D ≠ 0 μ D > 0 μ D < 0 μ_D\neq 0 \\ μ_D> 0 \\ μ_D<0 μD=0μD>0μD<0 | ∣ t ∣ ⩾ t α / 2 ( n − 1 ) t ⩾ t α ( n − 1 ) t ⩽ − t α ( n − 1 ) \vert t\vert ⩾ t_{α/2}(n-1) \\ t ⩾ t_{α}(n-1) \\ t ⩽ -t_{α}(n-1) ∣t∣⩾tα/2(n−1)t⩾tα(n−1)t⩽−tα(n−1) |
置信区间和假设检验的关系(The relationship between confidence intervals and hypothesis testing)
置信区间和假设检验之间有明显的关联,先考察置信区间和双边检验的关系
设总体
X
∼
F
(
x
;
θ
)
,
θ
∈
Θ
X∼ F(x;θ),θ\inΘ
X∼F(x;θ),θ∈Θ,
X
1
,
X
2
,
⋯
,
X
n
X_1,X_2,\cdots,X_n
X1,X2,⋯,Xn是来自总体的样本,
x
1
,
x
2
,
⋯
,
x
n
x_1,x_2,\cdots,x_n
x1,x2,⋯,xn是样本值
(1) 设区间
(
θ
^
L
,
θ
^
U
)
(\hatθ_L,\hatθ_U)
(θ^L,θ^U)为θ的置信水平为
1
−
α
1-α
1−α 的置信区间。则对于
∀
θ
∈
Θ
∀ θ\inΘ
∀θ∈Θ,满足
P
{
θ
^
L
<
θ
<
θ
^
U
}
⩾
1
−
α
(1)
P\{\hatθ_L<θ<\hatθ_U\}⩾ 1-α \tag{1}
P{θ^L<θ<θ^U}⩾1−α(1)其中
θ
^
L
=
θ
^
L
(
X
1
,
X
2
,
…
,
X
n
)
,
θ
^
U
=
θ
^
U
(
X
1
,
X
2
,
…
,
X
n
)
(
θ
^
L
<
θ
^
U
)
\hat θ_L=\hat θ_L(X_1,X_2,…,X_n),\hat θ_U=\hat θ_U(X_1,X_2,…,X_n)(\hatθ_L<\hatθ_U)
θ^L=θ^L(X1,X2,…,Xn),θ^U=θ^U(X1,X2,…,Xn)(θ^L<θ^U)
(2) 考虑显著性水平为
α
α
α 的双边检验
H
0
:
θ
=
θ
0
,
H
1
:
θ
≠
θ
0
H_0: θ=θ_0,\quad H_1: θ\neq θ_0
H0:θ=θ0,H1:θ=θ0,由(1)式得
P
{
θ
0
⩽
θ
^
L
∪
θ
0
⩾
θ
^
U
}
⩽
α
P\{θ_0⩽ \hatθ_L∪ θ_0⩾ \hatθ_U\}⩽ α
P{θ0⩽θ^L∪θ0⩾θ^U}⩽α,则假设检验的拒绝域为
θ
0
⩽
θ
^
L
∪
θ
0
⩾
θ
^
U
θ_0⩽ \hatθ_L∪ θ_0⩾ \hatθ_U
θ0⩽θ^L∪θ0⩾θ^U,接受域为
θ
^
L
<
θ
0
<
θ
^
U
\hatθ_L<θ_0<\hatθ_U
θ^L<θ0<θ^U
这就是说当我们做上述假设检验时,先求出θ的置信区间,若
θ
0
∈
(
θ
^
L
,
θ
^
U
)
θ_0\in(\hatθ_L,\hatθ_U)
θ0∈(θ^L,θ^U)则接受原假设
H
0
H_0
H0,若
θ
0
∉
(
θ
^
L
,
θ
^
U
)
θ_0\not\in(\hatθ_L,\hatθ_U)
θ0∈(θ^L,θ^U)则拒绝原假设
H
0
H_0
H0
反之,显著水平为
α
α
α 的双边检验也是置信水平为
1
−
α
1-α
1−α 的置信区间。
结论
(1) 置信水平为
1
−
α
1-α
1−α 的置信区间
(
θ
^
L
,
θ
^
U
)
(\hatθ_L,\hatθ_U)
(θ^L,θ^U) 与显著水平为
α
α
α 的双边检验
H
0
:
θ
=
θ
0
,
H
1
:
θ
≠
θ
0
H_0: θ=θ_0,\quad H_1: θ\neq θ_0
H0:θ=θ0,H1:θ=θ0 有对应关系
(2) 置信水平为
1
−
α
1-α
1−α 的单侧置信区间
(
−
∞
,
θ
^
U
)
(-∞,\hatθ_U)
(−∞,θ^U) 与显著水平为
α
α
α 的左边检验
H
0
:
θ
⩾
θ
0
,
H
1
:
θ
<
θ
0
H_0: θ⩾ θ_0,\quad H_1: θ< θ_0
H0:θ⩾θ0,H1:θ<θ0 有对应关系
(3) 置信水平为
1
−
α
1-α
1−α 的单侧置信区间
(
θ
^
L
,
+
∞
)
(\hatθ_L,+∞)
(θ^L,+∞) 与显著水平为
α
α
α 的右边检验
H
0
:
θ
⩽
θ
0
,
H
1
:
θ
>
θ
0
H_0: θ⩽ θ_0,\quad H_1: θ> θ_0
H0:θ⩽θ0,H1:θ>θ0 有对应关系
非参数检验(Nonparametric tests)
前面介绍的各种检验都是在总体服从正态分布前提下,对参数进行假设检验的。但是,在数据分析过程中,由于种种原因,总体服从何种理论分布并不知道,此时参数检验的方法就不再适用了。
百度百科:非参数检验
Pearson
χ
2
χ^2
χ2 拟合优度检验(Goodness of Fit Test)
拟合优度检验是用卡方统计量根据样本数据,推断总体分布与期望分布或某一理论分布是否存在显著差异,是一种吻合性检验,通常适于对有多项分类值的总体分布的分析。它的原假设是:样本来自的总体分布与期望分布或某一理论分布无差异。
设总体 X ∼ F ( x ) , F ( x ) X∼ F(x), F(x) X∼F(x),F(x)未知,检验假设 H 0 : F ( x ) = F 0 ( x ) ∀ x ∈ R H_0:F(x)=F_0(x)\quad ∀ x\in\R H0:F(x)=F0(x)∀x∈R,其中 F 0 ( x ) F_0(x) F0(x)为已知分布函数。
拟合优度检验的基本原理和步骤:
(1) 在
H
0
H_0
H0 下,总体
X
X
X 取值的全体分成 k 个两两不相交的子集
A
1
,
⋯
,
A
k
A_1,\cdots,A_k
A1,⋯,Ak
(2) 以
n
i
(
i
=
1
,
⋯
,
k
)
n_i(i=1,\cdots,k)
ni(i=1,⋯,k) 记样本观察值
x
1
,
⋯
,
x
n
x_1,\cdots,x_n
x1,⋯,xn 中落在
A
i
A_i
Ai的个数(实际频数)
(3) 当
H
0
H_0
H0 为真且
F
0
(
x
)
F_0(x)
F0(x) 完全已知时,计算事件
A
i
A_i
Ai发生的概率
p
i
=
P
{
X
=
a
i
}
p_i=P\{X=a_i\}
pi=P{X=ai}
当
F
0
(
x
)
F_0(x)
F0(x)含有
r
r
r 个未知参数时,先利用极大似然法估计
r
r
r 个未知参数,然后求得的
p
i
p_i
pi 估计
p
^
i
\hat p_i
p^i
此时称
n
p
i
np_i
npi(或
n
p
^
i
n\hat p_i
np^i )为理论频数。
(4) 以
F
0
(
x
)
F_0(x)
F0(x) 完全已知为例,由大数定理得
n
i
n
→
p
i
(
n
→
∞
)
\dfrac{n_i}{n}\to p_i\quad(n\to∞)
nni→pi(n→∞)
⟹
∣
n
i
−
n
p
i
∣
\implies |n_i-np_i|
⟹∣ni−npi∣ 应偏小
⟹
∑
i
=
1
k
C
i
(
n
i
/
n
−
p
i
)
2
\implies \displaystyle\sum^{k}_{i=1}C_i(n_i/n-p_i)^2
⟹i=1∑kCi(ni/n−pi)2 应偏小
其中 C i C_i Ci为给定常数,Pearson证明,如果 C i = n / p i C_i=n/p_i Ci=n/pi,则统计量 χ 2 = ∑ i = 1 k n p i ( n i n − p i ) 2 = ∑ i = 1 k ( n i 2 n p i − n ) χ^2=\displaystyle\sum^{k}_{i=1}\frac{n}{p_i}(\frac{n_i}{n}-p_i)^2=\displaystyle\sum^{k}_{i=1}(\frac{n_i^2}{np_i}-n) χ2=i=1∑kpin(nni−pi)2=i=1∑k(npini2−n) 具有以下定理中的性质
定理:不论总体服从什么分布,若 n充分大,当
H
0
H_0
H0 为真时
统计量
χ
2
=
∑
i
=
1
k
(
n
i
2
n
p
i
−
n
)
χ^2=\displaystyle\sum^{k}_{i=1}(\frac{n_i^2}{np_i}-n)
χ2=i=1∑k(npini2−n)近似服从
χ
2
(
k
−
1
)
χ^2(k-1)
χ2(k−1),拒绝域为
χ
2
⩾
χ
α
2
(
k
−
1
)
χ^2⩾ χ^2_{α}(k-1)
χ2⩾χα2(k−1)
统计量
χ
2
=
∑
i
=
1
k
(
n
i
2
n
p
^
i
−
n
)
χ^2=\displaystyle\sum^{k}_{i=1}(\frac{n_i^2}{n\hat p_i}-n)
χ2=i=1∑k(np^ini2−n)近似服从
χ
2
(
k
−
r
−
1
)
χ^2(k-r-1)
χ2(k−r−1),拒绝域为
χ
2
⩾
χ
α
2
(
k
−
r
−
1
)
χ^2⩾ χ^2_{α}(k-r-1)
χ2⩾χα2(k−r−1)
其中 k 为分类数,
r
r
r 为
F
0
(
x
)
F_0(x)
F0(x) 中被估未知参数的个数,
α
α
α为显著性水平。
注:根据实践,要求
n
⩾
50
,
n
p
i
(
或
n
p
^
i
)
⩾
5
n⩾ 50,np_i(或n\hat p_i)⩾ 5
n⩾50,npi(或np^i)⩾5。否则应适当合并相邻的类,以满足要求。
Pearson
χ
2
χ^2
χ2 独立性检验(test of independence)
考虑二维总体
(
X
,
Y
)
∼
F
(
x
,
y
)
(X,Y)∼ F(x,y)
(X,Y)∼F(x,y),且
X
∼
F
1
(
x
)
,
Y
∼
F
2
(
y
)
X∼ F_1(x),Y∼ F_2(y)
X∼F1(x),Y∼F2(y),
(
X
1
,
Y
1
)
,
(
X
2
,
Y
2
)
,
⋯
,
(
X
n
,
Y
n
)
(X_1,Y_1),(X_2,Y_2),\cdots,(X_n,Y_n)
(X1,Y1),(X2,Y2),⋯,(Xn,Yn)是来自总体的样本
检验假设
H
0
:
X
,
Y
独
立
H
1
:
X
,
Y
不
独
立
H_0:X,Y独立\quad H_1:X,Y不独立
H0:X,Y独立H1:X,Y不独立
即
H
0
:
F
(
x
,
y
)
=
F
1
(
x
)
F
2
(
y
)
H
1
:
F
(
x
,
y
)
≠
F
1
(
x
)
F
2
(
y
)
H_0:F(x,y)=F_1(x)F_2(y)\quad H_1:F(x,y)\neq F_1(x)F_2(y)
H0:F(x,y)=F1(x)F2(y)H1:F(x,y)=F1(x)F2(y)
(1) 将两个指标 X , Y X , Y X,Y 分别划分成两两互不相交的子集 A 1 , A 2 , ⋯ , A r A_1,A_2,\cdots,A_r A1,A2,⋯,Ar 和 B 1 , B 2 , ⋯ , B s B_1,B_2,\cdots,B_s B1,B2,⋯,Bs,样本值落入小区域的个数为
X / Y B 1 B 2 ⋯ B j ⋯ B s { X ∈ A i } A 1 n 11 n 12 ⋯ n 1 j ⋯ n 1 s n 1 ∙ A 2 n 21 n 22 ⋯ n 2 j ⋯ n 2 s n 2 ∙ ⋮ ⋯ ⋯ ⋯ ⋮ A i n i 1 n i 2 ⋯ n i j ⋯ n i s n i ∙ ⋮ ⋯ ⋯ ⋯ ⋮ A r n r 1 n r 2 ⋯ n r j ⋯ n r s n r ∙ { Y ∈ B j } n ∙ 1 n ∙ 2 ⋯ p ∙ j ⋯ n ∙ s n \begin{array}{c|cccccc|c} X/Y &B_1&B_2&\cdots&B_j&\cdots &B_s&\{X\in A_i\} \\ \hline A_1 &n_{11}&n_{12}&\cdots&n_{1j}&\cdots &n_{1s}&n_{1 \bullet} \\ A_2 &n_{21}&n_{22}&\cdots&n_{2j}&\cdots &n_{2s}&n_{2 \bullet} \\ \vdots &\cdots&&\cdots&&\cdots& &\vdots \\ A_i&n_{i1}&n_{i2}&\cdots&n_{ij}&\cdots &n_{is}&n_{i \bullet} \\ \vdots &\cdots&&\cdots&&\cdots &&\vdots \\ A_r&n_{r1}&n_{r2}&\cdots&n_{rj}&\cdots &n_{rs}&n_{r \bullet} \\ \hline \{Y\in B_j\} &n_{\bullet 1}&n_{\bullet 2}&\cdots&p_{\bullet j}&\cdots &n_{\bullet s}&n \end{array} X/YA1A2⋮Ai⋮Ar{Y∈Bj}B1n11n21⋯ni1⋯nr1n∙1B2n12n22ni2nr2n∙2⋯⋯⋯⋯⋯⋯⋯⋯Bjn1jn2jnijnrjp∙j⋯⋯⋯⋯⋯⋯⋯⋯Bsn1sn2snisnrsn∙s{X∈Ai}n1∙n2∙⋮ni∙⋮nr∙n
(2) 记
p
i
j
p_{ij}
pij 为
(
X
,
Y
)
( X , Y )
(X,Y) 落入小区域
(
i
,
j
)
( i , j )
(i,j) 中概率
边缘分布
p
i
∙
=
∑
j
=
1
s
p
i
j
,
p
∙
j
=
∑
i
=
1
r
p
i
j
p_{i\bullet}=\displaystyle\sum^{s}_{j=1}p_{ij},\quad p_{\bullet j}=\displaystyle\sum^{r}_{i=1}p_{ij}
pi∙=j=1∑spij,p∙j=i=1∑rpij
则检验转化为
H
0
:
p
i
j
=
p
i
∙
⋅
p
∙
j
H_0:p_{ij}=p_{i\bullet}\cdot p_{\bullet j}
H0:pij=pi∙⋅p∙j
由于
p
i
∙
,
p
∙
j
p_{i\bullet}, p_{\bullet j}
pi∙,p∙j未知,用估计量代替
p
^
∙
j
=
n
i
∙
n
,
p
^
∙
j
=
n
∙
j
n
\hat p_{\bullet j}=\dfrac{n_{i\bullet}}{n},\quad \hat p_{\bullet j}=\dfrac{n_{\bullet j}}{n}
p^∙j=nni∙,p^∙j=nn∙j
(3) 当
H
0
H_0
H0为真时,
n
i
j
n
≈
p
^
i
∙
⋅
p
^
∙
j
\dfrac{n_{ij}}{n}\approx \hat p_{i\bullet}\cdot \hat p_{\bullet j}
nnij≈p^i∙⋅p^∙j
统计量
χ
2
=
∑
i
=
1
r
∑
j
=
1
s
(
n
i
2
n
p
^
i
∙
⋅
p
^
∙
j
−
n
)
∼
χ
2
(
(
r
−
1
)
(
s
−
1
)
)
χ^2=\displaystyle\sum^{r}_{i=1}\sum^{s}_{j=1}(\frac{n_i^2}{n\hat p_{i\bullet}\cdot \hat p_{\bullet j}}-n)∼ χ^2((r-1)(s-1))
χ2=i=1∑rj=1∑s(np^i∙⋅p^∙jni2−n)∼χ2((r−1)(s−1))
拒绝域为
χ
2
⩾
χ
α
2
(
(
r
−
1
)
(
s
−
1
)
)
χ^2⩾ χ^2_{α}((r-1)(s-1))
χ2⩾χα2((r−1)(s−1))
偏度、峰度检验
Pearson
χ
2
χ^2
χ2 检验方法适合于一般的总体分布,但若总体分布为正态分布,其犯第
I
I
\mathrm{II}
II 类错误的概率往往比较高。为此统计学家们对检验正态总体的方法进行了比较,通过大量模拟运算的结果,认为检验正态分布以偏度和峰度检验法和夏皮罗-威尔克检验(Shapiro–Wilk test)最为有效。
随机变量
X
X
X 的偏度和峰度为标准变化量
Y
=
[
X
−
E
(
X
)
]
/
D
(
X
)
Y=[X-E(X)]/\sqrt{D(X)}
Y=[X−E(X)]/D(X) 的三阶矩和四阶矩。
ν
1
=
E
(
Y
3
)
=
B
3
B
2
3
/
2
,
ν
2
=
E
(
Y
4
)
=
B
4
B
2
2
ν_1=E(Y^3)=\dfrac{B_3}{B_2^{3/2}},\quad ν_2=E(Y^4)=\dfrac{B_4}{B_2^2}
ν1=E(Y3)=B23/2B3,ν2=E(Y4)=B22B4
其中
B
k
=
E
[
(
X
−
E
(
X
)
)
k
]
B_k=E[(X-E(X))^k]
Bk=E[(X−E(X))k]为
X
X
X的 k 阶中心距。当
X
∼
N
(
μ
,
σ
2
)
X∼ N(μ,σ^2)
X∼N(μ,σ2) 时,
ν
1
=
0
,
ν
2
=
3
ν_1=0,ν_2=3
ν1=0,ν2=3
(1) 设
X
1
,
X
2
,
⋯
,
X
n
X_1,X_2,\cdots,X_n
X1,X2,⋯,Xn是来自总体的样本,偏度和峰度
ν
1
,
ν
2
ν_1,ν_2
ν1,ν2 的矩估计量分别是
G
1
,
G
2
G_1,G_2
G1,G2
当
X
∼
N
(
μ
,
σ
2
)
X∼ N(μ,σ^2)
X∼N(μ,σ2) ,可证明
n
n
n 充分大时, 近似有
G
1
∼
N
(
0
,
σ
1
2
)
,
G
2
∼
N
(
μ
2
,
σ
2
2
)
G_1∼ N(0,σ_1^2) , G_2∼ N(μ_2,σ_2^2)
G1∼N(0,σ12),G2∼N(μ2,σ22)
其中
σ
1
2
=
6
(
n
−
2
)
(
n
+
1
)
(
n
+
3
)
,
μ
2
=
3
−
6
n
+
1
,
σ
2
2
=
24
(
n
−
2
)
(
n
−
3
)
(
n
+
1
)
2
(
n
+
3
)
(
n
+
5
)
σ_1^2=\dfrac{6(n-2)}{(n+1)(n+3)},μ_2=3-\dfrac{6}{n+1},σ_2^2=\dfrac{24(n-2)(n-3)}{(n+1)^2(n+3)(n+5)}
σ12=(n+1)(n+3)6(n−2),μ2=3−n+16,σ22=(n+1)2(n+3)(n+5)24(n−2)(n−3)
(2) 现在来检验假设
H
0
:
X
∼
N
(
μ
,
σ
2
)
H_0:X∼ N(μ,σ^2)
H0:X∼N(μ,σ2)
将
G
1
,
G
2
G_1,G_2
G1,G2标准化
U
1
=
G
1
/
σ
1
,
U
2
=
(
G
2
−
μ
2
)
/
σ
2
U_1=G_1/σ_1,U_2=(G_2-μ_2)/σ_2
U1=G1/σ1,U2=(G2−μ2)/σ2
由大数定律得
G
1
→
P
ν
1
(
n
→
∞
)
G_1\xrightarrow{P}ν_1\quad (n\to∞)
G1Pν1(n→∞)
G
2
→
P
ν
2
(
n
→
∞
)
G_2\xrightarrow{P}ν_2\quad (n\to∞)
G2Pν2(n→∞)
当
H
0
H_0
H0 为真,
n
n
n 充分大时(一般大于100为宜),有
U
1
∼
N
(
0
,
1
)
,
U
2
∼
N
(
0
,
1
)
U_1∼ N(0,1),U_2∼ N(0,1)
U1∼N(0,1),U2∼N(0,1)
∣
U
1
∣
,
∣
U
2
∣
|U_1|,|U_2|
∣U1∣,∣U2∣的观察值
∣
u
1
∣
,
∣
u
2
∣
|u_1|,|u_2|
∣u1∣,∣u2∣不应与
ν
1
=
0
,
ν
2
=
3
ν_1=0,ν_2=3
ν1=0,ν2=3 偏离过大,取显著性水平为
α
α
α,可求得拒绝域为
∣
u
1
∣
⩾
z
α
/
4
或
∣
u
2
∣
⩾
z
α
/
4
|u_1|⩾ z_{α/4}或|u_2|⩾ z_{α/4}
∣u1∣⩾zα/4或∣u2∣⩾zα/4
K-S 检验
单样本K-S检验方法能够利用样本数据推断样本来自的总体是否服从某一理论分布,是一种拟合优度的检验方法,适用于探索连续型随机变量的分布。原假设是:样本来自的总体与指定的理论分布无显著差异。
K-S检验不仅能够检验单个总体是否服从某一理论分布,还能够检验两总体分布是否存在显著差异。其原假设是:两组独立样本来自的两总体的分布无显著差异。
这里是以变量值的秩作为分析对象,而非变量值本身。
秩和检验(Rank Sum Test)
如果两个样本来自两个独立的但非正态或形态不清的两总体,要检验两样本之间的差异是否显著,不应运用参数检验中的T检验,而需采用秩和检验。
秩和检验方法最早是由维尔克松(Wilcoxon)提出,叫维尔克松两样本检验法。后来曼—惠特尼将其应用到两样本容量不等(
n
1
≠
n
2
n_1\ne n_2
n1=n2)的情况,因而又称为Mann-Whitney U test。这种方法主要用于比较两个独立样本的差异。
(1) 设两个总体
X
1
∼
F
1
(
x
)
,
X
2
∼
F
2
(
x
)
X_1∼ F_1(x),X_2∼ F_2(x)
X1∼F1(x),X2∼F2(x) 均为未知分布。
检验假设
H
0
:
F
1
(
x
)
=
F
2
(
x
)
,
H
1
:
F
1
(
x
)
≠
F
2
(
x
)
H_0:F_1(x)=F_2(x),\quad H_1:F_1(x)\neq F_2(x)
H0:F1(x)=F2(x),H1:F1(x)=F2(x)
(2) 秩的定义:设从总体
X
X
X 抽出容量为
n
n
n的样本,把样本观察值排列成
x
(
1
)
<
x
(
2
)
<
⋯
<
x
(
n
)
x_{(1)}<x_{(2)}<\cdots<x_{(n)}
x(1)<x(2)<⋯<x(n),称
x
(
i
)
x_{(i)}
x(i)的足标
i
(
i
=
1
,
2
,
⋯
,
n
)
i(i=1,2,\cdots,n)
i(i=1,2,⋯,n) 为秩。
现从两总体
X
1
,
X
2
X_1,X_2
X1,X2 中分别抽出容量为
n
1
,
n
2
n_1,n_2
n1,n2的样本(这里总假定
n
1
⩽
n
2
n_1⩽ n_2
n1⩽n2),把
n
1
+
n
2
n_1+n_2
n1+n2 个观察值混合排列求出每个观察值的秩。
计算
X
1
,
X
2
X_1,X_2
X1,X2 的样本观察值的秩相加,其和为
R
1
,
R
2
R_1,R_2
R1,R2,称为秩和。显然
R
1
,
R
2
R_1,R_2
R1,R2为离散型随机变量,且
R
1
+
R
2
=
1
2
(
n
1
+
n
2
)
(
n
1
+
n
2
+
1
)
R_1+R_2=\dfrac{1}{2}(n_1+n_2)(n_1+n_2+1)
R1+R2=21(n1+n2)(n1+n2+1)已知。
(3) 当
H
0
H_0
H0 为真时,两个样本其实来自同一整体,
X
1
X_1
X1 的样本观察值的秩应随机分散的在自然数
[
1
,
n
1
+
n
2
]
[1,n_1+n_2]
[1,n1+n2] 之间,不应过渡集中,即
R
1
∈
[
1
2
n
1
(
n
1
+
1
)
,
1
2
n
1
(
n
1
+
2
n
2
+
1
)
]
R_1\in [\dfrac{1}{2}n_1(n_1+1),\dfrac{1}{2}n_1(n_1+2n_2+1)]
R1∈[21n1(n1+1),21n1(n1+2n2+1)]且不应过度靠近区间两端。
据以上分析,对于双边检验,在显著性水平
α
α
α 下,
H
0
H_0
H0 的拒绝域为
r
1
⩽
C
U
(
α
/
2
)
或
r
1
⩾
C
L
(
α
/
2
)
r_1⩽ C_U(α/2)或r_1⩾ C_L(α/2)
r1⩽CU(α/2)或r1⩾CL(α/2),其中
r
1
r_1
r1为
R
1
R_1
R1的观察值,临界点
C
U
(
α
/
2
)
C_U(α/2)
CU(α/2) 是满足
P
{
R
1
⩽
C
U
(
α
/
2
)
}
⩽
α
/
2
P\{R_1⩽ C_U(α/2)\}⩽ α/2
P{R1⩽CU(α/2)}⩽α/2 的最大整数,
C
L
(
α
/
2
)
C_L(α/2)
CL(α/2) 是满足
P
{
R
1
⩾
C
L
(
α
/
2
)
}
⩽
α
/
2
P\{R_1⩾ C_L(α/2)\}⩽ α/2
P{R1⩾CL(α/2)}⩽α/2 的最小整数,而犯第
I
\mathrm{I}
I 类错误的概率为
α
α
α
如果知道
R
1
R_1
R1的分布,临界点是不难求得的。
(4) 当
H
0
H_0
H0 为真时,可知
X
1
X_1
X1的样本秩的不同取法共有
∁
n
1
+
n
2
n
1
∁_{n_1+n_2}^{n_1}
∁n1+n2n1种,且是等可能的。
可证明
μ
R
1
=
E
(
R
1
)
=
1
2
n
1
(
n
1
+
n
2
+
1
)
,
σ
R
1
2
=
D
(
R
1
)
=
1
12
n
1
n
2
(
n
1
+
n
2
+
1
)
μ_{R_1}=E(R_1)=\dfrac{1}{2}n_1(n_1+n_2+1),σ^2_{R_1}=D(R_1)=\dfrac{1}{12}n_1n_2(n_1+n_2+1)
μR1=E(R1)=21n1(n1+n2+1),σR12=D(R1)=121n1n2(n1+n2+1)
而当
n
1
,
n
2
⩾
10
n_1,n_2⩾ 10
n1,n2⩾10,
H
0
H_0
H0 为真时,近似有
R
1
∼
N
(
μ
R
1
,
σ
R
1
2
)
R_1∼ N(μ_{R_1},σ^2_{R_1})
R1∼N(μR1,σR12),此时我们可以采用
Z
=
(
R
1
−
μ
R
1
)
/
σ
R
1
2
Z=(R_1-μ_{R_1})/σ^2_{R_1}
Z=(R1−μR1)/σR12 作为统计量
在显著性水平
α
α
α 下的双边检验、右边检验、左边检验分别为
∣
z
∣
⩾
z
α
/
2
,
z
⩾
z
α
,
z
⩽
−
z
α
|z|⩾ z_{α/2},\quad z⩾ z_{α},\quad z⩽ -z_{α}
∣z∣⩾zα/2,z⩾zα,z⩽−zα
这里
z
z
z是
Z
Z
Z 的观察值。
假设检验问题中的 p-value 检验法(The p-value test method in hypothesis testing problems)
以上讨论的问题都是临界值法,在假设检验中常见到P值(P-Value,Probability,Pr),P值是进行检验决策的另一个依据。
为理解 p-value 从一个例子引入
设总体
X
∼
N
(
μ
,
σ
2
)
,
μ
X∼ N(μ,σ^2),μ
X∼N(μ,σ2),μ未知,现有样本
X
1
,
X
2
,
⋯
,
X
n
X_1,X_2,\cdots,X_n
X1,X2,⋯,Xn,观察值为
x
1
,
x
2
,
⋯
,
x
n
x_1,x_2,\cdots,x_n
x1,x2,⋯,xn
检验假设
H
0
:
μ
=
μ
0
,
H
1
:
μ
≠
μ
0
H_0:μ=μ_0,\quad H_1:μ\neq μ_0
H0:μ=μ0,H1:μ=μ0
采用 Z检验法,统计量为
Z
=
X
ˉ
−
μ
0
σ
/
n
Z=\dfrac{\bar X-μ_0}{σ/\sqrt{n}}
Z=σ/nXˉ−μ0,将样本数据代入得Z的观察值为
z
0
=
x
ˉ
−
μ
0
σ
/
n
z_0=\dfrac{\bar x-μ_0}{σ/\sqrt{n}}
z0=σ/nxˉ−μ0
概率
P
{
Z
⩾
∣
z
0
∣
}
=
1
−
Φ
(
∣
z
0
∣
)
P\{Z⩾ |z_0|\}=1-Φ( |z_0|)
P{Z⩾∣z0∣}=1−Φ(∣z0∣),此即为标准正态曲线下位于
∣
z
0
∣
|z_0|
∣z0∣ 右端的尾部面积。
此概率称为 Z 检验法右边检验的 p-value,若显著性水平
α
⩾
p
α⩾ p
α⩾p,则对应的临界值
z
α
⩽
z
0
z_α⩽ z_0
zα⩽z0,因而拒绝原假设。
p-value 的定义:假设检验的 p-value 由检验统计量的样本观察值得出的原假设可被拒绝的最小显著性水平。
(1) 若
p
⩽
α
p⩽ α
p⩽α,则在显著性水平 α 下拒绝
H
0
H_0
H0
(2) 若
p
>
α
p> α
p>α,则在显著性水平 α 下接受
H
0
H_0
H0
如果
p
⩽
0.01
p⩽ 0.01
p⩽0.01,说明是较强的判定结果。
如果
0.01
<
p
⩽
0.05
0.01<p⩽ 0.05
0.01<p⩽0.05,说明较弱的判定结果。
如果
0.05
<
p
⩽
0.1
0.05<p⩽ 0.1
0.05<p⩽0.1,说明拒绝的理由较弱的,检验不显著。
如果
p
>
0.1
p>0.1
p>0.1,没有理由拒绝。
方差分析(Analysis of Variance, ANOVA)
方差分析(ANOVA),是由英国统计学家费歇尔(Fisher)在20世纪20年代提出的,可用于推断两个或两个以上总体均值是否有差异的显著性检验。
在方差分析中, 通常把研究对象的特征值,即所考察的试验结果称为试验指标(Test Indicators)。
对试验指标产生影响的原因称为因素(factor)。因素中各个不同状态称为水平(level)。
如果在一项试验中,只有一个因素改变,称为单因素试验(Single-factor experiment),如果有多个因素改变称为多因素试验(Multiple-factor experiment)。
单因素方差分析(One-factor analysis of variance)
设因素 A A A 有 r r r 个水平 A 1 , A 2 , ⋯ , A r A_1,A_2,\cdots,A_r A1,A2,⋯,Ar,在水平 A i A_i Ai下进行 n i ( n i ⩾ 2 ) n_i(n_i⩾ 2) ni(ni⩾2)次试验,得到如下结果,我们假定各个水平下的样本符合正态分布,且方差相等(均值和方差均未知),且不同水平下样本之间相互独立。
水平 | 观察结果 | 样本均值 |
---|---|---|
A 1 : N ( μ 1 , σ 2 ) A_1: N(μ_1,σ^2) A1:N(μ1,σ2) | X 11 , X 12 , ⋯ , X 1 n 1 X_{11},X_{12},\cdots,X_{1n_1} X11,X12,⋯,X1n1 | X ˉ 1 ∙ \bar X_{1\bullet} Xˉ1∙ |
A 2 : N ( μ 2 , σ 2 ) A_2: N(μ_2,σ^2) A2:N(μ2,σ2) | X 21 , X 22 , ⋯ , X 2 n 2 X_{21},X_{22},\cdots,X_{2n_2} X21,X22,⋯,X2n2 | X ˉ 2 ∙ \bar X_{2\bullet} Xˉ2∙ |
⋯ \cdots ⋯ | ⋯ \cdots ⋯ | ⋯ \cdots ⋯ |
A r : N ( μ r , σ 2 ) A_r: N(μ_r,σ^2) Ar:N(μr,σ2) | X r 1 , X r 2 , ⋯ , X r n r X_{r1},X_{r2},\cdots,X_{rn_r} Xr1,Xr2,⋯,Xrnr | X ˉ r ∙ \bar X_{r\bullet} Xˉr∙ |
由于
X
i
j
∼
N
(
μ
i
,
σ
2
)
X_{ij}∼ N(μ_i,σ^2)
Xij∼N(μi,σ2),即有
X
i
j
−
μ
i
∼
N
(
0
,
σ
2
)
X_{ij}-μ_i∼ N(0,σ^2)
Xij−μi∼N(0,σ2),故
ϵ
i
j
=
X
i
j
−
μ
i
ϵ_{ij}=X_{ij}-μ_i
ϵij=Xij−μi可看做随机误差。各总体间相互独立.。因此, 有如下的数学模型
{
X
i
j
=
μ
i
+
ϵ
i
j
ϵ
i
j
∼
N
(
0
,
σ
2
)
,
各
ϵ
i
j
独
立
i
=
1
,
2
,
⋯
,
r
j
=
1
,
2
,
⋯
,
n
i
\begin{cases} X_{ij}=μ_i+ϵ_{ij} \\ ϵ_{ij}∼ N(0,σ^2),\quad 各ϵ_{ij}独立 \\ i=1,2,\cdots,r \\ j=1,2,\cdots,n_i \end{cases}
⎩⎪⎪⎪⎨⎪⎪⎪⎧Xij=μi+ϵijϵij∼N(0,σ2),各ϵij独立i=1,2,⋯,rj=1,2,⋯,ni
记
n
=
∑
i
=
1
r
n
i
n=\displaystyle\sum_{i=1}^{r}n_i
n=i=1∑rni
方差分析的任务是对上述模型:
- 比较这 r r r 个总体的均值差异,即检验假设 H 0 : μ 1 = μ 2 = ⋯ = μ r , H 1 : μ 1 , μ 2 , ⋯ , μ r H_0:μ_1=μ_2=\cdots=μ_r,\quad H_1:μ_1,μ_2,\cdots,μ_r H0:μ1=μ2=⋯=μr,H1:μ1,μ2,⋯,μr不全相等
- 作出未知参数 μ 1 , μ 2 , ⋯ , μ r , σ 2 μ_1,μ_2,\cdots,μ_r,σ^2 μ1,μ2,⋯,μr,σ2 的估计
(1) 为了便于讨论,我们引入加权平均值
μ
=
1
n
∑
i
=
1
r
n
i
μ
i
μ=\displaystyle\frac{1}{n}\sum_{i=1}^{r}n_iμ_i
μ=n1i=1∑rniμi,称为总平均(total mean)。
再引入
δ
i
=
μ
i
−
μ
δ_i=μ_i-μ
δi=μi−μ,称为水平
A
i
A_i
Ai 的效应(effect)。此时
∑
i
=
1
r
n
i
δ
i
=
0
\displaystyle\sum_{i=1}^{r}n_iδ_i=0
i=1∑rniδi=0
于是上述模型可改写为
{
X
i
j
=
μ
+
δ
i
+
ϵ
i
j
ϵ
i
j
∼
N
(
0
,
σ
2
)
,
各
ϵ
i
j
独
立
i
=
1
,
2
,
⋯
,
r
j
=
1
,
2
,
⋯
,
n
i
∑
i
=
1
r
n
i
δ
i
=
0
\begin{cases} X_{ij}=μ+δ_i+ϵ_{ij} \\ ϵ_{ij}∼ N(0,σ^2),\quad 各ϵ_{ij}独立 \\ i=1,2,\cdots,r \\ j=1,2,\cdots,n_i \\ \displaystyle\sum_{i=1}^{r}n_iδ_i=0 \end{cases}
⎩⎪⎪⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎪⎪⎧Xij=μ+δi+ϵijϵij∼N(0,σ2),各ϵij独立i=1,2,⋯,rj=1,2,⋯,nii=1∑rniδi=0
假设等价于
H
0
:
δ
1
=
δ
2
=
⋯
=
δ
r
=
0
,
H
1
:
δ
1
,
δ
2
,
⋯
,
δ
r
H_0:δ_1=δ_2=\cdots=δ_r=0,\quad H_1:δ_1,δ_2,\cdots,δ_r
H0:δ1=δ2=⋯=δr=0,H1:δ1,δ2,⋯,δr不全为零
(2) 平方和分解
引入总偏差平方和
S
T
=
∑
i
=
1
r
∑
j
=
1
n
i
(
X
i
j
−
X
ˉ
)
2
\displaystyle S_T=\sum_{i=1}^{r}\sum_{j=1}^{n_i}(X_{ij}-\bar X)^2
ST=i=1∑rj=1∑ni(Xij−Xˉ)2,能反应全部数据间的差异,称为总离差平方和(total sum of deviation square)。其中
X
ˉ
=
1
n
∑
i
=
1
r
∑
j
=
1
n
i
X
i
j
\displaystyle\bar X=\frac{1}{n}\sum_{i=1}^{r}\sum_{j=1}^{n_i}X_{ij}
Xˉ=n1i=1∑rj=1∑niXij是数据的总平均。
S
T
S_T
ST 可分解为两部分
一部分是有因素
A
A
A 引起的差异,称为效应平方和(sum of squares of effects)
S
A
=
∑
i
=
1
r
(
X
i
∙
−
X
ˉ
)
2
S_A=\displaystyle\sum_{i=1}^{r}(X_{i\bullet}-\bar X)^2
SA=i=1∑r(Xi∙−Xˉ)2
另一部分则由随机误差所引起的差异,称为误差平方和(sum of squares of errors)
S
E
=
∑
i
=
1
r
∑
j
=
1
n
i
(
X
i
j
−
X
i
∙
)
2
S_E=\displaystyle\sum_{i=1}^{r}\sum_{j=1}^{n_i}(X_{ij}-X_{i\bullet})^2
SE=i=1∑rj=1∑ni(Xij−Xi∙)2
其中水平
A
i
A_i
Ai 下的样本均值
X
ˉ
i
∙
=
1
n
i
∑
j
=
1
n
i
X
i
j
\bar X_{i\bullet}=\displaystyle\frac{1}{n_i}\sum_{j=1}^{n_i}X_{ij}
Xˉi∙=ni1j=1∑niXij
(3) 拒绝域
S
T
=
S
A
+
S
E
S_T=S_A+S_E
ST=SA+SE
S
E
/
σ
2
∼
χ
2
(
n
−
r
)
S_E/σ^2∼ χ^2(n-r)
SE/σ2∼χ2(n−r)
S
A
S_A
SA与
S
E
S_E
SE相互独立
当
H
0
H_0
H0为真时,
S
A
/
σ
2
∼
χ
2
(
r
−
1
)
,
S
A
/
(
r
−
1
)
S
E
/
(
n
−
r
)
∼
F
(
r
−
1
,
n
−
r
)
S_A/σ^2∼ χ^2(r-1),\dfrac{S_A/(r-1)}{S_E/(n-r)}∼ F(r-1,n-r)
SA/σ2∼χ2(r−1),SE/(n−r)SA/(r−1)∼F(r−1,n−r)
单因素试验方差分析表
方差来源 | 平方和(SumSquare) | 自由度(df) | 均方(MS) | F |
---|---|---|---|---|
因素A(组间) | S A S_A SA | r − 1 r-1 r−1 | M S A = S A / ( r − 1 ) MS_A=S_A/(r-1) MSA=SA/(r−1) | M S A / M S E MS_A/MS_E MSA/MSE |
误差(组内) | S E S_E SE | n − r n-r n−r | M S E = S E / ( n − r ) MS_E=S_E/(n-r) MSE=SE/(n−r) | |
总和 | S T S_T ST | n − 1 n-1 n−1 |
当 F ⩾ F α ( r − 1 , n − r ) F⩾ F_α(r-1,n-r) F⩾Fα(r−1,n−r)时,拒绝原假设
(4) 未知参数估计
无论
H
0
H_0
H0 是否为真,
E
(
S
E
n
−
r
)
=
σ
2
E(\dfrac{S_E}{n-r})=σ^2
E(n−rSE)=σ2,
σ
2
σ^2
σ2的无偏估计
σ
^
2
=
S
E
/
(
n
−
r
)
=
M
S
E
\hat σ^2=S_E/(n-r)=MS_E
σ^2=SE/(n−r)=MSE
μ
i
μ_i
μi 的无偏估计
μ
^
i
=
X
ˉ
i
∙
(
i
=
1
,
2
,
⋯
,
r
)
\hat μ_i=\bar X_{i\bullet}\ (i=1,2,\cdots,r)
μ^i=Xˉi∙ (i=1,2,⋯,r)
(5) 当拒绝
H
0
H_0
H0时,只能说明均值不全相等。那么,它们中有没有部分是相等的?
比较
N
(
μ
i
,
σ
2
)
,
N
(
μ
j
,
σ
2
)
(
i
≠
j
)
N(μ_i,σ^2),N(μ_j,σ^2)\ (i\neq j)
N(μi,σ2),N(μj,σ2) (i=j) 差异的两个方法
a. 作
μ
i
−
μ
j
μ_i-μ_j
μi−μj 的区间估计
由于
E
(
X
ˉ
i
∙
−
X
ˉ
j
∙
)
=
μ
i
−
μ
j
,
D
(
X
ˉ
i
∙
−
X
ˉ
j
∙
)
=
σ
2
(
1
n
i
+
1
n
j
)
E(\bar X_{i\bullet}-\bar X_{j\bullet})=μ_i-μ_j,\ D(\bar X_{i\bullet}-\bar X_{j\bullet})=σ^2(\dfrac{1}{n_i}+\dfrac{1}{n_j})
E(Xˉi∙−Xˉj∙)=μi−μj, D(Xˉi∙−Xˉj∙)=σ2(ni1+nj1),且
X
ˉ
i
∙
−
X
ˉ
j
∙
\bar X_{i\bullet}-\bar X_{j\bullet}
Xˉi∙−Xˉj∙与
σ
^
2
=
S
E
/
(
n
−
r
)
=
M
S
E
\hat σ^2=S_E/(n-r)=MS_E
σ^2=SE/(n−r)=MSE相互独立
则
(
X
ˉ
i
∙
−
X
ˉ
j
∙
)
−
(
μ
i
−
μ
j
)
M
S
E
(
1
/
n
i
+
1
/
n
j
)
∼
t
(
n
−
r
)
\dfrac{(\bar X_{i\bullet}-\bar X_{j\bullet})-(μ_i-μ_j)}{\sqrt{MS_E(1/n_i+1/n_j)}}∼ t(n-r)
MSE(1/ni+1/nj)(Xˉi∙−Xˉj∙)−(μi−μj)∼t(n−r),由此得
μ
i
−
μ
j
μ_i-μ_j
μi−μj在
1
−
α
1-α
1−α置信水平下的置信区间为
(
(
X
ˉ
i
∙
−
X
ˉ
j
∙
)
±
t
α
/
2
(
n
−
r
)
M
S
E
(
1
/
n
i
+
1
/
n
j
)
)
((\bar X_{i\bullet}-\bar X_{j\bullet})± t_{α/2}(n-r)\sqrt{MS_E(1/n_i+1/n_j)})
((Xˉi∙−Xˉj∙)±tα/2(n−r)MSE(1/ni+1/nj))
b. 检验假设
H
0
:
μ
i
=
μ
j
,
H
1
:
μ
i
≠
μ
j
H_0:μ_i=μ_j,\quad H_1:μ_i\neqμ_j
H0:μi=μj,H1:μi=μj
检验统计量
t
=
X
ˉ
i
∙
−
X
ˉ
j
∙
M
S
E
(
1
/
n
i
+
1
/
n
j
)
t=\dfrac{\bar X_{i\bullet}-\bar X_{j\bullet}}{\sqrt{MS_E(1/n_i+1/n_j)}}
t=MSE(1/ni+1/nj)Xˉi∙−Xˉj∙
当
H
0
H_0
H0成立时,
t
∼
t
(
n
−
r
)
t∼ t(n-r)
t∼t(n−r)
拒绝域
W
=
{
∣
t
∣
⩾
t
α
/
2
(
n
−
r
)
}
W=\{|t|⩾ t_{α/2}(n-r)\}
W={∣t∣⩾tα/2(n−r)}
方差分析的前提
(1) 独立性:数据是来自r 个独立总体的简单随机样本.
(2) 正态性: r 个独立总体均为正态总体.
(3) 方差齐性: r 个正态总体的方差是相同
方差分析和其它统计推断一样,样本的独立性对方差分析是非常重要的,在实际应用中会经常遇到非随机样本的情况,这时使用方差分析得出的结论不可靠。
在实际中,几乎没有一个总体真正服从正态分布,而方差分析却依赖于正态性假设,不过根据经验,方差分析F检验对正态性的假设并不是非常敏感,即实际所得到的数据,若没有异常值和偏性, 或者说, 数据显示的分布比较对称的话, 即使样本容量比较小(如每个水平下的样本容量仅为5左右), 方差分析的结果仍是值得依赖的。
方差齐性对于方差分析是非常重要的, 方差齐性的检验可采用如下的经验准则:当最大样本标准差不超过最小样本标准差的两倍时, 方差分析F检验结果近似正确。
双因素方差分析(Two-factor analysis of variance)
双因素等重复试验的方差分析
设有两个因素作用于试验指标,因素
A
A
A 有
r
r
r 个水平
A
1
,
A
2
,
⋯
,
A
r
A_1,A_2,\cdots,A_r
A1,A2,⋯,Ar,因素
B
B
B 有
s
s
s 个水平
B
1
,
B
2
,
⋯
,
B
s
B_1,B_2,\cdots,B_s
B1,B2,⋯,Bs,在水平
(
A
i
,
B
j
)
(A_i,B_j)
(Ai,Bj)下都进行
t
(
t
⩾
2
)
t(t ⩾ 2)
t(t⩾2)次试验,得到如下结果
因素A / 因素B | B 1 B_1 B1 | B 2 B_2 B2 | ⋯ \cdots ⋯ | B s B_s Bs |
---|---|---|---|---|
A 1 A_1 A1 | X 111 , X 112 , ⋯ , X 11 t X_{111},X_{112},\cdots,X_{11t} X111,X112,⋯,X11t | X 121 , X 122 , ⋯ , X 12 t X_{121},X_{122},\cdots,X_{12t} X121,X122,⋯,X12t | ⋯ \cdots ⋯ | X 1 s 1 , X 1 s 2 , ⋯ , X 1 s t X_{1s1},X_{1s2},\cdots,X_{1st} X1s1,X1s2,⋯,X1st |
A 2 A_2 A2 | X 211 , X 212 , ⋯ , X 21 t X_{211},X_{212},\cdots,X_{21t} X211,X212,⋯,X21t | X 221 , X 222 , ⋯ , X 22 t X_{221},X_{222},\cdots,X_{22t} X221,X222,⋯,X22t | ⋯ \cdots ⋯ | X 2 s 1 , X 2 s 2 , ⋯ , X 2 s t X_{2s1},X_{2s2},\cdots,X_{2st} X2s1,X2s2,⋯,X2st |
⋮ \vdots ⋮ | ⋮ \vdots ⋮ | ⋮ \vdots ⋮ | ⋮ \vdots ⋮ | |
A r A_r Ar | X r 11 , X r 12 , ⋯ , X r 1 t X_{r11},X_{r12},\cdots,X_{r1t} Xr11,Xr12,⋯,Xr1t | X r 21 , X r 22 , ⋯ , X r 2 t X_{r21},X_{r22},\cdots,X_{r2t} Xr21,Xr22,⋯,Xr2t | ⋯ \cdots ⋯ | X r s 1 , X r s 2 , ⋯ , X r s t X_{rs1},X_{rs2},\cdots,X_{rst} Xrs1,Xrs2,⋯,Xrst |
并设
X
i
j
k
∼
N
(
μ
i
j
,
σ
2
)
X_{ijk}∼ N(μ_{ij},σ^2)
Xijk∼N(μij,σ2),且方差相等(均值和方差均未知),且不同水平下样本之间相互独立。或写成如下的数学模型
{
X
i
j
k
=
μ
i
j
+
ϵ
i
j
k
ϵ
i
j
k
∼
N
(
0
,
σ
2
)
,
各
ϵ
i
j
k
独
立
i
=
1
,
2
,
⋯
,
r
;
j
=
1
,
2
,
⋯
,
s
;
k
=
1
,
2
,
⋯
,
t
\begin{cases} X_{ijk}=μ_{ij}+ϵ_{ijk} \\ ϵ_{ijk}∼ N(0,σ^2),\quad 各ϵ_{ijk}独立 \\ i=1,2,\cdots,r ;\quad j=1,2,\cdots,s ;\quad k=1,2,\cdots,t \\ \end{cases}
⎩⎪⎨⎪⎧Xijk=μij+ϵijkϵijk∼N(0,σ2),各ϵijk独立i=1,2,⋯,r;j=1,2,⋯,s;k=1,2,⋯,t
引入记号
μ
=
1
r
s
∑
i
=
1
r
∑
j
=
1
s
μ
i
j
μ
i
∙
=
1
s
∑
j
=
1
s
μ
i
j
(
i
=
1
,
2
,
⋯
,
r
)
μ
∙
j
=
1
r
∑
i
=
1
r
μ
i
j
(
j
=
1
,
2
,
⋯
,
s
)
α
i
=
μ
i
∙
−
μ
(
i
=
1
,
2
,
⋯
,
r
)
β
j
=
μ
∙
j
−
μ
(
j
=
1
,
2
,
⋯
,
s
)
\displaystyle μ=\frac{1}{rs}\sum_{i=1}^{r}\sum_{j=1}^{s}μ_{ij} \\ μ_{i\bullet}=\frac{1}{s}\sum_{j=1}^{s}μ_{ij} \quad(i=1,2,\cdots,r)\\ μ_{\bullet j}=\frac{1}{r}\sum_{i=1}^{r}μ_{ij} \quad(j=1,2,\cdots,s)\\ α_i=μ_{i\bullet}-μ \quad(i=1,2,\cdots,r)\\ β_j=μ_{\bullet j}-μ \quad(j=1,2,\cdots,s)
μ=rs1i=1∑rj=1∑sμijμi∙=s1j=1∑sμij(i=1,2,⋯,r)μ∙j=r1i=1∑rμij(j=1,2,⋯,s)αi=μi∙−μ(i=1,2,⋯,r)βj=μ∙j−μ(j=1,2,⋯,s)
易见
∑
i
=
1
r
α
i
=
0
,
∑
j
=
1
s
β
j
=
0
\displaystyle\sum_{i=1}^{r}α_i=0,\quad\sum_{j=1}^{s}β_j=0
i=1∑rαi=0,j=1∑sβj=0
称
μ
μ
μ 为总平均,
α
i
α_i
αi为水平
A
i
A_i
Ai 的效应,
β
j
β_j
βj为水平
B
j
B_j
Bj 的效应,这样可将
μ
i
j
μ_{ij}
μij 表示为
μ
i
j
=
μ
+
α
i
+
β
j
+
γ
i
j
μ_{ij}=μ+α_i+β_j+γ_{ij}
μij=μ+αi+βj+γij
γ
i
j
=
μ
i
j
−
μ
i
∙
−
μ
∙
j
+
μ
γ_{ij}=μ_{ij}-μ_{i\bullet}-μ_{\bullet j}+μ
γij=μij−μi∙−μ∙j+μ 称为水平
A
i
A_i
Ai和水平
β
j
β_j
βj的交互效应(interaction)。这是有两个因素联合起来搭配作用而引起的。
易见
∑
i
=
1
r
γ
i
j
=
0
,
∑
j
=
1
s
γ
i
j
=
0
\displaystyle\sum_{i=1}^{r}γ_{ij}=0,\quad\sum_{j=1}^{s}γ_{ij}=0
i=1∑rγij=0,j=1∑sγij=0
于是上述模型可改写为
{
X
i
j
=
μ
+
α
i
+
β
j
+
γ
i
j
+
ϵ
i
j
k
ϵ
i
j
∼
N
(
0
,
σ
2
)
,
各
ϵ
i
j
独
立
i
=
1
,
2
,
⋯
,
r
;
j
=
1
,
2
,
⋯
,
s
;
k
=
1
,
2
,
⋯
,
t
∑
i
=
1
r
α
i
=
0
,
∑
j
=
1
s
β
j
=
0
,
∑
i
=
1
r
γ
i
j
=
0
,
∑
j
=
1
s
γ
i
j
=
0
\begin{cases} X_{ij}=μ+α_i+β_j+γ_{ij}+ϵ_{ijk} \\ ϵ_{ij}∼ N(0,σ^2),\quad 各ϵ_{ij}独立 \\ i=1,2,\cdots,r ;\quad j=1,2,\cdots,s ;\quad k=1,2,\cdots,t \\ \displaystyle\sum_{i=1}^{r}α_i=0,\quad\sum_{j=1}^{s}β_j=0 ,\quad\sum_{i=1}^{r}γ_{ij}=0,\quad\sum_{j=1}^{s}γ_{ij}=0 \end{cases}
⎩⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎧Xij=μ+αi+βj+γij+ϵijkϵij∼N(0,σ2),各ϵij独立i=1,2,⋯,r;j=1,2,⋯,s;k=1,2,⋯,ti=1∑rαi=0,j=1∑sβj=0,i=1∑rγij=0,j=1∑sγij=0
其中
μ
,
α
i
,
β
j
,
γ
i
j
μ,α_i,β_j,γ_{ij}
μ,αi,βj,γij 都是未知参数。
对于这个模型我们要检验三个假设:
H
01
:
α
1
=
α
2
=
⋯
=
α
r
=
0
,
H
11
:
α
1
,
α
2
,
⋯
,
α
r
H_{01}:α_1=α_2=\cdots=α_r=0,\quad H_{11}:α_1,α_2,\cdots,α_r
H01:α1=α2=⋯=αr=0,H11:α1,α2,⋯,αr不全为零
H
02
:
β
1
=
β
2
=
⋯
=
β
s
=
0
,
H
12
:
β
1
,
β
2
,
⋯
,
β
s
H_{02}:β_1=β_2=\cdots=β_s=0,\quad H_{12}:β_1,β_2,\cdots,β_s
H02:β1=β2=⋯=βs=0,H12:β1,β2,⋯,βs不全为零
H
03
:
γ
i
j
=
0
,
H
13
:
γ
i
j
H_{03}:γ_{ij}=0,\quad H_{13}:γ_{ij}
H03:γij=0,H13:γij不全为零
与单因素类似,这些问题也建立在平方和分解上,先引入
X
ˉ
=
1
r
s
t
∑
i
=
1
r
∑
j
=
1
s
∑
k
=
1
t
X
i
j
k
X
ˉ
i
j
∙
=
1
t
∑
k
=
1
t
X
i
j
k
(
i
=
1
,
2
,
⋯
,
r
;
j
=
1
,
2
,
⋯
,
s
)
X
ˉ
i
∙
∙
=
1
s
t
∑
j
=
1
s
∑
k
=
1
t
X
i
j
k
(
i
=
1
,
2
,
⋯
,
r
)
X
ˉ
∙
j
∙
=
1
r
t
∑
i
=
1
r
∑
k
=
1
t
X
i
j
k
(
j
=
1
,
2
,
⋯
,
s
)
\displaystyle \bar X=\frac{1}{rst}\sum_{i=1}^{r}\sum_{j=1}^{s}\sum_{k=1}^{t}X_{ijk} \\ \bar X_{ij\bullet}=\frac{1}{t}\sum_{k=1}^{t}X_{ijk} \quad(i=1,2,\cdots,r ;\quad j=1,2,\cdots,s)\\ \bar X_{i\bullet\bullet}=\frac{1}{st}\sum_{j=1}^{s}\sum_{k=1}^{t}X_{ijk} \quad(i=1,2,\cdots,r)\\ \bar X_{\bullet j\bullet}=\frac{1}{rt}\sum_{i=1}^{r}\sum_{k=1}^{t}X_{ijk} \quad(j=1,2,\cdots,s)
Xˉ=rst1i=1∑rj=1∑sk=1∑tXijkXˉij∙=t1k=1∑tXijk(i=1,2,⋯,r;j=1,2,⋯,s)Xˉi∙∙=st1j=1∑sk=1∑tXijk(i=1,2,⋯,r)Xˉ∙j∙=rt1i=1∑rk=1∑tXijk(j=1,2,⋯,s)
总偏差平方和
S
T
=
∑
i
=
1
r
∑
j
=
1
s
∑
k
=
1
t
(
X
i
j
k
−
X
ˉ
)
2
\displaystyle S_T=\sum_{i=1}^{r}\sum_{j=1}^{s}\sum_{k=1}^{t}(X_{ijk}-\bar X)^2
ST=i=1∑rj=1∑sk=1∑t(Xijk−Xˉ)2
平方和的分解式
S
T
=
S
E
+
S
A
+
S
B
+
S
A
×
B
S_T=S_E+S_A+S_B+S_{A× B}
ST=SE+SA+SB+SA×B
S
E
=
∑
i
=
1
r
∑
j
=
1
s
∑
k
=
1
t
(
X
i
j
k
−
X
ˉ
i
j
∙
)
2
S
A
=
s
t
∑
i
=
1
r
(
X
ˉ
i
∙
∙
−
X
ˉ
)
2
S
B
=
r
t
∑
j
=
1
s
(
X
ˉ
∙
j
∙
−
X
ˉ
)
2
S
A
×
B
=
t
∑
i
=
1
r
∑
j
=
1
s
(
X
ˉ
i
j
∙
−
X
ˉ
i
∙
∙
−
X
ˉ
∙
j
∙
+
X
ˉ
)
2
\displaystyle S_E=\displaystyle\sum_{i=1}^{r}\sum_{j=1}^{s}\sum_{k=1}^{t}( X_{ijk}-\bar X_{ij\bullet})^2 \\ S_A=st\sum_{i=1}^{r}(\bar X_{i\bullet\bullet}-\bar X)^2 \\ S_B=rt\sum_{j=1}^{s}(\bar X_{\bullet j\bullet}-\bar X)^2 \\ S_{A× B}=t\sum_{i=1}^{r}\sum_{j=1}^{s}(\bar X_{ij\bullet}-\bar X_{i\bullet\bullet}-\bar X_{\bullet j\bullet}+\bar X)^2
SE=i=1∑rj=1∑sk=1∑t(Xijk−Xˉij∙)2SA=sti=1∑r(Xˉi∙∙−Xˉ)2SB=rtj=1∑s(Xˉ∙j∙−Xˉ)2SA×B=ti=1∑rj=1∑s(Xˉij∙−Xˉi∙∙−Xˉ∙j∙+Xˉ)2
S
E
S_E
SE称为误差平方和,
S
A
,
S
B
S_A,S_B
SA,SB称为因素
A
,
B
A,B
A,B 的效应平方和,
S
A
×
B
S_{A× B}
SA×B称为
A
,
B
A,B
A,B 交互效应平方和。
且可以证明下列
E
(
S
E
r
s
(
t
−
1
)
)
=
σ
2
E
(
S
A
r
−
1
)
=
σ
2
+
s
t
r
−
1
∑
i
=
1
r
α
i
2
E
(
S
B
s
−
1
)
=
σ
2
+
r
t
s
−
1
∑
j
=
1
s
β
j
2
E
(
S
A
×
B
(
r
−
1
)
(
s
−
1
)
)
=
σ
2
+
t
(
r
−
1
)
(
s
−
1
)
∑
i
=
1
r
∑
j
=
1
s
γ
i
j
2
\displaystyle E\left(\dfrac{S_E}{rs(t-1)}\right)=σ^2\\ E\left(\dfrac{S_A}{r-1}\right)=σ^2+\dfrac{st}{r-1}\sum_{i=1}^{r}α_i^2 \\ E\left(\dfrac{S_B}{s-1}\right)=σ^2+\dfrac{rt}{s-1}\sum_{j=1}^{s}β_j^2 \\ E\left(\dfrac{S_{A× B}}{(r-1)(s-1)}\right)=σ^2+\dfrac{t}{(r-1)(s-1)}\sum_{i=1}^{r}\sum_{j=1}^{s}γ_{ij}^2
E(rs(t−1)SE)=σ2E(r−1SA)=σ2+r−1sti=1∑rαi2E(s−1SB)=σ2+s−1rtj=1∑sβj2E((r−1)(s−1)SA×B)=σ2+(r−1)(s−1)ti=1∑rj=1∑sγij2
双因素试验方差分析表
方差来源 | 平方和 SumSquare | 自由度 df | 均方 MS | F |
---|---|---|---|---|
因素A(组间) | S A S_A SA | r − 1 r-1 r−1 | M S A = S A ( r − 1 ) MS_A=\dfrac{S_A}{(r-1)} MSA=(r−1)SA | F A = M S A / M S E F_A=MS_A/MS_E FA=MSA/MSE |
因素B(组间) | S B S_B SB | s − 1 s-1 s−1 | M S B = S B ( s − 1 ) MS_B=\dfrac{S_B}{(s-1)} MSB=(s−1)SB | F B = M S B / M S E F_B=MS_B/MS_E FB=MSB/MSE |
交互作用 | S A × B S_{A× B} SA×B | ( r − 1 ) ( s − 1 ) (r-1)(s-1) (r−1)(s−1) | M S A × B = S A ( r − 1 ) ( s − 1 ) MS_{A× B}=\dfrac{S_A}{(r-1)(s-1)} MSA×B=(r−1)(s−1)SA | F A × B = M S A × B / M S E F_{A× B}=MS_{A× B}/MS_E FA×B=MSA×B/MSE |
误差 (组内) | S E S_E SE | r s ( t − 1 ) rs(t-1) rs(t−1) | M S E = S E r s ( t − 1 ) MS_E=\dfrac{S_E}{rs(t-1)} MSE=rs(t−1)SE | |
总和 | S T S_T ST | r s t − 1 rst-1 rst−1 |
当
H
01
H_{01}
H01为真时,可以证明
F
A
∼
F
(
r
−
1
,
r
s
(
t
−
1
)
)
F
B
∼
F
(
s
−
1
,
r
s
(
t
−
1
)
)
F
A
×
B
∼
F
(
(
r
−
1
)
(
s
−
1
)
,
r
s
(
t
−
1
)
)
F_A∼ F(r-1,rs(t-1)) \\ F_B∼ F(s-1,rs(t-1)) \\ F_{A× B}∼ F((r-1)(s-1),rs(t-1))
FA∼F(r−1,rs(t−1))FB∼F(s−1,rs(t−1))FA×B∼F((r−1)(s−1),rs(t−1))
假设
H
01
H_{01}
H01的拒绝域为
F
A
⩾
F
α
(
r
−
1
,
r
s
(
t
−
1
)
)
F_A ⩾ F_α(r-1,rs(t-1))
FA⩾Fα(r−1,rs(t−1))
假设
H
02
H_{02}
H02的拒绝域为
F
B
⩾
F
α
(
s
−
1
,
r
s
(
t
−
1
)
)
F_B ⩾ F_α(s-1,rs(t-1))
FB⩾Fα(s−1,rs(t−1))
假设
H
03
H_{03}
H03的拒绝域为
F
A
×
B
⩾
F
α
(
(
r
−
1
)
(
s
−
1
)
,
r
s
(
t
−
1
)
)
F_{A× B} ⩾ F_α((r-1)(s-1),rs(t-1))
FA×B⩾Fα((r−1)(s−1),rs(t−1))
双因素无重复试验的方差分析
在以上的讨论中,我们考虑了双因素实验中两个因素的交互作用,为了检验交互作用的效应是否显著,对于两个因素的每个组合
(
A
i
,
B
j
)
(A_i,B_j)
(Ai,Bj) 至少要做两次试验,如果在实际问题中已经知道不存在交互作用,或效应很小,此时可以不考虑交互作用。每组制作一次试验
k
=
1
k=1
k=1 也能分析因素
A
,
B
A,B
A,B 的效应,此时称为无重复实验的方差分析。
回归分析(Regression Analysis)
在客观世界中,普遍存在着变量之间的额关系。变量之间的关系一般可分为确定性关系和相关性关系。
确定性关系(deterministic):
当自变量给定一个值时,就确定应变量的值与之对应。如:在自由落体中,物体下落的高度
h
h
h与下落时间
t
t
t之间有函数关系:
h
=
1
2
g
t
2
h=\frac{1}{2}gt^2
h=21gt2
相关性关系(correlation):
变量之间的关系并不确定,而是表现为具有随机性的一种“趋势”。即对自变量
x
x
x 的同一值,在不同的观测中,因变量
Y
Y
Y可以取不同的值,而且取值是随机的,但对应
x
x
x在一定范围的不同值,对
Y
Y
Y进行观测时,可以观察到
Y
Y
Y随
x
x
x的变化而呈现有一定趋势的变化。回归分析是研究相关关系得一种数学工具。
一元线性回归(Simple Linear Regression)
设随机变量
Y
Y
Y(因变量)与普通变量
x
x
x(自变量)存在着相关关系,对于
x
x
x的确定值,
Y
Y
Y有他的分布(图中
C
1
,
C
2
C_1,C_2
C1,C2分别是
x
1
,
x
2
x_1,x_2
x1,x2处
Y
Y
Y 的概率密度曲线),用
F
(
y
∣
x
)
F(y|x)
F(y∣x) 表示
Y
Y
Y的分布函数,然而此函数讨论过于复杂,我们转而讨论
E
(
Y
)
=
μ
(
x
)
E(Y)=μ(x)
E(Y)=μ(x)与
x
x
x的关系,称为
Y
Y
Y 关于
x
x
x 的回归函数(regression function)。同时以回归函数作为
Y
Y
Y 的近似,其均方误差
E
[
Y
−
μ
(
x
)
]
2
E[Y-μ(x)]^2
E[Y−μ(x)]2 最小。
对从整体
(
x
,
Y
)
(x,Y)
(x,Y) 中抽取的一个样本
(
x
1
,
Y
1
)
,
(
x
2
,
Y
2
)
,
⋯
,
(
x
n
,
Y
n
)
(x_1,Y_1),(x_2,Y_2),\cdots,(x_n,Y_n)
(x1,Y1),(x2,Y2),⋯,(xn,Yn),样本值为
(
x
1
,
y
1
)
,
(
x
2
,
y
2
)
,
⋯
,
(
x
n
,
y
n
)
(x_1,y_1),(x_2,y_2),\cdots,(x_n,y_n)
(x1,y1),(x2,y2),⋯,(xn,yn),若
μ
(
x
)
μ(x)
μ(x) 为线性函数
μ
(
x
)
=
a
+
b
x
μ(x)=a+bx
μ(x)=a+bx,此时的问题称为一元线性回归问题。
Y
∼
N
(
a
+
b
x
,
σ
2
)
Y∼ N(a+bx,σ^2)
Y∼N(a+bx,σ2),设
ϵ
=
Y
−
(
a
+
b
x
)
ϵ=Y-(a+bx)
ϵ=Y−(a+bx)则
Y
=
a
+
b
x
+
ϵ
,
ϵ
∼
N
(
0
,
σ
2
)
Y=a+bx+ϵ,\quad ϵ∼ N(0,σ^2)
Y=a+bx+ϵ,ϵ∼N(0,σ2)其中
a
,
b
,
σ
2
a,b,σ^2
a,b,σ2 都是不依赖于
x
x
x 的未知参数,因变量
Y
Y
Y 由线性函数和随机误差
ϵ
ϵ
ϵ 组成。
一元线性回归要解决的问题:
参数估计:
{
a
,
b
的
估
计
σ
2
的
估
计
\begin{cases} a,b 的估计 \\ σ^2 的估计 \end{cases}
{a,b的估计σ2的估计
参数检验及模型应用:
{
线
性
假
设
的
显
著
性
检
验
回
归
系
数
b
的
置
信
区
间
Y
的
点
预
测
\begin{cases} 线性假设的显著性检验 \\ 回归系数b的置信区间 \\ Y 的点预测 \end{cases}
⎩⎪⎨⎪⎧线性假设的显著性检验回归系数b的置信区间Y的点预测
(1)
a
,
b
a,b
a,b 的估计:最小二乘法
由一元线性回归模型得
Y
i
=
a
+
b
x
i
+
ϵ
,
ϵ
i
∼
N
(
0
,
σ
2
)
Y_i=a+bx_i+ϵ,\quad ϵ_i∼ N(0,σ^2)
Yi=a+bxi+ϵ,ϵi∼N(0,σ2)且相互独立,用最大似然法来估计未知参数。只需函数取
Q
(
a
,
b
)
=
∑
i
=
1
n
(
y
i
−
a
−
b
x
i
)
2
Q(a,b)=\displaystyle\sum_{i=1}^{n}(y_i-a-bx_i)^2
Q(a,b)=i=1∑n(yi−a−bxi)2取最小值,若
Y
Y
Y 不是正态分布,使
Y
Y
Y 的观察值
y
i
y_i
yi 与
a
+
b
x
i
a+bx_i
a+bxi 偏差的平方和为最小,这种方法叫最小二乘法(least square method)。
求估计值
a
^
,
b
^
\hat a,\hat b
a^,b^ ,使
Q
(
a
^
,
b
^
)
=
min
a
,
b
Q
(
a
,
b
)
\displaystyle Q(\hat a,\hat b)=\min_{a,b}Q(a,b)
Q(a^,b^)=a,bminQ(a,b)
∂
Q
∂
a
=
0
,
∂
Q
∂
b
=
0
\dfrac{∂ Q}{∂ a}=0,\dfrac{∂ Q}{∂ b}=0
∂a∂Q=0,∂b∂Q=0得方程组
{
n
a
+
(
∑
i
=
1
n
x
i
)
b
=
∑
i
=
1
n
y
i
(
∑
i
=
1
n
x
i
)
a
+
(
∑
i
=
1
n
x
i
2
)
b
=
∑
i
=
1
n
x
i
y
i
\begin{cases} \displaystyle na+(\sum_{i=1}^{n}x_i)b=\sum_{i=1}^{n}y_i \\ \displaystyle (\sum_{i=1}^{n}x_i)a+(\sum_{i=1}^{n}x_i^2)b=\sum_{i=1}^{n}x_iy_i \end{cases}
⎩⎪⎪⎪⎨⎪⎪⎪⎧na+(i=1∑nxi)b=i=1∑nyi(i=1∑nxi)a+(i=1∑nxi2)b=i=1∑nxiyi称为正规方程组(normal system of equations)。由于
x
i
x_i
xi不全相同,系数行列式不为零,故上述方程有唯一解。则估计值为
{
b
^
=
S
x
y
S
x
x
a
^
=
y
ˉ
−
b
^
x
ˉ
\begin{cases} \hat b=\dfrac{S_{xy}}{S_{xx}} \\ \hat a=\bar y-\hat b\bar x \end{cases}
⎩⎨⎧b^=SxxSxya^=yˉ−b^xˉ其中
S
x
x
=
∑
i
=
1
n
(
x
i
−
x
ˉ
)
2
S
y
y
=
∑
i
=
1
n
(
y
i
−
y
ˉ
)
2
S
x
y
=
∑
i
=
1
n
(
x
i
−
x
ˉ
)
(
y
i
−
y
ˉ
)
x
ˉ
=
1
n
∑
i
=
1
n
x
i
,
y
ˉ
=
1
n
∑
i
=
1
n
y
i
\displaystyle S_{xx}=\sum_{i=1}^{n}(x_i-\bar x)^2 \\ S_{yy}=\sum_{i=1}^{n}(y_i-\bar y)^2 \\ S_{xy}=\sum_{i=1}^{n}(x_i-\bar x)(y_i-\bar y) \\ \bar x=\dfrac{1}{n}\sum_{i=1}^{n}x_i,\quad\bar y=\dfrac{1}{n}\sum_{i=1}^{n}y_i
Sxx=i=1∑n(xi−xˉ)2Syy=i=1∑n(yi−yˉ)2Sxy=i=1∑n(xi−xˉ)(yi−yˉ)xˉ=n1i=1∑nxi,yˉ=n1i=1∑nyi
由此得经验回归方程(Empirical regression equation)
y
^
=
a
^
+
b
^
x
\hat y=\hat a+\hat bx
y^=a^+b^x
(2)
σ
2
σ^2
σ2 的估计
记
e
i
=
y
i
−
y
^
e_i=y_i-\hat y
ei=yi−y^为残差(residual),注意到
σ
2
=
D
(
ϵ
)
+
E
(
ϵ
2
)
=
E
(
ϵ
2
)
σ^2=D(ϵ)+E(ϵ^2)=E(ϵ^2)
σ2=D(ϵ)+E(ϵ2)=E(ϵ2)
用残差平方和(residual sum of squares)
Q
e
=
∑
i
=
1
n
(
y
i
−
y
^
i
)
2
\displaystyle Q_e=\sum_{i=1}^{n}(y_i-\hat y_i)^2
Qe=i=1∑n(yi−y^i)2估计
σ
2
σ^2
σ2
可以证明
Q
e
/
σ
2
∼
χ
2
(
n
−
2
)
Q_e/σ^2∼ χ^2(n-2)
Qe/σ2∼χ2(n−2),于是
E
(
Q
e
)
=
(
n
−
2
)
σ
2
E(Q_e)=(n-2)σ^2
E(Qe)=(n−2)σ2
因此
σ
^
2
=
Q
e
/
(
n
−
2
)
\hat σ^2=Q_e/(n-2)
σ^2=Qe/(n−2) 是
σ
2
σ^2
σ2的无偏估计。
(3) 线性假设的显著性检验
检验假设
H
0
:
b
=
0
,
H
1
:
b
≠
0
H_0:b=0,\quad H_1:b\neq0
H0:b=0,H1:b=0
若原假设被拒绝,说明回归效果是显著的。否则,若接受原假设,说明 Y与x不是线性关系,回归方程无意义。回归效果不显著的原因可能有以下几种:
a. 影响Y取值的,除了x,还有其他不可忽略的因素;
b. E(Y)与x的关系不是线性关系,而是其他关系;
c. Y与x不存在关系。
我们用 t 检验法来检验,可证明
b
^
∼
N
(
b
,
σ
2
/
S
x
x
)
⟹
(
n
−
2
)
σ
^
2
/
σ
2
=
Q
e
/
σ
2
∼
χ
2
(
n
−
2
)
\hat b∼ N(b,σ^2/S_{xx}) \implies (n-2)\hat σ^2/σ^2=Q_e/σ^2∼ χ^2(n-2)
b^∼N(b,σ2/Sxx)⟹(n−2)σ^2/σ2=Qe/σ2∼χ2(n−2)且
b
^
\hat b
b^与
Q
e
Q_e
Qe独立。故可推得
(
b
^
−
b
)
S
x
x
/
σ
^
∼
t
(
n
−
2
)
(\hat b-b)\sqrt{S_{xx}}/\hat σ∼ t(n-2)
(b^−b)Sxx/σ^∼t(n−2)
若
H
0
H_0
H0 为真,
t
=
b
^
S
x
x
/
σ
^
∼
t
(
n
−
2
)
t=\hat b\sqrt{S_{xx}}/\hat σ∼ t(n-2)
t=b^Sxx/σ^∼t(n−2),且
E
(
b
^
)
=
b
=
0
E(\hat b)=b=0
E(b^)=b=0既得
H
0
H_0
H0的拒绝域为
∣
t
∣
⩾
t
α
/
2
(
n
−
2
)
|t|⩾ t_{α/2}(n-2)
∣t∣⩾tα/2(n−2)
(4) 回归系数(coefficient)
b
b
b 的置信区间
由上面推论,
b
b
b 的置信水平为
1
−
α
1-α
1−α 的置信区间为:
b
^
±
t
α
/
2
(
n
−
2
)
σ
^
/
S
x
x
\hat b± t_{α/2}(n-2)\hat σ/\sqrt{S_{xx}}
b^±tα/2(n−2)σ^/Sxx
(5) 回归函数
Y
=
a
+
b
x
Y=a+bx
Y=a+bx 的点估计和置信区间
考虑估计量
Y
^
0
=
a
^
+
b
^
x
0
\hat Y_0=\hat a+\hat bx_0
Y^0=a^+b^x0,可求得
E
(
Y
^
0
)
=
a
+
b
x
0
E(\hat Y_0)=a+bx_0
E(Y^0)=a+bx0为无偏估计。
并且
Y
^
0
−
(
a
+
b
x
0
)
σ
^
1
/
n
+
(
x
0
−
x
ˉ
)
2
/
S
x
x
∼
t
(
n
−
2
)
\dfrac{\hat Y_0-(a+bx_0)}{\hat σ\sqrt{1/n+(x_0-\bar x)^2/S_{xx}}}∼ t(n-2)
σ^1/n+(x0−xˉ)2/SxxY^0−(a+bx0)∼t(n−2),于是得到
a
+
b
x
0
a+bx_0
a+bx0 的置信区间为
(
Y
^
0
±
t
α
/
2
(
n
−
2
)
σ
^
1
n
+
(
x
0
−
x
ˉ
)
2
S
x
x
)
\left(\hat Y_0± t_{α/2}(n-2)\hat σ\sqrt{\dfrac{1}{n}+\dfrac{(x_0-\bar x)^2}{S_{xx}}}\right)
(Y^0±tα/2(n−2)σ^n1+Sxx(x0−xˉ)2)
这一置信区间的长度是
x
0
x_0
x0的函数,他随
∣
x
0
−
x
ˉ
∣
|x_0-\bar x|
∣x0−xˉ∣的增加而增加。
(6)
Y
Y
Y 的观察值
Y
0
Y_0
Y0 的点估计和置信区间
观察值
Y
0
=
a
+
b
x
0
+
ϵ
0
,
ϵ
0
∼
N
(
0
,
σ
2
)
Y_0=a+bx_0+ϵ_0,\quad ϵ_0∼ N(0,σ^2)
Y0=a+bx0+ϵ0,ϵ0∼N(0,σ2)
⟹
Y
^
0
=
Y
ˉ
+
b
^
(
x
0
−
x
ˉ
)
\implies \hat Y_0=\bar Y+\hat b(x_0-\bar x)
⟹Y^0=Yˉ+b^(x0−xˉ)
⟹
Y
^
0
−
(
a
+
b
x
0
)
σ
^
1
+
1
/
n
+
(
x
0
−
x
ˉ
)
2
/
S
x
x
∼
t
(
n
−
2
)
\implies\dfrac{\hat Y_0-(a+bx_0)}{\hat σ\sqrt{1+1/n+(x_0-\bar x)^2/S_{xx}}}∼ t(n-2)
⟹σ^1+1/n+(x0−xˉ)2/SxxY^0−(a+bx0)∼t(n−2)
于是得
Y
0
Y_0
Y0 的预测区间(prediction interval)
(
Y
^
0
±
t
α
/
2
(
n
−
2
)
σ
^
1
+
1
n
+
(
x
0
−
x
ˉ
)
2
S
x
x
)
\left(\hat Y_0± t_{α/2}(n-2)\hat σ\sqrt{1+\dfrac{1}{n}+\dfrac{(x_0-\bar x)^2}{S_{xx}}}\right)
(Y^0±tα/2(n−2)σ^1+n1+Sxx(x0−xˉ)2)
回归分析表
Coef. | Std. error | t Stat | P value | Lower 95% | Upper 95% | |
---|---|---|---|---|---|---|
I n t e r c e p t \rm Intercept Intercept | a ^ \hat a a^ | |||||
X \rm X X | b ^ \hat b b^ |
可化为一元线性回归的例子
方程 | 条件 | 转化 |
---|---|---|
Y = a e b x ⋅ ϵ Y=ae^{bx}\cdot ϵ Y=aebx⋅ϵ | ln ϵ ∼ N ( 0 , σ 2 ) \ln ϵ∼ N(0,σ^2) lnϵ∼N(0,σ2) | ln Y = ln a + ln b x + ln ϵ \ln Y=\ln a+ \ln bx+\ln ϵ lnY=lna+lnbx+lnϵ |
Y = a x b ⋅ ϵ Y=ax^{b}\cdot ϵ Y=axb⋅ϵ | ln ϵ ∼ N ( 0 , σ 2 ) \ln ϵ∼ N(0,σ^2) lnϵ∼N(0,σ2) | ln Y = ln a + b ln x + ln ϵ \ln Y=\ln a+ b\ln x+\ln ϵ lnY=lna+blnx+lnϵ |
Y = a + b h ( x ) + ϵ Y=a+bh(x)+ ϵ Y=a+bh(x)+ϵ | ϵ ∼ N ( 0 , σ 2 ) ϵ∼ N(0,σ^2) ϵ∼N(0,σ2) | Y = a + b x ′ + ϵ , x ′ = h ( x ) Y=a+bx'+ ϵ,\quad x'=h(x) Y=a+bx′+ϵ,x′=h(x) |
如果回归函数
μ
(
x
;
θ
1
,
θ
2
,
⋯
,
θ
p
)
μ(x;θ_1,θ_2,\cdots,θ_p)
μ(x;θ1,θ2,⋯,θp) 是参数
θ
1
,
θ
2
,
⋯
,
θ
p
θ_1,θ_2,\cdots,θ_p
θ1,θ2,⋯,θp 的线性函数(不必是
x
x
x 的线性函数),则称线性回归模型(linear regressive model),若为非线性函数则称非线性回归模型(non-linear regressive model),但是他们都能转化为线性回归模型。
又如
Y
=
θ
1
e
θ
1
x
⋅
ϵ
,
ϵ
∼
N
(
0
,
σ
2
)
Y=θ_1e^{θ_1x}\cdot ϵ,\quad ϵ∼ N(0,σ^2)
Y=θ1eθ1x⋅ϵ,ϵ∼N(0,σ2),他不能转化为线性回归模型,称为本质的非线性回归模型。又如
Y
=
(
θ
1
+
θ
2
x
+
θ
3
x
2
)
−
1
+
ϵ
(
H
o
l
i
d
a
y
模
型
)
Y
=
θ
1
1
+
e
x
p
(
θ
2
+
θ
3
x
)
+
ϵ
(
L
o
g
i
s
t
i
c
模
型
)
Y=(θ_1+θ_2x+θ_3x^2)^{-1}+ϵ\quad(Holiday 模型)\\ Y=\dfrac{θ_1}{1+exp(θ_2+θ_3^x)}+ϵ\quad(Logistic 模型)
Y=(θ1+θ2x+θ3x2)−1+ϵ(Holiday模型)Y=1+exp(θ2+θ3x)θ1+ϵ(Logistic模型)
都是本质的非线性回归模型。
多元线性回归(Multiple Linear Regression)
在实际问题中,随机变量
Y
Y
Y 往往与多个变量
x
1
,
x
2
,
⋯
,
x
p
(
p
>
1
)
x_1,x_2,\cdots,x_p(p>1)
x1,x2,⋯,xp(p>1) 有关,若
Y
Y
Y 的数学期望存在,记为
μ
(
x
1
,
x
2
,
⋯
,
x
p
)
μ(x_1,x_2,\cdots,x_p)
μ(x1,x2,⋯,xp)它就是
Y
Y
Y 关于
x
x
x的回归函数为,再次,我们只讨论线性回归模型
Y
=
b
0
+
b
1
x
1
+
b
2
x
2
=
⋯
+
b
p
x
p
+
ϵ
,
ϵ
∼
N
(
0
,
σ
2
)
Y=b_0+b_1x_1+b_2x_2=\cdots+b_px_p+ϵ,\quad ϵ∼ N(0,σ^2)
Y=b0+b1x1+b2x2=⋯+bpxp+ϵ,ϵ∼N(0,σ2) 其中
b
0
,
b
1
,
b
2
,
⋯
,
b
p
b_0,b_1,b_2,\cdots,b_p
b0,b1,b2,⋯,bp 都是与
x
1
,
x
2
,
⋯
,
x
p
x_1,x_2,\cdots,x_p
x1,x2,⋯,xp 无关的未知参数。
设
(
x
11
,
x
12
,
⋯
,
x
1
p
)
,
⋯
,
(
x
n
1
,
x
n
2
,
⋯
,
x
n
p
)
(x_{11},x_{12},\cdots,x_{1p}),\cdots,(x_{n1},x_{n2},\cdots,x_{np})
(x11,x12,⋯,x1p),⋯,(xn1,xn2,⋯,xnp) 是一个样本,我们用最小二乘法来估计未知参数,只需函数
Q
(
b
0
,
b
1
,
b
2
,
⋯
,
b
p
)
=
∑
i
=
1
n
[
y
i
−
(
b
0
+
b
1
x
1
+
b
2
x
2
=
⋯
+
b
p
x
p
)
]
2
Q(b_0,b_1,b_2,\cdots,b_p)=\displaystyle\sum_{i=1}^{n}[y_i-(b_0+b_1x_1+b_2x_2=\cdots+b_px_p)]^2
Q(b0,b1,b2,⋯,bp)=i=1∑n[yi−(b0+b1x1+b2x2=⋯+bpxp)]2取最小值,即估计值满足
Q
(
b
^
0
,
b
^
1
,
b
^
2
,
⋯
,
b
^
p
)
=
min
Q
(
b
0
,
b
1
,
b
2
,
⋯
,
b
p
)
\displaystyle Q(\hat b_0,\hat b_1,\hat b_2,\cdots,\hat b_p)=\min Q(b_0,b_1,b_2,\cdots,b_p )
Q(b^0,b^1,b^2,⋯,b^p)=minQ(b0,b1,b2,⋯,bp)求偏导数
∂
Q
∂
b
0
=
0
,
∂
Q
∂
b
j
=
0
(
j
=
1
,
2
,
⋯
,
p
)
\dfrac{∂ Q}{∂ b_0}=0,\dfrac{∂ Q}{∂ b_j}=0(j=1,2,\cdots,p)
∂b0∂Q=0,∂bj∂Q=0(j=1,2,⋯,p)得方程组
{
b
0
n
+
b
1
∑
i
=
1
n
x
i
1
+
b
2
∑
i
=
1
n
x
i
2
+
⋯
+
b
p
∑
i
=
1
n
x
i
p
=
∑
i
=
1
n
y
i
b
0
∑
i
=
1
n
x
i
1
+
b
1
∑
i
=
1
n
x
i
1
2
+
b
2
∑
i
=
1
n
x
i
1
∑
i
=
1
n
x
i
2
+
⋯
+
b
p
∑
i
=
1
n
x
i
1
∑
i
=
1
n
x
i
p
=
∑
i
=
1
n
x
i
1
y
i
⋯
b
0
∑
i
=
1
n
x
i
p
+
b
1
∑
i
=
1
n
x
i
1
∑
i
=
1
n
x
i
p
+
b
2
∑
i
=
1
n
x
i
p
∑
i
=
1
n
x
i
2
+
⋯
+
b
p
∑
i
=
1
n
x
i
p
2
=
∑
i
=
1
n
x
i
p
y
i
\begin{cases} \displaystyle b_0n+b_1\sum_{i=1}^{n}x_{i1}+b_2\sum_{i=1}^{n}x_{i2}+\cdots+b_p\sum_{i=1}^{n}x_{ip}=\sum_{i=1}^{n}y_i \\ \displaystyle b_0\sum_{i=1}^{n}x_{i1}+b_1\sum_{i=1}^{n}x_{i1}^2+b_2\sum_{i=1}^{n}x_{i1}\sum_{i=1}^{n}x_{i2}+\cdots+b_p\sum_{i=1}^{n}x_{i1}\sum_{i=1}^{n}x_{ip}=\sum_{i=1}^{n}x_{i1}y_i \\ \cdots \\ \displaystyle b_0\sum_{i=1}^{n}x_{ip}+b_1\sum_{i=1}^{n}x_{i1}\sum_{i=1}^{n}x_{ip}+b_2\sum_{i=1}^{n}x_{ip}\sum_{i=1}^{n}x_{i2}+\cdots+b_p\sum_{i=1}^{n}x_{ip}^2=\sum_{i=1}^{n}x_{ip}y_i \end{cases}
⎩⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎧b0n+b1i=1∑nxi1+b2i=1∑nxi2+⋯+bpi=1∑nxip=i=1∑nyib0i=1∑nxi1+b1i=1∑nxi12+b2i=1∑nxi1i=1∑nxi2+⋯+bpi=1∑nxi1i=1∑nxip=i=1∑nxi1yi⋯b0i=1∑nxip+b1i=1∑nxi1i=1∑nxip+b2i=1∑nxipi=1∑nxi2+⋯+bpi=1∑nxip2=i=1∑nxipyi称为正规方程组,引入矩阵
X
=
(
1
x
11
x
12
⋯
x
1
p
1
x
21
x
22
⋯
x
2
p
⋮
⋮
⋱
⋮
1
x
n
1
x
n
2
⋯
x
n
p
)
,
Y
=
(
y
1
y
2
⋮
y
p
)
,
B
=
(
b
0
b
1
⋮
b
p
)
,
\bf{X}=\begin{pmatrix} 1&x_{11}&x_{12}&\cdots&x_{1p} \\ 1&x_{21}&x_{22}&\cdots&x_{2p} \\ \vdots&\vdots&\ddots&\vdots \\ 1&x_{n1}&x_{n2}&\cdots&x_{np} \\ \end{pmatrix}, \bf{Y}=\begin{pmatrix} y_1 \\ y_2 \\ \vdots \\ y_p \\ \end{pmatrix}, \bf{B}=\begin{pmatrix} b_0 \\ b_1 \\ \vdots \\ b_p \\ \end{pmatrix},
X=⎝⎜⎜⎜⎛11⋮1x11x21⋮xn1x12x22⋱xn2⋯⋯⋮⋯x1px2pxnp⎠⎟⎟⎟⎞,Y=⎝⎜⎜⎜⎛y1y2⋮yp⎠⎟⎟⎟⎞,B=⎝⎜⎜⎜⎛b0b1⋮bp⎠⎟⎟⎟⎞,因
X
T
X
=
(
n
∑
i
=
1
n
x
i
1
⋯
∑
i
=
1
n
x
i
p
∑
i
=
1
n
x
i
1
∑
i
=
1
n
x
i
1
2
⋯
∑
i
=
1
n
x
i
1
x
i
p
⋮
⋮
⋱
⋮
∑
i
=
1
n
x
i
p
∑
i
=
1
n
x
i
1
x
i
p
⋯
∑
i
=
1
n
x
i
p
2
)
,
X
T
Y
=
(
∑
i
=
1
n
y
i
⋮
∑
i
=
1
n
x
i
1
y
i
∑
i
=
1
n
x
i
p
y
i
)
\bf{X^TX}=\begin{pmatrix} n&\displaystyle \sum_{i=1}^{n}x_{i1}&\cdots&\displaystyle \sum_{i=1}^{n}x_{ip} \\ \displaystyle \sum_{i=1}^{n}x_{i1} &\displaystyle \sum_{i=1}^{n}x_{i1}^2&\cdots&\displaystyle \sum_{i=1}^{n}x_{i1}x_{ip} \\ \vdots&\vdots&\ddots&\vdots \\ \displaystyle \sum_{i=1}^{n}x_{ip} &\displaystyle \sum_{i=1}^{n}x_{i1}x_{ip}&\cdots&\displaystyle \sum_{i=1}^{n}x_{ip}^2 \\ \end{pmatrix}, \bf{X^TY}=\begin{pmatrix} \displaystyle\sum_{i=1}^{n}y_i \\ \vdots \\ \displaystyle\sum_{i=1}^{n}x_{i1}y_i \\ \displaystyle\sum_{i=1}^{n}x_{ip}y_i \end{pmatrix}
XTX=⎝⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎛ni=1∑nxi1⋮i=1∑nxipi=1∑nxi1i=1∑nxi12⋮i=1∑nxi1xip⋯⋯⋱⋯i=1∑nxipi=1∑nxi1xip⋮i=1∑nxip2⎠⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎞,XTY=⎝⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎛i=1∑nyi⋮i=1∑nxi1yii=1∑nxipyi⎠⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎞上式化为矩阵方程
X
T
X
B
=
X
T
Y
\bf{X^TXB=X^TY}
XTXB=XTY可求得其解为
B
^
=
(
b
^
0
b
^
1
⋮
b
^
p
)
=
(
X
T
X
)
−
1
X
T
Y
\bf{\hat B}=\begin{pmatrix} \hat b_0 \\ \hat b_1 \\ \vdots \\ \hat b_p \\ \end{pmatrix}=\bf{(X^TX)^{-1}X^TY}
B^=⎝⎜⎜⎜⎛b^0b^1⋮b^p⎠⎟⎟⎟⎞=(XTX)−1XTY 于是我们求得经验回归方程为
y
^
=
b
^
0
+
b
^
1
x
1
+
b
^
2
x
2
=
⋯
+
b
^
p
x
p
\hat y=\hat b_0+\hat b_1x_1+\hat b_2x_2=\cdots+\hat b_px_p
y^=b^0+b^1x1+b^2x2=⋯+b^pxp