Этот пример показывает, как оценить, имеют ли многомерные временные ряды несколько cointegrating отношений с помощью теста Йохансена.
Загрузите Data_Canada в MATLAB® Workspace. Набор данных содержит термин структура канадских процентных ставок [81]. Извлеките краткосрочный, среднесрочный, и долгосрочный ряд процентной ставки.
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 является размерностью данных. Сводные данные показывают, что первый тест отклоняет ранг коинтеграции 0 (никакая коинтеграция) и едва-едва отклоняет ранг коинтеграции 1, но не удается отклонить ранг коинтеграции 2. Вывод - то, что данные показывают 1 или 2 cointegrating отношения. С дополнительной задержкой в модели второму тесту не удается отклонить любой из рангов коинтеграции, обеспечивая мало посредством вывода. Важно определить разумную продолжительность задержки для модели 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