十分钟学习极大自然似估计

EndlessLethe原创文章,转载请注明: 转载自小楼吹彻玉笙寒
原文链接地址: 十分钟学习极大似然估计

前言

参数估计是机器学习里面的一个重要主题,而极大似然估计是最传统、使用最广泛的估计方法之一。

本文主要介绍了极大似然估计,简单说明了其和矩估计、贝叶斯估计的异同,其他估计(如MAP)并不涉及。

为什么要用极大似然估计

对于一系列观察数据,我们常常可以找到一个具体分布来描述,但不清楚分布的参数。这时候我们就需要用极大似然估计来求解这个分布的参数。换句话说,极大似然估计提供了一种给定观察数据来评估模型参数的方法,即:“模型已定,参数未知”。

极大似然估计概述

下面结合一个例子介绍极大似然估计法的思想和方法:

设一个袋子中有黑、白两种球,摸到白球的概率为p,现在要估计p的值。
我们令总体X为

X = \left. \begin{cases} 0,\quad 从袋中取得一白球,\\ 1,\quad 从袋中取得一黑球.\\ \end{cases} \right.

则X服从01分布B(1,p)

我们先进行有放回地摸球10次,其结果可用随机变量xi表示,则x1,x2,⋯,x10是来自总体X的一个样本。其值=(1,0,1,0,0,0,1,0,0,0),则似然函数为L(p)=p^3 (1−p)^7

极大似然估计其实是理想地认为,对于极少的样本观测,我们观测到的样本很可能就是发生概率最大的。

似然函数L(p)是每个样本出现概率的乘积\prod_{i=1}^N {P({x_i})},因为显然样本是独立同分布的。
根据极大似然估计的思想,我们需要让L(p)最大,把这时对应的\hat p作为我们的估计值。

求解L(p)的最大值点\hat p,可由一阶导数

\frac{dL(p)}{dp}=0

确定。更一般的,我们通常可以假设白球出现次数为k,可以解得

\hat p = \frac{k}{N}

这里带入k=3\hat p,所以我们把0.3作为摸到白球的概率

值得注意的是,根据似然函数来求解参数的过程,与样本数量是无关的。我们可以使用变量x_i来描述样本观测值,并将模型参数θ用来x_i表示。当样本较少时,极大似然估计偏差较大。但随着样本的增多(样本逐渐靠近总体分布),偏差慢慢减少为0。这意味着,极大似然估计是非常普适的。

实际上,即使直观上“极大似然估计”似乎是非常自然的想法,但它能在统计学中拥有堪比牛顿力学在物理学中的地位,是因为这种朴素的想法背后蕴含了估计量的泛函不变性、相合性、渐近有效性和渐进正态等诸多逆天的性质。

Note:极大似然估计暗合了切比雪夫大数定律。比如在这个例子中,如果放回次数变得极大,那么根据大数定律也有\hat p = \frac{k}{N} 。所以在用“局部估计整体”时,可以说使用了极大似然估计法,也可以说根据大数定律。

极大似然估计的具体步骤

我们需要做四步:表示似然函数、假设样本观测值、求解方程和代入数据。

似然函数

对于离散型和连续型随机变量,极大似然估计值\hat \theta都满足:

L(\hat \theta)=\max{L(\theta)}

只不过似然函数L(θ)的表示方式略有不同:
离散型随机变量的似然函数是L(\theta) = \prod_{i=1}^N P({x_i}),而连续型是L(\theta) = \prod_{i=1}^N f({x_i})

样本假设

假设样本观测值为x_i

求解方程
由定义可知,估计值可由一阶导数

\frac{dL(p)}{dp}=0

解得。但由于lnL和L在同一位置取得最大值,所以极大似然估计值也可以由对数似然方程

\frac{d(\ln{L(p)})}{dp}=0

解得。

深入的数理统计理论可以证明:当总体分布服从单峰分布时,如果上两式有解,则其解就是θ的极大似然估计值。
Note:当方程无解时,应从定义出发,考虑L(θ)的单调性,找到maxL(θ)对应的估计值。

带入数据

x_i用真实数据替换

多参数极大似然估计

当总体X的分布中含有多个未知参数,即\theta=(\theta_1,\theta_2,⋯,\theta_k)时,似然函数为L(\theta)=L(\theta_1,\theta_2,⋯,\theta_k)。有对数似然方程组:

\left\{ \begin{aligned} \frac{\partial(\ln L(\theta))}{\partial{\theta_1}} & =0\\ \frac{\partial(\ln L(\theta))}{\partial{\theta_2}} & =0\\ \cdots \cdots \cdots\\ \frac{\partial(\ln L(\theta))}{\partial{\theta_k}} & =0\\ \end{aligned} \right.

方程的解是其对应的极大似然估计值。

总结

最大似然估计的特点

  1. 比其他估计方法更加简单
  2. 收敛性:无偏或者渐近无偏
  3. 如果假设的类条件概率模型正确,则通常能获得较好的结果。但如果假设模型出现偏差,将导致非常差的估计结果。

常用分布的估计值

任意概率分布:\hat p = \frac{k}{N}
正态分布:\hat \mu = \overline x = \frac{1}{n}\sum_{i=1}^n{x_i}\hat{\sigma^2}=\frac{1}{n}\sum_{i=1}^N(x_i – \overline x)^2
泊松分布:\hat \lambda = \overline x = \frac{1}{n}\sum_{i=1}^n{x_i}
均匀分布:\hat a = \max{x_i},\hat b = \max{x_i}

最大似然估计和矩估计

对于正态分布和泊松分布,两种估计都是一致的(均匀分布不一致)。矩估计的优点是简单,只需知道总体的矩,总体的分布形式不必知道。而最大似然估计则必须知道总体分布形式,并且在一般情况下,似然方程组的求解较复杂,往往需要在计算机上通过迭代运算才能计算出其近似解。

最大似然估计和贝叶斯估计

最大似然估计是传统频率派估计参数的方法,而贝叶斯估计是贝叶斯统计派估计参数的方法。前者认为θ是一个固定的值,而后者认为θ是一个满足某个分布的随机变量。后者的泛化能力更好,但计算更复杂。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 219,635评论 6 508
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,628评论 3 396
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 165,971评论 0 356
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,986评论 1 295
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 68,006评论 6 394
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,784评论 1 307
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,475评论 3 420
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,364评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,860评论 1 317
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,008评论 3 338
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,152评论 1 351
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,829评论 5 346
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,490评论 3 331
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,035评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,156评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,428评论 3 373
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,127评论 2 356

推荐阅读更多精彩内容

  • 文章作者:Tyan博客:noahsnail.com[http://noahsnail.com] | CSDN[ht...
    SnailTyan阅读 58,165评论 6 99
  • 极大似然估计 以前多次接触过极大似然估计,但一直都不太明白到底什么原理,最近在看贝叶斯分类,对极大似然估计...
    MiracleJQ阅读 1,221评论 0 8
  • 》拟合训练数据和寻找能够泛化到新数据的参数有哪些不同。 》如何使用额外的数据设置超参数。 》机器学习本质上属于应用...
    N89757阅读 2,139评论 0 0
  • 循环的世界 今天你 明天我
    多年以后DNYH阅读 171评论 0 0
  • 叶子想去见那人,可是终究忍住了。 她已快忘了他长什么样。只是,隔一些时日,他就会出现在她的梦里。 在他那里,叶子始...
    亦缘君阅读 359评论 0 0