0

我正在编写一些代码来确定一组财务和公共利​​益数据的格兰杰因果关系。我在包中的causality()函数语法上遇到了一些问题。VAR这是代码示例及其潜在结果：

data = cbind(x, y, z, price, vol)
data_VAR = VAR(data, type="both", lag.max=30, ic="AIC")
causality(data_VAR, cause="x")$Granger

Granger causality H0: x do not Granger-cause y z price vol

data:  VAR object data_VAR
F-Test = 1.6696, df1 = 120, df2 = 185, p-value = 0.0008476

这将给出与 x 不会导致 、 和 变化的假设y相反z的price结果vol。

如果我想测试x并y作为格兰杰导致其他变量的变量，那么语法是什么？根据我在网上找到的文档，如果愿意的话，可以使用多个变量作为“原因”来运行它，但是根据函数的代码，我似乎无法准确地弄清楚如何读取多个变量。

提前感谢您的帮助！

4

1 回答 1

1

您需要将所有原因放入一个向量中。

> library(vars)
> data(Canada)
> var.2c <- VAR(Canada, p = 2, type = "const")
> causality(var.2c, cause = c("e", "prod"))$Granger

    Granger causality H0: e prod do not Granger-cause rw U

data:  VAR object var.2c
F-Test = 6.8545, df1 = 8, df2 = 292, p-value = 2.919e-08
于 2018-04-10T07:55:37.360 回答