1.最大似然估计 (MLE):

 什么是最大似然估计?
     问题:给定一组观察数据还有一个参数待定的模型,如何来估计这个未知参数呢?

观察数据(x1,y1)......(xn,yn)   待定模型参数为θ,模型为f(x;θ)。这时候可以借助观察数据来估计这个θ。这就是最大似然函数估计。

     举个例子:

        假设我们有一个袋子,里面装着白球和黑球,但是我们不知道他们分别有多少个,这时候需要我们估计每次取出一个球是白球的概率是多少?如何估计呢? 可以通过连续有放回的从袋子里面取一百次,看看是白球还是黑球。假设取100次里面 白球占70次,黑球30次。设抽取是白球的概率为P。 那么一百次抽取的总概率为 p(x;p)

         p(x;p)=p(x1,x2.......x100;θ)=p(x1;θ)*p(x2;θ)........p(x100;θ

  =p70*(1-p)30

那么这时候我们希望可以使这个概率最大。

求导:logp(x;p)=logp70*(1-p)30   另导数为0则可以求出p=0.7(同理可以用到连续变量里面,这时候就是概率密度函数的乘积so easy)

       是不是很简单,对!就是这么简单!其实最大似然估计就是在给定一组数据和一个待定参数模型,如何确定这个模型未知参数,使得这个确定后的参数模型产生的已知数据概率最大。当然这里我只是举了一个只有一个未知参数的估计方法,多个未知参数的做法是一样的,就是求似然函数求导取最大值。其实并不是所有似然函数都可以求导的,当似然函数无法求导时就需要根据定义求使得L(θ)最大的θ。

举个例子,以抛筛子为例:

2.最大后验估计 (MAP):

最大后验估计就是在原来的MLE的基础上加上了先验知识:

EM算法 :

相信大家对似然函数已经手到擒来了。那么我们就来看看高深的。

一个概率模型有时候既含有观察变量,有含有隐变量。如果只有观察变量那么我们可以用最大似然法(或者贝叶斯)估计未知参数,但是如果还含有隐变量就不能如此简单解决了。这时候就需要EM算法。

大家可能对这种问题不是很明白,也不太明白隐变量是什么意思。我举个例子(引用统计学习方法的例子):

有3枚硬币分别记为A,B,C,并且出现正面概率分别为p ,q ,k.规则如下:先抛硬币A,如果为正面就选择B,否则选择C,然后再将选择的硬币(B或者C抛),然后观测结果。正面为1 反面为0.独立重复实验10次结果如下:1,1,1,0,0,0,1,1,1,0。我们并不知道抛A硬币时为正面还是方面,只知道最后的结果,问如何估计p,q,k的值?

如果我们知道抛的是哪个硬币就可以使用最大似然估计来估计这些参数,但是我们不知道。因为有p的原因,所以无法估计,这个p就是隐变量

       log(Θ)=Σlogp(x;Θ)=Σlogp(x,p;Θ),Θ就是要求的q,k 待定参数,x为观测数据,因为这个p导致我们无法求解MaxΣlogp(x;Θ)。

还比如说调查 男生 女生身高的问题。身高肯定是服从高斯分布。以往我们可以通过对男生抽样进而求出高斯分布的参数,女生也是,但是如果我们只能知道某个人的高度,却不能知道他是男生或者女生(隐含变量),这时候就无法使用似然函数估计了。这个时候就可以使用EM方法。

分为E和M两步:

E步:

首先通过随机赋值一个我们要求的参数,然后求出另外一个隐含参数的后验概率。这是期望计算过程,我们首先通过随便赋予模型参数的初始值p,q,k,求出各个数据到模型的结果。

M步

用求出来的隐含参数的后验概率进行对传统的似然函数估计,对要求参数进行修正。迭代直到前后两次要求的参数一样为止

其实可以这么简单理解:就是在无监督聚类的时候,我们不知道模型的参数(比如为高斯分布),这时候我们就随便赋值给模型的待定参数(u和ó)。然后我们就可以计算出各个数据分别属于那一类。然后我们用这些分类好的数据重新估计u和ó。

   参考 : http://www.cnblogs.com/jerrylead/archive/2011/04/06/2006936.html

      

转载于:https://www.cnblogs.com/GuoJiaSheng/p/3966812.html

MLE MAP EM相关推荐

  1. 线性回归中 LSE MLE MAP之间的关系

    记目标方程为f(X)=w⊤X{f(X) = w^\top X}f(X)=w⊤X,其中w{w}w为P{P}P维向量,X{X}X为大小为N×P{N\times P}N×P的矩阵 xi{x_i}xi​对应的 ...

  2. MLE, MAP and LSM

    文章目录 1 最大似然估计法 1.1 似然函数 1.2 最大似然估计 1.3 求解步骤 1.4 Demo(高斯函数的均值和方差) 2 最大似然估计和最小二乘法的联系 3 最大后验概率估计 [补充] 最 ...

  3. 频率学派还是贝叶斯学派?聊一聊机器学习中的MLE和MAP

    作者:夏飞 Google | 软件工程师 量子位 已获授权编辑发布 转载请联系原作者 本文作者夏飞,清华大学计算机软件学士,卡内基梅隆大学人工智能硕士,现为谷歌软件工程师. 在这篇文章中,他探讨了机器 ...

  4. (九)GD、L1和L2、MAP、MLE、Lasso

    文章目录 一.GD 二.从LR看L1.L2 2.1 模型复杂度 2.2 正则 三.K折交叉验证 四.MAP与MLE 五.Lasso Regression 一.GD GD: 梯度稳定,所有样本的梯度 S ...

  5. Kaldi 对说话人识别GMM-UBM的MAP 参数更新和对数似然概率解读

    写博客=写日记,为自己记录工作进度和理论知识,如果有恰好路过的大牛经过,可以驻足看看我的理解 本人刚接触说话人识别不到一个月,因工作需求研究了kaldi.大致弄懂了GMM-UBM,正在研究Ivecto ...

  6. mybaits八:select查询返回map集合

    返回一条纪录的map package com.atChina.dao;import java.util.List; import java.util.Map;import org.apache.iba ...

  7. 极大似然估计(MLE)相关总结

    一.极大似然估计概述 极大似然估计是频率学派的进行参数估计的法宝,基于以下两种假设前提: ①某一事件发生是因为该事件发生概率最大. ②事件发生与模型参数θ有关,模型参数θ是一个定值.          ...

  8. Java SortedMap接口总结 有序Map要实现的接口规范

    SortedMap接口源码重点 SortedMap接口继承自Map接口,自然具有Map接口定义的其他方法,Map源码可以看我这篇文章 Map 实现SortedMap接口表明该Map的键是根据排序规则有 ...

  9. 文本语言模型的参数估计-最大似然估计、MAP及贝叶斯估计

    以PLSA和LDA为代表的文本语言模型是当今统计自然语言处理研究的热点问题.这类语言模型一般都是对文本的生成过程提出自己的概率图模型,然后利用观察到的语料数据对模型参数做估计.有了语言模型和相应的模型 ...

  10. 面试AI Lab能力测评

    几年前如果熟练使用TensorFlow,同时掌握基本的AI算法就可以很容易找到一份高薪的工作,但现在不一样了,AI岗位的要求越来越高,对知识的深度也提出了更高的要求. 如果现在一个面试官让你从零推导S ...

最新文章

  1. 浏览Github必备的5款神器级别的Chrome插件
  2. boost::mp11::mp_replace_at相关用法的测试程序
  3. python字符串截取方法_如何使用python语言中的字符串方法截取字符串
  4. 修改器内置脚本编写_Node.js 中实践 Redis Lua 脚本
  5. android半透明闪退,(Android)react-native-splash-screen实践-解决react-native打包好后启动白屏的问题...
  6. [编织消息框架][网络IO模型]AIO
  7. MFC中让自定义的类能响应消息
  8. 跨考中科院计算机 学硕,【20082P】中国科学院大学心理学考研学硕认知神经科学学姐:放弃保研,跨考二战上岸...
  9. 555定时器的工作原理
  10. 2019JAVA面试题附答案-精心整理-很全面面试大大节省时间
  11. linux病毒扫描工具,linux病毒扫描工具ClamAV使用
  12. linux下玩三国志游戏,三国志威力无双手游官网版
  13. 次世代教程——美队盾牌雕刻
  14. connect() failed (111: Connection refused) while connecting to upstream 解决方法
  15. 考研中国科学院大学计算机分数线,中国科学院大学2017年考研分数线
  16. 导弹发射各项参数计算涉及计算机应用,按计算机用途分类
  17. 新浪微博和腾讯微博图标
  18. UVA 12576 - Simply Loopy
  19. IOS证书制作,最简单的在线制作流程
  20. 全局最优和局部最优的理解

热门文章

  1. [30期] 第一个项目总结-—数码黑店在这一年开张
  2. NetApp收购Data Domain 当上冤大头?
  3. List of X$ Tables and how the names are derived
  4. CCDA认证的详细综合叙述
  5. Android中使用ShareSDK功能学习
  6. wxcharts.js结合java后台在小程序中的使用(柱状图,其它同理)
  7. 判断手机上是否安装某个APP(iOS)
  8. 基于 display 和 javaScript 封装一个页面布局小插件
  9. Mysql时间戳转日期-日期转时间戳
  10. OFFICE技术讲座:WP/SS文字渲染不开启的条件