最大似然估计提供了一种给定观察数据来评估模型参数的方法,即:“模型已定,参数未知”。简单而言,假设我们要统计全国人口的身高,首先假设这个身高服从服从正态分布,但是该分布的均值与方差未知。我们没有人力与物力去统计全国每个人的身高,但是可以通过采样,获取部分人的身高,然后通过最大似然估计来获取上述假设中的正态分布的均值与方差。

最大似然估计中采样需满足一个很重要的假设,就是所有的采样都是独立同分布的。下面我们具体描述一下最大似然估计:

首先,假设

为独立同分布的采样,θ为模型参数,f为我们所使用的模型,遵循我们上述的独立同分布假设。参数为θ的模型f产生上述采样可表示为

回到上面的“模型已定,参数未知”的说法,此时,我们已知的为

,未知为θ,故似然定义为:

在实际应用中常用的是两边取对数,得到公式如下:

其中

称为对数似然,而

称为平均对数似然。而我们平时所称的最大似然为最大的对数平均似然,即:

举个别人博客中的例子,假如有一个罐子,里面有黑白两种颜色的球,数目多少不知,两种颜色的比例也不知。我 们想知道罐中白球和黑球的比例,但我们不能把罐中的球全部拿出来数。现在我们可以每次任意从已经摇匀的罐中拿一个球出来,记录球的颜色,然后把拿出来的球 再放回罐中。这个过程可以重复,我们可以用记录的球的颜色来估计罐中黑白球的比例。假如在前面的一百次重复记录中,有七十次是白球,请问罐中白球所占的比例最有可能是多少?很多人马上就有答案了:70%。而其后的理论支撑是什么呢?

我们假设罐中白球的比例是p,那么黑球的比例就是1-p。因为每抽一个球出来,在记录颜色之后,我们把抽出的球放回了罐中并摇匀,所以每次抽出来的球的颜 色服从同一独立分布。这里我们把一次抽出来球的颜色称为一次抽样。题目中在一百次抽样中,七十次是白球的概率是P(Data | M),这里Data是所有的数据,M是所给出的模型,表示每次抽出来的球是白色的概率为p。如果第一抽样的结果记为x1,第二抽样的结果记为x2... 那么Data = (x1,x2,…,x100)。这样,

P(Data | M)

= P(x1,x2,…,x100|M)

= P(x1|M)P(x2|M)…P(x100|M)

= p^70(1-p)^30.

那么p在取什么值的时候,P(Data |M)的值最大呢?将p^70(1-p)^30对p求导,并其等于零。

70p^69(1-p)^30-p^70*30(1-p)^29=0。

解方程可以得到p=0.7。

在边界点p=0,1,P(Data|M)=0。所以当p=0.7时,P(Data|M)的值最大。这和我们常识中按抽样中的比例来计算的结果是一样的。

假如我们有一组连续变量的采样值(x1,x2,…,xn),我们知道这组数据服从正态分布,标准差已知。请问这个正态分布的期望值为多少时,产生这个已有数据的概率最大?

P(Data | M) = ?

根据公式

可得:

对μ求导可得

,则最大似然估计的结果为μ=(x1+x2+…+xn)/n

由上可知最大似然估计的一般求解过程:

(1) 写出似然函数;

(2) 对似然函数取对数,并整理;

(3) 求导数 ;

(4) 解似然方程

注意:最大似然估计只考虑某个模型能产生某个给定观察序列的概率。而未考虑该模型本身的概率。这点与贝叶斯估计区别。贝叶斯估计方法将在以后的博文中描述

本文参考

http://en.wikipedia.org/wiki/Maximum_likelihood

http://www.shamoxia.com/html/y2010/1520.html

最大似然估计_什么是最大似然估计相关推荐

  1. 二维正态分布的最大似然估计_机器学习系列(二)多元正态分布

    一元正态分布回顾 如果随机变量 服从均值为 方差为 的正态分布 (Univariate normal distribution), ,则其概率密度函数为: 整个分布可以仅用均值及方差来刻画 如果变量之 ...

  2. 最大似然估计_机器学习最大似然估计

    很早之前接触了费曼学习法,费曼学习法是说"如果你真的想确保你对知识的理解没什么问题,就把它教给非专业的另一个人,如果他能听懂,就说明你真的理解了".前段时间看卷积神经网络时特别想用 ...

  3. 最大似然估计_状态估计的基本概念(2)最大似然估计和最大后验估计

    (1)最大似然估计ML和最大后验估计MAP 最大似然估计量 非贝叶斯方法通常是最大化似然函数: 其中 被称为 的最大似然估计量,它是 的函数. 最大后验估计量 估计随机参数的通常方法是最大化后验分布函 ...

  4. 二维正态分布的最大似然估计_最大似然估计-高斯分布

    前言:介绍了最简单的 问题 (这里都是玩具数据,为了方便理解才列出) 0123456789101112 X 1 2 3 4 4.2 4.4 4.6 4.8 5 6 7 8 y 0 0 0 0 1 1 ...

  5. logit模型应用实例_最大似然估计(上)——离散选择模型之十

    原创文章,如需转载请联系作者! 希望这篇文章能讲清楚什么是"最大似然估计". 通过前文的推理,我们已经得到了二项Probit和二项Logit的模型表达式.在二项Probit模型中, ...

  6. 参数估计—最大似然估计与贝叶斯估计

    文章目录 一 参数估计 二 最大似然估计 2.1 参数分量 2.2 基本原理 2.3 高斯情况 2.3.1 协方差矩阵Σ\SigmaΣ已知,而均值μ\muμ未知 2.3.2 协方差矩阵Σ\SigmaΣ ...

  7. 独家 | 一文读懂最大似然估计(附R代码)

    作者:阿尼·辛格 翻译: 陈之炎 校对:丁楠雅 本文约4200字,建议阅读10+分钟. 本文将研究MLE是如何工作的,以及它如何用于确定具有任何分布的模型的系数. 简介 解释模型如何工作是数据科学中最 ...

  8. 专栏 | 机器学习 - 最大似然估计

    https://www.toutiao.com/a6672959716013900301/ 生活实例 我们在生活中就经常应用到最大似然估计的思想. 比如你高中的班主任上课时从教室门缝进行扒头观测,10 ...

  9. 最大似然估计(MLE:样本观测总体参数)是如何工作的?

    1. MLE的意义:样本估计总体分布参数 假定一个事件的观测样本服从如下分布,我们如何确定总体数据的分布模型? 首先应该想到是建立线性回归模型,然而由于该变量不是正态分布的,而且是不对称的,因此不符合 ...

  10. 最小二乘与最大似然估计之间的关系

    最小二乘与最大似然估计之间的关系 1.结论 测量误差(测量)服从高斯分布的情况下, 最小二乘法等价于极大似然估计. 2.最大似然估计概念       最大似然估计就是通过求解最大的(1)式得到可能性最 ...

最新文章

  1. exit的用法python_python 中exit,sys.exit,os._exit用法
  2. OEA 框架演示 - 快过原型的开发
  3. xen虚拟化实战系列(六)之xen虚拟机破解密码
  4. 信号通路怎么读_奥利给~~代谢通路富集分析原来还可以这么做?
  5. java接收的文件转换成临时文件_Android中怎样使用createTempFile实现将字节数据创建到临时文件并转换成FileOutputStream和FileInputStream...
  6. 2021“MINIEYE杯”中国大学生算法设计超级联赛(2)I love exam(背包)
  7. 小米12比我的小米10还便宜
  8. CoInitialize\CoCreateGuid\CoCreateInstance
  9. 美国已批准马斯克的SpaceX发射1.2万颗卫星 以打造一个天基互联网
  10. js如何在当前页面加载springmvc返回的页面_手写SpringMVC学习
  11. 老罗Android开发视频教程( android解析json数据 )4集集合
  12. NetworkManager概述
  13. 又回来了~工作告一段落了,终于有时间看书写点东西了
  14. vue调用接口修改密码_vue开发前后端分离前端如何调用后端接口?
  15. OPPO推送:怎样开通?
  16. Python实现鼠标自动在屏幕上随机移动
  17. 字节码增强之Javassist
  18. linux消息队列优缺点,linux消息队列概念
  19. Linux线程优先级
  20. entity framework 新手入门篇(2)-entity framework基本的增删改查

热门文章

  1. RemotingTooMuchRequestException: sendDefaultImpl call timeout
  2. QueryWrapper、LambdaQueryWrapper以及LambdaQueryChainWrapper用法
  3. mysql表格字放大_删除MySQL表中内容,表大小反而变大了
  4. Centos7重置Mysql 8.0.1 root 密码
  5. Chladni Figure CodeForces - 1162D (暴力,真香啊~)
  6. Eclipse @override报错解决
  7. Compiz Check测试Linux桌面3D兼容性
  8. docker网络、bridge、host
  9. mariadb 存储引擎mysql_MySQL/MariaDB---查询缓存与存储引擎
  10. Java中的类加载器ClassLoader