拓端tecdat|R语言用Backfitting MCMC抽样算法进行贝叶斯推理案例
原文链接:http://tecdat.cn/?p=3429
BART是贝叶斯非参数模型,可以使用Backfitting MCMC进行拟合 。
我不使用任何软件包...... MCMC是从头开始实现的。
考虑协变量数据和成果为主题,。在这个 示例中,数据看起来像这样:
我们可能会考虑以下概率模型
基本上我们使用三次多项式对条件均值进行建模。请注意,这是更一般的 模型的特例
在这种情况下和。该模型在参数矢量的每个元素上具有平坦的先验和在方差参数上具有形状和速率的反伽马先验。
每个条件后验都是高斯(因为共轭)。我们可以使用共轭Gibbs或Metropolis从中进行采样。我们也可以将整个参数矢量作为一个块进行采样,但是在这篇文章中我们将坚持反向拟合 - 这本身就是一个Gibbs采样器。我们仍然从其他条件的每个条件的条件后验中进行抽样。然而,我们利用关键的洞察力,每个条件后验取决于其他beta ,仅由残差表示
直观地,是在减去其他项(非)所解释的部分之后的左手平均值的部分。它也是正常分布的,
在正常之前,后验可以通过共轭来计算
Backfitting MCMC如下进行。首先,初始化所有测试版除外。这完全是任意的 - 您可以从任何参数开始。然后,在每个Gibbs迭代中,
- 计算与值在当前迭代。来自后验的样本 以电流抽取为条件。
- 计算与值在当前迭代。注意,使用步骤1中的值。来自后部的样本 。
- 对所有beta参数继续此过程。
- 绘制完所有参数后,进行采样。这个后验是另一个反伽马。
术语反向拟合似乎是合适的,因为在每次迭代中,我们都“退出” 我们想要使用其他测试版进行采样的分布。
为了获得拟合的回归线,我们需要从后验预测分布中进行采样。我们在每个Gibbs迭代中的步骤4之后通过绘制值来执行此操作
上标表示使用来自Gibbs迭代的值的参数。
非常感谢您阅读本文,有任何问题请在下面留言!
拓端tecdat|R语言用Backfitting MCMC抽样算法进行贝叶斯推理案例相关推荐
- R语言用Backfitting MCMC抽样算法进行贝叶斯推理案例
BART是贝叶斯非参数模型,可以使用Backfitting MCMC进行拟合 . 我不使用任何软件包...... MCMC是从头开始实现的. 考虑协变量数据和成果为主题,.在这个玩具示例中,数 ...
- 拓端tecdat|R语言逻辑回归(Logistic回归)模型分类预测病人冠心病风险
最近我们被客户要求撰写关于冠心病风险的研究报告,包括一些图形和统计输出. 相关视频:R语言逻辑回归(Logistic回归)模型分类预测病人冠心病风险 逻辑回归Logistic模型原理和R语言分类预测冠 ...
- 拓端tecdat|R语言用LOESS(局部加权回归)季节趋势分解(STL)进行时间序列异常检测
最近我们被客户要求撰写关于LOESS(局部加权回归)的研究报告,包括一些图形和统计输出. 这篇文章描述了一种对涉及季节性和趋势成分的时间序列的中点进行建模的方法.我们将对一种叫做STL的算法进行研究, ...
- 拓端tecdat|R语言向量误差修正模型 (VECMs)分析长期利率和通胀率影响关系
最近我们被客户要求撰写关于向量误差修正模型的研究报告,包括一些图形和统计输出. 向量自回归模型估计的先决条件之一是被分析的时间序列是平稳的.但是,经济理论认为,经济变量之间在水平上存在着均衡关系,可以 ...
- 拓端tecdat|R语言线性回归和时间序列分析北京房价影响因素可视化案例
最近我们被客户要求撰写关于北京房价影响因素的研究报告,包括一些图形和统计输出. 目的 房价有关的数据可能反映了中国近年来的变化: 人们得到更多的资源(薪水),期望有更好的房子 人口众多 独生子女政策: ...
- 拓端tecdat荣获掘金社区入驻新人奖
2021年7月,由掘金发起了"入驻成长礼"颁奖活动.本次活动邀请到知名开发者.服务机构代表等业界人士. 据了解,掘金社区"新入驻创作者礼"主要对已经积累了一定历 ...
- 拓端tecdat荣获2022年度51CTO博主之星
相信技术,传递价值,这是51CTO每一个技术创作者的动力与信念,2022 年度,拓端tecdat 作为新锐的数据分析咨询公司,在51CTO平台上,不断的输出优质的技术文章,分享前沿创新技术,输出最佳生 ...
- 基于R语言混合效应模型(mixed model)案例研究
全文链接: http://tecdat.cn/?p=2596 在本文中,我们描述了灵活的竞争风险回归模型.回归模型被指定为转移概率,也就是竞争性风险设置中的累积发生率(点击文末"阅读原文&q ...
- R语言诊断试验数据处理与ROC分析实战案例2
R语言诊断试验数据处理与ROC分析实战案例2 目录 R语言诊断试验数据处理与ROC分析实战案例2 #ROC指标 #样例数据
- R语言诊断试验数据处理与ROC分析实战案例1
R语言诊断试验数据处理与ROC分析实战案例1 目录 R语言诊断试验数据处理与ROC分析实战案例1 #ROC指标 #样例数据
最新文章
- vmoptions默认配置_AndroidStudio基本配置
- Problem A: 平面上的点——Point类 (I)
- nofollow标签_如何Nofollow外链
- 虚拟服务器的磁盘,虚拟服务器磁盘扩容
- EventBus设计与实现分析——事件的发布
- 代码注释掉还能执行_日志消息是可执行代码和注释
- IntelliJ IDEA 单行注释调整
- C语言(CED)01背包——动态规划第二题
- (转)Java atomic原子类的使用方法和原理(一)
- Linux-Windows-Mac-RabbitMQ安装教程
- Object-C中的非正式协议与正式协议
- position absolute定位之所属的containing box
- 设A和B是n阶矩阵 ,A,B,A+B均可逆,证明:(A^-1+B^-1)也可逆,并求其逆。
- 【线性代数】20 基变换,基变换公式,坐标变换公式
- 国外一些DICOM资源下载网址
- Java生成文本水印
- html设置word页脚,Word页码从任意页开始如何设置?
- 获取文件图片后缀名和图片转base64
- 【渝粤教育】电大中专计算机网络基础 (2)作业 题库
- 淘宝技术发展(Oracle/支付宝/旺旺)