全文链接：http://tecdat.cn/?p=22458最近咱们被客户要求撰写对于动静模型均匀的钻研报告，包含一些图形和统计输入。
本文提供了一个经济案例。着重于原油市场的例子。简要地提供了在经济学中应用模型均匀和贝叶斯办法的论据，应用了动静模型平均法（DMA），并与ARIMA、TVP等办法进行比拟
简介心愿对经济和金融畛域的从业人员和钻研人员有用。
动机事实上，DMA将计量经济学建模的几个特点联合在一起。首先，最终预测是通过模型平均化从几个回归模型中产生的。其次，该办法是贝叶斯办法，也就是说，概率是以置信水平的形式解释的。例如，对工夫t的DMA预测只基于截至工夫t-1的数据。此外，新数据的取得间接导致参数的更新。因而，在DMA中，回归系数和赋予模型的权重都随工夫变动。
贝叶斯办法不是古代计量经济学的支流。然而，这些办法最近正取得越来越多的关注。这其中有各种起因。首先，咱们能够将其与钻研中日益增多的数据量分割起来。因为技术提高，人们通常面临着许多潜在的解释变量的状况。只管大多数变量可能并不重要，但研究者通常不晓得哪些变量应该被剔除。
当然，到某种程度上依然能够应用惯例办法。但因为不足足够的信息，通常无奈对参数进行准确预计。最简略的例子是当解释变量的数量大于工夫序列中的察看值的数量时。例如，即便在线性回归的状况下，规范的一般最小二乘法预计也会呈现一个奇怪矩阵，导致不可能取其倒数。在贝叶斯框架下，依然能够得出一个有意义的公式。贝叶斯办法仿佛也能更好地解决适度参数化和适度拟合问题。
在最近的预测趋势中能够发现各种办法。以原油价格为例，预测办法通常能够分为工夫序列模型、构造模型和其余一些办法，如机器学习、神经网络等。一般来说，工夫序列模型的重点是对稳定的建模，而不是对现货价格的建模。构造模型顾名思义包含因果关系，但它们通常在某些期间有很好的预测能力，而在其余期间则很差。另外，基于小波合成、神经网络等的其余办法通常疏忽了其余因素的影响，只关注繁多工夫序列。这些使得DMA成为从业者的一个乏味的办法。
DMA的下一个方面是，它容许回归系数是随工夫变动的。事实上，在经济呈现迟缓和疾速（结构性中断）变动的状况下，计量经济学模型的这种属性是十分可取的。当然，这样的办法也存在于传统的方法论中，例如，递归或滚动窗口回归。
实践框架咱们将简短地形容fDMA的实践框架。特地是，动静模型平均化（DMA）、动静模型抉择（DMS）、中位概率模型。
动静模型均匀（DMA）DMA在[1]的原始论文中失去了十分具体的介绍。然而，上面是一个简短的阐述，对于了解fDMA中每个函数的作用是必要的。
假如yt是预测的工夫序列（因变量），让x（k）t是第k个回归模型中独立变量的列向量。例如，有10个潜在的原油价格驱动因素。如果它们中的每一个都由一个适合的工夫序列来示意，那么就能够构建2^10个可能的线性回归模型。每个变量都能够包含或不包含在一个模型中。因而，每个变量有两种抉择，形成了2^10种可能性。这包含一个只有常数的模型。因而，一般来说，有潜在的有用的m个独立变量，最多能够构建K=2^m个模型。换句话说，状态空间模型是由以下几个局部组成的
其中k = 1, ... . ，K，t是回归系数的列向量。假如误差遵循正态分布，即e（k）t∼N（0，V（k）t）和（k）t∼N（0，W（k）t）。
在此请留神，有m个潜在的解释变量，2m是构建模型的下限。然而，本文形容的所有办法（如果没有特地阐明的话）都实用于这些2m模型的任何子集，即K≤2m。
动静模型抉择(DMS)动静模型抉择（DMS）是基于雷同的理念，与DMA的理念雷同。惟一的区别是，在DMA中进行的是模型平均化，而在DMS中是模型抉择。换句话说，对于每个期间t，抉择具备最高后验概率的模型。这意味着，只需将公式批改为
其中HT示意k模型。
一个例子：原油市场咱们举一个原油市场的例子。据此能够说，在哪些工夫序列能够作为预测现货原油价格的有用解释变量方面，存在着不确定性。
xts对象crudeoil蕴含来自原油市场的选定数据，即。
-WTI代表WTI（西德克萨斯中质油）现货价格，以每桶计。
MSCI代表MSCI世界指数。TB3MS代表3个月国库券二级市场利率（%）。CSP代表粗钢产量，单位是千吨（能够作为掂量寰球经济流动的一种形式）。TWEXM代表贸易加权的指数（1973年3月=100）。PROD代表原油产品供应量，单位为千桶。CONS代表经合组织的原油产品总消费量。VXO代表规范普尔100指数的隐含稳定率（即股票市场稳定率）。这些数据的频率为每月一次。它们涵盖了1990年1月至2016年12月的期间。
xts对象的趋势蕴含来自谷歌的对于选定搜索词的互联网数量的数据。
stock_markets代表Google Trends的 "股票市场"。interest_rate代表Google Trends的 "利率"。economic_activity示意 "经济流动 "的Google趋势。exchange_rate代表 "汇率 "的谷歌趋势。oil_production示意 "石油生产 "的Google趋势。oil_consumption代表 "石油生产 "的谷歌趋势。market_stress代表Google Trends的 "市场压力"。这些数据也是以月度为频率的。它们涵盖了2004年1月至2016年12月这段时间，因为谷歌趋势没有涵盖更早的期间。从经济角度来看，思考这些工夫序列的对数差分是正当的
R> drivers <- (lag(crudeoil[ , -1], k = 1))[-1, ]R> l.wti <- (diff(log(wti)))[-1, ]R> l.drivers <- (diff(log(driv )))[-1, ]
R> archtest(ld.wti)
R> descstat((ld))
除了PROD的一些问题，所有的工夫序列都能够在5%的显著性程度上被认为是安稳的。对于WTI差分也存在ARCH效应。因而，在DMA中思考指数加权挪动均匀（EWMA）预计方差仿佛是正当的。此外，还能够测试一些忘记因子。依据倡议，对月度工夫序列采取=0.97。所有的方差都小于1。因而，仿佛没有必要对工夫序列进行从新标准化。在DMA的预计中，采取initvar=1仿佛也足够了。
DMA(y = lwti, x = ldrivers,+ alpha = ra, lambda = rl, meth = "ewma" )
依据最小化RMSE，最佳DMA模型是=0.99和=0.97的模型。因而，对这个模型稍作钻研。
plot(x$y, type="l", ylim=c(min(x$y,x$y.hat),max(x$y,x$y.hat)), xlab="", ylab="", main="理论值和预测值", axes = F)比拟图1和图2能够看出，在市场的动荡期间，DMA迅速适应，对有更多变量的模型赋予更高的权重。事实上，这与图3统一。在这一时期，所有解释变量的绝对变量重要性都在回升。咱们还能够看到，自2007年以来，发达的股票市场的作用有所增加。然而，在2013年之后，这种作用变得越来越小；而其余变量的作用开始减少。这一点非常明显，特地是对于汇率。
图3应与图4能够看出。尽管，绝对变量的重要性可能很高，但这个变量的回归系数的预期值可能在0左右。事实上，高的绝对变量重要性同时察看到MSCI、CSP和TWEXM的预期回归系数不为零。所以，这个剖析当初证实了这三个因素在2007年和2013年之间对原油价格起到了重要的预测作用。自2013年以来，股票市场的作用缩小了，被汇率所取代。在2013年前后，最重要的作用是由发达股票市场施展的。
图1
for (i in 1:7) { inc[i+1] <- floor(i * nrow( post.incl)/7) } plot( exp.var, type="l" ylim=c(0,ncol(x$models)) main="变量数量期望值 ", axes = F)
图2
for (i in 1:(ncol( post.incl)-1))plot( post.incl[,i+1], type="l", col=col[i+1], ylim=c(0,1), xlab="", ylab="", main="后蕴含概率", axes = F)
图3
点击题目查阅往期内容
[](http://mp.weixin.qq.com/s?__b...)R语言：EM算法和高斯混合模型的实现
左右滑动查看更多
01
02
03
...