密度估计公式

  1. 极大似然估计:

y = p ( x 1 , x 2 , x 3 , . . . , x n ) = 1 2 π σ e − ( x 1 − μ ) 2 2 σ 2 1 2 π σ e − ( x 2 − μ ) 2 2 σ 2 . . . 1 2 π σ e − ( x n − μ ) 2 2 σ 2 y = p(x_1,x_2,x_3,...,x_n) = \frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_1-\mu)^2}{2\sigma ^ 2}} \frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_2-\mu)^2}{2\sigma ^ 2}}...\frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_n-\mu)^2}{2\sigma ^ 2}} y=p(x1,x2,x3,...,xn)=2π σ1e−2σ2(x1−μ)22π σ1e−2σ2(x2−μ)2...2π σ1e−2σ2(xn−μ)2

l n y = l n p ( x 1 , x 2 , x 3 , . . . , x n ) = l n ( 1 2 π σ e − ( x 1 − μ ) 2 2 σ 2 1 2 π σ e − ( x 2 − μ ) 2 2 σ 2 . . . 1 2 π σ e − ( x n − μ ) 2 2 σ 2 ) = − n l n ( 2 π σ ) − ∑ i = 1 n ( x i − μ ) 2 2 σ 2 lny = ln p(x_1,x_2,x_3,...,x_n) =ln( \frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_1-\mu)^2}{2\sigma ^ 2}} \frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_2-\mu)^2}{2\sigma ^ 2}}...\frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_n-\mu)^2}{2\sigma ^ 2}} ) =\\ -nln(\sqrt{2 \pi} \sigma ) - \sum_{i=1}^{n}\frac{(x_i-\mu)^2}{2\sigma ^ 2} lny=lnp(x1,x2,x3,...,xn)=ln(2π σ1e−2σ2(x1−μ)22π σ1e−2σ2(x2−μ)2...2π σ1e−2σ2(xn−μ)2)=−nln(2π σ)−i=1∑n2σ2(xi−μ)2

要求y的极限值(将 μ 和 σ \mu 和 \sigma μ和σ视为变量,x视为常量),只需要对上述等式两边对x求导并令导数为0:

∂ ln ⁡ y ∂ μ = − ∑ i = 1 n ( x i − μ ) σ 2 = 0 \frac{\partial \ln y}{\partial \mu} = -\sum_{i=1}^{n} \frac{(x_i - \mu)}{\sigma ^ 2} = 0 ∂μ∂lny=−i=1∑nσ2(xi−μ)=0

即: μ = 1 n ∑ i = 1 n ( x i ) \mu = \frac{1}{n}\sum_{i=1}^{n} (x_i) μ=n1∑i=1n(xi)

∂ ln ⁡ y ∂ σ = − n 1 σ + ∑ i = 1 n ( x i − μ ) 2 σ 3 = 0 \frac{\partial \ln y}{\partial \sigma} =-n\frac{1}{\sigma} +\sum_{i=1}^{n} \frac{(x_i - \mu)^2}{\sigma ^ 3} = 0 ∂σ∂lny=−nσ1+i=1∑nσ3(xi−μ)2=0

即: σ 2 = ∑ i = 1 n ( x i − μ ) 2 n \sigma^2 = \sum_{i=1}^{n} \frac{(x_i - \mu)^2}{n} σ2=∑i=1nn(xi−μ)2

  1. 先验估计:

y = p ( x 1 , x 2 , x 3 , . . . , x n ; θ 0 ) = 1 2 π σ e − ( x 1 − μ ) 2 2 σ 2 1 2 π σ e − ( x 2 − μ ) 2 2 σ 2 . . . 1 2 π σ e − ( x n − μ ) 2 2 σ 2 1 2 π σ 0 e − ( μ 0 − μ ) 2 2 σ 0 2 y = p(x_1,x_2,x_3,...,x_n;\theta_0) = \frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_1-\mu)^2}{2\sigma ^ 2}} \frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_2-\mu)^2}{2\sigma ^ 2}}...\frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_n-\mu)^2}{2\sigma ^ 2}} \frac{1}{\sqrt{2\pi} \sigma_0} e ^{-\frac{(\mu_0-\mu)^2}{2\sigma_0 ^ 2}} y=p(x1,x2,x3,...,xn;θ0)=2π σ1e−2σ2(x1−μ)22π σ1e−2σ2(x2−μ)2...2π σ1e−2σ2(xn−μ)22π σ01e−2σ02(μ0−μ)2

l n y = l n p ( x 1 , x 2 , x 3 , . . . , x n ; θ 0 ) = l n ( 1 2 π σ e − ( x 1 − μ ) 2 2 σ 2 1 2 π σ e − ( x 2 − μ ) 2 2 σ 2 . . . 1 2 π σ e − ( x n − μ ) 2 2 σ 2 1 2 π σ 0 e − ( μ 0 − μ ) 2 2 σ 0 2 ) = − n l n ( 2 π σ ) − ∑ i = 1 n ( x i − μ ) 2 2 σ 2 − l n ( 2 π σ 0 ) − ( μ 0 − μ ) 2 2 σ 0 2 lny = ln p(x_1,x_2,x_3,...,x_n;\theta_0) =ln( \frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_1-\mu)^2}{2\sigma ^ 2}} \frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_2-\mu)^2}{2\sigma ^ 2}}...\frac{1}{\sqrt{2\pi} \sigma} e ^{-\frac{(x_n-\mu)^2}{2\sigma ^ 2}} \frac{1}{\sqrt{2\pi} \sigma_0}e ^{-\frac{(\mu_0-\mu)^2}{2\sigma_0 ^ 2}} )=\\ -nln(\sqrt{2 \pi} \sigma ) - \sum_{i=1}^{n}\frac{(x_i-\mu)^2}{2\sigma ^ 2} -ln(\sqrt{2 \pi} \sigma _0) -\frac{(\mu_0-\mu)^2}{2\sigma_0 ^ 2} lny=lnp(x1,x2,x3,...,xn;θ0)=ln(2π σ1e−2σ2(x1−μ)22π σ1e−2σ2(x2−μ)2...2π σ1e−2σ2(xn−μ)22π σ01e−2σ02(μ0−μ)2)=−nln(2π σ)−i=1∑n2σ2(xi−μ)2−ln(2π σ0)−2σ02(μ0−μ)2

要求y的极限值(将 μ 和 σ \mu 和 \sigma μ和σ视为变量,x视为常量),只需要对上述等式两边对x求导并令导数为0:

∂ ln ⁡ y ∂ μ = − ∑ i = 1 n ( x i − μ ) σ 2 − ( μ 0 − μ ) σ 0 2 = 0 \frac{\partial \ln y}{\partial \mu} = -\sum_{i=1}^{n} \frac{(x_i - \mu)}{\sigma ^ 2} -\frac{(\mu_0- \mu)}{\sigma_0 ^ 2} = 0 ∂μ∂lny=−i=1∑nσ2(xi−μ)−σ02(μ0−μ)=0

这里要注意的是,贝叶斯估计在 θ 0 \theta_0 θ0处的先验概率的计算方式,此时要将 u 0 和 σ 0 当作先验参数 u_0和\sigma_0当作先验参数 u0和σ0当作先验参数

相关推荐
2601_961845157 小时前
新高考一卷真题2025|真题PDF全科整理
线性代数·矩阵·pdf·动态规划·概率论·高考
changjh12 天前
线性回归的似然函数推导
机器学习·线性回归·概率论
AI科技星3 天前
精细结构常数α的多维度物理比值特性及空间螺旋模型研究
人工智能·线性代数·架构·概率论·学习方法
AI科技星3 天前
基于奇合数边界的离散解析数论与双螺旋宇大统一体系(中英文双语纯净终稿)
人工智能·线性代数·架构·概率论·学习方法
chsmiao4 天前
深度学习之概率论
人工智能·深度学习·概率论
词元Max4 天前
3.4 概率论基础:贝叶斯思维与不确定性
概率论
Eloudy4 天前
光子的单缝衍射模型
线性代数·机器学习·概率论
AI科技星5 天前
依托Gε₀ = e²/(4παmₚ²)核心方程:全新公式推导+原创理论提炼+全维度精算验证
人工智能·线性代数·架构·概率论·学习方法
AI科技星8 天前
国家重点研发计划项目申报书
人工智能·线性代数·架构·概率论·学习方法
AI科技星9 天前
基于**v=c(空间光速螺旋运动)唯一第一性原理**重新完整求导证明
人工智能·线性代数·算法·机器学习·架构·概率论·学习方法