楼主: 我的素质低
22663 67

[R] BSTS贝叶斯结构时间序列模型(Bayesian Structural Time Series)   [推广有奖]

学术权威

83%

还不是VIP/贵宾

-

TA的文库  其他...

〖素质文库〗

结构方程模型

考研资料库

威望
8
论坛币
23388 个
通用积分
28302.3504
学术水平
2705 点
热心指数
2881 点
信用等级
2398 点
经验
223563 点
帖子
2977
精华
52
在线时间
2175 小时
注册时间
2012-11-24
最后登录
2024-1-13

一级伯乐勋章 初级学术勋章 初级热心勋章 初级信用勋章 中级热心勋章 中级学术勋章 中级信用勋章 高级学术勋章 高级热心勋章 高级信用勋章 特级学术勋章

相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币

BSTS贝叶斯结构时间序列模型(Bayesian Structural Time Series)





一、引言——“胖”回归下变量选择+即时预测

    人们判断影响变量都是通过自己的主观看法,但是主观看法下,就会出现"fat regression"问题(变量极多,但是观测值不够),下面BSTS方法就是针对这一问题提出使用的。


    当预测变量比观测变量还要多的时候,当在短期时间序列预测过程中,我们的方法是结合三种贝叶斯方法:卡尔曼滤波法、spike-and-slab回归、模型平均。


    Choi and Varian展示了如何应用搜索引擎的数据进行宏观经济指标的预测,这种即时预测的方式,非常得中央银行的喜欢,并进行了很多相关研究。Arola and Galan [2012], McLaren and Shanbhoge [2011], Hellerstein and Middeldorp [2012],Suhoy [2009], Carri`ere-Swallow and Labb′e [2011]. Choi and Varian [2012] 就进行了相关的研究。


二、变量选择问题



    Castle et al. [2009, 2010] 描绘了21中变量选择的技术,时间序列的预测重要归类为四种方法:显著性检验(向前、向后逐步回归法);信息准则(AIC BIC);主成分因子分析模型;lasso 岭回归和其他惩罚性模型。



三、BSTS模型的子模型族




    1、BSTS是一个模型集合体,BSTS结合了三种统计方法进一个集体系统中(两个部分:趋势部分+回归部分)。

    a、用基础的状态空间结构模型做趋势和季节拟合(模型的趋势部分),卡尔曼滤波用来估计先验p的分布;

    b、用spike-slab 回归做变量选择(模型的回归部分);

    c、用贝叶斯模型平均法(Bayesian model averaging)做预测;

    d、估计方法,用MCMC模拟后验分布生成的样本集,来进行估计。


    2、状态空间时间序列与卡尔曼滤波法

    很多下述学者应用了卡尔曼滤波法对时间序列数据进行预测,空间状态结构模型将时间序列进行划分,成水平项,常规趋势,季节效应和误差项。这个模型是一个泛随机化的模型,这个模型加入了季节效应的作为一个回归组成。

    卡尔曼滤波法是很好的估计状态空间模型的方法,他能够计算p(x/x)预测分布滤波与平滑是常见的空间状态模型的计算步骤


    3、spike-and-slab回归部分

    George and McCulloch[1997]) and Madigan and Raftery [1994].论文研究了spike-slab对于变量的选择。

    spike-and-slab在回归估计之前的变量稀疏问题,这个系统从庞大的模型集中寻找有贡献的模型,我们用他来进行即期预测,但这一方法可以用在很多短期预测


    4、贝叶斯模型平均法

    5、模型的估计

    Scott and Varian【2012】,用MCMC技术通过后验分布来生成模拟样本这些技术生成了参数后验分布的样本集,这些样本集可以作为预测用。



    四、结果分析




    1、后验包含分布

    等式8强加给系数为0的变量一个正向概率(与概率密度不同)(给系数施加一个密度约束条件的方程),因此这个等式是全后验分布的一个共同特征,并不是简单的模型中的值。

    这个的全概率无法计算,但是边缘包含概率可以被蒙特卡罗模拟的样本估计出来。


    2、initial claims data模型后验分布

    在最初数据中对状态的贡献。图4.5是每一个时点潜在状态的动态后验分布

    图4代表联合状态,代表没有数据噪声下序列的平稳值。

    图5代表,每一个组成部分地贡献(脉冲图类似)。图5是一个最初变量被趋势、季节和回归部分所解释的多少。趋势部分两次触底,但是季节和回归部分,展示了更显著的变化。

黑色实线是状态均值,1%的中值置信是浅颜色的,一直到99%置信,颜色逐渐变淡。


    3、预测误差图

    可以看到前期预测的效果不是很好,波动幅度很大,所以造成前期预测误差。后期预测误差趋于稳定,直到很小。


    4、y的增量图

    增量显示图,季节,趋势,某两个个体的趋势图下面还有误差的图。

    消费者信心使用趋势预测数据的分解,MAE是平均绝对误差,这个与随机森林变量选择中的(重要性评分)原理相同,越大越好。



(需要Code的,加我好友,私信我,帖子下面的回复不太能够看到...


image-0001.jpg
image-0002.jpg
image-0003.jpg
image-0004.jpg
image-0005.jpg
image-0006.jpg


二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:Structural Bayesian Series 时间序列模型 struct Series 模型

已有 5 人评分经验 论坛币 学术水平 热心指数 信用等级 收起 理由
saihuatuo + 5 + 2 + 2 + 1 精彩帖子
crystal8832 + 30 + 2 + 2 + 2 精彩帖子
xddlovejiao1314 + 100 + 100 + 5 + 5 + 5 精彩帖子
Nicolle + 100 + 1 + 1 + 1 精彩帖子
niuniuyiwan + 60 + 60 + 5 + 5 + 5 精彩帖子

总评分: 经验 + 160  论坛币 + 295  学术水平 + 15  热心指数 + 15  信用等级 + 14   查看全部评分

本帖被以下文库推荐

心晴的时候,雨也是晴;心雨的时候,晴也是雨!
扣扣:407117636,欢迎一块儿吐槽!!
沙发
niuniuyiwan 在职认证  发表于 2015-8-8 22:22:55 |只看作者 |坛友微信交流群

回帖奖励 +1

研究的好深啊,以前只知道Matlab来完成卡尔曼滤波法,好帖,感谢分享。
已有 2 人评分经验 论坛币 收起 理由
xddlovejiao1314 + 10 + 3 鼓励积极发帖讨论
我的素质低 + 10 恩,恩 这个方法我也是这次才认真看...

总评分: 经验 + 10  论坛币 + 13   查看全部评分

使用道具

藤椅
liuzhen123123 发表于 2015-8-8 22:27:56 |只看作者 |坛友微信交流群

回帖奖励 +1

支持!如果再配上code就更好了!
已有 2 人评分经验 论坛币 热心指数 收起 理由
xddlovejiao1314 + 10 + 3 鼓励积极发帖讨论
我的素质低 + 30 + 3 + 5 精彩帖子

总评分: 经验 + 40  论坛币 + 6  热心指数 + 5   查看全部评分

使用道具

板凳
我的素质低 学生认证  发表于 2015-8-8 22:46:43 |只看作者 |坛友微信交流群
liuzhen123123 发表于 2015-8-8 22:27
支持!如果再配上code就更好了!
  有需要可以单独法,code,太长了
已有 1 人评分经验 论坛币 收起 理由
xddlovejiao1314 + 10 + 3 鼓励积极发帖讨论

总评分: 经验 + 10  论坛币 + 3   查看全部评分

使用道具

报纸
liuzhen123123 发表于 2015-8-9 14:55:35 |只看作者 |坛友微信交流群
我的素质低 发表于 2015-8-8 22:46
有需要可以单独法,code,太长了
求code,谢谢,麻烦你了,408813997@qq.com
已有 1 人评分经验 论坛币 收起 理由
xddlovejiao1314 + 10 + 3 鼓励积极发帖讨论

总评分: 经验 + 10  论坛币 + 3   查看全部评分

使用道具

地板
xddlovejiao1314 学生认证  发表于 2015-8-9 22:05:34 |只看作者 |坛友微信交流群

回帖奖励 +1

好贴,谢谢分享。
已有 1 人评分论坛币 热心指数 收起 理由
niuniuyiwan + 10 + 3 精彩帖子

总评分: 论坛币 + 10  热心指数 + 3   查看全部评分

使用道具

7
condmn 发表于 2015-8-11 17:27:02 |只看作者 |坛友微信交流群

回帖奖励 +1

where is the code
已有 1 人评分经验 论坛币 收起 理由
xddlovejiao1314 + 10 + 3 鼓励积极发帖讨论

总评分: 经验 + 10  论坛币 + 3   查看全部评分

使用道具

8
lzt341 发表于 2015-8-12 15:31:33 |只看作者 |坛友微信交流群

回帖奖励 +1

求code   245405108@qq.com
已有 1 人评分经验 论坛币 收起 理由
xddlovejiao1314 + 10 + 3 鼓励积极发帖讨论

总评分: 经验 + 10  论坛币 + 3   查看全部评分

使用道具

9
Alfred_G 学生认证  发表于 2015-9-1 19:10:58 |只看作者 |坛友微信交流群

回帖奖励 +1

太强了。。!
已有 1 人评分经验 论坛币 收起 理由
xddlovejiao1314 + 10 + 3 鼓励积极发帖讨论

总评分: 经验 + 10  论坛币 + 3   查看全部评分

使用道具

10
pilibaofeng 发表于 2015-10-20 23:04:43 |只看作者 |坛友微信交流群

回帖奖励 +1

已有 1 人评分经验 论坛币 收起 理由
xddlovejiao1314 + 10 + 3 鼓励积极发帖讨论

总评分: 经验 + 10  论坛币 + 3   查看全部评分

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-20 05:01