1. 总体、个体和样本
在使用数理统计方法研究某个实际问题时,往往把所研究对象的全体称为总体,而把组成总体的每个基本单位称为个体。
例如:某一批灯泡组成一个总体,而每个灯泡都称为个体;又如某个学校的所有学生被称为总体,每个学生都称为个体。
在数理统计中,总体或者是个体本身并不是我们所关心的,我们关心的是它们的某个特性或者说某个指标,如灯泡的寿命、学生的身高等,寿命和身高都是随机变量。
当我们需要对总体的某个特征进行研究时,往往不能对所有个体进行试验,一种可行的方法就是对总体进行抽样调查。数理统计的基本任务就是通过抽样结果来推断总体的统计规律。
对总体,把在相同条件下对随即变量X进行的n次重复独立观察,成为n次简单随机抽样,简称简单抽样。用Xi表示第i次观察结果,那么对总体X的n次重复独立观察所得结果记为X1,X2,…,Xn。每个Xi都是一个与总体X有相同分布的随机变量,且X1,X2,…,Xn相互独立,称随机变量X1,X2,…,Xn为来自总体X的简单随机样本,简称简单样本。
抽样观察的结果是n个具体数值x1,x2,…,xn,称其为简单样本X1,X2,…,Xn的一组观察值,成为样本观察值。
一般来说,不同的抽样得到样本X1,X2,…,Xn的不同观察值,样本所有可能取值的全体称为样本空间,用Ω表示,样本观察值实际就是样本空间中的一个点,称其为样本点。
例如:把一枚硬币重复独立的掷三次,正面记为1,反面记为0。用X1,X2,X3表示三次可能的结果,可知样本空间Ω=
{(0,0,0),(0,0,1),(0,1,0),(0,1,1),(1,0,0),(1,0,1),(1,1,0),(1,1,1)}。样本观察值可能为(0,1,1),它是样本空间中的一个样本点。
此外,样本的数量称为样本容量。
2. 统计量与充分统计量
统计量
我们知道,样本是总体的代表和反映,是对总体进行统计分析和推断的依据。因此,为了提取出样本中有关总体的信息,我们需要引入统计量的概念。
统计量是样本的函数,它的作用是汇集样本中有关总体的信息。
值得特别注意的是:统计量仅为样本 X 1 , X 2 , . . . X n X_1, X_2, ...X_n X1,X2,...Xn的函数,它不包含总体的任何未知量。这一点非常重要,是理解、求解和利用统计量的重点。
下面给出统计量的定义:
定义1: 设
X
1
,
X
2
,
.
.
.
X
n
X_1, X_2, ... X_n
X1,X2,...Xn为来自总体X的简单样本,若样本的函数
T
(
X
1
,
X
2
,
.
.
.
X
n
T(X_1, X_2, ... X_n
T(X1,X2,...Xn中不包含任何未知参数,则称此函数为统计量。
例如, X ˉ = 1 n ∑ i = 1 n X i \displaystyle\bar{X}=\frac{1}{n}\sum_{i=1}^{n}X_i Xˉ=n1i=1∑nXi仅为样本 X 1 , X 2 , . . . X n X_1, X_2, ... X_n X1,X2,...Xn的函数,因此它是统计量;而 1 n ∑ i = 1 n ( X i − μ ) 2 \displaystyle\frac{1}{n}\sum_{i=1}^{n}(X_i-\mu)^2 n1i=1∑n(Xi−μ)2中包含未知参数 μ \mu μ,因此它不是统计量。
常用统计量
下面介绍几个常用的统计量。设 X 1 , X 2 , . . . X n X_1, X_2, ... X_n X1,X2,...Xn为来自总体X的简单样本。
(1) 统计量 X ˉ = 1 n ∑ i = 1 n X i \displaystyle\bar{X}=\frac{1}{n}\sum_{i=1}^{n}X_i Xˉ=n1i=1∑nXi称为样本均值,它反映了样本的平均取值,描述了样本的集中取值趋势。
(2) 统计量 S 2 = 1 n − 1 ∑ i = 1 n ( X i − X ˉ ) 2 \displaystyle S^2=\frac{1}{n-1}\sum_{i=1}^{n}{(X_i-\bar{X})^2} S2=n−11i=1∑n(Xi−Xˉ)2称为为样本方差,它描述了样本取值的分散程度大小,称 S S S为样本标准差。
(3) 统计量 A k = 1 n ∑ i = 1 n X i k \displaystyle A_k=\frac{1}{n}\sum_{i=1}^{n}{X_i^k} Ak=n1i=1∑nXik称为样本k阶原点矩。
(4) 统计量 B k = 1 n ∑ i = 1 n ( X i − X ˉ ) k \displaystyle B_k=\frac{1}{n}\sum_{i=1}^{n}{(X_i-\bar{X})^k} Bk=n1i=1∑n(Xi−Xˉ)k称为样本k阶中心矩。
顺序统计量
除了样本矩之外,还有一类十分重要的统计量,称为样本顺序统计量,它在数理统计中有着广泛的应用。下面给出样本顺序统计量的定义:
定义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)}
x(1),x(2),...x(n)。其中,
x
(
1
)
,
x
(
2
)
,
.
.
.
x
(
n
)
x_{(1)}, x_{(2)}, ... x_{(n)}
x(1),x(2),...x(n)满足
x
(
1
)
≤
x
(
2
)
≤
.
.
.
≤
x
(
n
)
x_{(1)}≤x_{(2)}≤ ... ≤x_{(n)}
x(1)≤x(2)≤...≤x(n)
定义排在
k
k
k位置的数
x
(
k
)
x{(k)}
x(k)为随机变量
X
(
k
)
X_{(k)}
X(k)的观察值,显然:
X
(
1
)
≤
X
(
2
)
≤
.
.
.
≤
X
(
n
)
X_{(1)}≤X_{(2)}≤ ... ≤X_{(n)}
X(1)≤X(2)≤...≤X(n)
由于
X
(
k
)
X_{(k)}
X(k)是样本
X
1
,
X
2
,
.
.
.
X
n
X_1, X_2, ... X_n
X1,X2,...Xn的函数,且不含任何未知数,因此
X
(
k
)
X_{(k)}
X(k)是统计量。称
X
(
1
)
,
X
(
2
)
,
.
.
.
X
(
n
)
X_{(1)}, X_{(2)}, ... X_{(n)}
X(1),X(2),...X(n)为顺序统计量。
特别的, X ( 1 ) X_{(1)} X(1)=min { 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 ( n ) X_{(n)} X(n)=max { 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(1)和 X ( n ) X_{(n)} X(n)的分布函数。
设
X
1
,
X
2
,
.
.
.
X
n
X_1, X_2, ... X_n
X1,X2,...Xn为总体X的简单样本,总体X的分布函数为F(x),根据分布函数的定义:
F
(
X
(
1
)
<
x
)
=
1
−
F
(
X
(
1
)
≥
x
)
=
1
−
F
(
X
(
n
)
≥
.
.
.
≥
X
(
2
)
≥
X
(
1
)
≥
x
)
=
1
−
F
(
X
(
n
)
≥
x
)
.
.
.
F
(
X
(
2
)
≥
x
)
⋅
F
(
X
(
1
)
≥
x
)
=
1
−
[
1
−
F
(
X
(
n
)
<
x
)
]
.
.
.
[
1
−
F
(
X
(
2
)
<
x
)
]
⋅
[
1
−
F
(
X
(
1
)
<
x
)
]
=
1
−
[
1
−
F
(
x
)
]
n
\begin{aligned} F(X_{(1)}<x) &=1-F(X_{(1)}≥x)\\ &=1-F(X_{(n)}≥...≥X_{(2)}≥X_{(1)}≥x) \\&=1-F(X_{(n)}≥x)...F(X_{(2)}≥x)·F(X_{(1)}≥x) \\&=1-[1-F(X_{(n)}<x)]...[1-F(X_{(2)}<x)]·[1-F(X_{(1)}<x)] \\&=1-[1-F(x)]^n \end{aligned}
F(X(1)<x)=1−F(X(1)≥x)=1−F(X(n)≥...≥X(2)≥X(1)≥x)=1−F(X(n)≥x)...F(X(2)≥x)⋅F(X(1)≥x)=1−[1−F(X(n)<x)]...[1−F(X(2)<x)]⋅[1−F(X(1)<x)]=1−[1−F(x)]n
记作:
F
X
(
n
)
(
x
)
=
1
−
[
1
−
F
(
x
)
]
n
F_{X_{(n)}}(x)=1-[1-F(x)]^n
FX(n)(x)=1−[1−F(x)]n
同理可求得:
F
X
(
1
)
(
x
)
=
F
n
(
x
)
F_{X_{(1)}}(x)=F^n(x)
FX(1)(x)=Fn(x)
基于顺序统计量,我们可以构造一些常用的统计量,如样本极差,样本分位数等。
(1) 样本极差:样本极差定义为R=
X
(
n
)
−
X
(
1
)
X_{(n)}-X_{(1)}
X(n)−X(1),同样本方差类似,它也可以用来描述样本的分散程度。
(2) 样本分位数:样本的p(0<p<1)分位数
m
p
m_p
mp定义为:
m
p
=
{
X
(
[
n
p
+
1
]
)
n
p
不
为
整
数
1
2
(
X
(
n
p
)
+
X
(
n
p
+
1
)
)
n
p
为
整
数
m_p=\begin{cases} X_{([np+1])} & np不为整数 \\\frac{1}{2}(X_{(np)}+X_{(np+1)}) & np为整数 \end{cases}
mp={X([np+1])21(X(np)+X(np+1))np不为整数np为整数
其中, [x]表示取不超过数x的最大整数。
注1:若将样本观察值从小到大进行排序,p分位数
m
p
m_p
mp的意义是排在它前的观察值个数约占观察值总数的p%。
注2:一般情况下,要给出
m
p
m_p
mp的精确分布是很困难的,但当样本容量n很大时,可以证明
m
p
m_p
mp为渐近正态。
注3:样本极差和样本分位数都是随机变量,理解这一点非常重要。
充分统计量
下面将介绍充分统计量的概念,为了便于理解,先看下例:
例1:掷硬币试验,设正面出现的概率为 θ \theta θ,总体X的分布族为 { P θ : θ ∈ ( 0 , 1 ) } \{P_{\theta}:\theta∈(0, 1)\} {Pθ:θ∈(0,1)},其中 P θ { X = x } = θ x ( 1 − θ ) 1 − x , x i = 0 , 1 P_{\theta}\{X=x\}=\theta^x(1-\theta)^{1-x},x_i=0,1 Pθ{X=x}=θx(1−θ)1−x,xi=0,1。设 X 1 , X 2 , . . . X n X_1, X_2, ... X_n X1,X2,...Xn是来自总体X的简单样本,样本空间 S = { ( x 1 , x 2 , . . . , x n ) : x i = 0 , 1 , i = 1 , 2 , . . . , n } S=\{(x_1, x_2, ..., x_n):x_i=0,1,i=1, 2, ..., n\} S={(x1,x2,...,xn):xi=0,1,i=1,2,...,n},它包含了 2 n 2^n 2n个样本点。
对于上述问题,总体X服从两点分布,因此 θ \theta θ实际上反映的是正面出现的概率。当样本数量足够多时,我们也可以说 θ \theta θ决定了正面出现的次数,这可以由样本点 ( x 1 , x 2 , . . . , x n ) (x_1, x_2, ..., x_n) (x1,x2,...,xn)中1的个数来体现。同样的,样本点 ( x 1 , x 2 , . . . , x n ) (x_1, x_2, ..., x_n) (x1,x2,...,xn)中1的个数也很好的体现了 θ \theta θ的信息。
下面考虑统计量 T = ∑ i = 1 n X i \displaystyle T=\sum_{i=1}^nX_i T=i=1∑nXi,它是样本 X 1 , X 2 , . . . X n X_1, X_2, ... X_n X1,X2,...Xn的函数,提取了样本点中"1"的个数信息,剔除了样本带你中"1"的位置信息,起到了很好的信息压缩作用。同时,样本点中"1"的个数信息正好是参数 θ \theta θ的全部信息,因此T既起到了提取参数 θ \theta θ的全部信息的作用,也起到了剔除了无关信息的作用。
为了理解一个统计量包含了样本中有关某参数的全部信息这一个说法,可以通过下面的例子来说明:
例2:当我们用样本均值 X ˉ \bar{X} Xˉ估计总体均值 μ = E X \mu=EX μ=EX时,用 X ˉ \bar{X} Xˉ的观察值 x ˉ \bar{x} xˉ作为 μ \mu μ的估计值,这时统计量 X ˉ \bar{X} Xˉ已经包含了样本中有关参数 μ \mu μ的所有信息,当我们进一步知道 ( X 1 , X 2 , . . . X n ) (X_1,X_2, ...X_n) (X1,X2,...Xn)中任何一个变量的取值时并不能帮助我们获取参数 μ \mu μ的更多信息。换句话说,对于样本中关于参数 μ \mu μ的信息,我们不需要记录所有样本的取值,只需要记录 X ˉ \bar{X} Xˉ即可。
也就是说,当我们用样本来估计总体的某个参数时,不需要分析样本的所有信息,只需要找到一个关于该参数的充分统计量即可,因为该统计量包含了样本中关于该参数的全部信息。
下面举一个更加具体的例子:
设总体X服从正态分布 N ( μ , 1 ) , μ N(\mu, 1),\mu N(μ,1),μ是需要估计的未知参数, X 1 , X 2 X_1, X_2 X1,X2是来自总体的样本,用样本均值 X ˉ = 1 2 ( X 1 + X 2 ) \displaystyle \bar{X}=\frac{1}{2}(X_1+X_2) Xˉ=21(X1+X2)来估计 μ \mu μ。这时样本空间 S S S为平面区域,当给定 X ˉ \bar{X} Xˉ=a时, 1 2 ( X 1 + X 2 ) = a \displaystyle \frac{1}{2}(X_1+X_2)=a 21(X1+X2)=a表示样本空间 S S S中的一条直线:
( X 1 , X 2 ) (X_1, X_2) (X1,X2)服从二维正态分布 N ( μ , 1 ; μ , 1 ; 0 ) N(\mu, 1;\mu, 1;0) N(μ,1;μ,1;0),其密度函数为:
p ( x 1 , x 2 ; μ ) = 1 2 π e x p { − 1 2 ∑ i = 1 2 ( x i − μ ) 2 } p(x_1, x_2;\mu)=\frac{1}{2\pi}exp\{-\frac{1}{2}\sum_{i=1}^{2}(x_i-\mu)^2\} p(x1,x2;μ)=2π1exp{−21i=1∑2(xi−μ)2}
样本均值 X ˉ = 1 2 ( X 1 + X 2 ) \displaystyle \bar{X}=\frac{1}{2}(X_1+X_2) Xˉ=21(X1+X2)服从正态分布 N { μ , 1 2 } \displaystyle N\{\mu, \frac{1}{2}\} N{μ,21},其密度函数为:
p ( a ; μ ) = 1 π e x p { − ( a − μ ) 2 } p(a;\mu)=\frac{1}{\sqrt\pi}exp\{-(a-\mu)^2\} p(a;μ)=π1exp{−(a−μ)2}
在直线 X ˉ = a \bar{X}=a Xˉ=a上, ( X 1 , X 2 ) (X_1, X_2) (X1,X2)的分布可以用条件密度函数
p ( x 1 , x 2 ∣ x ˉ = a ) = { p ( x 1 , x 2 ; μ ) p ( a ; μ ) x ˉ = a 0 x ˉ ≠ a \displaystyle p(x_1, x_2|\bar{x}=a)=\begin{cases} \frac{p(x_1, x_2;\mu)}{p(a;\mu)} &\bar{x}=a \\0 & \bar{x}≠a \end{cases} p(x1,x2∣xˉ=a)={p(a;μ)p(x1,x2;μ)0xˉ=axˉ̸=a
来描述,其中 x ˉ \bar{x} xˉ为 X ˉ \bar{X} Xˉ的观察值,而a为 X ˉ \bar{X} Xˉ的给定值。所以
p ( x 1 , x 2 ∣ x ˉ = y ) = { 1 2 π e x p { − 1 2 ∑ i = 1 2 ( x i − a ) 2 x ˉ = a 0 x ˉ ≠ a \displaystyle p(x_1, x_2|\bar{x}=y)=\begin{cases} \displaystyle \frac{1}{2\sqrt\pi}exp\{-\frac{1}{2}\sum_{i=1}^{2}(x_i-a)^2 & \bar{x}=a \\0 & \bar{x}≠a \end{cases} p(x1,x2∣xˉ=y)=⎩⎪⎨⎪⎧2π1exp{−21i=1∑2(xi−a)20xˉ=axˉ̸=a
由上面的条件密度函数的表达式可知,它与参数 μ \mu μ无关,这表明在已知 x ˉ = y \bar{x}=y xˉ=y时,如果进一步知道样本的观察值 ( x 1 , x 2 ) (x_1, x_2) (x1,x2)在这条直线上的具体位置,也并没有增加关于参数 μ \mu μ的任何新信息。换句话说, x ˉ \bar{x} xˉ已充分提取了样本 X 1 , X 2 X_1, X_2 X1,X2关于参数 μ \mu μ的全部信息,因此 X ˉ \bar{X} Xˉ是参数 μ \mu μ的充分统计量。
根据上面的解释,根据概率论中的条件分布给出充分统计量的定义:
定义: 设总体分布族为 { P θ : θ ∈ Θ } , X 1 , X 2 , . . . , X n \{P_\theta:\theta\in\Theta\},X_1, X_2, ..., X_n {Pθ:θ∈Θ},X1,X2,...,Xn是来自总体的简单样本,样本的函数 T { X 1 , X 2 , . . . X n } T\{X_1, X_2, ...X_n\} T{X1,X2,...Xn}是统计量。如果在给定 T { X 1 , X 2 , . . . X n } T\{X_1, X_2, ...X_n\} T{X1,X2,...Xn}=t的条件下,样本X_1, X_2, …X_n的条件分布函数 F θ ( x 1 , x 2 , . . . x n ∣ t ) F_\theta(x_1, x_2, ... x_n|t) Fθ(x1,x2,...xn∣t)与参数 θ \theta θ无关,则称统计量 T { X 1 , X 2 , . . . X n } T\{X_1, X_2, ...X_n\} T{X1,X2,...Xn}为参数 θ \theta θ的充分统计量。
注: 由于条件分布函数与条件分布列或条件密度函数是一一对应的,因此可以用条件分布列或者条件密度函数来代替上面的条件分布函数。
例3:掷硬币试验,设正面出现的概率为 θ \theta θ,总体X的分布族为 { P θ : θ ∈ ( 0 , 1 ) } \{P_{\theta}:\theta∈(0, 1)\} {Pθ:θ∈(0,1)},其中 P θ { X = x } = θ x ( 1 − θ ) 1 − x , x i = 0 , 1 P_{\theta}\{X=x\}=\theta^x(1-\theta)^{1-x},x_i=0,1 Pθ{X=x}=θx(1−θ)1−x,xi=0,1。设 X 1 , X 2 , . . . X n X_1, X_2, ... X_n X1,X2,...Xn是来自总体X的简单样本,求 θ \theta θ的充分统计量
解:总体X的分布律为: { P θ : θ ∈ ( 0 , 1 ) } \{P_{\theta}:\theta∈(0, 1)\} {Pθ:θ∈(0,1)}
样本的联合分布列为:
P ( x ; θ ) = θ ∑ i = 1 n x i ( 1 − θ ) n − ∑ i − 1 n x i \displaystyle P(x;\theta )=\theta^{\sum_{i=1}^{n}x_i}(1-\theta)^{n-\sum_{i-1}^{n}x_i} P(x;θ)=θ∑i=1nxi(1−θ)n−∑i−1nxi
其中 x i = 0 或 1 x_i=0或1 xi=0或1.给定 T = ∑ i − 1 n X i = t T=\sum_{i-1}^{n}X_i=t T=∑i−1nXi=t的条件下,样本 X 1 , X 2 , . . . X n X_1, X_2, ... X_n X1,X2,...Xn的条件分布列为:
P ( x ; θ ) = P { X 1 = x 1 , X 2 = x 2 , . . . , X n = x n , ∑ i = 1 n X i = t } P { ∑ i = 1 n X i = t } = θ t ( 1 − θ ) n − t C n t θ t ( 1 − θ ) n − t = 1 C n t \begin{aligned} P(x;\theta )&=\frac{P\{X_1=x_1,X_2=x_2,..., X_n=x_n, \sum_{i=1}^{n}X_i=t\}}{P\{\sum_{i=1}^{n}X_i=t\}} \\&=\frac{\theta^t(1-\theta)^{n-t}}{C_n^t\theta^t(1-\theta)^{n-t}} \\&=\frac{1}{C_n^t} \end{aligned} P(x;θ)=P{∑i=1nXi=t}P{X1=x1,X2=x2,...,Xn=xn,∑i=1nXi=t}=Cntθt(1−θ)n−tθt(1−θ)n−t=Cnt1
因此,条件分布列与参数 θ \theta θ无关,由定义,统计量 T = ∑ i − 1 n X i T=\sum_{i-1}^{n}X_i T=∑i−1nXi是充分统计量。
一般情况下,利用条件分布直接证明给定的统计量是否是充分统计量是比较困难的。下面给出一个定理来求解样本的充分统计量。
定理(Fisher-Neyman准则):设总体分布族为
{
P
θ
:
θ
∈
Θ
}
\{ P_{\theta}:\theta\in\Theta\}
{Pθ:θ∈Θ},则统计量
T
(
x
)
T(x)
T(x)是充分统计量的,当且仅当存在一个定义在
I
×
Θ
I\times\Theta
I×Θ上的实值函数
g
(
t
,
θ
)
g(t, \theta)
g(t,θ)及定义在样本空间
S
S
S上的不依赖于参数
θ
\theta
θ的实值函数
h
(
x
)
h(x)
h(x)使样本
x
1
,
x
2
,
.
.
.
x
n
x_1, x_2, ...x_n
x1,x2,...xn的联合分布
p
(
x
;
θ
)
p(x;\theta)
p(x;θ)的分解式
p
(
x
;
θ
)
=
g
(
T
(
x
)
,
θ
)
h
(
x
)
p(x;\theta)=g(T(x), \theta)h(x)
p(x;θ)=g(T(x),θ)h(x)
对所有的
x
∈
S
x\in S
x∈S都成立,其中
I
I
I是统计量
T
(
x
)
T(x)
T(x)的值域,
p
(
x
,
θ
)
p(x, \theta)
p(x,θ)是样本的联合分布列或密度函数。
例4. 设 x 1 , x 2 , . . . x n x_1, x_2, ...x_n x1,x2,...xn是来自Possion分布总体 P ( x ; λ ) P(x;\lambda) P(x;λ)的简单样本,求参数 λ \lambda λ的充分统计量。
解:泊松分布的分布律为:
P { X = k } = λ k k ! e − λ , k = 0 , 1 , . . . P\{X=k\}=\frac{\lambda^k}{k!}e^{-\lambda},k=0, 1, ... P{X=k}=k!λke−λ,k=0,1,...
则样本的联合分布列为:
p ( x 1 , x 2 , . . . x n ; λ ) = e − n λ x 1 ! x 2 ! . . . x n ! λ ∑ i = 1 n x i = ( e − n λ λ ∑ i = 1 n x i ) e − n λ x 1 ! x 2 ! . . . x n ! \begin{aligned} p(x_1, x_2, ...x_n;\lambda)&=\frac{e^{-n\lambda}}{x_1!x_2!...x_n!}\lambda^{\sum^n_{i=1}x_i} \\&=(e^{-n\lambda}\lambda^{\sum^n_{i=1}x_i})\frac{e^{-n\lambda}}{x_1!x_2!...x_n!} \end{aligned} p(x1,x2,...xn;λ)=x1!x2!...xn!e−nλλ∑i=1nxi=(e−nλλ∑i=1nxi)x1!x2!...xn!e−nλ
取 T ( x ) = ∑ i = 1 n x i = t , g ( t , λ ) = e − n λ λ t , h ( x ) = e − n λ x 1 ! x 2 ! . . . x n ! T(x)=\sum_{i=1}^nx_i=t,g(t, \lambda)=e^{-n\lambda}\lambda^t,h(x)=\frac{e^{-n\lambda}}{x_1!x_2!...x_n!} T(x)=∑i=1nxi=t,g(t,λ)=e−nλλt,h(x)=x1!x2!...xn!e−nλ
由定理得 T ( x ) = ∑ i = 1 n x i T(x)=\sum_{i=1}^nx_i T(x)=∑i=1nxi为参数 λ \lambda λ的充分统计量。
下面是一个多参数的例子:
例5. 设 x 1 , x 2 , . . . x 3 x_1, x_2, ...x_3 x1,x2,...x3是来自正态总体 N ( μ , σ 2 ) N(\mu, \sigma^2) N(μ,σ2)的简单样本,其中 μ 和 σ 2 \mu和\sigma^2 μ和σ2都是未知的,试求参数 θ = ( μ , σ 2 ) \theta=(\mu, \sigma^2) θ=(μ,σ2)的完全统计量。
解:样本的联合密度函数为:
p ( x 1 , x 2 . . . x n ; θ ) = ( 1 2 π σ 2 ) n e x p { − 1 2 σ 2 ∑ i = 1 n ( x i − μ ) 2 } = ( 1 2 π σ 2 ) n e x p { − 1 2 σ 2 ( ∑ i = 1 n x i 2 − 2 μ ∑ i = 1 n x i + n μ 2 ) } \begin{aligned} p(x_1, x_2...x_n;\theta)&=(\frac{1}{\sqrt{2\pi\sigma^2}})^nexp\{-\frac{1}{2\sigma^2}\sum_{i=1}^n(x_i-\mu)^2\} \\&=(\frac{1}{\sqrt{2\pi\sigma^2}})^nexp\{-\frac{1}{2\sigma^2}(\sum_{i=1}^nx_i^2-2\mu\sum_{i=1}^nx_i+n\mu^2)\} \end{aligned} p(x1,x2...xn;θ)=(2πσ21)nexp{−2σ21i=1∑n(xi−μ)2}=(2πσ21)nexp{−2σ21(i=1∑nxi2−2μi=1∑nxi+nμ2)}
取 T ( x ) = ( T 1 ( x ) , T 2 ( x ) ) = ( ∑ i = 1 n x i , ∑ i = 1 n x i 2 ) , h ( x ) ≡ 1 T(x)=(T_1(x), T_2(x))=(\sum_{i=1}^nx_i, \sum_{i=1}^nx_i^2),h(x)\equiv1 T(x)=(T1(x),T2(x))=(∑i=1nxi,∑i=1nxi2),h(x)≡1及
g ( t , θ ) = ( 1 2 π σ 2 ) 2 e x p { − 1 2 σ 2 ) ( t 2 − 2 μ t 1 + n μ 2 ) } g(t, \theta)=(\frac{1}{\sqrt{2\pi\sigma^2}})^2exp\{-\frac{1}{2\sigma^2})(t_2-2\mu t_1+n\mu^2)\} g(t,θ)=(2πσ21)2exp{−2σ21)(t2−2μt1+nμ2)}
其中, t = ( t 1 , t 2 ) t=(t1, t2) t=(t1,t2),则有
p ( x 1 , x 2 , . . . x n ; θ ) = g ( t , θ ) h ( x ) p(x_1, x_2, ...x_n;\theta)=g(t, \theta)h(x) p(x1,x2,...xn;θ)=g(t,θ)h(x)
由定理得, T ( x ) = ( ∑ i = 1 n x i , ∑ i = 1 n x i 2 ) T(x)=(\sum_{i=1}^nx_i, \sum_{i=1}^nx_i^2) T(x)=(∑i=1nxi,∑i=1nxi2)是 θ \theta θ的充分统计量。
又因为样本的联合密度函数可表示为:
p ( x 1 , x 2 . . . x n ; θ ) = ( 1 2 π σ 2 ) n e x p { − 1 2 σ 2 ∑ i = 1 n ( x i + X ˉ − X ˉ − μ ) 2 } = ( 1 2 π σ 2 ) n e x p { − 1 2 σ 2 [ ∑ i = 1 n ( x i − x ˉ ) 2 + n ( x ˉ − μ ) 2 ] } \begin{aligned} p(x_1, x_2...x_n;\theta)&=(\frac{1}{\sqrt{2\pi\sigma^2}})^nexp\{-\frac{1}{2\sigma^2}\sum_{i=1}^n(x_i+\bar{X}-\bar{X}-\mu)^2\} \\&=(\frac{1}{\sqrt{2\pi\sigma^2}})^nexp\{-\frac{1}{2\sigma^2}[\sum_{i=1}^n(x_i-\bar{x})^2+n(\bar{x}-\mu)^2]\} \end{aligned} p(x1,x2...xn;θ)=(2πσ21)nexp{−2σ21i=1∑n(xi+Xˉ−Xˉ−μ)2}=(2πσ21)nexp{−2σ21[i=1∑n(xi−xˉ)2+n(xˉ−μ)2]}
可取 S ( x ) = ( S 1 ( x ) , S 2 ( x ) ) = ( ∑ i − 1 n x i , ∑ i = 1 n ( x i − x ˉ ) , h ( x ) ≡ 1 S(x)=(S_1(x), S_2(x))=(\sum_{i-1}^{n}x_i, \sum_{i=1}^n(x_i-\bar{x}), h(x)\equiv1 S(x)=(S1(x),S2(x))=(∑i−1nxi,∑i=1n(xi−xˉ),h(x)≡1及
g ( s , θ ) = ( 1 2 π σ 2 ) n e x p { − 1 2 σ 2 [ s 2 + n ( s 1 n − μ ) 2 ] } g(s, \theta)=(\frac{1}{\sqrt{2\pi\sigma^2}})^nexp\{-\frac{1}{2\sigma^2}[s_2+n(\frac{s_1}{n}-\mu)^2]\} g(s,θ)=(2πσ21)nexp{−2σ21[s2+n(ns1−μ)2]}
由定理得, S ( x ) = ( ∑ i − 1 n x i , ∑ i = 1 n ( x i − x ˉ ) , h ( x ) ≡ 1 S(x)=(\sum_{i-1}^{n}x_i, \sum_{i=1}^n(x_i-\bar{x}), h(x)\equiv1 S(x)=(∑i−1nxi,∑i=1n(xi−xˉ),h(x)≡1是 θ \theta θ的充分统计量。
上例说明,参数 θ \theta θ的充分统计量是不唯一的。另一方面,当参数的维数和充分统计量的维数相等时,应当将充分统计量看作是一个有机的整体,不能拆分,认为充分统计量的分量是相应参数向量的分量的充分统计量是错误的。
下面再举一个例子,通过下面这个例子我们可以看到,顺序统计量有时候也是充分统计量。
例6. 设 x 1 , x 2 , . . . x n x_1, x_2, ...x_n x1,x2,...xn是来自均匀分布总体U[a, b]的简单样本,其中a和b都是未知的,且a<b,令参数 θ = ( a , b ) \theta=(a, b) θ=(a,b),求 θ \theta θ的充分统计量。
解:样本的联合密度函数为:
p ( x 1 , x 2 , . . . x n ) = 1 ( b − a ) n I { a ≤ x ( 1 ) ≤ x ( n ) ≤ b } ( x ( 1 ) , x ( 2 ) , . . . x ( n ) ) = 1 ( b − a ) n I { a ≤ x ( 1 ) ≤ x ( n ) ≤ b } ( x ( 1 ) , x ( n ) ) \begin{aligned} p(x_1, x_2, ...x_n)&=\frac{1}{(b-a)^n}I_{\{a\le x_{(1)}\le x_{(n)}\le b\}}(x_{(1)}, x_{(2)}, ...x_{(n)}) \\&=\frac{1}{(b-a)^n}I_{\{a\le x_{(1)}\le x_{(n)}\le b\}}(x_{(1)}, x_{(n)}) \end{aligned} p(x1,x2,...xn)=(b−a)n1I{a≤x(1)≤x(n)≤b}(x(1),x(2),...x(n))=(b−a)n1I{a≤x(1)≤x(n)≤b}(x(1),x(n))
其中
I A ( u ) = { 1 , u ∈ A 0 , u ∉ A I_A(u)=\begin{cases}1,&u\in A \\ 0, &u\notin A \end{cases} IA(u)={1,0,u∈Au∈/A
称其为示性函数。
令 T ( x ) = ( x ( 1 ) , x ( n ) ) , g ( t , θ ) = 1 ( b − a ) n I { a ≤ t 1 ≤ t 2 ≤ b } ( t 1 , t 2 ) , h ( x ) ≡ 1 T(x)=(x_{(1)},x_{(n)}),g(t, \theta)=\frac{1}{(b-a)^n}I_{\{a\le t_1\le t_2\le b\}}(t_1, t_2),h(x)\equiv1 T(x)=(x(1),x(n)),g(t,θ)=(b−a)n1I{a≤t1≤t2≤b}(t1,t2),h(x)≡1
其中, t = ( t 1 , t 2 ) t=(t_1, t_2) t=(t1,t2),则联合密度函数有分解式
p ( x 1 , x 2 , . . . x n ; θ ) = g ( t , θ ) h ( x ) p(x_1, x_2, ...x_n;\theta)=g(t, \theta)h(x) p(x1,x2,...xn;θ)=g(t,θ)h(x)
由定理知, T ( x ) = ( x ( 1 ) , x ( n ) ) T(x)=(x_{(1)},x_{(n)}) T(x)=(x(1),x(n))是参数 θ \theta θ的充分统计量。
最后需要提及的是,对一个分布族 { P θ : θ ∈ Θ } \{P_\theta:\theta\in\Theta\} {Pθ:θ∈Θ}而言,充分统计量往往不是唯一的,都包含样本中关于参数的全部信息,但它们在压缩数据能力方面是不同的。例如:对于两点分布族,统计量 ( x 1 , x 2 , . . . x n ) (x_1, x_2, ...x_n) (x1,x2,...xn)和统计量 ∑ i = 1 n x i \sum_{i=1}^nx_i ∑i=1nxi都是充分统计量,都保留了与参数 θ \theta θ有关的“1”的个数的信息,但是在剔除无关信息方面, ∑ i = 1 n x i \sum_{i=1}^nx_i ∑i=1nxi好得多。因此,我们往往需要求的是最小充分统计量。