最大似然估计

  最大似然估计,貌似最没用一种数学方法。一般情况下,其计算结果凭直觉就能猜出来。或许特殊情况下,当问题超出人类直觉能力时,这种方法会有用武之地。
  
  我们先来看一个例子。

例子

[例] 从一批产品中随机抽取5件,经过测试,其中1件不合格。问这批产品不合格率是多少?

  
  所谓不合格率,就是不合格的总数量除以这一批产品的总数量。换句话讲,就是任意抽取一件,产品不合格的概率。
  
  这个问题猪脑袋也能想明白呀!不合格率:
  

p=15

质疑

  面对上面的做法,必然会有以下疑问:
  

问:用抽样产品的合格率估计整体合格率,为什么可以这样做?   
答:不这样做,难道还有更好的选择吗?

问:为什么没有更好地选择了?   
答:如果有,你给我找出来呀!

  
  难道说,暂时找不到更好的方法,就意味着不存在更好的办法了吗?最大似然估计,能够给上述的做法做出一个合理的解释。

最大似然估计给出的解释

  利用最大似然估计方法,可以证明,当整体合格率是 15 时,抽取5个样品,出现1个不合格这样一种抽样结果的概率最大。证明方法如下:
  
  假设不合格率为 p ,则当前抽样结果出现的概率为:

L(p)=p(1p)4

  两边取对数,得

l(p)=ln(L(p))=ln(p)+4ln(1p)
 
  接下来看看 p 取何值,L(p) 取得最大值。上式对 p 求导
l(p)p=1p41p=0

解得,
p=15

 
  也就是说,整体不合格率为 15 时,产生上述抽样结果的概率最大。说不上多有道理,但是最大似然估计给人类的直觉找到了一个貌似合理的解释。
  
  最大似然估计是建立在这样的思想上: 已知某个参数能使这个样本出现的概率最大,而其它参数使这个样本出现的概率减小,所以干脆就把这个参数作为估计的真实值。
  

最大似然估计的一般模式

  
  假设学生考试成绩符合正态分布,
  

p(x)=12πσexp((xμ)22σ2)

  我们随机抽取3个成绩 x1x2x3 。用最大似然估计的方法判断一下这次考试的平均成绩。哈哈,你可能已经猜到答案了!不过,建议还是耐心看一下计算过程,了解一下这种方法。

  1、构造似然函数:似然函数实际上就是当前抽样结果出现的概率。本问题要估值的参数是均值 μ ,因此,似然函数的自变量就是 μ ,

L(μ)=p(x1)p(x2)p(x3)

  2、两边取对数,

l(μ)=ln(L(μ))=ln(p(x1))+ln(p(x2))+ln(p(x3))

把正态分布表达式代入,化简,
l(μ)=(x1μ)22σ2(x2μ)22σ2(x3μ)22σ2+C

  3、求导,求极大值
l(μ)μ=K(x1+x2+x33μ)=0

于是, μ 取下面的值时,本抽样结果出现的概率最大,
μ=x1+x2+x33

  前面推导过程中, CK 表示常数项。

  虽然经过一番高等数学的推导,结果和用小学数学知识计算的结果完全一样。这一结果再次表明,最大似然估计真是毫无用途的统计方法呀!

  • 3
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

许野平

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值