1 / 11
文档名称:

MCMC案例学习.pdf

格式:pdf   页数:11页
下载后只包含 1 个 PDF 格式的文档,没有任何的图纸或源代码,查看文件列表

如果您已付费下载过本站文档,您可以点这里二次下载

分享

预览

MCMC案例学习.pdf

上传人:1954311756 2014/6/10 文件大小:0 KB

下载得到文件列表

MCMC案例学习.pdf

文档介绍

文档介绍:MCMC 案例学习
Charles
这个案例,我们不关心题目的具体意义,重点放在利用贝叶斯的观点来解决问题
时,MCMC 在后续的计算中所发挥的巨大作用。我们知道,贝叶斯的结果往往是一
个后验分布。这个后验分布往往很复杂,我们难以用经典的方法求解其期望与方差等
一系列的数据特征,这时 MCMC 来了,将这一系列问题通过模拟来解决。从这个意
义上说,MCMC 是一种计算手段。依频率学派看来,题目利用广义线性模型可以解
决,在贝叶斯看来同样可以解决,但是遇到了一个问题,就是我们得到的非标准后验
分布很复杂。我们正是利用 MCMC 来解决了这个分布的处理问题。本文的重点也在
于此。
在使用 MCMC 时作者遵循了这样的思路,首先依照贝叶斯解决问题的套路,构
建了非标准后验分布函数。然后初步运行 MCMC,确定合适的 scale。继而,确定适
当的模拟批次和每批长度(以克服模拟取样的相关性)。最后,估计参数并利用 delta
方法估计标准误。
1 问题的提出
这是一个关于 R 软件中 mcmc 包的应用案例。问题出自明尼苏达大学统计系博
士入学考试试题。这个问题所需要的数据存放在 logit 数据集中。在这个数据集中
有五个变量,其中四个自变量 x1、x2、x3、x4,一个响应变量 y。
对于这个问题,频率学派的处理方法是利用广义线性模型进行参数估计,下面是
相应的 R 代码以及结果:
> library(mcmc)
> data(logit)
> out <- glm(y ~ x1 + x2 + x3 + x4, data = logit, family = binomial(), x = T)
> summary(out)
Call:
glm(formula = y ~ x1 + x2 + x3 + x4, family = binomial(), data = logit,
1
1 问题的提出 2
x = T)
Deviance Residuals:
Min 1Q Median 3Q Max
- -
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) *
x1 *
x2 **
x3
x4 .
---
Signif. codes: 0 '***' '**' '*' '.' ' ' 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: on 99 degrees of freedom
Residual deviance: on 95 degrees of freedom
AIC:
Number of Fisher Scoring iterations: 6
但是,使用频率学派的分析方法解决这个问题并不是我们想要的,我们希望使用
Bayesian 分析方法。对于 Bayesian 分析而言,我们假定数据模型(即广义线性模型)
与频率学派一致。同时假定,五个参数(回归系数) 相互独立,并服从均值为 0、标
准差为 2 的先验正态分布。
定义下面的 R 函数来计算非标准的对数后验分布概率密度(先验密度与似然函
数相乘)。我们为什么要定义成密度函数的对数形式?因为虽然我们是从分布中采样,
但是 MCMC 算法的执行函数 metrop() 需要的一个参数正是这个分布的密度函数的
对数形式。
> x <- out$x
> y <- out$y
> lupost <- function(beta, x, y) {
+ eta <- (x %*% beta)
+ logp <- ifelse(eta < 0, eta - log1p(exp(eta)), -log1p(exp(-eta)))
+ logq <- ifelse(eta < 0, -log1p(exp(eta)), -eta - log1p(exp(-eta)))
+ logl <- sum(logp[y == 1]) + sum(logq[y == 0])
2 开始

最近更新

2025年郑州商贸旅游职业学院单招综合素质考试.. 74页

2025年郑州城建职业学院单招职业技能测试题库.. 74页

2025年智能生活与科技实践探索与畅想作文七年.. 9页

2025年郑州智能科技职业学院单招职业适应性测.. 73页

2025年重庆海联职业技术学院单招职业倾向性考.. 76页

2025年郑州轨道工程职业学院单招职业技能测试.. 75页

2025年年终工作总结ppt报告2025 30页

2025年重庆能源职业学院单招职业适应性考试题.. 74页

2025年重庆三峡医药高等专科学校单招职业技能.. 73页

2025年重庆人文科技学院单招职业适应性考试题.. 75页

2025年金山职业技术学院单招职业适应性考试题.. 72页

2025年无线传感器网络中的目标定位技术教程 3页

2025年重庆工商职业学院单招职业适应性考试题.. 75页

2025年无人机遥感技术的应用研究 5页

2025年长春师范高等专科学校单招职业技能测试.. 73页

2025年重庆应用技术职业学院单招职业技能测试.. 74页

2025年大学个人自我鉴定总结最新篇 18页

2025年长沙民政职业技术学院单招职业适应性考.. 75页

2025年长沙电力职业技术学院单招职业适应性考.. 74页

2025年重庆艺术工程职业学院单招职业适应性测.. 75页

2025年闽北职业技术学院单招职业倾向性测试题.. 73页

2025年阜新高等专科学校单招职业倾向性测试题.. 77页

2025年锡林郭勒职业学院单招职业适应性考试题.. 75页

2025年阳光学院单招职业适应性考试题库汇编 75页

菜鸟驿站员工工作合同 5页

学徒工聘用劳动协议汇编3篇 10页

2025年农村别墅可行性方案模板 31页

钢结构厂房施工合同范本 11页

员工入职健康体检表健康证明 2页

基地种植主管岗位职责(共8篇) 14页