第七章 参数估计
关键词:
矩法估计
极大似然估计
置信区间
置信水平(置信度)
枢轴量
参数:反映总体某方面特征的量
I 90
90 1 (90 )
X X
p P X
2
设浙江大学大一学生某学年的《微积分》成绩 服从正态分布,当时为优秀,则优秀率
也是一个参数,它是和 例:
的函数。
当总体的参数未知时,需利用样本资料对其 给出估计——参数估计。
3
两类参数估计方法:点估计和区间估计
1 n
X
X
X X
设 设 设 设 设 设 设 设 设 , ,设 设 设 设 设 设 设 设 设
4
1
1
1
ˆˆ ( ˆ
ˆ(
n
n n
X X
x x
x x
设设设设设设设设设设设设设设, , )
设设设设设设设设设设设设设设
设设设设设 设设设设设设设设, ,设设
设, , )设设设设
设设设设设常用的点估计方法:
矩法、极大似然法
7.1
参数的点估计:以样本矩估计总体矩,以样本矩的 函数估计总体矩
统计思想
的函数。
5
:辛钦大数定律和依概率收
理论根据 敛的性质
1 1
( ;
, ,
k k
k
k
F x
k
1 1
设总体的分布函数为, , ),其中, ,是 待估的未知参数,假定总体的前阶原点矩
或中心矩, ,存在。
(一)矩估计法
(
i) ( , , ),
11, ,
i
E X h
i ki
k
k k
设1设设设设设设设设设设设设设设设
6
( , ,
1), 1, ,
i
g
i ki k
k
(2)设设设设设设设设设设设设1
ˆ 1 1,
, ,
, k
i i k
A A
kg A A i k
1
(3)以样本各阶矩, ,代替总体各阶矩,
得各参数的矩估
, 计
( , ),
基本步骤
.
i i i
i
B
注:在实际应用时,为求解方便,也可以用 中心矩代替原点矩,相应地以样本中心矩
估计
7
1
( , ) , ,
nX U a b a b
X X a b
设设设设设设设设设设设设 设设设设设设设设设设设
设1设
设设设设
8
2
1 2
( )
( ) ,
2 12
a b a b
E X 解(1)求矩关于参数的函数
9
1 3 ,2 1 3 2
a b
(2) 求参数关于矩的反函数
2
1 2
1
2 2
, 1 (
ˆ
)
ˆ 3 3
n i i
A X B X
a X B b X B
n X a b
1
2
=
(3) 以样本阶矩代替总体矩 代替,得参数和
, =
的矩估计
10
1
1 2
2
0 1
0 0
, , n
X
x x
f x
X X X X
例:设总体的密度为:
为未知参数,
其他
,为取自的样本,求的矩估计。
11
1 1
E X
xf x dx
解:()
1
1
0
x dx
23 ˆ 1 X
X
()
2 1
1
2 1
()
2
2
~ ( , ,
X N
:某工厂生产的零件长度),
其中未知,
测定当长度落在[46, 50]时,产品合格,
并以参数代表该厂生产零件的合格率.
从中随机抽取10个,测得长度为
46, 51, 48, 47, 50, 44, 48, 49, 50, 47, 使用矩法给出合格率的
例3
估计值.
50
46
解:由正态分布的性质知13
2 2
50 46
ˆ X X
B B
因此的矩估计为
48,
24
ˆx B
由样本资料计算得代入
=2(1)- 1
得
=0. 6826
极(最)大似然估计的原理介绍
考察以下例子:
假设在一个罐中放着许多白球和黑球,并假定 已经知道两种球的数目之比是 1:3 ,但不知道哪种 颜色的球多。如果用放回抽样方法从罐中取 5 个球,
观察结果为:黑、白、黑、黑、黑,估计取到黑球 的概率 p.
二 大似然估计法:极1 3
, .
4 4
p p
解:设抽到黑球的概率为则本例中,或
15
3 31 1 .
4 4 4 1024
p 4
当时,出现本次观察结果的概率为
3 3 1 81 .
4 4 4 1024
p 4
当时,出现本次观察结果的概率为
3 81 3 1
1024 1024 4 4
ˆ 3 .
4
p p
p
由于,因此认为比更有可能,
于是取为更合理
16
~ ( ; )
X p x
一般地,设离散型总体,,未知。1, , n 1, , n
X X X x x
从总体中取得样本,其观察值为,
1 1
1 1 1
1
, ,
( ) , , ( ; )... ( ; ) ( ; ).
n n
n
n n n i
i
X x X x
L P X x X x p x p x p x
设设设设设设设设设
ˆ 1
( ( , , )) max ( ).n
L x x L
设设设设设设设
似然 函数
1
1 ML
ˆ , ,
ˆ , , E)
n
n
x x
X X
设设设设设设
设设设设设
设
设设设设设设设设 设设
设设设设 设设
设设
设
,
X
f x
若总体为连续型的,概率密度为,,为未知参数。
1 2 1 2
1
, , , , , ,
,
n n
n
i i
X X X x x x
L
f x
设设设设设设设设设设 设设设设 设
ˆ 1
( ( , , )) max ( ).n
L x x L
设设设设设设设
18
1 2
1. , , ,k
设设 设设设设设设设设设设设设设设设 设
2.
0, 1, 2,..., . ˆi 1, 2,..., .
i
L lnL
lnL
lnL i k i k
在求的最大值时,通常转换为求:的最大值,
称为对数似然函数.
利用解得,
3. i
i
L
若关于某个是单调增减函数,
此时的极大似然估计在其边界取得;
4.若是的极大似然估计,则的极大似然估计为。 ˆˆ g g
11
0 1 0
, , n
X
x x
f x
X X X
设设设设设设设设设设
设设 设
设4设
设设设设设设设设设设设设设设设设
20
2
2
1
ˆ
n i i
n lnX
的极大似然估计量为:
2
1
1 1
1 1
,
n
n
i i
n n
i i
i i
L f x
x x
解:似然函数
1
1 ln 2
n i i
lnL n ln x
1
1 1 ln 0 2 2
n i i
dlnL n x
d
令
1
ln
n i i
n x
即:
1
1
, 0
0, , , , ,
,
x
n
X
e x
f x
X X X
设 设 设 设 设 设 设 设 设 设
设 设 设 设 设 设 设 设 设 设
设 设 设 设 设 设 设 设 设 设 设 设 设 设 设
设
设 设
5设
22
1 解:矩估计
1 E X xf x dx
2
1 2
v
v
得
22 ( )
v
D X
E X
2 1
2 1
ˆ 1 ( )
ˆ 1 ( )
n i i
n i i
X X n
X X X
n
1 x
x e dx
2
1
(
x)
e x dx
2 2
0
( ) 1
t x
t
e
t dt
23
2 极大似然估计
1
, n 1 xi
i
L
e
此处不能通过求偏导数获得的极大似然估计量,
(1) 1 2
,
, , ,
i
n
x
x min x x x
设设设设设设设设设设设
1
1 1 , 1, 2,..., .
n i i
x
n
e
x
i i n
24
2 极大似然估计
,
1 1 1n i i
x n
L
ne
L
注意到,是的增函数,
取到最大值时,达到最大。
2
1
1
1 n i 0
i
dlnL n X X
d
令
1
1, 2, , n
,X min X X X
设
1
ˆ X X
1
1 n i ˆ
i
lnL nln
X 又
1 2
0, 0
, , ,
n Xx x x
设设设设设设设设设设设设设设设 设设设设设设设设设设设设设设
设6设
设设设设
26
1
解:极大似然估计
;
1 00
X f x x
因的概率密度为:
其它
1 0 1, , ,20
n x x xn
L
设 设 设 设 设 设 设 设 设 设
设 设
0, ˆL
dln n
d
由于 不能用微分法求
ˆ :L 从义发
以下定出求
1 2
0 xi , x n max x x, , , xn 设 设 设 设 设 设 设 设 设 设 设
1
ˆ
n n
L n
L x
L x L
又对的是减函数,
越小,越大,故时,最大;
1 2
ˆL X n max X X, , , Xn
设设设设设设设设设设设
2 矩估计 0 1 2
E X xdx X
由ˆ 2X
2 ,
0 2 ,
3
X
1 2 3
设总体的概率分布律为:
2 1- 3 其中,未知
现得到样本观测值2,3,2,1,3,
求的矩估计值与极大似然估计
7值。
例:
30
1 解:矩估计
1
E X x p
k k
3 5 2
2 2 3 (1 3 2)
2.2
X ˆ 0.32
1
2 (3 )
5
ˆ 2 (3 )
5 X
2 极大似然估计
( ) ( 2)(1 3 2)( 2) (1 3 2)
L
3 2
1
16 (2 3 )
ln ( )
L
ln16 3ln 2ln(2 3 ) ln ( ) 3 6 2 3 0 d L
d
ˆ 0.4
7.2 估计量的评选准则
从前一节看到,对总体的未知参数可用不 同方法求得不同的估计量,如何评价好坏?
四条评价准则:
无偏性准则,有效性准则,均方误差准则和 相合性准则
32
1. 无偏性准则
1 2ˆˆ
,
, , ,
ˆ
ˆ
, n
X X X
E
设设设设设设设设设设设设
设设设
设设
设设设设设设设
ˆˆ ,, ˆ
n
E li
E m E
若那么称为估计量的 若则
偏差 渐近
称是的 无偏估计量
22
2
, ,
X
E X D X
X S
设总体的一阶和二阶矩存在,分布是 任意的,记
证明:样本均 例
值和样本方差分别是 和
1:
的无偏估计。
35
1, 2, , n
X X
X X
设 设 设 设 设 设 设 设 设 设 设
X
故是的无偏估计.
11 n
i i
E X E X
n
1
1
ni i
n
E X
1 n n
36
2 21
1 ( )
1
n i i
E S E X X
n
22 1
1 ( )
1
n i i
E X n X
n
2 2
21 1
n
n
2 2
S
故是的无偏估计.
21
1 ( )
1
n i i
E X X
n
1
1 1
n
i i
D X nD X
n
.1 0,
ˆ 2
ˆ
L nX X
X
检验7节例6(即总体服从上的均匀分布)
的矩估计量与
例
极大似
2然估计量的
:
无偏性。
38
0, , ,
X U
E X
2 设 设 1, , n
X
X X
设 设 设 设 设 设
ˆ 2
E
E X
1
2 n
i i
n
E X
2 n n 2
ˆ 2X
因此是的无偏估计
39
ˆL
X
nX
n
为考察 的无偏性,先求 的分布,
由第三章第 节知:5
F
X n x
F x
n ,
1
0
0
n
n n X
nx x
f x
于是
其它
1 0
n
x nx
ndx
ˆL
nE
E X
因此有:n
1n
ˆ
L
X
n
所以是有偏的。
纠偏方法
1
ˆ , , , 0
1 ˆ
2 1 ,
ˆ , ,
ˆ
n n n
n
E a b a b a
a b
X n X X
n
X X
设设 设设设设设设设 设设设设设设设设
设设设设设设设设设设设设
设设设设设设设设设设设设设设设设设设设 设设设设设设设设设设设设设
设设设设设设设设设设设设设设设
40
2.
有效性准则
22
1 1 2
1
ˆˆ,
ˆ
ˆ
,ˆ
D ˆ D
定义:设是的两个估计,
如果对一切成立,
且不等号至少对某一成立,则称 比
有效。
无偏
41
1
1 2
1 2
~ 0, , , ,
ˆˆ 2 , 1
ˆˆ 2
n
n
X U X X X
X n X
n n
设设设设设设设设设设
设设设设设设设设设设设 设设设设设设设
设3设
设设
43
ˆ1
2 4 122 32D D X
n
n
解:
1 0 0
n
n n X
nx x
f x
由
其它
ˆ
2 n 21
2 2n
n 2
D E X E X
n 于是
ˆˆˆˆ1 32
2 2
2 2 1D D
n n n
因为比更有效
2n 0 nxnn 1 n 2 2E X dx
n
2
2
n n
ˆ ( ˆ )2
( ).ˆ
E Mse
定义:设是参数的点估计,方差存在,则称 是估计量的均方误差,记为
44
ˆ
Mse( ˆˆ) D( )
若是的无偏估计,则有
在实际应用中,均方误差准则比无偏性准则更重要.
3. 均方误差准则
45
2
2
S
B
2:试利用均方误差准则,对用样本方差 和样本二阶中心矩分别估计正态总体方差时 进
例4
行评价.
46
2 2 2
(n 1)S / ~ (n 1).
解:根据第六章抽样分布定理,在正态总体下,
2 2
2 2 2 4
( ) ( ) .
1 S
Mse S D S
n
又因是的无偏估计,因此
47
2 2
2 2
( ) [( ) Mse B E B 而
2 2
2 2
( ) [ ( ) ] D B E B
2 2 2 2
1 1
(n ) [ (n ) ]
D S E S
n n
4
2n 2 1 n
2
2 2
2 1 2
1 1
. n n
n n
B S
当时,有,
因此在均方误差准则下,优于
4.
相合性准则
1, ,
0, 0
n n
n
n
n
lim X
n P
X
设设设设设设设设设
设设设设设设设设设设设设设设设设 设设设设设设
设 设设设
设设设设设设
设设设 设设设设设
48
1
1
1
2 2 2
2 1
( ) ( 2) , , ,
(1) ( )
2 1 , 2,..., , 2,...,
3 1 ( ) , ( )
k
k n
n l
l i l
i n i i
X k E X k X X
X
X E X
A X l k l k
n
B X X S D X
n S
设设设设设设设5 设设设 设设设设设设设设设设
设设设设设设设 设设设设设设设设设
设设设设设设设设设 (4)设设设设设设设
50
1
1
( ), 1, 2,..., . (1) 2
n l
i i
l l
n X
E X l k
证明:由辛钦大数定律知,依概率收敛到
因此,()成立。
1 1
1
1 1
,..., ,..., ( ,..., )
( ,..., ) ( ,..., )
k k
k
k k
A A
g g A A g
根据依概率收敛的性质,由是 的相合估计,若是连续函数,
则是的相合估计。
2 2
2 1
2 2 2
2 2
1
2 2 2
2 2
( )
1 ( )
1
n i i
D X
B X X A X
n
S n B S
n
S S
因为,所以
是的相合估计,
注意到,因此也是的相合估计;
是的相合估计。
1
1 2
0, , , , ˆˆ 2 1
n
n
X U X X X
X n X
n
设设设设设设设设设设设 设设设设设
设6
设设设设设设
52
ˆˆ1 2 ,E E 证:
0, n
由切比雪夫不等式,当时,
ˆ1
D
2ˆ1P
有: 2 2 0
3n
1 2
ˆˆ
所以和都是的相合估计。
ˆ1 32 ,D n D
ˆ2 n n 2 2
ˆ2
D
2ˆ2P
同理:
2
2 0
2 n n
7.3 区间估计
53
1
1 1 1 2 2 1
1 2
, ,
ˆˆˆˆ , , , , ,
ˆˆ,
n
n n
X X X
X X X X
设设设设设设设设设设设
设设设设设设设设设设设设设设设
设设设设设设设设设设设设设设
54
1
;
, , n 0 1 ,
X F x
X X X
设设7. 3. 1设设设设设设设设设设设设设设设设设设
设设设设设设设设设设设设设设
L L U
U
U
1
1 1
L 1, , , , 1
, , , , ,
n n 7 1
n n
P
X
X X X
X X
X X
设设设设设设设设设设设设
L U
L U
, 1
双侧置信区间置 则称随机区间是的;称为;
和分别
信度 双侧置信下限双侧置
称为和 信上限。
( 一)置信区间的定义
55
置信区间的含义:
),
, , ,
100(1 )%.
n
若反复抽样多次(各次得到的样本容量相等,都为每个 样本值确定一个区间每个这样的区间或者包含的真值
或者不包含的真值。按伯努里大数定律,在这些区间中,
包含真值的约占
0.05, 95%
0.01, 99%
如反复抽样10000次,当即置信水平为时,
10000个区间中不包含的真值的约为500个;当
即置信水平为时, 10000个区间中不包含的真值的约 为100个。
单侧置信限
1
L 1
L L
7 1
, , 1 , 7 2 , ,
, 1
n
n
P X
X X
X
设设7. 3. 2 设设设设设设设设设设设设设
设设设设设
设设设设设设设设设设设
设设设设设设
设设设设设设设
57
U
U 1
U
1, , 1
7 2
7 3 ,
, 1
,
,
n
P X n
X X
X
设设设设设设设
设设设设设
设设设设设设设设设设设
设设设设设设
设设设设设设设
单侧置信限和双侧置信区间的关系:
L L 1
U U 1
L U
, , , ,
n
n
X X
X X
1
2
1 2
设设设设设设设设设1-设设设设设设设设
设设设设设设设设1-设设设设设设设设 设设设设设设设设设设设设1-设设设设设设设设
60
L U U L
[ ,
ˆˆˆˆ ]E
(
) 定义:称置信区间的平均长度为区间的,并称二分之一区间的 平均长度为置信区间的
精确度
误差限。
说明:在给定的样本容量下,置信水平和精确度 是相互制约的。