"> "> 数理统计 - 点估计 | Yufei Luo's Blog

数理统计 - 点估计

定义

定义在统计结构分布族 P={Pθ,θΘ} 上的一个实值泛函 g(P) 称为参数,而统计结构 (X,B,P) 上的用来估计 g(P) 的实值统计量称为 g(P) 的点估计量,简称估计。为简单起见,常用 θ 表示参数,以 θ^=θ^(X) 表示其点估计。

一般来说,任何定义在 Θ 上的实值函数都可以称为参数,例如均值、方差等特征数。

估计的优良性

均方误差

为了排除样本随机性的影响,最常用的评价估计的标准是均方误差: MSE=E(θ^(X)θ)2 自然我们希望 MSE 越小越好,但是使得均方误差一致达到最小的最优估计是不存在的(可以用反证法证明,假设存在,那么对于任意的 θ0Θ 应取 θ^(X)=θ0。但是由于 θ0 具有任意性,因此无法找到这样的最优估计)。因此通常先对估计提出一些合理要求,然后在满足合理要求的范围内选取好的估计。

无偏性

由于 MSE=E(θ^(X)θ)2=E(θ^(X)22θ^(X)θ+θ2)=E(θ^(X)2)2E(θ^(X))E(θ)+E(θ2)=E(θ^(X)2)E(θ^(X))2+E(θ^(X))22E(θ^(X))θ+θ2=Var(θ^(X))+(E(θ^(X))θ)2 也就是说 MSE 可以分解为 θ^(X) 的偏差和方差之和。如果偏差 E(θ^(X))θ 等于 0,就是所谓的无偏估计。

无偏估计的正式定义如下:

(X,B,P={Pθ,θΘ}) 为可控参数统计结构,g(θ) 是未知数,X=(X1,,Xn) 为来自该统计结构的一个样本,如果用 g^(X) 估计 g(θ),且 Eθ(g^(X))=g(θ),则称 g^(X) g(θ) 的无偏估计。

无偏估计体现了一种频率思想,只有在大量重复使用时,无偏性才有意义。例如某一工厂每天对其生产的样品进行抽检,如果假定其生产过程相对稳定,则估计的无偏性要求便是合理的。

例:设 X1,,Xn 是来自于 N(μ,σ2) 的一个样本,μ σ 都是未知参数。其常用的估计分别是样本均值 X¯=1ni=1nXi 和样本方差 sn2=1ni=1n(XiX¯)2。由于 E(X¯)=μ,而 nsn2/σχn12E(sn2)=n1nσ2,因此 X¯ μ 的无偏估计,而 s2n 不是 σ2 的无偏估计。但是如果将 sn2 修正为 s2=nn1sn2,则可以得到 σ2 的无偏估计。

当样本容量较大时,sn2 s2 将会很接近,此时称 sn2 σ2 的渐进无偏估计。

g^n=g^n(X1,,Xn) g(θ) 的估计量,如果 limnEθ(g^n)=g(θ),  θΘ 则称 g^n g(θ) 的渐进无偏估计。

相合性

估计量是与样本容量有关的,假设用 θ^n=θ^n(X1,,Xn) 估计 θ,不可能做到对某一个 nMSE(θn^) 对所有的 θΘ 任意小,但是当 n 时往往可以做到这一点,这对应于相合性的概念。需要注意的是,相合性只是指出了当 n 时估计量的性质,但是对于任意有限的 n,相合性是没有意义的。而且相合估计可以不止一个。

θ^n=θ^n(X1,,Xn) θ 的估计量,如果当 n 时有 θ^nPθ,则称 θ^n θ 的弱相合估计。进一步,如果 θ^nθ,a.s.,则称 θ^n θ 的强相合估计。在统计研究中,一般弱相合性就已经足够。

相合性被认为是对估计的一个最基本要求,如果无论有多少观测值都无法把要估计的参数估计到任意指定的精度,那么这个估计是很值得怀疑的。

要证明相合性,需要使用描述极限的 ϵ 语言,通常使用大数定律或者是依概率分布的定义来构造。

例:设 X1,,Xn 是来自于 U(0,θ) 的一个样本,最大次序统计量 X(n) θ 的常用估计。由于 X(n) 的密度函数为 p(t,θ)=ntn1θn,0<t<θ,可求得 E(X(n))=nθ/(n+1)。因此 X(n) 只是 θ 的渐进无偏估计而不是无偏估计。此外,对于任意的 ϵ>0Pθ(|X(n)θ|ϵ)=Pθ(X(n)θϵ)=0θϵntn1θndt=(θϵθ)n0(n) 因此 X(n) θ 的相合估计。

对于相合估计有如下性质,可以方便计算:

Tjn=Tjn(X1,,Xn) gj(θ) 的相合估计,j=1,,k。函数 h() (g1(θ),,gk(θ)) 处连续,则 h(T1n,,Tkn) h(g1(θ),,gk(θ)) 的相合估计。

例:设 X1,,Xn 是来自于密度函数为 p(x;θ)=θ(θ+1)xθ1(1x),0<x<1,θ>0 的一个样本,由于 E(X1)=θ/(θ+2),为了方便将其记为 η

由强大数定律,有:X¯η, a.s.,而 θ=2η/(1η) η 的连续函数,因此根据上述定理可直接计算得到 2X¯/(1X¯) θ 的相合估计。

渐进正态性

θ^n=θ^n(X1,,Xn) θ 的估计量,如果存在 σn2(n) 满足 (θ^nθ)/σn(θ)LN(0,1) 则称 θ^n θ 的渐进正态估计,σn2(n) 称为 θ^n 的渐进方差,记作 θ^nAN(θ,σn2(n))。根据中心极限定律,一般 σn2(n) 可以取为 θn^ 的方差。

渐进正态性只是反映了当 n 时估计的性质,无法说明达到所需要的精度所需的样本容量。

此外,g(θ) 的渐进正态估计可以有很多,它们之间的优劣可以使用相对渐进效来对比,即两个不同估计 σn2(n) 的比值。

有效性

信息不等式:设 {Pθ,θΘ} 是 C-R 正则族,ΘRk,其 Fisher 信息 I(θ) 为非奇矩阵。并设 g(θ)=(g1(θ),,gs(θ)),sk,且 gi(θ)/θj 都存在。假设 T(X) g(θ) 的模平方可积的无偏估计,记 Δ=E(TSθ)=ddθg(θ),则有: Varθ(T(X))ΔI1(θ)Δ 其中 ΔI1(θ)Δ 称为 g(θ) 的无偏估计协差阵的下界,或者称为 g(θ) 的无偏估计的 C-R 下界。在一维的情况下,不等式简化为:Varθ(T(X))(g(θ))2I1(θ)

如果某个无偏估计的方差可以达到这个下界,那么它也自然就是 UMVUE。

估计的效:设 {Pθ,θΘ} 是 C-R 正则族,g(θ) 是可估参数,T(X) g(θ) 的一个无偏估计,则称 (g(θ))2I1(θ)/Varθ(T(X)) 为估计 T(X) 的效,如果效等于 1,则称 T(X) g(θ) 的有效无偏估计。

无偏估计

无偏性是统计问题中应用很广的一个准则,对于无偏估计,需要注意以下三点:

  1. 无偏估计不一定存在。在统计中,一般将存在无偏估计的参数函数称为可估参数。
  2. 对于可估参数,无偏估计一般不唯一。
  3. 无偏估计不一定是好估计

一致最小方差无偏估计

g(θ) 为可估参数,如果 T(X) g(θ) 的无偏估计,且对于所有无偏估计中的任意一个估计 φ(X),都有 Varθ(T(X))Varθ(φ(X)),  θΘ 成立,则称 T(X) g(θ) 的一致最小方差无偏估计(Uniformly Minimum Variance Unbiased Estimate,UMVUE)。

根据如下引理,在构造 UMVUE 时,通常需要使用到充分统计量,它可以降低无偏估计的方差:

S(X) 是分布族 {pθ,θΘ} 的充分统计量,φ(X) g(θ) 的无偏估计,令 T(X)=E(φ(X)|S(X)),则 T(X) 也是 g(θ) 的无偏估计,且 Var(T(X))Var(φ(X))

这一引理的证明如下:

因为 S(X) 为充分统计量,故 T(X) θ 无关,它也是统计量。

根据重期望公式,E(T(X))=E[E(φ(X)|S(X))]=E(φ(X))=g(θ)

Var(φ(X))=E(φ(X)g(θ))2=E(φ(X)T(X))2+E(T(X)g(θ))2Var(T(X))

上式中的交叉乘积项 E[(φ(X)T(X))(T(X)g(θ))]=E{E[(φ(X)T(X))(T(X)g(θ))|S(X)]}=E{[T(X)g(θ)]E[(φ(X)T(X))|S(X)]}(因 T (X)) 是 S(X) 的函数,故在条件期望的情况下可以做常数处理)=0E[(φ(X)T(X))|S(X)]=0

由上述引理,可以得到如下的定理:设 S(X) 为分布族 {pθ,θΘ} 的完备充分统计量,g(θ) 是可估参数,则 g(θ) 的 UMVUE 存在,它是 S(X) 的函数而且在几乎处处意义下是唯一的。证明过程如下:

在所有的无偏估计中任取一个 φ(X),令 T(X)=E(φ(X)|S(X)),根据上述引理,T(X) S(X) 的函数,且是 UMVUE。

在无偏估计中可以在另外任取一个 g^(X),令 T(X)=E(g^(X)|S(X)),它也是 S(X) 的函数,并且有 E(T(X)T(X))=0 成立。即 (T(X)T(X))p(s)ds=0

根据 S(X) 的完备性,可得 T(X)=T(X)。这也就是说,从任意无偏估计出发,均可以得到一个相同的 T(X),它是完备充分统计量 S(X) 的函数,是 g(θ) 几乎处处唯一的 UMVUE。

因此,UMVUE 的求法有两种:

  1. 寻找完备充分统计量的函数使之称为 g(θ) 的无偏估计。通常可以通过计算完备充分统计量的期望、方差等的含有参数 θ 的表达式,来反推出 g(θ)

    例:设 X1,,Xn 是来自 b(1,θ),0<θ<1 的一个样本,由指数型分布族的性质可知,S(X)=Xi 是完备充分统计量。

    对于 θ,由于 E(S(X))=nθ,因此 X¯=1nS(X) θ 的 UMVUE。

  2. 任取 g(θ) 的一个无偏估计,然后对完备充分统计量求条件期望。由于无偏估计可以任取,因此简便起见,通常是将无偏估计构造成指示函数的形式。

    例:接上面的例子,对于 g(θ)=θk+(1θ)nk,如果要直接找一个 S(X) 的函数 h(S(X)) 使之成为 g(θ) 的无偏估计是很困难的。因此考虑先使用指示函数构造一个无偏估计,令: φ1(X)={1,  i=1kXi=k0,   elseφ2(X)={1,  i=k+1nXi=00,   else φ(X)=φ1(X)+φ2(X),则 φ(X) g(θ) 的无偏估计。因此,可以通过计算 E(φ(X)|S(X)) 来得到 UMVUE。记 S(X)=s,此时要分三种情况讨论:

    1. s<k,此时 φ1(X)=0 一定成立,因此有: E(φ(X)|S(X)=s)=E(φ2(X)|S(X)=s)=P(φ2(X)=1|S(X)=s)=P(i=k+1nXi=0,S(X)=s)P(S(X)=s)=(1θ)nk(ks)θs(1θ)ks(ns)θs(1θ)ns=(1θ)nk(ks)θs(1θ)ks(ns)θs(1θ)ns=(ks)(ns)

    2. s>k,此时 φ2(X)=0 一定成立,此时有: E(φ(X)|S(X)=s)=E(φ1(X)|S(X)=s)=P(φ1(X)=1|S(X)=s)=P(i=0kXi=1,S(X)=s)P(S(X)=s)=(θ)k(nksk)θsk(1θ)ns(ns)θs(1θ)ns=(nksk)(ns)

    3. s=k,此时有: E(φ(X)|S(X)=s)=P(φ1(X)=1|S(X)=s)+P(φ2(X)=1|S(X)=s)=P(i=0kXi=1,S(X)=s)P(S(X)=s)+P(i=k+1nXi=0,S(X)=s)P(S(X)=s)=(θ)s(1θ)ns(ns)θs(1θ)ns+(θ)s(1θ)ns(ns)θs(1θ)ns=2(ns)

    综上,记 T(X) φ(X) 的 UMVUE,可得: T(X)={(ks)(ns), s<k2(ns), k=s(nksk)(ns), s>k

    例 2:某厂生产的产品其废品率为 θ,现将该产品包装成盒,每盒抽 n 个产品进行检验(设盒中产品数远大于 n),得到废品数 X,可以认为 Xb(n,θ)。当 X2 时接收该盒产品,其它情况拒收。通过概率 g(θ)=(1θ)n+n(1θ)n1θ+n(n1)2(1θ)n2θ2 厂方自然很关心 g(θ) 的估计。

    假设抽了 r 盒产品进行检验,第 i 盒的废品数记为 XiXib(n,θ)。则 S(X)=i=1rXr 为完备充分统计量,且 S(X)b(rn,θ)

    φ(X)={1, X12(第一盒被接收)0, 其它 E(φ(X))=g(θ) g(θ) 的无偏估计量。下面计算 E(φ(X)|S(X)) 的值。为了叙述方便,记 B0B1B2 分别对应于事件第 1 盒样品中的废品数为 0,1,2,这三个事件不相容,因此有: E(φ(X)|S(X)=s)=P(B0|S(X)=s)+P(B1|S(X)=s)+P(B2|S(X)=s)=P(B0,S(X)=s)+P(B1,S(X)=s)+P(B2,S(X)=s)P(S(X)=s)=(rnns)+n(rnns1)+(n1)n/2(rnns2)(rns) 上式即为 g(θ) 的 UMVUE。

估计方法

矩估计

定义

矩估计指的是用样本矩及其相应的函数估计相应的总体矩及其函数。具体表述如下:

X1,,Xn 是来自于某总体的一个样本,以 μr 表示总体的 r 阶原点矩,mr 表示从样本 X1,,Xn 得到的 r 阶样本原点矩,即 μr=EX1rmr=1nXir。如果某个参数 θ 可以表示为总体前 k 阶矩的函数,即 θ=g(μ1,,μk),则我们可以用 θ^(X)=g(m1,,mk) 去估计 θθ^(X) 即称为 θ 的矩估计。

需要注意的是,矩估计不唯一。例如来自泊松分布 P(λ) 的样本 X1,,Xn,由于 EX1=λVar(X1)=λ,因此样本均值和样本方差都可以作为 λ 的矩估计。在矩估计不唯一时,通常根据如下两个原则来选择矩估计:1. 涉及到矩的阶数尽可能小,常用的矩估计通常只涉及一二阶矩;2. 所用的估计最好是充分统计量的函数。

特点与性质

矩估计的特点为:

  1. 矩估计基于经验分布函数,而经验分布函数逼近真实分布函数的前提条件是样本容量较大。因此矩估计是以大样本为应用对象的。
  2. 矩估计没有使用到总体分布的任何信息,本质上讲它是一种非参数方法。

对于矩估计,有如下两个大样本性质:

相合性:设 X1,,Xn 为独立同分布变量序列,E|X1|k<θ=g(μ1,,μk),并且设 g 是连续的,则矩估计 θ^=g(m1,,mk) θ 的相合估计。

渐近正态性:设 X1,,Xn 为独立同分布变量序列,E|X1|2k<,记 θ=(θ1,,θs),其中 θi=gi(μ1,,μk),并且设 θi μj 有连续偏导数。记 G=(gi/μj)s×k,则对 θ 的矩估计 θ^=(g1(m1,,mk),,gs(m1,,mk)),有: n(θ^θ)LNs(0,GΣG) 其中 Σ k×k 阶矩阵,其 (i,j) 元素为 μi+jμiμj

极大似然估计

定义

p(x;θ)θΘ (Rn,BRn) 上的一族联合概率密度函数,对于给定的 x,称 L(θ;x)=p(x;θ) θ 的似然函数。进一步,如果存在 (Rn,BRn) (Θ,BΘ) 的统计量 θ^(x) 使得 L(θ^(x);x)=supθL(θ;x) θ^(x) 称为 θ 的一个极大似然估计(Maximum Likelihood Estimate,MLE)。

由于概率密度大多具有指数函数形式,因此通常使用似然函数的对数。由于对数变换严格单调增,因此求极大值时等价。寻找 MLE 的常用办法是求导数。

性质

不变性:如果 θ^ θ 的 MLE,g() 是可测函数,则 g^(θ) 也是 g(θ) 的 MLE。

相合性:如果 lnp(x;θ) Θ 上可微,并设 p(x;θ) 是可识别的(θθ{x:p(x;θ)p(x;θ)} 不是零测集),则似然方程在 n 时以概率 1 有解,且此解关于 θ 是相合的。

渐近正态性:假设 Θ 为开区间,概率密度函数 p(x;θ),θΘ 满足:1. 在参数真值 θ0 的邻域内,lnp/θ2lnp/θ23lnp/θ3 对所有的 x 都存在;2. 在参数真值 θ0 的邻域内, |3lnp/θ3|H(x),且 E(H(x))<;3. 在参数真值 θ0 处,Eθ0[p(X,θ0)p(X,θ0)]=0Eθ0[p(X,θ0)p(X,θ0)]=0I(θ0)=Eθ0[p(X,θ0)p(X,θ0)]2>0(Fisher 信息量)。记 θ^n n 时似然方程的相合解,则有: n(θ^nθ0)LN(0,I1(θ0)) 在 C-R 正则族场合,定理条件一般是满足的,因此 MLE 是渐近正态的。

渐近有效性:在一定的正则条件下,MLE 的渐近方差为 [nI(θ)]1,它正是由容量为 n 的样本得到的 θ 的无偏估计的方差下界。因此它具有渐近有效性。

最小二乘估计

定义

最小二乘法是一种常用的估计方法,最常见于线性模型。考虑 Gauss-Markov 模型:EY=XβVar(Y)=σ2In,其中 Y n×1 维观测向量,X 为已知的 n×p 维矩阵,β p×1 维未知参数,σ2 未知。这一模型也被称为独立线性观测模型,通常记为 (Y,Xβ,σ2In)。如果 (YXβ^)(YXβ^)=minβ(YXβ)(YXβ) 则称 β^ β 的最小二乘估计(Least Squares Estimate,LSE)。其中, β^=(XX)1XY

性质

  1. β^ β 的无偏估计,它的协方差矩阵为 σ2(XX)1
  2. 如果 X 是列满秩矩阵,那么 β^ β 唯一的最好线性无偏估计(Best Linear Unbiased Estimate,BLUE),也是 UMVUE。
  3. 如果 X 是列满秩矩阵,那么 σ^2=(YXβ^)(YXβ^)/(np) σ2 的 UMVUE。
0 条评论
未登录用户
支持 Markdown 语法

来做第一个留言的人吧!