深入理解交叉熵损失 CrossEntropyLoss - 最大似然估计

2024-06-10 07:44

本文主要是介绍深入理解交叉熵损失 CrossEntropyLoss - 最大似然估计,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

深入理解交叉熵损失 CrossEntropyLoss - 最大似然估计

flyfish

下面有详细的例子和公式的说明。

最大似然估计的概念

最大似然估计是一种统计方法,用来估计模型参数,使得在这些参数下观测到的数据出现的概率(即似然)最大。

具体步骤

  1. 定义似然函数
  • 给定一个参数化的概率模型 P ( X ∣ θ ) P(X|\theta) P(Xθ),其中 θ \theta θ 是模型的参数, X X X 是观测数据。
  • 似然函数 L ( θ ∣ X ) L(\theta|X) L(θX) 表示在参数 θ \theta θ 下,观测数据 X X X 出现的概率。
  1. 计算似然函数
  • 对于独立同分布的数据集 { x 1 , x 2 , … , x n } \{x_1, x_2, \ldots, x_n\} {x1,x2,,xn},似然函数是各数据点概率的乘积:
    L ( θ ∣ X ) = P ( X ∣ θ ) = ∏ i = 1 n P ( x i ∣ θ ) L(\theta | X) = P(X|\theta) = \prod_{i=1}^{n} P(x_i|\theta) L(θX)=P(Xθ)=i=1nP(xiθ)
  1. 取对数得到对数似然函数
  • 为了简化计算,通常取似然函数的对数,即对数似然函数:
    log ⁡ L ( θ ∣ X ) = ∑ i = 1 n log ⁡ P ( x i ∣ θ ) \log L(\theta | X) = \sum_{i=1}^{n} \log P(x_i|\theta) logL(θX)=i=1nlogP(xiθ)
  1. 最大化对数似然函数
  • 找到使对数似然函数最大的参数 θ \theta θ
    θ ^ = arg ⁡ max ⁡ θ log ⁡ L ( θ ∣ X ) \hat{\theta} = \arg\max_{\theta} \log L(\theta | X) θ^=argmaxθlogL(θX)

似然函数的定义

假设我们有一个概率模型 P ( X ∣ θ ) P(X|\theta) P(Xθ),其中 θ \theta θ 是模型的参数, X X X 是观测数据。似然函数 L ( θ ∣ X ) L(\theta | X) L(θX) 表示在参数 θ \theta θ 下,观测数据 X X X 出现的概率。

对于独立同分布的数据

如果我们有独立同分布的数据集 { x 1 , x 2 , … , x n } \{x_1, x_2, \ldots, x_n\} {x1,x2,,xn},似然函数是各数据点概率的乘积:

L ( θ ∣ X ) = P ( X ∣ θ ) = ∏ i = 1 n P ( x i ∣ θ ) L(\theta | X) = P(X|\theta) = \prod_{i=1}^{n} P(x_i|\theta) L(θX)=P(Xθ)=i=1nP(xiθ)

公式拆解
  • L ( θ ∣ X ) L(\theta | X) L(θX):似然函数,表示参数 θ \theta θ 给定的情况下,观测数据 X X X 出现的概率。
  • θ \theta θ:模型参数,我们希望估计的未知量。
  • X X X:观测数据的集合。
  • { x 1 , x 2 , … , x n } \{x_1, x_2, \ldots, x_n\} {x1,x2,,xn}:独立同分布的观测数据点。
  • P ( X ∣ θ ) P(X|\theta) P(Xθ):观测数据 X X X 在参数 θ \theta θ 下的联合概率。
  • ∏ i = 1 n \prod_{i=1}^{n} i=1n:从 1 到 n n n 的乘积符号,表示对所有数据点的概率进行乘积。
  • P ( x i ∣ θ ) P(x_i|\theta) P(xiθ):单个数据点 x i x_i xi 在参数 θ \theta θ 下的概率。

对数似然函数

为了简化计算,通常我们对似然函数取对数,得到对数似然函数:

log ⁡ L ( θ ∣ X ) = ∑ i = 1 n log ⁡ P ( x i ∣ θ ) \log L(\theta | X) = \sum_{i=1}^{n} \log P(x_i|\theta) logL(θX)=i=1nlogP(xiθ)

公式拆解
  • log ⁡ L ( θ ∣ X ) \log L(\theta | X) logL(θX):对数似然函数。
  • ∑ i = 1 n \sum_{i=1}^{n} i=1n:从 1 到 n n n 的求和符号,表示对所有数据点的对数概率求和。
  • log ⁡ P ( x i ∣ θ ) \log P(x_i|\theta) logP(xiθ):单个数据点 x i x_i xi 在参数 θ \theta θ 下的对数概率。

举例说明:投掷硬币

假设我们投掷硬币10次,结果是6次正面朝上,我们希望估计硬币正面朝上的概率 p p p

定义似然函数

对于二项分布,似然函数为:

L ( p ∣ X ) = ( 10 6 ) p 6 ( 1 − p ) 4 L(p | X) = \binom{10}{6} p^6 (1-p)^4 L(pX)=(610)p6(1p)4

公式拆解
  • L ( p ∣ X ) L(p | X) L(pX):似然函数,表示在正面概率为 p p p 的情况下,观测到6次正面和4次反面的概率。
  • ( 10 6 ) \binom{10}{6} (610):组合数,表示从10次投掷中选择6次正面的组合数。
  • p 6 p^6 p6:正面出现6次的概率。
  • ( 1 − p ) 4 (1-p)^4 (1p)4:反面出现4次的概率。
对数似然函数

对似然函数取对数:

log ⁡ L ( p ∣ X ) = log ⁡ ( ( 10 6 ) ) + 6 log ⁡ ( p ) + 4 log ⁡ ( 1 − p ) \log L(p | X) = \log \left( \binom{10}{6} \right) + 6 \log(p) + 4 \log(1-p) logL(pX)=log((610))+6log(p)+4log(1p)

最大化对数似然函数

通过求导数并设为0,可以找到使对数似然函数最大的参数 p p p

d d p log ⁡ L ( p ∣ X ) = 6 p − 4 1 − p = 0 \frac{d}{dp} \log L(p | X) = \frac{6}{p} - \frac{4}{1-p} = 0 dpdlogL(pX)=p61p4=0

解这个方程得到:

6 p = 4 1 − p \frac{6}{p} = \frac{4}{1-p} p6=1p4
6 ( 1 − p ) = 4 p 6(1-p) = 4p 6(1p)=4p
6 − 6 p = 4 p 6 - 6p = 4p 66p=4p
6 = 10 p 6 = 10p 6=10p
p = 6 10 = 0.6 p = \frac{6}{10} = 0.6 p=106=0.6

代码

import numpy as np
from scipy.optimize import minimize# 定义对数似然函数,加入小偏移量避免除零错误
def log_likelihood(p, data, epsilon=1e-10):n = len(data)k = np.sum(data)p = np.clip(p, epsilon, 1 - epsilon)  # 确保 p 在 (epsilon, 1 - epsilon) 之间return -(k * np.log(p) + (n - k) * np.log(1 - p))# 模拟数据:10次投掷,6次正面朝上
data = [1] * 6 + [0] * 4# 最大化对数似然函数
result = minimize(log_likelihood, x0=[0.5], args=(data), bounds=[(0, 1)])
p_hat = result.x[0]
print(f'Estimated probability of heads: {p_hat}')
Estimated probability of heads: 0.5999999961321424

最大化对数似然函数与最小化负对数似然函数在本质上是一样的。它们都是为了找到模型参数,使得观测数据在模型下的概率最大化。让我们详细解释一下这个关系。

对数似然函数

首先,我们有似然函数 L ( θ ∣ X ) L(\theta | X) L(θX),表示在参数 θ \theta θ 下,观测数据 X X X 出现的概率。为了简化计算,通常我们对似然函数取对数,得到对数似然函数:

log ⁡ L ( θ ∣ X ) \log L(\theta | X) logL(θX)

最大化对数似然函数就是找到参数 θ \theta θ,使得 log ⁡ L ( θ ∣ X ) \log L(\theta | X) logL(θX) 最大化:

θ ^ = arg ⁡ max ⁡ θ log ⁡ L ( θ ∣ X ) \hat{\theta} = \arg\max_{\theta} \log L(\theta | X) θ^=argθmaxlogL(θX)

负对数似然函数

负对数似然函数是对数似然函数取负号:

− log ⁡ L ( θ ∣ X ) -\log L(\theta | X) logL(θX)

最小化负对数似然函数就是找到参数 θ \theta θ,使得 − log ⁡ L ( θ ∣ X ) -\log L(\theta | X) logL(θX) 最小化:

θ ^ = arg ⁡ min ⁡ θ − log ⁡ L ( θ ∣ X ) \hat{\theta} = \arg\min_{\theta} -\log L(\theta | X) θ^=argθminlogL(θX)

等价关系

最大化对数似然函数和最小化负对数似然函数在数学上是等价的。因为一个数的负数和这个数的大小关系相反,所以在求极值时:

arg ⁡ max ⁡ θ log ⁡ L ( θ ∣ X ) = arg ⁡ min ⁡ θ − log ⁡ L ( θ ∣ X ) \arg\max_{\theta} \log L(\theta | X) = \arg\min_{\theta} -\log L(\theta | X) argmaxθlogL(θX)=argminθlogL(θX)

例子:投掷硬币

假设我们有10次投掷硬币的结果,6次正面朝上,我们希望估计正面朝上的概率 p p p

  1. 对数似然函数
    log ⁡ L ( p ∣ X ) = log ⁡ ( ( 10 6 ) ) + 6 log ⁡ ( p ) + 4 log ⁡ ( 1 − p ) \log L(p | X) = \log \left( \binom{10}{6} \right) + 6 \log(p) + 4 \log(1-p) logL(pX)=log((610))+6log(p)+4log(1p)

  2. 最大化对数似然函数

p ^ = arg ⁡ max ⁡ p [ log ⁡ ( ( 10 6 ) ) + 6 log ⁡ ( p ) + 4 log ⁡ ( 1 − p ) ] \hat{p} = \arg\max_{p} \left[ \log \left( \binom{10}{6} \right) + 6 \log(p) + 4 \log(1-p) \right] p^=argpmax[log((610))+6log(p)+4log(1p)]

  1. 负对数似然函数
    − log ⁡ L ( p ∣ X ) = − log ⁡ ( ( 10 6 ) ) − 6 log ⁡ ( p ) − 4 log ⁡ ( 1 − p ) -\log L(p | X) = -\log \left( \binom{10}{6} \right) - 6 \log(p) - 4 \log(1-p) logL(pX)=log((610))6log(p)4log(1p)
  2. 最小化负对数似然函数

p ^ = arg ⁡ min ⁡ p [ − log ⁡ ( ( 10 6 ) ) − 6 log ⁡ ( p ) − 4 log ⁡ ( 1 − p ) ] \hat{p} = \arg\min_{p} \left[ -\log \left( \binom{10}{6} \right) - 6 \log(p) - 4 \log(1-p) \right] p^=argpmin[log((610))6log(p)4log(1p)]

这篇关于深入理解交叉熵损失 CrossEntropyLoss - 最大似然估计的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1047562

相关文章

回调的简单理解

之前一直不太明白回调的用法,现在简单的理解下 就按这张slidingmenu来说,主界面为Activity界面,而旁边的菜单为fragment界面。1.现在通过主界面的slidingmenu按钮来点开旁边的菜单功能并且选中”区县“选项(到这里就可以理解为A类调用B类里面的c方法)。2.通过触发“区县”的选项使得主界面跳转到“区县”相关的新闻列表界面中(到这里就可以理解为B类调用A类中的d方法

如何理解redis是单线程的

写在文章开头 在面试时我们经常会问到这样一道题 你刚刚说redis是单线程的,那你能不能告诉我它是如何基于单个线程完成指令接收与连接接入的? 这时候我们经常会得到沉默,所以对于这道题,笔者会直接通过3.0.0源码分析的角度来剖析一下redis单线程的设计与实现。 Hi,我是 sharkChili ,是个不断在硬核技术上作死的 java coder ,是 CSDN的博客专家 ,也是开源

MySQL理解-下载-安装

MySQL理解: mysql:是一种关系型数据库管理系统。 下载: 进入官网MySQLhttps://www.mysql.com/  找到download 滑动到最下方:有一个开源社区版的链接地址: 然后就下载完成了 安装: 双击: 一直next 一直next这一步: 一直next到这里: 等待加载完成: 一直下一步到这里

PyTorch模型_trace实战:深入理解与应用

pytorch使用trace模型 1、使用trace生成torchscript模型2、使用trace的模型预测 1、使用trace生成torchscript模型 def save_trace(model, input, save_path):traced_script_model = torch.jit.trace(model, input)<

isa指针的理解

D3实例isa指向D3类对象。D3类的话isa指向D3元类对象。D3元类保存类中的方法调度列表,包括类方法和对象方法

从《深入设计模式》一书中学到的编程智慧

软件设计原则   优秀设计的特征   在开始学习实际的模式前,让我们来看看软件架构的设计过程,了解一下需要达成目标与需要尽量避免的陷阱。 代码复用 无论是开发何种软件产品,成本和时间都最重要的两个维度。较短的开发时间意味着可比竞争对手更早进入市场; 较低的开发成本意味着能够留出更多营销资金,因此能更广泛地覆盖潜在客户。 代码复用是减少开发成本时最常用的方式之一。其意图

[大师C语言(第三十六篇)]C语言信号处理:深入解析与实战

引言 在计算机科学中,信号是一种软件中断,它允许进程之间或进程与内核之间进行通信。信号处理是操作系统中的一个重要概念,它允许程序对各种事件做出响应,例如用户中断、硬件异常和系统调用。C语言作为一门接近硬件的编程语言,提供了强大的信号处理能力。本文将深入探讨C语言信号处理的技术和方法,帮助读者掌握C语言处理信号的高级技巧。 第一部分:C语言信号处理基础 1.1 信号的概念 在Unix-lik

WeakHashMap深入理解

这一章,我们对WeakHashMap进行学习。 我们先对WeakHashMap有个整体认识,然后再学习它的源码,最后再通过实例来学会使用WeakHashMap。 第1部分 WeakHashMap介绍 第2部分 WeakHashMap数据结构 第3部分 WeakHashMap源码解析(基于JDK1.6.0_45) 第4部分 WeakHashMap遍历方式 第5部分 WeakHashMap示例

netty中常用概念的理解

目录   目录ChannelHandler ChannelHandler功能介绍通过ChannelHandlerAdapter自定义拦截器ChannelHandlerContext接口ChannelPipeline ChannelPipeline介绍ChannelPipeline工作原理ChannelHandler的执行顺序   在《Netty权威指南》(第二版)中,ChannelP

安全科普:理解SSL(https)中的对称加密与非对称加密

今天刚好为站点的后台弄了下https,就来分享我了解的吧。 密码学最早可以追溯到古希腊罗马时代,那时的加密方法很简单:替换字母。 早期的密码学:   古希腊人用一种叫 Scytale 的工具加密。更快的工具是 transposition cipher—:只是把羊皮纸卷在一根圆木上,写下信息,羊皮纸展开后,这些信息就加密完成了。 虽然很容易被解密,但它确实是第一个在现实中应用加密的