var对象指定了p阶平稳的多变量向量自回归模型(VAR(p))模型的函数形式并存储了参数值。
varm
对象的关键组成部分 包括时间序列的数量和多元自回归多项式 ( p )的阶数,因为它们完全指定了模型结构。
其他模型组件包括将相同的外生预测变量与每个序列相关联的回归成分,以及常数和时间趋势项。
可下载资源
例子
创建和修改默认模型
创建一个由一个序列组成的零阶 VAR 模型。
![](https://img-blog.csdnimg.cn/e076fee9da83487fac4ccfbb56efca17.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_14,color_FFFFFF,t_70,g_se,x_16)
Mdl
是一个 varm
模型对象。它包含一个序列、一个未知常数和一个未知方差。模型的属性出现在命令行中。
假设您的问题在滞后 1 处有一个自回归系数。要创建这样的模型,请将自回归系数属性 ( AR
) 设置为包含NaN
使用点表示法的值的单元格 。
![](https://img-blog.csdnimg.cn/27aaf1d09d034d82ba9b4a415294b787.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_13,color_FFFFFF,t_70,g_se,x_16)
如果您的问题包含多个序列,则使用不同的语法来创建模型。
为参数估计创建 VAR(4) 模型
为消费者价格指数 (CPI) 和失业率创建 VAR(4) 模型。
声明 CPI和失业率变量。
cpi DCP; ura = aaTeUAE;
创建默认的 VAR(4) 模型。
var(2,4)
![](https://img-blog.csdnimg.cn/91e56d49a2334265a0e922e3da573783.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_19,color_FFFFFF,t_70,g_se,x_16)
Mdl
是一个 varm
模型对象。例如,该 Constant
属性是一个 2×1 的NaN
值向量 。因此,模型常数是要估计的活动模型参数。
通过将Trend
属性设置为NaN,
使用点表示法来 包含未知的线性时间趋势项 。
![](https://img-blog.csdnimg.cn/91c082d37b304ab6a8f7e74819d50e60.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_20,color_FFFFFF,t_70,g_se,x_16)
扩展 NaN
到适当的长度,即一个 2×1 的NaN
值向量 。
指定 VAR 模型的所有参数值
为三个任意序列创建一个 VAR 模型。指定此方程组中的参数值。
![](https://img-blog.csdnimg.cn/8e6c7c60470d40178cb460363a1148ea.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_11,color_FFFFFF,t_70,g_se,x_16)
假设是多元高斯分布,均值为 0,协方差矩阵
![](https://img-blog.csdnimg.cn/c3541a8063f34157bb1b355aa9fc82c9.png)
为参数值创建变量。
使用适当的名称-值对参数创建一个 VAR(1) 模型对象,表示动态方程组。
var('Coan',cAR',i1're,dta,ovaice'Sa)
![](https://img-blog.csdnimg.cn/d012146af0a441bc862bd221e0b4b594.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_20,color_FFFFFF,t_70,g_se,x_16)
Mdl
是一个完全指定的 varm
模型对象。默认情况下, varm
将自回归系数归因于第一个滞后。
您可以使用圆点表示法调整模型属性。例如,考虑另一个 VAR 模型,该模型将自回归系数矩阵归因于 Phi1
第二个滞后项,为第一个滞后系数指定一个零矩阵,并将所有其他项视为等于 Mdl
。创建此 VAR(2) 模型。
M2R= Phi
![](https://img-blog.csdnimg.cn/10dd8eb58db54b36b19c4ebf164c94d0.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_20,color_FFFFFF,t_70,g_se,x_16)
或者,您可以使用varm
与 for 相同的语法 创建另一个模型对象 Mdl
,但另外指定 'Lags',2
.
估计的 VAR (4) 模型
将 VAR(4) 模型拟合到消费者价格指数 (CPI) 和失业率数据。
在不同的图上绘制两个序列。
figure; plot(atal.Te,DaTa.);
![](https://img-blog.csdnimg.cn/7f1344f817ad4ea582fdc9a58894b5f9.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_17,color_FFFFFF,t_70,g_se,x_16)
figure; plot(DaTTie,DatTE);
![](https://img-blog.csdnimg.cn/39033260e9d54cc6a27a7eada293e39d.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_17,color_FFFFFF,t_70,g_se,x_16)
通过将 CPI 转换为一系列增长率来稳定 CPI。通过从失业率序列中删除第一个观测值来同步这两个序列。
随时关注您喜欢的主题
prce2rt(DaTlL);
创建默认的 VAR(4) 模型。
![](https://img-blog.csdnimg.cn/b62ded394c124fd5a4dc668038afabc9.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_20,color_FFFFFF,t_70,g_se,x_16)
Mdl
是一个 var
模型对象。所有包含NaN
值的属性都 对应于给定数据要估计的参数。
使用整个数据集估计模型。
estate(Mdl)
![](https://img-blog.csdnimg.cn/bdd8cc332c434613880ba457c6e4187a.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_17,color_FFFFFF,t_70,g_se,x_16)
EstMdl
是一个估计的 varm
模型对象。它是完全指定的,因为所有参数都有已知值。说明表明自回归多项式是平稳的。
显示估计的汇总统计信息。
summari
![](https://img-blog.csdnimg.cn/20e8dbffcbfd4a4885e290524b2b6f2a.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_20,color_FFFFFF,t_70,g_se,x_16)
创建并估计 CPI 增长率和失业率的 VAR(4) 模型。将最后十个时期视为预测范围。
![](https://img-blog.csdnimg.cn/0a8d8f10ffad4354adfb4094c27fdf6f.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_20,color_FFFFFF,t_70,g_se,x_16)
VAR(4) 模型的预测
cp = pre2rt(ci); EMl = estme(dl,Y(1(end-10),:));
使用估计模型和样本内数据作为样本前观察预测 10 个数据。
freca(Estl);
在单独的图上绘制带有预测值的序列部分。
plot(Tie(ed - 50:ed),ci(nd - 50:ed));
![](https://img-blog.csdnimg.cn/d2c30e08039249e9ad24045c56518fc7.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_17,color_FFFFFF,t_70,g_se,x_16)
plot(Time(nd - 50:ed),ue(ed - 50:ed));
![](https://img-blog.csdnimg.cn/9d968259b73a4365b9062ef9358c8278.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5ouT56uv56CU56m25a6k,size_17,color_FFFFFF,t_70,g_se,x_16)
可下载资源
关于作者
Kaizong Ye是拓端研究室(TRL)的研究员。Kaizong Ye是拓端研究室(TRL)的研究员。在此对他对本文所作的贡献表示诚挚感谢,他在上海财经大学完成了统计学专业的硕士学位,专注人工智能领域。擅长Python.Matlab仿真、视觉处理、神经网络、数据分析。
本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。
非常感谢您阅读本文,如需帮助请联系我们!