Inquiry on Johansen method (jcitest function)
2 次查看(过去 30 天)
显示 更早的评论
Dear all,
I've got a quick question with regard to the code segment in
>>
load Data_Canada Y = Data(:,3:end); % Interest rate data [~,~,~,~,reg] = egcitest(Y,'test','t2'); c0 = reg.coeff(1); b = reg.coeff(2:3); beta = [1; -b];
[~,~,~,~,mles] = jcitest(Y,'model','H1*'); BJ2 = mles.r2.paramVals.B; c0J2 = mles.r2.paramVals.c0;
% Normalize the 2nd cointegrating relation with respect to % the 1st variable, to make it comparable to Engle-Granger: BJ2n = BJ2(:,2)/BJ2(1,2); c0J2n = c0J2(2)/BJ2(1,2);
% Plot the normalized Johansen cointegrating relation together % with the original Engle-Granger cointegrating relation:
COrd = get(gca,'ColorOrder');
plot(dates,Y*beta-c0,'LineWidth',2,'Color',COrd(4,:)) hold on plot(dates,Y*BJ2n+c0J2n,'--','LineWidth',2,'Color',COrd(5,:)) legend('Engle-Granger OLS','Johansen MLE','Location','NW') title('{\bf Cointegrating Relation}') axis tight grid on hold off
>>
********************** Results Summary (Test 1)
Data: Y Effective sample size: 40 Model: H1* Lags: 0 Statistic: trace Significance level: 0.05
r h stat cValue pValue eigVal ---------------------------------------- 0 1 38.8360 35.1929 0.0194 0.4159 1 0 17.3256 20.2619 0.1211 0.2881 2 0 3.7325 9.1644 0.5229 0.0891
>>
According to the test result, it implies that there is only one cointegrating relation. Yet there exist two columns in BJ2 which means that only one column of two is a candidate for the single cointegrating relation.
At this point, the code simply assumes that it is the second column for the relation, which justifies the following part.
>>
% Normalize the 2nd cointegrating relation with respect to % the 1st variable, to make it comparable to Engle-Granger: BJ2n = BJ2(:,2)/BJ2(1,2); c0J2n = c0J2(2)/BJ2(1,2);
>>
Now my question is how do you know the cointegrating relation is the second column, not the first one?
Please let me know.
Thanks in advance!
0 个评论
回答(1 个)
Hang Qian
2014-4-24
I think the number of cointegrations inferred from the data is suggestive rather than conclusive. If we have some theory that backs up the belief that there are 2 cointegrations among these three variables, then trust that belief. My understanding is that the demo illustrates that the EG cointegration matches one of the two cointegrations (or some linear combinations of the two) estimated by Johansen’s FIML, though I am not sure if we can always find the same on all other datasets.
2 个评论
Hang Qian
2014-5-13
For the first question: no, I do not expect the estimated single integrating relation fits any of the columns, since it could be a linear combination of the multiple cointegrations, if the sample size were infinity. Actually, if we invert the dependent and explanatory variables in the EG regression, we could obtain different estimation on the cointegration coefficients, for a given dataset. I tend to take the estimated cointegration results as a reference, not a golden rule.
For the second question: yes, I think you are right. Each column of the matrix refers to a cointegration relation produced by Johansen’s full information maximum likelihood. The JCITEST is both a hypothesis testing tool and an estimation tool.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Cointegration Analysis 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!