格兰杰检验的基本步骤_Toda-Yamamoto 格兰杰因果检验 TY-Granger方法

本文详细介绍了格兰杰因果检验的定义、检验步骤及特殊情况,重点阐述了Toda-Yamamoto方法在处理非平稳数据和协整关系时的适用性,确保了检验的有效性。通过E-Views软件演示了TY-Granger方法的具体操作,包括单位根检验、确定最佳滞后阶数、构建VAR模型和进行Wald检验。同时,提到了Granger非因果关系的测试注意事项,以及Bauer & Maynard提出的surplus-lag检验方法作为进一步的发展。
摘要由CSDN通过智能技术生成

一、格兰杰因果关系定义

对于因变量,找到有助于预测的协变量。"X is said to Granger-causeY if Y can be better predicted using the histories of bothX and Y than it can by using the history of Y alone."

二、格兰杰因果检验

格兰杰因果检验本质是对VAR模型的参数进行线性约束的检验(一般为检验系数是否为0),它使用Wald检验。Wald检验有效,建立在统计量服从渐进卡方分布的假设下。如果该假设被破坏,则Wald检验非有效,格兰杰检验也非有效。

那么在什么条件下假设会被破坏?比如某些变量是非平稳的;出现非线性约束时;预检验技术效力低...

在这之前,我们回到一般步骤:数据预检验+建模格兰杰检验

首先,数据预检验:单位根检验(ADF PP检验),协整检验(Johansen检验)

其次,建模和格兰杰检验:以下有三种情况

第一种情况:变量们都存在单位根(经济数据一般是一阶单位根,或者是在0~1之间的分数积整),且不存在协整关系,那么做一阶差分处理后,差分数据应用VAR建模,这样在对VAR系数进行检验时,传统渐进理论是有效的。

第二种情况:都存在单位根,存在协整关系,在水平数据(没有经过差分的数据)上应用ECM建模,再进行系数检验。

第三种情况:不管变量是否平稳,不管变量间是否存在协整关系,我们可以直接在水平数据上应用Wald检验来检验线性或非线性约束。这就是Toda Yamam

格兰杰因果检验是一种用于分析时间序列数据中因果关系的方法。R语言中可以使用quantile-on-quantile approach进行格兰杰因果检验。该方法基于分位数回归分析,通过对两个时间序列的分位数进行回归,来判断其中一个序列是否对另一个序列有因果影响。以下是一个使用quantile-on-quantile approach进行格兰杰因果检验的例子[^3]: ```R library(quantreg) library(ggplot2) # 生成两个时间序列数据 set.seed(123) x <- arima.sim(model = list(ar = 0.7), n = 100) y <- arima.sim(model = list(ar = 0.5), n = 100) # 将数据分为训练集和测试集 train.x <- x[1:80] test.x <- x[81:100] train.y <- y[1:80] test.y <- y[81:100] # 使用quantile-on-quantile approach进行格兰杰因果检验 g.fit <- rq(test.y ~ test.x, tau = 0.9) g.coef <- coef(g.fit) # 绘制分位数回归图像 ggplot(data.frame(x = test.x, y = test.y), aes(x, y)) + geom_point() + geom_quantreg(mapping = aes(color = "quantile regression"), formula = y ~ x, tau = 0.9, size = 1.2) + geom_smooth(mapping = aes(color = "lowess"), method = "loess", size = 1.2) + geom_abline(mapping = aes(intercept = g.coef, slope = g.coef, color = "quantile-on-quantile"), size = 1.2) + scale_color_manual(name = "Regression", values = c("quantile regression" = "black", "lowess" = "gray40", "quantile-on-quantile" = "blue"), labels = c("Quantile regression", "Lowess", "Quantile-on-quantile")) + theme_minimal() ``` 上述代码中,首先生成了两个时间序列x和y,然后将数据分为训练集和测试集。接着,使用`rq()`函数进行quantile-on-quantile approach格兰杰因果检验,其中`tau`参数为分位数回归的分位数。最后,通过`ggplot2`包绘制了分位数回归图像,并用蓝色直线表示quantile-on-quantile approach的结果。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值