《贝叶斯统计》最大后验密度可信集(HPD)Python实现

贝叶斯统计 最大后验密度可信集 91,95页

import numpy as np
import sympy as sy
from matplotlib import pyplot as pltplt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = Falsealpha_bar=1-0.05
x=np.arange(80,140,0.01)
x=np.around(x,decimals=2)
#pdf
y=np.around(1/(np.sqrt(2*np.pi)*8.32)*np.exp(-1*((x-110.38)**2)/(2*8.32**2)),decimals=9)#N(110.38,8.32^2)
plt.plot(x,y)#符号变量
sy_x=sy.symbols('x')
sy_y=1/(sy.sqrt(2*sy.pi)*8.32)*sy.exp(-1*((sy_x-110.38)**2)/(2*8.32**2))#N(110.38,8.32^2)#初始赋值
k=np.median(y)
theta_1_k=x[np.argwhere(y==k)[0]]
theta_2_k=x[np.argwhere(y==k)[1]]
s=sy.integrate(sy_y,(sy_x,theta_1_k,theta_2_k))
n=1
while True:if s-alpha_bar>0.0005:k=min(y[y>k])print(n,":",s-alpha_bar)n=n+1elif alpha_bar-s>0.0005:k=max(y[y<k])print(n,":",alpha_bar-s)n=n+1else:breakprint(theta_1_k,theta_2_k)theta_1_k=x[np.argwhere(y==k)[0]]theta_2_k=x[np.argwhere(y==k)[1]]s=sy.integrate(sy_y,(sy_x,theta_1_k,theta_2_k))print(theta_1_k,theta_2_k)

结果：

迭代128次 可信区间[94.1，126.66]

《贝叶斯统计》最大后验密度可信集(HPD)Python实现相关推荐

每个后验密度用计算机画出图形,遥感总复习题库(含答案)
第一章电磁波及遥感物理基础名词解释: 1.电磁波 (变化的电场能够在其周围引起变化的磁场,这一变化的磁场又在较远的区域内引起新的变化电场,并在更远的区域内引起新的变化磁场.) 变化电场和磁场的交替产 ...
MAP 最大后验——利用经验数据获得对未观测量的点态估计
Map (最大后验) 在贝叶斯统计学中,最大后验(Maximum A Posteriori,MAP)估计可以利用经验数据获得对未观测量的点态估计.它与Fisher的最大似然估计(Maximum Lik ...
最大后验（英文为Maximum a posteriori，
在统计学中,最大后验(英文为Maximum a posteriori,缩写为MAP)估计方法根据经验数据获得对难以观察的量的点估计.它与最大似然估计中的 Fisher 方法有密切关系,但是它使用了一个 ...
高密度ARM集群服务器的应用场景
1 什么是高密度ARM集群服务器 ARM集群服务器,是在一个机箱内设计N个高性能ARM处理器,每个ARM处理器配置独立的内存及存储,所有的处理器通过高速网络(千兆或万兆网络)进行相互之间以及对外通信的 ...
我们能从后验分布中学到什么?贝叶斯后验的频率解释
假设我们从未知分布 q 中观察到 N 个独立且同分布的 (iid) 样本 X = (x1, - , xN).统计学中的一个典型问题是"样本集 X 能告诉我们关于分布 q 的什么信息?&quo ...
最小二乘,最大似然,最大后验,查准率与查全率
最小二乘,最大似然,都带有最字,那他们之间有什么关系呢? 一.最大似然定义: 最大似然估计(maximum likelihood estimation, MLE)一种重要而普遍的求估计量的方法. 最 ...
最大后验（Maximum a Posteriori，MAP）概率估计详解
最大后验(Maximum A Posteriori,MAP)概率估计注:阅读本文需要贝叶斯定理与最大似然估计的部分基础最大后验(Maximum A Posteriori,MAP)估计可以利用经 ...
arm集群服务器_高密度ARM集群服务器的应用场景
1 什么是高密度ARM集群服务器 ARM集群服务器,是在一个机箱内设计N个高性能ARM处理器,每个ARM处理器配置独立的内存及存储,所有的处理器通过高速网络(千兆或万兆网络)进行相互之间以及对外通信的 ...
语音后验图特征PPG(Phonetic Posteriorgram)特征简介
最近有个数字人的项目,接触了下后验图PPG这个特征. 简介 PPG的全称是 phonetic posteriorgrams,即语音后验概率,PPG是一个时间对类别的矩阵,其表示对于一个话语的每个特定时 ...
贝叶斯学习--极大后验假设学习
我们假定学习器考虑的是定义在实例空间X上的有限的假设空间H,任务是学习某个目标概念c:X→{0,1}.如通常那样,假定给予学习器某训练样例序列〈〈x1,d1,〉-〈xm,dm〉〉,其中xi为X中的某实 ...

《贝叶斯统计》最大后验密度可信集(HPD)Python实现

《贝叶斯统计》最大后验密度可信集(HPD)Python实现相关推荐

最新文章

热门文章