深入理解交叉熵损失 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

相关文章

深入理解go中interface机制

《深入理解go中interface机制》本文主要介绍了深入理解go中interface机制,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学... 目录前言interface使用类型判断总结前言go的interface是一组method的集合,不

深入解析Java NIO在高并发场景下的性能优化实践指南

《深入解析JavaNIO在高并发场景下的性能优化实践指南》随着互联网业务不断演进,对高并发、低延时网络服务的需求日益增长,本文将深入解析JavaNIO在高并发场景下的性能优化方法,希望对大家有所帮助... 目录简介一、技术背景与应用场景二、核心原理深入分析2.1 Selector多路复用2.2 Buffer

Java Spring的依赖注入理解及@Autowired用法示例详解

《JavaSpring的依赖注入理解及@Autowired用法示例详解》文章介绍了Spring依赖注入(DI)的概念、三种实现方式(构造器、Setter、字段注入),区分了@Autowired(注入... 目录一、什么是依赖注入(DI)?1. 定义2. 举个例子二、依赖注入的几种方式1. 构造器注入(Con

深入理解Go语言中二维切片的使用

《深入理解Go语言中二维切片的使用》本文深入讲解了Go语言中二维切片的概念与应用,用于表示矩阵、表格等二维数据结构,文中通过示例代码介绍的非常详细,需要的朋友们下面随着小编来一起学习学习吧... 目录引言二维切片的基本概念定义创建二维切片二维切片的操作访问元素修改元素遍历二维切片二维切片的动态调整追加行动态

从原理到实战深入理解Java 断言assert

《从原理到实战深入理解Java断言assert》本文深入解析Java断言机制,涵盖语法、工作原理、启用方式及与异常的区别,推荐用于开发阶段的条件检查与状态验证,并强调生产环境应使用参数验证工具类替代... 目录深入理解 Java 断言(assert):从原理到实战引言:为什么需要断言?一、断言基础1.1 语

一文深入详解Python的secrets模块

《一文深入详解Python的secrets模块》在构建涉及用户身份认证、权限管理、加密通信等系统时,开发者最不能忽视的一个问题就是“安全性”,Python在3.6版本中引入了专门面向安全用途的secr... 目录引言一、背景与动机:为什么需要 secrets 模块?二、secrets 模块的核心功能1. 基

Go学习记录之runtime包深入解析

《Go学习记录之runtime包深入解析》Go语言runtime包管理运行时环境,涵盖goroutine调度、内存分配、垃圾回收、类型信息等核心功能,:本文主要介绍Go学习记录之runtime包的... 目录前言:一、runtime包内容学习1、作用:① Goroutine和并发控制:② 垃圾回收:③ 栈和

深入解析 Java Future 类及代码示例

《深入解析JavaFuture类及代码示例》JavaFuture是java.util.concurrent包中用于表示异步计算结果的核心接口,下面给大家介绍JavaFuture类及实例代码,感兴... 目录一、Future 类概述二、核心工作机制代码示例执行流程2. 状态机模型3. 核心方法解析行为总结:三

Android NDK版本迭代与FFmpeg交叉编译完全指南

《AndroidNDK版本迭代与FFmpeg交叉编译完全指南》在Android开发中,使用NDK进行原生代码开发是一项常见需求,特别是当我们需要集成FFmpeg这样的多媒体处理库时,本文将深入分析A... 目录一、android NDK版本迭代分界线二、FFmpeg交叉编译关键注意事项三、完整编译脚本示例四

spring IOC的理解之原理和实现过程

《springIOC的理解之原理和实现过程》:本文主要介绍springIOC的理解之原理和实现过程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、IoC 核心概念二、核心原理1. 容器架构2. 核心组件3. 工作流程三、关键实现机制1. Bean生命周期2.