第七章参数估计讲解
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第七章 参数估计
参数估计是数理统计研究的主要问题之一. 假设总体X ~N (μ,σ2),μ,σ2是未知参数,X 1,X 2,…,X n 是来自X 的样本,样本值是x 1,x 2,…,x n ,我们要由样本值来确定μ和σ2的估计值,这就是参数估计问题,参数估计分为点估计(Point estimation )和区间估计(Interval estimation).
第一节 点估计
所谓点估计是指把总体的未知参数估计为某个确定的值或在某个确定的点上,故点估计又称为定值估计.
定义7.1 设总体X 的分布函数为F (x ,θ),θ是未知参数,X 1,X 2,…,X n 是X 的一样本,样本值为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 )的方法很多,下面仅介绍矩法和极大似然估计法. 1.矩法
矩法(Moment method of estimation )是一种古老的估计方法.它是由英国统计学家皮尔逊(K .Pearson )于1894年首创的.它虽然古老,但目前仍常用.
矩法估计的一般原则是:用样本矩作为总体矩的估计,若不够良好,再作适当调整. 矩法的一般作法:设总体X ~F (X ;θ1,θ2,…,θl )其中θ1,θ2,…,θl 均未知. (1) 如果总体X 的k 阶矩μk =E (X k ) (1≤k ≤l)均存在,则
μk =μk (θ1,θ2,…,θl ),(1≤k ≤l ).
(2) 令⎪⎪
⎩⎪⎪⎨⎧.
),,,(,),,,(,
),,,(212
2121211l l l l l A A A θθθμθθθμθθθμ
其中A k (1≤k ≤l )为样本k 阶矩.
求出方程组的解,ˆ,,ˆ,ˆ21l θθθ 我们称),,,(ˆˆ21n k k X X X θθ=为参数θk (1≤k ≤l )的矩估计量, ),,,(ˆˆ21n
k k x x x θθ=为参数θk 的矩估计值. 例7.1 设总体X 的密度函数为:
f (x )=⎩⎨⎧-><<+.,
0),
1(,10,)1(其他αααx x
其中α未知,样本为(X 1,X 2,…,X n ),求参数α的矩法估计.
解 A 1=X .由μ1=A 1及
μ1=E (X )=
2
1
)1()(1
++=
+=⎰
⎰+∞
∞
-ααααx x x x x xf d d , 有21++=
ααX ,得1
21ˆ--=X X
α.
例7.2 设X ~N (μ,σ2),μ,σ2未知,试用矩法对μ,σ2进行估计. 解
⎪⎪⎩
⎪⎪⎨⎧======∑∑==.
1)(,1)(122221
11n
i i n
i i X n A X E X n A X E μμ 又 E (X )=μ, E (X 2)=D (X )+(EX )2=σ2+μ2,
那么 .1ˆˆ,ˆ2
222S n
n A X -=-==μσμ
. 例7.3 在某班期末数学考试成绩中随机抽取9人的成绩.结果如下:
试求该班数学成绩的平均分数、标准差的矩估计值.
解 设X 为该班数学成绩,μ=E (X ),σ2=D (X )
)558994(9
1
9191+++==∑= i i x x =75;
2
/19
122)(819898⎥⎦
⎤⎢⎣⎡-⋅=∑=i i x x s =12.14.
⎪⎪⎩
⎪⎪⎨⎧======∑∑==.91)(,91)(9
12
2229
1
11i i i i X A X E X A X E μμ 由于E (X 2)=D (X )+(EX )2=σ2+μ2,那么,
2222228
ˆˆˆ,().9
X A A x S μ
σμ==-=-= 所以,该班数学成绩的平均分数的矩估计值x =μ
ˆ=75分,标准差的矩估计值2
9
8ˆs =σ
=12.14. 作矩法估计时无需知道总体的概率分布,只要知道总体矩即可.但矩法估计量有时不惟
一,如总体X 服从参数为λ的泊松分布时,X 和B 2都是参数λ的矩法估计.
2.极(最)大似然估计法
极大似然估计法(Maximum likelihood estimation)只能在已知总体分布的前提下进行,为了对
它的思想有所了解,我们先看一个例子.
例7.4 假定一个盒子里装有许多大小相同的黑球和白球,并且假定它们的数目之比为3∶1,但不知是白球多还是黑球多,现在有放回地从盒中抽了3个球,试根据所抽3个球中黑球的数目确定是白球多还是黑球多.
解 设所抽3个球中黑球数为X ,摸到黑球的概率为p ,则X 服从二项分布
P {X =k }=k 3C p k
(1-p )3-k , k =0,1,2,3.
问题是p =1/4还是p =3/4?现根据样本中黑球数,对未知参数p 进行估计.抽样后,共有4种
可能结果,其概率如表7-1所示.
假如某次抽样中,只出现一个黑球,即X =1,p =1/4时,P {X =1}=27/64;p =3/4时,P {X =1}=9/64,这时我们就会选择p =1/4,即黑球数比白球数为1∶3.因为在一次试验中,事件“1个黑球”发生了.我们认为它应有较大的概率27/64(27/64>9/64),而27/64对应着参数p =1/4,同样可以考虑X =0,2,3的情形,最后可得
p =⎪⎩⎪⎨⎧==.3,2,4
3
,1,0,41
时当时当x x
(1) 似然函数
在极大似然估计法中,最关键的问题是如何求得似然函数(定义下文给出),有了似然函数,问题就简单了,下面分两种情形来介绍似然函数. (a ) 离散型总体
设总体X 为离散型,P {X =x }=p (x ,θ),其中θ为待估计的未知参数,假定x 1,x 2,…,x n 为样本X 1,X 2,…,X n 的一组观测值.
P {X 1=x 1,X 2=x 2,…,X n =x n }=P {X 1=x 1}P {X 2=x 2}…P {X n =x n }
=p (x 1,θ)p (x 2,θ)…p (x n ,θ)
=
∏=n
i i
x p 1
),(θ.
将
∏=n
i i
x p 1
),(θ看作是参数θ的函数,记为L (θ)
,即 L (θ)=
∏=n
i i
x p 1
),(θ. (7.1)
(b ) 连续型总体
设总体X 为连续型,已知其分布密度函数为f (x ,θ),θ为待估计的未知参数,则样本(X 1,X 2,…,X n )的联合密度为:
f (x 1,θ)f (x 2,θ)…f (x n ,θ)=
∏=n
i i
x f 1
),(θ.