jcitest

Тест коинтеграции Йохансена

Синтаксис

[h,pValue,stat,cValue,mles] = jcitest(Y)
[h,pValue,stat,cValue,mles] = jcitest(Y,Name,Value)

Описание

Тесты Йохансена оценивают нулевую гипотезу, которую H (r) коинтеграции оценивает меньше чем или равный r среди numDims- размерные временные ряды в Y против альтернатив H (numDimsтрассировка протестируйте) или H (r +1) (maxeig тест. Тесты также производят оценки наибольшего правдоподобия параметров в модели векторного исправления ошибок (VEC) cointegrated ряда.

[h,pValue,stat,cValue,mles] = jcitest(Y) выполняет тест коинтеграции Йохансена на матрице данных Y.

[h,pValue,stat,cValue,mles] = jcitest(Y,Name,Value) выполняет тест коинтеграции Йохансена на матрице данных Y с дополнительными опциями, заданными одним или несколькими Name,Value парные аргументы.

Входные параметры

Y

numObs- numDims матрица, представляющая numObs наблюдения за numDims- размерные временные ряды yt, с последним наблюдением новое. Y не может иметь больше чем 12 столбцов. Наблюдения, содержащие NaN значения удалены. Начальные значения для изолированных переменных по оценке модели VEC приняты с начала данных.

Аргументы в виде пар имя-значение

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

'model'

Вектор символов, такой как 'H2', или вектор ячейки векторов символов, задающих форму детерминированных компонентов модели VEC(q) yt:

Δyt=Cyt1+B1Δyt1++BqΔytq+DX+εt

Если r <numDims ранг коинтеграции, затем C = AB, где A является numDims- r матрица скоростей исправления ошибок и B numDims- r матрица базисных векторов для пробела cointegrating отношений. X содержит любые внешние условия, представляющие детерминированные тренды в данных. Для оценки наибольшего правдоподобия это принято, что ε t ~ NID (0, Q), где Q является инновационной ковариационной матрицей.

Значения model рассмотренные Йохансеном [3].

ЗначениеФорма C y t −1 + D X
'H2'

AB ´yt−1. Нет никаких прерываний или трендов в cointegrated ряду и нет никаких детерминированных трендов на уровнях данных.

'H1*'

A (B 'yt−1+c0). Существуют прерывания в cointegrated ряду и нет никаких детерминированных трендов на уровнях данных.

'H1'

A (B 'yt−1+c0) +c1. Существуют прерывания в cointegrated ряду и существуют детерминированные линейные тренды на уровнях данных. Это - значение по умолчанию.

'H*'A (B 'yt−1+c0+d0t) +c1. Существуют прерывания и линейные тренды в cointegrated ряду и существуют детерминированные линейные тренды на уровнях данных.
'H'A (B 'yt−1+c0+d0t) +c1+d1t. Существуют прерывания и линейные тренды в cointegrated ряду и существуют детерминированные квадратичные тренды на уровнях данных.

Детерминированные условия за пределами cointegrating отношений, c 1 и d 1, идентифицированы путем проектирования коэффициентов постоянной и линейной регрессии, соответственно, на ортогональное дополнение A.

'lags'

Скаляр или вектор неотрицательных целых чисел, указывающих на номер q изолированных различий в модели VEC(q) yt.

Отставание и дифференцирование временные ряды уменьшает объем выборки. Отсутствующий любые преддемонстрационные значения, если yt задан для t = 1:N, то изолированная серия y tk задана для t = k + 1:N. Дифференцирование уменьшает основу времени до k +2:N. С изолированными различиями q общей основой времени является q +2:N, и эффективным объемом выборки является T = N − (q +1).

Значение по умолчанию: 0

'test'

Вектор символов, такой как 'trace', или вектор ячейки векторов символов, указывающих на тип теста, который будет выполняться. Значениями является 'trace' или 'maxeig'. Значением по умолчанию является 'trace'. Оба теста оценивают нулевую гипотезу, которую H (r) коинтеграции оценивает меньше чем или равный r. Статистические данные вычисляются с помощью эффективного объема выборки T и упорядоченные оценки собственных значений C = A B ′, λ 1>...> λd, где d = numDims.

  • Когда значением является 'trace', альтернативной гипотезой является H (numDims). Статистические данные:

    T[log(1λr+1)++log(1λnumDims)]

  • Когда значением является 'maxeig', альтернативной гипотезой является H (r +1). Статистические данные:

    Tlog(1λr+1)

'alpha'

Скаляр или вектор номинальных уровней значения для тестов. Значения должны быть между 0,001 и 0.999.

Значение по умолчанию: 0.05

'display'

Вектор символов, такой как 'off', или вектор ячейки векторов символов, указывающих, отобразить ли сводные данные результатов испытаний и параметра, оценивает в Командном окне.

ЗначениеОтображение
'off'Никакое отображение к командному окну. Это - значение по умолчанию если jcitest вызван только одним выходным аргументом (h).
'summary'Отобразите сводные данные результатов испытаний. Пустой указатель оценивает r = 0:numDims − 1 отображены в первом столбце каждых сводных данных. Несколько тестов отображены в отдельных сводных данных. Это - значение по умолчанию если jcitest вызван больше чем одним выходным аргументом (то есть, если pValue вычисляется), и недоступно если jcitest вызван только одним выходным аргументом (h).
'params'Отобразите оценки наибольшего правдоподобия значений параметров, сопоставленных с моделью VEC(q) уменьшаемого ранга yt. Это отображение только доступно если jcitest вызван пятью выходными аргументами (то есть, если mles вычисляется). Отображенные значения параметров возвращены в mles.rnM). paramVals для пустого ранга r = n и протестируйте m.
'full'Отобразите оба summary и params.

Значения векторов символов расширены до продолжительности любого векторного значения (количество тестов). Векторные значения должны иметь равную длину.

Выходные аргументы

h

numTests- numDims табличный массив булевых решений для тестов.

Строки h соответствуйте тестам, заданным входными параметрами, и программное обеспечение помечает строки t1t2 ..., tu, где u = numTests. Переменные h соответствуйте различным, обеспеченным рангам коинтеграции r = 0..., numDims – 1, и программное обеспечение помечает переменные r0R1 ..., rR, где R = numDims – 1. К результатам доступа, сохраненным в h, например, результат для теста m из пустого ранга n, используйте h.rnM).

Значения h равняйтесь 1 TRUE) укажите, что отклонение пустого указателя коинтеграции оценивает r в пользу альтернативы. Значения h равняйтесь 0 ложь) укажите на отказ отклонить пустой указатель.

pValue

numTests- numDims табличный массив вероятностей правильного хвоста тестовой статистики.

Строки pValue соответствуйте тестам, заданным входными параметрами, и программное обеспечение помечает строки t1t2 ..., tu, где u = numTests. Переменные pValue соответствуйте различным, обеспеченным рангам коинтеграции r = 0..., numDims – 1, и программное обеспечение помечает переменные r0R1 ..., rR, где R = numDims – 1. К результатам доступа, сохраненным в pValue, например, результат для теста m из пустого ранга n, используйте pValue.rnM).

stat

numTests- numDims табличный массив тестовой статистики, определенной test аргумент пары "имя-значение".

Строки stat соответствуйте тестам, заданным входными параметрами, и программное обеспечение помечает строки t1t2 ..., tu, где u = numTests. Переменные stat соответствуйте различным, обеспеченным рангам коинтеграции r = 0..., numDims – 1, и программное обеспечение помечает переменные r0R1 ..., rR, где R = numDims – 1. К результатам доступа, сохраненным в stat, например, результат для теста m из пустого ранга n, используйте статистику rnM).

cValue

numTests- numDims табличный массив критических значений для вероятностей правильного хвоста, определенных alpha аргумент пары "имя-значение". jcitest таблицы загрузок критических значений из файла Data_JCITest.mat, затем линейно интерполирует тестовые критические значения из таблиц. Сведенные в таблицу значения были вычислены с помощью методов, описанных в [4].

Строки cValue соответствуйте тестам, заданным входными параметрами, и программное обеспечение помечает строки t1t2 ..., tu, где u = numTests. Переменные cValue соответствуйте различным, обеспеченным рангам коинтеграции r = 0..., numDims – 1, и программное обеспечение помечает переменные r0R1 ..., rR, где R = numDims – 1. К результатам доступа, сохраненным в cValue, например, результат для теста m из пустого ранга n, используйте cValue.rnM).

mles

numTests- numDims табличный массив структур оценок наибольшего правдоподобия сопоставлен с моделью VEC(q) yt. Каждая структура содержит эти поля.

Поле Описание
paramNames

Вектор ячейки названий параметра, формы:

A BB1..., Bqc0d0 C1 d1}

Элементы зависят от значений lags и model.

paramValsСтруктура параметра оценивает с именами полей, соответствующими названиям параметра в paramNames.
res T-by-numDims матрица остаточных значений, где T является эффективным объемом выборки, полученным, подбирая модель VEC(q) yt к входным данным.
EstCovПредполагаемая ковариация Q инновационного процесса εt.
eigValСобственное значение сопоставило с H (r).
eigVecСобственный вектор сопоставлен с собственным значением в eigVal. v собственных векторов нормирован так, чтобы v S11v = 1, где S 11 задан как в [3].
rLL Ограниченная логарифмическая правдоподобность Y под пустым указателем.
uLLНеограниченная логарифмическая правдоподобность Y под альтернативой.

Строки mles соответствуйте тестам, заданным входными параметрами, и программное обеспечение помечает строки t1t2 ..., tu, где u = numTests. Переменные mles соответствуйте различным, обеспеченным рангам коинтеграции r = 0..., numDims – 1, и программное обеспечение помечает переменные r0R1 ..., rR, где R = numDims – 1. К результатам доступа, сохраненным в mles, например, результат для теста m из пустого ранга n, используйте mles.rnM). Можно далее получить доступ к полям структуры с помощью записи через точку, например, введите mles.rnM). paramNames для названий параметра.

Примеры

свернуть все

Загрузите данные по структуре термина процентных ставок в Канаде:

load Data_Canada
Y = Data(:,3:end);
names = series(3:end);
plot(dates,Y)
legend(names,'location','NW')
grid on

Тест для коинтеграции:

[h,pValue,stat,cValue,mles] = jcitest(Y,'model','H1');
************************
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  37.6886   29.7976  0.0050   0.4101  
1  1  16.5770   15.4948  0.0343   0.2842  
2  0  3.2003    3.8415   0.0737   0.0769  
h,pValue
h=1×3 table
           r0       r1       r2  
          _____    _____    _____

    t1    true     true     false

pValue=1×3 table
             r0           r1          r2   
          _________    ________    ________

    t1    0.0050497    0.034294    0.073661

График оценил cointegrating отношения Byt-1+c0:

YLag = Y(2:end,:);
T = size(YLag,1);
B = mles.r2.paramVals.B;
c0 = mles.r2.paramVals.c0;  
plot(dates(2:end),YLag*B+repmat(c0',T,1))
grid on

Алгоритмы

  • Если jcitest сбои, чтобы отклонить пустой указатель коинтеграции оценивают r = 0, вывод - то, что коэффициент исправления ошибок, C является нулем и моделью VEC(q), уменьшает до стандартной модели VAR (q) в первых различиях. Если jcitest отклонения вся коинтеграция оценивают r меньше, чем numDims, вывод - то, что C имеет полный ранг, и yt является стационарным на уровнях.

  • A параметров и B в модели VEC(q) уменьшаемого ранга не однозначно определяются, хотя их продукт C = A B. jcitest построения B = V(: 1:r) использование ортонормированных собственных векторов V возвращенный eig, затем повторно нормирует так, чтобы V'*S11*V = I, как в [3].

  • Чтобы протестировать линейные ограничения на скорости исправления ошибок A и пробел cointegrating отношений, заполненных B, используйте jcontest.

  • Временные ряды в Y может быть стационарным на уровнях или первых различиях (т.е. I (0) или I (1)). Вместо того, чтобы предварительно тестировать ряд на модульные корни (использование, например, adftest, pptest, kpsstest, или lmctest), процедура Йохансена формулирует вопрос в модели. Серия I (0) сопоставлена со стандартным единичным вектором в течение cointegrating отношений, и его присутствие может быть протестировано с помощью jcontest.

  • Преобразовывать VEC (q) параметры модели в mles выведите к VAR (q +1) параметры модели, используйте vec2var.

  • Deterministic cointegration, где cointegrating отношения, возможно, с прерыванием, производят стационарный ряд, является традиционным смыслом коинтеграции, введенной Энглом и Грейнджером [1] (см. egcitest). Stochastic cointegration, где cointegrating отношения производят стационарный трендом ряд (то есть, d0 является ненулевым), расширяет определение коинтеграции, чтобы разместить большее множество экономических рядов.

  • Если тренды высшего порядка на самом деле не присутствуют в данных, модели с меньшим количеством ограничений могут произвести хорошие подгонки в выборке, но плохие прогнозы из выборки.

Ссылки

[1] Энгл, R. F. и К. В. Дж. Грейнджер. "Коинтеграция и Исправление ошибок: Представление, Оценка и Тестирование". Econometrica. v. 55, 1987, стр 251–276.

[2] Гамильтон, J. D. Анализ Временных Рядов. Принстон, NJ: Издательство Принстонского университета, 1994.

[3] Йохансен, S. Основанный на вероятности вывод в векторных авторегрессивных моделях Cointegrated. Оксфорд: Издательство Оксфордского университета, 1995.

[4] Маккиннон, J. G. А. А. Хог и Л. Мичелис. “Числовые Функции распределения Тестов Отношения правдоподобия для Коинтеграции”. Журнал Прикладной Эконометрики. v. 14, 1999, стр 563–577.

[5] Токарь, пополудни “Тестирующий на Коинтеграцию Используя Подход Йохансена: Мы Используем Правильные Критические значения?” Журнал Прикладной Эконометрики. v. 24, 2009, стр 825–831.

Введенный в R2011a

Для просмотра документации необходимо авторизоваться на сайте