使用Engle-Granger测试进行协整测试
此示例显示了如何测试零假设,即组成多元模型的响应序列之间没有协整的关系。
加载data_canada
进入MATLAB®工作区。数据集包含加拿大利率的术语结构[141]。提取短期,中期和长期利率序列。
加载data_canaday = data(:,3:end);%多元响应系列
绘制响应系列。
图图(日期,Y,'行宽',2)Xlabel'年';ylabel'百分';名称=系列(3:end);图例(名称,'地点',,,,'NW') 标题'{\ bf加拿大利率,1954- 1994年}';轴紧的网格on
情节显示evidence of cointegration among the three series, which move together with a mean-reverting spread.
要测试协整,请计算
((T1
) 和
((T2
)Dickey-Fuller统计。egcitest
将测试统计量与Engle-Granger临界值的表值进行比较。
[h,pvalue,stat,cvalue] = egcitest(y,'测试',{'T1',,,,'t2'})
h =1x2逻辑数组0 1
PVALUE =1×20.0526 0.0202
Stat =1×2-3.9321 -25.4538
cvalue =1×2-3.9563 -22.1153
这 测试未能拒绝无协整的零,但几乎没有拒绝p-value only slightly above the default 5% significance level, and a statistic only slightly above the left-tail critical value. The 测试确实拒绝无协整的零。
测试回归y(:,1)
onY(:,2:结束)
(默认情况下)拦截C0
。残差系列是
[y(:,1)y(:,2:end)]*beta
-C0
=y(:,1)
-y(:,2:end)*b
-C0
。
第五个输出参数egcitest
包括回归统计数据,还包含回归系数C0
和b
。
检查回归系数以检查假设的协整载体beta
=[1;-b]
。
[~, ~, ~, ~, reg] = egcitest (Y,'测试',,,,'t2');c0 = reg.coeff(1);b = reg.coeff(2:3);beta = [1; -b];H = GCA;绳索= H.ColorOrder;H.NextPlot ='ReplaceChildren';H.ColorOrder = Circshift(CORD,3);
情节(日期,y*beta-c0,'行宽',2);标题'{\ bf协整关系}';轴紧的;传奇离开;网格on;
正如测试所证实的那样,该组合似乎相对固定。