Этот пример показывает, как оценить, имеет ли многомерные временные ряды несколько коинтегрирующих отношений, используя тест Йохансена.
Загрузка Data_Canada в рабочую область MATLAB ®. Набор данных содержит структуру терминов канадских процентных ставок [137]. Извлечение серии краткосрочных, среднесрочных и долгосрочных процентных ставок.
load Data_Canada Y = Data(:,3:end); % Interest rate data
Чтобы проиллюстрировать вход и структуру output 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
Тест «trace» по умолчанию оценивает нулевые гипотезы коинтеграции меньше или равен r против альтернативы , где 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, соответственно, , , и против . Заголовки строк t1 и t2 указать два отдельных теста (две отдельные структуры задержки), заданные входными параметрами.
Получите доступ к результату второго теста с нулевым рангом использование табличной индексации массива.
h20 = h.r0(2)
h20 = logical
0