estimate

Подбирайте модель векторного исправления ошибок (VEC) к данным

Описание

пример

EstMdl = estimate(Mdl,Y) возвращает полностью заданную модель VEC (p - 1). Эта модель хранит предполагаемые значения параметров, следующие из подбора кривой модели VEC (p - 1) Mdl к наблюдаемой многомерной серии Y ответа использование наибольшего правдоподобия.

пример

EstMdl = estimate(Mdl,Y,Name,Value) дополнительные опции использования заданы одним или несколькими аргументами пары "имя-значение". Например, 'Model',"H1*",'X',X задает H1* форма Йохансена детерминированных условий и X как внешние данные о предикторе для компонента регрессии.

пример

[EstMdl,EstSE] = estimate(___) возвращает предполагаемые, асимптотические стандартные погрешности предполагаемых параметров с помощью любого входного параметра в предыдущих синтаксисах.

пример

[EstMdl,EstSE,logL,E] = estimate(___) возвращает оптимизированное значение целевой функции логарифмической правдоподобности (logL) и многомерные остаточные значения (E).

Примеры

свернуть все

Рассмотрите модель VEC для следующих семи макроэкономических рядов, и затем подбирайте модель к данным.

  • Валовой внутренний продукт (ВВП)

  • GDP неявный ценовой дефлятор

  • Заплаченная компенсация сотрудников

  • Несельскохозяйственные часы делового сектора всех людей

  • Эффективная ставка по федеральным фондам

  • Частные потребительские расходы

  • Грубые частные внутренние инвестиции

Предположим, что cointegrating ранг 4 и один срок короткого промежутка времени является соответствующим, то есть, рассмотрите модель VEC(1).

Загрузите Data_USEconVECModel набор данных.

load Data_USEconVECModel

Для получения дополнительной информации о наборе данных и переменных, введите Description в командной строке.

Определите, должны ли данные быть предварительно обработаны путем графического вывода ряда на отдельных графиках.

figure;
subplot(2,2,1)
plot(FRED.Time,FRED.GDP);
title('Gross Domestic Product');
ylabel('Index');
xlabel('Date');
subplot(2,2,2)
plot(FRED.Time,FRED.GDPDEF);
title('GDP Deflator');
ylabel('Index');
xlabel('Date');
subplot(2,2,3)
plot(FRED.Time,FRED.COE);
title('Paid Compensation of Employees');
ylabel('Billions of $');
xlabel('Date');
subplot(2,2,4)
plot(FRED.Time,FRED.HOANBS);
title('Nonfarm Business Sector Hours');
ylabel('Index');
xlabel('Date');

figure;
subplot(2,2,1)
plot(FRED.Time,FRED.FEDFUNDS);
title('Federal Funds Rate');
ylabel('Percent');
xlabel('Date');
subplot(2,2,2)
plot(FRED.Time,FRED.PCEC);
title('Consumption Expenditures');
ylabel('Billions of $');
xlabel('Date');
subplot(2,2,3)
plot(FRED.Time,FRED.GPDI);
title('Gross Private Domestic Investment');
ylabel('Billions of $');
xlabel('Date');

Стабилизируйте весь ряд, кроме ставки по федеральным фондам, путем применяния логарифмического преобразования. Масштабируйте получившийся ряд 100 так, чтобы все ряды были по той же шкале.

FRED.GDP = 100*log(FRED.GDP);      
FRED.GDPDEF = 100*log(FRED.GDPDEF);
FRED.COE = 100*log(FRED.COE);       
FRED.HOANBS = 100*log(FRED.HOANBS); 
FRED.PCEC = 100*log(FRED.PCEC);     
FRED.GPDI = 100*log(FRED.GPDI);

Создайте модель VEC(1) с помощью краткого синтаксиса. Задайте имена переменных.

Mdl = vecm(7,4,1);
Mdl.SeriesNames = FRED.Properties.VariableNames
Mdl = 
  vecm with properties:

             Description: "7-Dimensional Rank = 4 VEC(1) Model with Linear Time Trend"
             SeriesNames: "GDP"  "GDPDEF"  "COE"  ... and 4 more
               NumSeries: 7
                    Rank: 4
                       P: 2
                Constant: [7×1 vector of NaNs]
              Adjustment: [7×4 matrix of NaNs]
           Cointegration: [7×4 matrix of NaNs]
                  Impact: [7×7 matrix of NaNs]
   CointegrationConstant: [4×1 vector of NaNs]
      CointegrationTrend: [4×1 vector of NaNs]
                ShortRun: {7×7 matrix of NaNs} at lag [1]
                   Trend: [7×1 vector of NaNs]
                    Beta: [7×0 matrix]
              Covariance: [7×7 matrix of NaNs]

Mdl vecm объект модели. Все свойства, содержащие NaN значения соответствуют параметрам, чтобы быть оцененными определенными данными.

Оцените модель с помощью целого набора данных и опций по умолчанию.

EstMdl = estimate(Mdl,FRED.Variables)
EstMdl = 
  vecm with properties:

             Description: "7-Dimensional Rank = 4 VEC(1) Model"
             SeriesNames: "GDP"  "GDPDEF"  "COE"  ... and 4 more
               NumSeries: 7
                    Rank: 4
                       P: 2
                Constant: [14.1329 8.77841 -7.20359 ... and 4 more]'
              Adjustment: [7×4 matrix]
           Cointegration: [7×4 matrix]
                  Impact: [7×7 matrix]
   CointegrationConstant: [-28.6082 109.555 -77.0912 ... and 1 more]'
      CointegrationTrend: [4×1 vector of zeros]
                ShortRun: {7×7 matrix} at lag [1]
                   Trend: [7×1 vector of zeros]
                    Beta: [7×0 matrix]
              Covariance: [7×7 matrix]

EstMdl предполагаемый vecm объект модели. Это полностью задано, потому что все параметры знали значения. По умолчанию, estimate налагает ограничения формы модели H1 Йохансен VEC путем удаления cointegrating тренда и линейных условий тренда из модели. Исключение параметра из оценки эквивалентно наложению ограничений равенства, чтобы обнулить.

Отобразите краткое изложение от оценки.

results = summarize(EstMdl)
results = struct with fields:
               Description: "7-Dimensional Rank = 4 VEC(1) Model"
                     Model: "H1"
                SampleSize: 238
    NumEstimatedParameters: 112
             LogLikelihood: -1.4939e+03
                       AIC: 3.2118e+03
                       BIC: 3.6007e+03
                     Table: [133x4 table]
                Covariance: [7x7 double]
               Correlation: [7x7 double]

Table поле results таблица оценок параметра и соответствующей статистики.

Считайте модель и данные в Оценке Моделью VEC, и предположите, что выборка оценки запускается в Q1 1 980.

Загрузите Data_USEconVECModel набор данных и предварительно обрабатывает данные.

load Data_USEconVECModel
FRED.GDP = 100*log(FRED.GDP);      
FRED.GDPDEF = 100*log(FRED.GDPDEF);
FRED.COE = 100*log(FRED.COE);       
FRED.HOANBS = 100*log(FRED.HOANBS); 
FRED.PCEC = 100*log(FRED.PCEC);     
FRED.GPDI = 100*log(FRED.GPDI);

Идентифицируйте индекс, соответствующий запуску выборки оценки.

estIdx = FRED.Time(2:end) > '1979-12-31';

Создайте модель VEC(1) по умолчанию с помощью краткого синтаксиса. Примите, что соответствующий ранг коинтеграции равняется 4. Задайте имена переменных.

Mdl = vecm(7,4,1);
Mdl.SeriesNames = FRED.Properties.VariableNames;

Оцените модель с помощью выборки оценки. Задайте все наблюдения перед выборкой оценки как преддемонстрационные данные. Кроме того, задайте оценку формы Х Йохансена модели VEC, которая включает все детерминированные параметры.

Y0 = FRED{~estIdx,:};
EstMdl = estimate(Mdl,FRED{estIdx,:},'Y0',Y0,'Model',"H")
EstMdl = 
  vecm with properties:

             Description: "7-Dimensional Rank = 4 VEC(1) Model with Linear Time Trend"
             SeriesNames: "GDP"  "GDPDEF"  "COE"  ... and 4 more
               NumSeries: 7
                    Rank: 4
                       P: 2
                Constant: [17.5698 3.74759 -20.1998 ... and 4 more]'
              Adjustment: [7×4 matrix]
           Cointegration: [7×4 matrix]
                  Impact: [7×7 matrix]
   CointegrationConstant: [85.4825 -57.3569 -81.7344 ... and 1 more]'
      CointegrationTrend: [-0.0264185 -0.00275396 -0.0249583 ... and 1 more]'
                ShortRun: {7×7 matrix} at lag [1]
                   Trend: [0.000514564 -0.000291183 0.00179965 ... and 4 more]'
                    Beta: [7×0 matrix]
              Covariance: [7×7 matrix]

Поскольку порядок p модели VEC равняется 2, estimate использование только последние два наблюдения (строки) в Y0 как предварительная выборка.

Считайте модель и данные в Оценке Моделью VEC.

Загрузите Data_USEconVECModel набор данных и предварительно обрабатывает данные.

load Data_USEconVECModel
FRED.GDP = 100*log(FRED.GDP);      
FRED.GDPDEF = 100*log(FRED.GDPDEF);
FRED.COE = 100*log(FRED.COE);       
FRED.HOANBS = 100*log(FRED.HOANBS); 
FRED.PCEC = 100*log(FRED.PCEC);     
FRED.GPDI = 100*log(FRED.GPDI);

Data_Recessions набор данных содержит начало и окончание последовательных дат рецессий. Загрузите этот набор данных. Преобразуйте матрицу порядковых номеров даты к массиву datetime.

load Data_Recessions
dtrec = datetime(Recessions,'ConvertFrom','datenum');

Создайте фиктивную переменную, которая идентифицирует периоды, в которые США были в рецессии или хуже. А именно, переменной должен быть 1 если FRED.Time происходит во время рецессии и 0 в противном случае.

isin = @(x)(any(dtrec(:,1) <= x & x <= dtrec(:,2)));
isrecession = double(arrayfun(isin,FRED.Time));

Создайте модель VEC(1) с помощью краткого синтаксиса. Примите, что соответствующий ранг коинтеграции равняется 4. Вы не должны задавать присутствие компонента регрессии при создании модели. Задайте имена переменных.

Mdl = vecm(7,4,1);
Mdl.SeriesNames = FRED.Properties.VariableNames;

Оцените модель с помощью целой выборки. Задайте предиктор, идентифицирующий, было ли наблюдение измерено во время рецессии. Возвратите стандартные погрешности.

[EstMdl,EstSE] = estimate(Mdl,FRED.Variables,'X',isrecession);

Отобразите коэффициент регрессии для каждого уравнения и соответствующих стандартных погрешностей.

EstMdl.Beta
ans = 7×1

   -1.1975
   -0.0187
   -0.7530
   -0.7094
   -0.5932
   -0.6835
   -4.4839

EstSE.Beta
ans = 7×1

    0.1547
    0.0581
    0.1507
    0.1278
    0.2471
    0.1311
    0.7150

EstMdl.Beta и EstSE.Beta 7 1 векторы. Строки соответствуют переменным отклика в EstMdl.SeriesNames и столбцы соответствуют предикторам.

Чтобы проверять, являются ли эффекты рецессий значительными, получите итоговую статистику из summarize, и затем отобразите результаты для Beta.

results = summarize(EstMdl);
isbeta = contains(results.Table.Properties.RowNames,'Beta');
betaresults = results.Table(isbeta,:)
betaresults=7×4 table
                   Value      StandardError    TStatistic      PValue  
                 _________    _____________    __________    __________

    Beta(1,1)      -1.1975       0.15469         -7.7411     9.8569e-15
    Beta(2,1)    -0.018738       0.05806        -0.32273         0.7469
    Beta(3,1)     -0.75305       0.15071         -4.9966     5.8341e-07
    Beta(4,1)     -0.70936       0.12776         -5.5521     2.8221e-08
    Beta(5,1)      -0.5932       0.24712         -2.4004       0.016377
    Beta(6,1)     -0.68353       0.13107         -5.2151      1.837e-07
    Beta(7,1)      -4.4839         0.715         -6.2712     3.5822e-10

whichsig = EstMdl.SeriesNames(betaresults.PValue < 0.05)
whichsig = 1x6 string
    "GDP"    "COE"    "HOANBS"    "FEDFUNDS"    "PCEC"    "GPDI"

Весь ряд кроме GDPDEF кажись, оказывать значительное влияние рецессий.

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

свернуть все

Модель VEC, содержащая неизвестные значения параметров в виде vecm объект модели возвращен vecm.

NaN- ценные элементы в свойствах указывают на неизвестные, допускающие оценку параметры. Указанные элементы указывают на ограничения равенства на параметры по оценке модели. Инновационная ковариационная матрица Mdl.Covariance не может содержать соединение NaN значения и вещественные числа; необходимо полностью задать ковариацию, или это должно быть абсолютно неизвестно (NaN(Mdl.NumSeries)).

Наблюдаемый многомерный ряд ответа, к который estimate подбирает модель в виде numobs- numseries числовая матрица.

numobs объем выборки. numseries количество переменных отклика (Mdl.NumSeries).

Строки соответствуют наблюдениям, и последняя строка содержит последнее наблюдение.

Столбцы соответствуют отдельным переменным отклика.

Y представляет продолжение преддемонстрационного ряда ответа в Y0.

Типы данных: double

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

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

Пример: 'Model',"H1*",'Y0',Y0,'X',X включает только детерминированные условия в форму H1 модели VEC, использует матричный Y0 как преддемонстрационные ответы, требуемые для оценки, и, включает компонент линейной регрессии, состоявший из данных о предикторе в X.

Преддемонстрационные ответы, чтобы инициировать оценку модели в виде разделенной запятой пары, состоящей из 'Y0' и numpreobs- numseries числовая матрица.

numpreobs количество преддемонстрационных наблюдений.

Строки соответствуют преддемонстрационным наблюдениям, и последняя строка содержит последнее наблюдение. Y0 должен иметь, по крайней мере, Mdl.P 'Строки' . Если вы предоставляете больше строк, чем необходимый, estimate использует последний Mdl.P наблюдения только.

Столбцы должны соответствовать ряду ответа в Y.

По умолчанию, estimate использование Y(1:Mdl.P,:) как преддемонстрационные наблюдения, и затем подбирает модель к Y((Mdl.P + 1):end,:). Это действие уменьшает эффективный объем выборки.

Типы данных: double

Данные о предикторе для компонента регрессии в модели в виде разделенной запятой пары, состоящей из 'X' и числовая матрица, содержащая numpreds столбцы.

numpreds количество переменных предикторов.

Строки соответствуют наблюдениям, и последняя строка содержит последнее наблюдение. estimate не использует компонент регрессии в преддемонстрационный период. X должен иметь, по крайней мере, столько наблюдений, сколько используются после преддемонстрационного периода.

  • Если вы задаете Y0, затем X должен иметь, по крайней мере, numobs строки (см. Y).

  • В противном случае, X должен иметь, по крайней мере, numobsMdl.P наблюдения с учетом преддемонстрационного удаления.

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

Столбцы соответствуют отдельным переменным предикторам. Все переменные предикторы присутствуют в компоненте регрессии каждого уравнения ответа.

По умолчанию, estimate исключает компонент регрессии, независимо от его присутствия в Mdl.

Типы данных: double

Форма Йохансена модели VEC (p - 1) детерминированные условия [2] в виде разделенной запятой пары, состоящей из 'Model' и значение в этой таблице (для определений переменной, см. Векторную Модель Исправления ошибок).

ЗначениеСрок исправления ошибокОписание
"H2"

AB 'yt − 1

Никакие прерывания или тренды не присутствуют в cointegrating отношениях, и никакие детерминированные тренды не присутствуют на уровнях данных.

Задайте эту модель только, когда все ряды ответа будут иметь среднее значение нуля.

"H1*"

A (B 'yt−1+c0)

Прерывания присутствуют в cointegrating отношениях, и никакие детерминированные тренды не присутствуют на уровнях данных.

"H1"

A (B 'yt−1+c0) +c1

Прерывания присутствуют в cointegrating отношениях, и детерминированные линейные тренды присутствуют на уровнях данных.

"H*"A (B 'yt−1+c0+d0t) +c1

Прерывания и линейные тренды присутствуют в cointegrating отношениях, и детерминированные линейные тренды присутствуют на уровнях данных.

"H"A (B 'yt−1+c0+d0t) +c1+d1t

Прерывания и линейные тренды присутствуют в cointegrating отношениях, и детерминированные квадратичные тренды присутствуют на уровнях данных.

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

Во время оценки, если полный постоянный, полный линейный тренд модели, cointegrating постоянные, или cointegrating линейные параметры тренда не находится в модели, то estimate ограничивает их обнулять. Если вы задаете различное ограничение равенства, то есть, если свойства, соответствующие тем детерминированным условиям, ограничиваемым обнулять, имеют значение кроме вектора из NaN значения или нули, затем estimate выдает ошибку. Чтобы осуществить поддерживаемые ограничения равенства, выберите модель Йохансена, содержащую детерминированный термин, который вы хотите ограничить.

Пример: 'Model',"H1*"

Типы данных: string | char

Отображение информации оценки вводит в виде разделенной запятой пары, состоящей из 'Display' и значение в этой таблице.

ЗначениеОписание
"off"estimate не отображает информацию об оценке в командной строке.
"table"estimate отображает таблицу информации об оценке. Строки соответствуют параметрам, и столбцы соответствуют оценкам, стандартным погрешностям, статистике t и значениям p.
"full"В дополнение к таблице итоговой статистики, estimate отображает предполагаемые инновации ковариационные и корреляционные матрицы, значение логарифмической правдоподобности, Критерий информации о Akaike (AIC), Байесов информационный критерий (BIC) и другая информация об оценке.

Пример: 'Display',"full"

Типы данных: string | char

Максимальное количество итераций решателя, позволенных в виде разделенной запятой пары, состоящей из 'MaxIterations' и положительный числовой скаляр.

estimate отправки MaxIterations к mvregress.

Типы данных: double

Примечание

NaN значения в Yy0 , и X укажите на отсутствующие значения. estimate удаляет отсутствующие значения из данных мудрым списком удалением.

  • Для предварительной выборки, estimate удаляет любую строку, содержащую по крайней мере один NaN.

  • Для выборки оценки, estimate удаляет любую строку конкатенированной матрицы данных [Y X] содержа по крайней мере один NaN.

Этот тип снижения объема данных уменьшает эффективный объем выборки.

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

свернуть все

Предполагаемая модель VEC (p - 1), возвращенная как vecm объект модели. EstMdl полностью заданный vecm модель.

estimate использование mvregress реализовывать многомерную нормальную оценку наибольшего правдоподобия. Для получения дополнительной информации смотрите Оценку Многомерных Моделей Регрессии.

Предполагаемые, асимптотические стандартные погрешности предполагаемых параметров, возвращенных как массив структур, содержащий поля в этой таблице.

Поле Описание
ConstantСтандартные погрешности полных констант модели (c), соответствующий оценкам в EstMdl.Constant, Mdl.NumSeries- 1 числовой вектор
AdjustmentСтандартные погрешности скоростей корректировки (A), соответствующий оценкам в EstMdl.Adjustment, Mdl.NumSeries- Mdl.Rank числовой вектор
ImpactСтандартные погрешности коэффициента удара (Π), соответствующий оценкам в EstMdl.Impact, Mdl.NumSeries- Mdl.NumSeries числовой вектор
ShortRunСтандартные погрешности коэффициентов короткого промежутка времени (Φ), соответствующий оценкам в EstMdl.ShortRun, вектор ячейки с элементами, соответствующими EstMdl.ShortRun
BetaСтандартные погрешности коэффициентов регрессии (β), соответствующий оценкам в EstMdl.Beta, Mdl.NumSeries- numpreds числовая матрица
TrendСтандартные погрешности полных линейных трендов времени (d), соответствующий оценкам в EstMdl.Trend, Mdl.NumSeries- 1 числовой вектор

Если estimate применяет ограничения равенства во время оценки путем фиксации любых параметров к значению, затем соответствующими стандартными погрешностями тех параметров является 0.

estimate извлечения все стандартные погрешности от инверсии ожидаемой матрицы информации о Фишере, возвращенной mvregress (см. стандартные погрешности).

Оптимизированное значение целевой функции логарифмической правдоподобности, возвращенное в виде числа.

Многомерные остаточные значения подобранной модели, возвращенной как числовая матрица, содержащая numseries столбцы.

  • Если вы задаете Y0, затем E имеет numobs строки (см. Y).

  • В противном случае, E имеет numobsMdl.P строки с учетом преддемонстрационного удаления.

Больше о

свернуть все

Векторная модель исправления ошибок

vector error-correction (VEC) model является многомерной, стохастической моделью временных рядов, состоящей из системы m = numseries уравнения отличного m, differenced переменные отклика. Уравнения в системе могут включать error-correction term, который является линейной функцией ответов на уровнях, используемых, чтобы стабилизировать систему. r cointegrating rank является количеством cointegrating relations, которые существуют в системе.

Каждое уравнение ответа может включать авторегрессивный полином, состоявший из первых различий ряда ответа (short-run polynomial степени p – 1), константа, тренд времени, внешние переменные предикторы, и постоянный тренд и тренд времени в сроке исправления ошибок.

Модель VEC (p - 1) в difference-equation notation и в reduced form может быть описана двумя способами:

  • Этим уравнением является component form модели VEC, где скорости корректировки коинтеграции и матрица коинтеграции являются явными, тогда как матрица удара подразумевается.

    Δyt=A(Byt1+c0+d0t)+c1+d1t+Φ1Δyt1+...+Φp1Δyt(p1)+βxt+εt=c+dt+AByt1+Φ1Δyt1+...+Φp1Δyt(p1)+βxt+εt.

    cointegrating отношениями является B' y t – 1 + c 0 +, d 0t и срок исправления ошибок является A (B' y t – 1 + c 0 + d 0t).

  • Этим уравнением является impact form модели VEC, где матрица удара является явной, тогда как скорости корректировки коинтеграции и матрица коинтеграции подразумеваются.

    Δyt=Πyt1+A(c0+d0t)+c1+d1t+Φ1Δyt1+...+Φp1Δyt(p1)+βxt+εt=c+dt+Πyt1+Φ1Δyt1+...+Φp1Δyt(p1)+βxt+εt.

В уравнениях:

  • yt является m-by-1 вектор из значений, соответствующих переменным отклика m во время t, где t = 1..., T.

  • Δyt = yty t – 1. Структурный коэффициент является единичной матрицей.

  • r является количеством cointegrating отношений и, в целом, 0 <r <m.

  • A является m-by-r матрица скоростей корректировки.

  • B является m-by-r матрица коинтеграции.

  • Π m-by-m матрица удара с рангом r.

  • c 0 является r-by-1 вектор из констант (прерывания) в cointegrating отношениях.

  • d 0 является r-by-1 вектор из линейных трендов времени в cointegrating отношениях.

  • c 1 является m-by-1 вектор из констант (deterministic linear trends в yt).

  • d 1 является m-by-1 вектор из линейных значений тренда времени (deterministic quadratic trends в yt).

  • c = A c 0 + c 1 и является полной константой.

  • d = A d 0 + d 1 и является полным коэффициентом тренда времени.

  • Φj является m-by-m матрица коэффициентов короткого промежутка времени, где j = 1..., p – 1 и Φp – 1 не является матрицей, содержащей только нули.

  • xt является k-by-1 вектор из значений, соответствующих k внешние переменные предикторы.

  • β является m-by-k матрица коэффициентов регрессии.

  • εt является m-by-1 вектор из случайных Гауссовых инноваций, каждого со средним значением 0 и коллективно m-by-m ковариационная матрица Σ. Для ts, εt и εs независимы.

Сжатый и в обозначении оператора задержки, система

Φ(L)(1L)yt=A(Byt1+c0+d0t)+c1+d1t+βxt+εt=c+dt+AByt1+βxt+εt

где Φ(L)=IΦ1Φ2...Φp1, I является m-by-m единичная матрица и L yt = y t – 1.

Если m = r, то модель VEC является устойчивой моделью VAR (p) на уровнях ответов. Если r = 0, то срок исправления ошибок является матрицей нулей и моделью VEC (p - 1), является устойчивой моделью VAR (p - 1) в первых различиях ответов.

Форма Йохансена

Johansen forms Модели VEC отличается относительно присутствия детерминированных условий. Как детализировано в [2], процедура оценки отличается среди форм. Следовательно, допустимые ограничения равенства на детерминированных условиях во время оценки отличаются среди форм. Для получения дополнительной информации смотрите Роль Детерминированных Условий.

Эта таблица описывает пять форм Йохансена и поддерживаемые ограничения равенства.

ФормаСрок исправления ошибокДетерминированные коэффициентыОграничения равенства
H2

AB 'yt − 1

c = 0 (Констант).

d = 0 (Тренд).

c 0 = 0 (CointegrationConstant).

d 0 = 0 (CointegrationTrend).

Можно полностью задать B.

Все детерминированные коэффициенты являются нулем.

H1*

A (B 'yt−1+c0)

c = A c 0.

d = 0.

d 0 = 0.

Если вы полностью задаете или B или c 0, то необходимо полностью задать другой.

MATLAB® получает значение c от c 0 и A.

Все детерминированные тренды являются нулем.

H1

A (B ´yt−1 + c 0) + c 1

c = A c 0 + c 1.

d = 0.

d 0 = 0.

Можно полностью задать B.

Можно задать смесь NaN и числовые значения для c.

MATLAB получает значение c 0 от c и A.

Все детерминированные тренды являются нулем.

H

A (B ´yt−1 + c 0 + d 0t) + c 1

c = A c 0 + c 1.

d = A d 0.

Если вы полностью задаете или B или d 0, то необходимо полностью задать другой.

Можно задать смесь NaN и числовые значения для c.

MATLAB получает значение c 0 от c и A.

MATLAB получает значение d от A и d 0.

H

A (B 'yt−1+c0+d0t) +c1+d1t

c = A c 0 + c 1.

d = A.d0 + d 1.

Можно полностью задать B.

Можно задать смесь NaN и числовые значения для c и d.

MATLAB получает значения c 0 и d 0 от c, d и A.

Алгоритмы

  • Если 1 ≤ Mdl.RankMdl.NumSeries– 1 , как с большинством моделей VEC, затем estimate выполняет оценку параметра два шага.

    1. estimate оценивает параметры cointegrating отношений, включая любые ограниченные прерывания и тренды времени, методом Йохансена [2].

      • Форма cointegrating отношений соответствует одной из пяти параметрических форм, рассмотренных Йохансеном в [2] (см. 'Model'). Для получения дополнительной информации смотрите jcitest и jcontest.

      • Параметр скорости корректировки (A) и матрица коинтеграции (B) в модели VEC (p - 1) не может быть однозначно определен. Однако продукт Π = A *Bʹ идентифицируется. На этом шаге оценки, B = V 1:r, где V 1:r является матрицей, состоявшей из всех строк и первых столбцов r матрицы собственного вектора V. V нормирован так, чтобы V ʹ *S11*V = I. Для получения дополнительной информации см. [2].

    2. estimate создает условия исправления ошибок из предполагаемых cointegrating отношений. То, estimate оценивает остающиеся условия в модели VEC путем построения векторной модели (VAR) авторегрессии в первых различиях и включая условия исправления ошибок с должности предикторов. Для моделей без cointegrating отношений (Mdl.Rank = 0) или с cointegrating матрицей полного ранга (Mdl.Rank = Mdl.Numseries), estimate выполняет этот шаг оценки VAR только.

  • Можно удалить стационарные ряды, которые сопоставлены со стандартными единичными векторами в течение cointegrating отношений от анализа коинтеграции. К предварительному отдельному ряду для стационарности использовать adftest, pptest, kpsstest, и lmctest. Как альтернатива, можно протестировать на стандартные единичные векторы в контексте полной модели при помощи jcontest.

  • Если 1Mdl.RankMdl.NumSeries– 1 , асимптотические ошибочные ковариации параметров в cointegrating отношениях (которые включают B, c 0 и d 0 соответствий Cointegration, CointegrationConstant, и CointegrationTrend свойства, соответственно), являются обычно негауссовыми. Поэтому estimate не оценивает или возвращает соответствующие стандартные погрешности.

    В отличие от этого ошибочные ковариации составной матрицы удара, которая задана как продукт A *Bʹ, являются асимптотически Гауссовыми. Поэтому estimate оценки и возвращают его стандартные погрешности. Подобные протесты содержат для стандартных погрешностей полного постоянного и линейного тренда (A *c0 и A *d0corresponding к Constant и Trend свойства, соответственно) H1* и H* формы Йохансена.

Ссылки

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

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

[3] Juselius, K. Модель VAR Cointegrated. Оксфорд: Издательство Оксфордского университета, 2006.

[4] Lütkepohl, H. Новое введение в несколько анализ временных рядов. Берлин: Спрингер, 2005.

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