随着计算机、互联网等信息技术的发展,马尔可夫链蒙特卡罗(MCMC)模拟技术使贝叶斯统计方法得以应用于许多领域的复杂问题.
本书在介绍常用MCMC 算法的基础上,着重介绍计算贝叶斯后验估计的MCMC 方法和新发展的贝叶斯随机搜索模型选择方法,特别是MCMC 方法在贝叶斯数据分析中的应用. 为了便于读者掌握 MCMC 方法,书中提供了大量的数据分析案例及相应的算法程序、图表和模拟分析结果.
样章试读
目录
- 目录
序
引言1
第1章贝时斯统计中的MCMC方法4
1.1Gibbs算法4
1.1.1基本Gibbs抽样4
1.1.2Gibbs抽样的收敛性10
1.1.3潜变模型的Gibbs抽样12
1.1.4分层模型的Gibbs抽样18
1.2M-H算法21
1.2.1基本M-H算法21
1.2.2M-H算法的收敛性23
1.2.3独立链和随机游走链M-H算法24
1.3混合算法32
1.4逃逸算法34
M可逆跳MCMC算法37
1.6MCMC算法的加速收敛38
1.6.1分组移动和多网格MC抽样39
1.6.2MCMC算法的协方差改进41
1.6.3评分数据例子42
1.6.4其他加速收敛方法45
1.7收敛性诊断49
1.7.1收敛到平稳分布49
1.7.2收敛到后验均值50
第2章贝时斯后验估计54
2.1计算后验均值54
2.1.1基本MC方法54
2.1.2标准误差的估计57
2.1.3MC估计的改进59
2.1.4控制模拟误差63
2.2计算后验众数65
2.2.1矩阵微分65
2.2.2Lindley-Smith最优化66
2.2.3随机近似方法68
2.3估计边缘后验密度71
2.3.1边缘后验密度71
2.3.2核估计方法73
2.3.3重要性加权密度估计74
2.3.4基于K-L散度的有效性80
第3章贝肘斯模型选择91
3.1贝叶斯模型比较91
3.1.1边缘似然91
3.1.2预测密度94
3.2随机搜索模型选择方法95
3.2.1线性回归模型的变量选择96
3.2.2自回归模型的滞后项选择100
3.3贝叶斯模型平均106
3.4模型选择的可逆跳MCMC方法109
3.4.1可逆跳MCMC方法110
3.4.2变阶数时间序列模型选择方法112
第4章多元贝肘斯回归模型118
4.1多元线性回归模型118
4.1.1共郭先验下的贝叶斯推断119
4.1.2广义先验下的贝叶斯推断122
4.1.3实例分析125
4.2向量自回归模型129
4.2.1VAR模型定义129
4.2.2似然函数和先验分布131
4.2.3VAR模型约束的随机搜索133
4.2.4模拟实验136
4.3SUR模型的Gibbs抽样137
4.3.1基本SUR模型138
4.3.2数值模拟139
4.4SUR模型的分层贝叶斯分析140
4.4.1分层模型的Gibbs抽样141
4.4.2具有相关序列误差的SUR模型142
4.4.3具有时变参数的SUR模型145
4.4.4实例分析147
4.5SUR模型的变系数MCMC方法149
4.5.1日UR模型的变换形式149
4.5.2模型参数变换152
4.5.3模拟实验153
4.5.4实例分析157
4.6非参数SUR模型的MCMC方法162
4.6.1基函数展开及分层分析163
4.6.2条件后验分布164
4.6.3MCMC抽样方法166
4.6.4模拟实验168
4.6.5实例分析171
4.7多元空间回归模型172
4.7.1公共分量模型173
4.7.2空间回归模型175
4.7.3先验分布178
4.7.4后验推断179
4.7.5实例分析182
第5章一些非线性时间序列模型185
5.1门限自回归模型185
5.1.1模型参数估计185
5.1.2模型选择的RJMCMC方法186
5.1.3抽样方案187
5.1.4模拟举例188
5.2门限自回归模型的变点分析l90
5.2.1模型的定义190
5.2.2模型参数的贝叶斯推断191
5.2.3模型后验概率的估计194
5.2.4模拟举例195
5.3多变点门限自回归模型198
5.3.1模型的参数估计l98
5.3.2门限估计的MCMC方法203
5.3.3模拟举例205
5.4门限自回归滑动平均模型208
5.4.1模型参数的贝叶斯推断208
5.4.2模型选择的RJMCMC方法210
5.4.3模拟举例210
5.5带有外生变量的门限自回归滑动平均模型212
5.5.1模型参数的贝叶斯推断213
5.5.2模拟举例216
第6章其他贝时斯模型218
6.1空气污染模型218
6.1.1贝叶斯模型分析218
6.1.2实际数据分析222
6.1.3结论229
6.2货币汇率模型230
6.2.1双门限GARCH模型230
6.2.2参数的贝叶斯推断231
6.2.3实证分析233
6.3大脑图像分割模型236
6.3.1功能性磁共振成像236
6.3.2分层混合模型的贝叶斯分析237
6.3.3分层混合模型的RJMCMC方法240
6.3.4大脑FMR1图像分割242
6.4洞流水文数据模型245
6.4.1变点水文时间序列概述245
6.4.2贝叶斯变点时序分析246
6.4.3贝叶斯变点回归分析254
6.4.4结论258
6.5随机波动模型259
6.5.1随机波动模型简介259
6.5.2随机波动模型的贝叶斯拉断260
6.5.3实证分析265
参考文献270
附录A常用统计分布284
A.1一元统计分布284
A.1.1二项分布284
A.1.2贝塔分布284
A.1.3正态分布285
A.1.4伽马分布285
A.1.5逆伽马分布286
A.1.6学生t分布286
A.1.7F分布287
A.2多元统计分布287
A.2.1多元正态分布287
A.2.2多元t分布288
A.3随机矩阵分布288
A.3.1矩阵正态分布288
A.3.2Wish缸t分布290
A.3.3逆Wishart分布290
A.3.4矩阵T分布291
附录B统计软件简介292
B.1R软件介绍292
B.l.l数据操作292
B.1.2统计分析的常见R函数及概率分布函数294
B.1.3可视化295
B.2WinBUGS软件介绍296
B.2.1软件筒介296
B.2.2编程要点297
B.2.3软件操作297
B.2.4举例298
B.3R与Winbugs接口299
B.3.1R2WinBUGS软件包介绍299
B.3.2coda软件包介绍300
B.3.3举例300
附录C4.6.4节模拟实验的R程序302