В этом примере показано, как оценить, имеет ли многомерный временной ряд множественные отношения коинтегрирования с помощью теста Йохансена.
Груз Data_Canada в рабочую область MATLAB ®. Набор данных содержит терминологическую структуру канадских процентных ставок [137]. Извлеките ряды краткосрочных, среднесрочных и долгосрочных процентных ставок.
load Data_Canada Y = Data(:,3:end); % Interest rate data
Для иллюстрации структуры ввода и вывода jcitest при проведении нескольких тестов проверьте коинтеграционный ранг с использованием модели H1 по умолчанию и двух различных структур запаздывания.
[h,pValue,stat,cValue] = jcitest(Y,'Model','H1','Lags',1:2);
************************ Results Summary (Test 1) Data: Y Effective sample size: 39 Model: H1 Lags: 1 Statistic: trace Significance level: 0.05 r h stat cValue pValue eigVal ---------------------------------------- 0 1 35.3442 29.7976 0.0104 0.3979 1 1 15.5568 15.4948 0.0490 0.2757 2 0 2.9796 3.8415 0.0843 0.0736 ************************ Results Summary (Test 2) Data: Y Effective sample size: 38 Model: H1 Lags: 2 Statistic: trace Significance level: 0.05 r h stat cValue pValue eigVal ---------------------------------------- 0 0 25.8188 29.7976 0.1346 0.2839 1 0 13.1267 15.4948 0.1109 0.2377 2 0 2.8108 3.8415 0.0937 0.0713
Тест «трассировки» по умолчанию оценивает нулевые гипотезы ) ранга коинтеграции, меньшего или равного r, по сравнению с альтернативой n), где n - размерность данных. Резюме показывают, что первый тест отклоняет ранг коинтеграции 0 (без коинтеграции) и просто едва отвергает ранг коинтеграции 1, но не отклоняет ранг коинтеграции 2. Вывод состоит в том, что данные показывают 1 или 2 отношения коинтегрирования. С дополнительным запаздыванием в модели, второй тест не может отвергнуть ни один из коинтеграционных рангов, обеспечивая мало в качестве вывода. Важно определить разумную длину задержки для модели VEC (а также общую форму модели) перед тестированием на коинтеграцию.
Поскольку метод Йохансена по своей природе проверяет множественные спецификации ранга для каждой спецификации остальных параметров модели, jcitest возвращает результаты в виде табличных массивов и индексов по нулевому рангу и номеру теста.
Просмотр результатов теста, h.
h
h=2×3 table
r0 r1 r2
_____ _____ _____
t1 true true false
t2 false false false
Заголовки столбцов указывают на тесты r0, r1, и r2соответственно Н ), 1) (2) (3). Заголовки строкt1 и t2 указать два отдельных теста (две отдельные структуры запаздывания), определяемые входными параметрами.
Получите доступ к результату второго теста с нулевым рангом 0 с помощью индексирования табличного массива.
h20 = h.r0(2)
h20 = logical
0