极大似然估计的意义_极大似然估计原理思想

极大似然估计的意义_极大似然估计原理思想极大似然估计是概率论中一个很常用的估计方法,在机器学习中的逻辑回归中就是基于它计算的损失函数,因此还是很有必要复习一下它的相关概念的。背景先来看看几个小例子:猎人师傅和徒弟一同去打猎,遇到一只兔


Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

极大似然估计是概率论中一个很常用的估计方法,在机器学习中的逻辑回归中就是基于它计算的损失函数,因此还是很有必要复习一下它的相关概念的。

背景

先来看看几个小例子:

  • 猎人师傅和徒弟一同去打猎,遇到一只兔子,师傅和徒弟同时放枪,兔子被击中一枪,那么是师傅打中的,还是徒弟打中的?
  • 一个袋子中总共有黑白两种颜色100个球,其中一种颜色90个,随机取出一个球,发现是黑球。那么是黑色球90个?还是白色球90个?

看着两个小故事,不知道有没有发现什么规律…由于师傅的枪法一般都高于徒弟,因此我们猜测兔子是被师傅打中的。随机抽取一个球,是黑色的,说明黑色抽中的概率最大,因此猜测90个的是黑色球。

他们有一个共同点,就是我们的猜测(估计),都是基于一个理论:概率最大的事件,最可能发生

其实我们生活中无时无刻不在使用这种方法,只是不知道它在数学中是如何确定或者推导的。而在数理统计中,它有一个专业的名词:

极大似然估计(maximum likelihood estimation, MLE),通俗的说就是 —— 最像估计法(最可能估计法)

数学过程

极大似然原理与数学表示

官方一点描述上面的过程,即:有n个实验结果,\({ A }_{ i }\)到\({ A }_{ n }\),如果\({ A }_{ j }\)发生了,则意味着\({ A }_{ j }\)发生的概率最大。

即,一次试验就发生的事件,这个事件本身发生概率最大

PS
极大似然估计的意义_极大似然估计原理思想
举个例子,我们在学校衡量学习成绩的标准就是考试成绩,高考更是一考定终身的感觉。高考成绩的好坏,则可以当做一个学生能力的体现,虽然有的人考试紧张考砸了,有的人超常发挥了,但是从概率上来说,高考的成绩基本可以判断这个人的(学习)能力。基于极大似然的解释就是,我们高考的成绩很大程度上反应了平时的学习能力,因此考得好的(当前发生的事件),可以认为是学习好的(所有事件发生概率最大的)。

再抽象一点,如果事件发生是关于 \(\theta\) 参数的,那么一次事件放生时,样本为\({x}_{1},…{x}_{k}\),那么\(\hat { \theta } ({x}_{1},…{x}_{k})\)就是\(\theta\)的估计值。当\(\theta=\hat { \theta } ({x}_{1},…{x}_{k})\)时,当前样本发生的概率最大。

PS
极大似然估计的意义_极大似然估计原理思想
再举个射箭的例子,在《权力的游戏》中有个场景,老徒利死的时候,尸体放在穿上,需要弓箭手在岸边发射火箭引燃。但是当时的艾德慕·徒利公爵射了三箭都没中,布林登·徒利实在看不下去了,通过旗帜判断风向,一箭命中!
因此箭能否射中靶心,不仅跟弓箭手的瞄准能力有关,还跟外界的风向有关系。假设不考虑人的因素,但看风向…同样的瞄准和力度,风太大不行、太小也不行….那我们给风的大小设置一个值为\(\theta\)。假设一名弓箭手射出了三只箭,分别是8环、6环、7环(即\({x}_{1}=8\),\({x}_{2}=6\),\({x}_{3}=7\)),当天风的大小为88。那么我们认为只有\(\theta=88\),发生上面事件的概率最大。

极大似然估计法

如果总体X为离散型

假设分布率为\(P=p(x;\theta )\),x是发生的样本,\(\theta\)是代估计的参数,\(p(x;\theta)\)表示估计参数为\(\theta\)时,发生x的的概率。

那么当我们的样本值为:\({x}_{1},{x}_{2},…,{x}_{n}\)时,

\[L(\theta )=L({ x }_{ 1 },{ x }_{ 2 },…,{ x }_{ n };\theta )=\prod _{ i=1 }^{ n }{ p({ x }_{ i };\theta ) } \]

其中\(L(\theta)\)成为样本的似然函数。

假设

\[L({ x }_{ 1 },{ x }_{ 2 },…,{ x }_{ n };\hat { \theta } )=\underset { \theta \in \Theta }{ max } L({ x }_{ 1 },{ x }_{ 2 },…,{ x }_{ n };\theta ) \]

有 \(\hat{\theta}\) 使得 \(L(\theta)\) 的取值最大,那么 \(\hat {\theta}\)就叫做参数 \(\theta\) 的极大似然估计值。

如果总体X为连续型

基本和上面类似,只是概率密度为\(f(x;\theta)\),替代p。

解法

  1. 构造似然函数\(L(\theta)\)
  2. 取对数:\(lnL(\theta)\)
  3. 求导,计算极值
  4. 解方程,得到\(\theta\)

解释一下,其他的步骤很好理解,第二步取对数是为什么呢?

因为根据前面你的似然函数公式,是一堆的数字相乘,这种算法求导会非常麻烦,而取对数是一种很方便的手段:

  • 由于ln对数属于单调递增函数,因此不会改变极值点
  • 由于对数的计算法则:\(ln{ a }^{ b }=blna\)、\(lnab=lna+lnb\) ,求导就很方便了

例子这里就不举了,感兴趣的话,可以看看参考的第二篇里面有好几个求解极大似然估计的例子。

参考

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/166716.html原文链接:https://javaforall.net

赞 (0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • db2 terminate作用_db2 truncate table immediate

    db2 terminate作用_db2 truncate table immediate SQLSTATE是按类代码进行分组的;对于子代码,请参阅相应的表。表2.SQLSTATE类代码类代码  含义要获得子代码,参阅…00完全成功完成表301警告表402无数据表507动态SQL错误表608连接异常表709触发操作异常表80A功能部件不受支持表90D目标类型规范无效表100F无效标记表11

    2022年10月29日
  • 名词复数变化_名词变复数专项训练题

    名词复数变化_名词变复数专项训练题以o结尾的名词正常情况直接加s;以下名词加es:英雄hero和黑人negro爱吃西红柿potato和土豆tomato。芒果mango加s和es都可以不规则变化child-childrenfoot-feettooth-teethmouse-miceman-men…

    2022年10月30日
  • 有关微信登录小程序接口40029报错

    有关微信登录小程序接口40029报错有关微信登录小程序接口40029报错微信的报错信息errcode:40029,errmsg:"invalidcode,hints:[req_id:HQd79a0747th31]我的可能是例外,但是应该不是个例!微信验签等代码是copy上一个项目的!而且我的项目用上一个的小程序/公众号写过测试,小程序/公众号中做过相应配置!所以以上错误很可能是你的APPID和app_…

    2022年5月5日
  • 实验室设备管理系统c语言程序设计_设备管理系统怎么实现

    实验室设备管理系统c语言程序设计_设备管理系统怎么实现#include#include#includetypedefstructequipment{intnum;//编号charname[20];//名称charmodel[20];//型号charfactory[20];//厂家charbuy_date[20];//购买日期struc

    2022年10月13日
  • [转]Win 10 的 Win 按键没反应

    [转]Win 10 的 Win 按键没反应来自:http://www.pc0359.cn/article/win10/66397.html方法步骤:  1、首先我们需要将任务管理器运行出来。方法有两种,第一种:使用鼠标右键单击屏幕下方的任务栏,在弹出的选框中单击“任务管理器”即可打开;第二种:使用快捷键:Ctrl+Shift+Esc将其调出来。  2、进入任务管理器界面之后,结束桌面窗口管理器任务。  3、会…

    2022年5月12日
  • PyTorch源码解读之torchvision.models「建议收藏」

    PyTorch源码解读之torchvision.models「建议收藏」PyTorch框架中有一个非常重要且好用的包:torchvision,该包主要由3个子包组成,分别是:torchvision.datasets、torchvision.models、torchvision.transforms。这3个子包的具体介绍可以参考官网:http://pytorch.org/docs/master/torchvision/index.html。具体代码可以参考github:

    2022年6月24日

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号