Exponenta Banner

irf

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

свернуть все на странице

Синтаксис

Response = irf(Mdl)
Response = irf(Mdl,Name,Value)
[Response,Lower,Upper] = irf(___)

Описание

irf функция возвращает динамический ответ или функцию импульсной характеристики (IRF), к шоку с одним стандартным отклонением для каждой переменной в модели VEC (p - 1). Полностью заданный vecm объект модели характеризует модель VEC.

IRFs прослеживают эффекты инновационного шока для одной переменной на ответе всех переменных в системе. В отличие от этого разложение отклонения ошибки прогноза (FEVD) предоставляет информацию об относительной важности каждых инноваций во влиянии на все переменные в системе. Оценить FEVD модели VEC, охарактеризованной vecm объект модели, смотрите fevd.

пример

Response = irf(Mdl) возвращает ортогонализируемый IRF с 20 периодами переменных временных рядов, которые составляют модель VEC (p - 1) Mdl охарактеризованный полностью заданным vecm объект модели. irf переменные шоков во время 0, и возвращают IRF в течение многих времен 0 до 19.

пример

Response = irf(Mdl,Name,Value) дополнительные опции использования заданы одним или несколькими аргументами пары "имя-значение". Например, 'NumObs',10,'Method',"generalized" задает оценку обобщенного IRF для 10 моментов времени, запускающихся во время 0, во время который irf применяет шок, и заканчивающийся в период 9.

пример

[Response,Lower,Upper] = irf(___) использование любая из комбинаций входных аргументов в предыдущих синтаксисах и возвращает более низкие и верхние доверительные границы, в течение каждого периода и переменной в IRF, 95%-го доверительного интервала на истинном IRF.

  • Если вы задаете серию остаточных значений при помощи E аргумент пары "имя-значение", затем irf оценивает доверительные границы путем начальной загрузки заданных остаточных значений.

  • В противном случае, irf оценочные доверительные границы путем проведения симуляции Монте-Карло.

Если Mdl пользовательский vecm объект модели (объект, не возвращенный estimate или измененный после оценки), irf может потребовать объема выборки для симуляции SampleSize или преддемонстрационные ответы Y0.

Примеры

свернуть все

Скрипт Open Live Script

Подбирайте 4-D модель VEC(2) к датским деньгам и поступите ряд уровня. Затем оцените ортогонализируемый IRF из предполагаемой модели.

Загрузите датские деньги и поступите набор данных.

load Data_JDanish

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

Создайте vecm объект модели, который представляет 4-D модель VEC(2) двумя cointegrating отношениями. Задайте имена переменных.

Mdl = vecm(4,2,2);
Mdl.SeriesNames = DataTable.Properties.VariableNames;

Mdl vecm объект модели, задающий структуру 4-D модели VEC(2); это - шаблон для оценки.

Подбирайте модель VEC(2) к набору данных.

Mdl = estimate(Mdl,DataTable.Series);

Mdl полностью заданный vecm объект модели, представляющий предполагаемую 4-D модель VEC(2).

Оцените ортогонализируемый IRF из предполагаемой модели VEC(2).

Response = irf(Mdl);

Response 20 массивом 4 на 4, представляющим IRF Mdl. Строки соответствуют последовательным моментам времени со времени от 0 до 19, столбцы соответствуют переменным, получающим инновационный шок с одним стандартным отклонением во время 0, и страницы соответствуют ответам переменных к потрясаемой переменной. Mdl.SeriesNames задает переменный порядок.

Отобразите IRF курса облигаций (переменные 3, IB) когда журнал реального дохода (переменные 2, Y) потрясен во время 0.

Response(:,2,3)
ans = 20×1

    0.0021
    0.0057
    0.0064
    0.0067
    0.0064
    0.0061
    0.0057
    0.0056
    0.0057
    0.0058
      ⋮

armairf графики функций IRF моделей VAR охарактеризованы содействующими матрицами AR. Постройте IRF модели VEC:

  1. Выражение модели VEC(2) как модель VAR (3) путем передачи Mdl к varm

  2. Передача содействующей и инновационной ковариационной матрицы AR модели VAR к armairf

Постройте модель VEC(2) IRF в течение 40 периодов.

VARMdl = varm(Mdl);
armairf(VARMdl.AR,[],"InnovCov",VARMdl.Covariance,...
    "NumObs",40);

Figure contains an axes object. The axes object with title Orthogonalized IRF of Variable 1 contains 4 objects of type line. These objects represent Shock to Variable 1, Shock to Variable 2, Shock to Variable 3, Shock to Variable 4.

Figure contains an axes object. The axes object with title Orthogonalized IRF of Variable 2 contains 4 objects of type line. These objects represent Shock to Variable 1, Shock to Variable 2, Shock to Variable 3, Shock to Variable 4.

Figure contains an axes object. The axes object with title Orthogonalized IRF of Variable 3 contains 4 objects of type line. These objects represent Shock to Variable 1, Shock to Variable 2, Shock to Variable 3, Shock to Variable 4.

Figure contains an axes object. The axes object with title Orthogonalized IRF of Variable 4 contains 4 objects of type line. These objects represent Shock to Variable 1, Shock to Variable 2, Shock to Variable 3, Shock to Variable 4.

Каждый график показывает четыре IRFs переменной, когда все другие переменные потрясены во время 0. Mdl.SeriesNames задает переменный порядок.

Скрипт Open Live Script

Считайте 4-D модель VEC(2) с двумя cointegrating отношениями в Оценке и Графике Моделью VEC IRF. Оцените обобщенный IRF системы в течение 50 периодов.

Загрузите датские деньги и поступите набор данных, затем оцените модель VEC(2).

load Data_JDanish

Mdl = vecm(4,2,2);
Mdl.SeriesNames = DataTable.Properties.VariableNames;
Mdl = estimate(Mdl,DataTable.Series);

Оцените обобщенный IRF из предполагаемой модели VEC(2).

Response = irf(Mdl,"Method","generalized","NumObs",50);

Response 50 массивом 4 на 4, представляющим обобщенный IRF Mdl.

Постройте обобщенный IRF курса облигаций, когда реальный доход будет потрясен во время 0.

figure;
plot(0:49,Response(:,2,3))
title("IRF of IB When Y Is Shocked")
xlabel("Observation Time")
ylabel("Response")
grid on

Figure contains an axes object. The axes object with title IRF of IB When Y Is Shocked contains an object of type line.

Когда реальный доход потрясен, курс облигаций реагирует и затем обосновывается приблизительно в 0,0032 после 15 периодов.

Скрипт Open Live Script

Считайте 4-D модель VEC(2) с двумя cointegrating отношениями в Оценке и Графике Моделью VEC IRF. Оцените и постройте его ортогонализируемый IRF и 95% доверительных интервалов Монте-Карло на истинном IRF.

Загрузите датские деньги и поступите набор данных, затем оцените модель VEC(2).

load Data_JDanish

Mdl = vecm(4,2,2);
Mdl.SeriesNames = DataTable.Properties.VariableNames;
Mdl = estimate(Mdl,DataTable.Series);

Оцените IRF и соответствующие 95% доверительных интервалов Монте-Карло из предполагаемой модели VEC(2).

rng(1); % For reproducibility
[Response,Lower,Upper] = irf(Mdl);

Responseниже, и Upper 20 массивами 4 на 4, представляющими ортогонализируемый IRF Mdl и соответствующие нижние и верхние границы доверительных интервалов. Для всех массивов строки соответствуют последовательным моментам времени со времени от 0 до 19, столбцы соответствуют переменным, получающим инновационный шок с одним стандартным отклонением во время 0, и страницы соответствуют ответам переменных к потрясаемой переменной. Mdl.SeriesNames задает переменный порядок.

Постройте ортогонализируемый IRF с его доверительными границами курса облигаций, когда реальный доход будет потрясен во время 0.

irfshock2resp3 = Response(:,2,3);
IRFCIShock2Resp3 = [Lower(:,2,3) Upper(:,2,3)];

figure;
h1 = plot(0:19,irfshock2resp3);
hold on
h2 = plot(0:19,IRFCIShock2Resp3,'r--');
legend([h1 h2(1)],["IRF" "95% Confidence Interval"])
xlabel("Time Index");
ylabel("Response");
title("IRF of IB When Y Is Shocked");
grid on
hold off

Figure contains an axes object. The axes object with title IRF of IB When Y Is Shocked contains 3 objects of type line. These objects represent IRF, 95% Confidence Interval.

Когда реальный доход потрясен, курс облигаций реагирует и затем обосновывается между –0.002 и 0.0095 с 95%-м доверием.

Скрипт Open Live Script

Считайте 4-D модель VEC(2) с двумя cointegrating отношениями в Оценке и Графике Моделью VEC IRF. Оцените и постройте его ортогонализируемый IRF и 90%-е доверительные интервалы начальной загрузки на истинном IRF.

Загрузите датские деньги и поступите набор данных, затем оцените модель VEC(2). Возвратите остаточные значения оценки модели.

load Data_JDanish

Mdl = vecm(4,2,2);
Mdl.SeriesNames = DataTable.Properties.VariableNames;
[Mdl,~,~,E] = estimate(Mdl,DataTable.Series);
T = size(DataTable,1) % Total sample size
T = 55

n = size(E,1)         % Effective sample size
n = 52

E 52 4 массив остаточных значений. Столбцы соответствуют переменным в Mdl.SeriesNames. estimate функция требует Mdl.P = 3 наблюдения, чтобы инициализировать модель VEC(2) для оценки. Поскольку преддемонстрационные данные (Y0) не задано, estimate берет первые три наблюдения в заданных данных об ответе, чтобы инициализировать модель. Поэтому получившимся эффективным объемом выборки является T – Mdl.P = 52, и строки E соответствуйте индексам наблюдения 4 через T.

Оцените ортогонализируемый IRF и соответствующие 90%-е доверительные интервалы начальной загрузки из предполагаемой модели VEC(2). Чертите 500 путей длины n от серии остаточных значений.

rng(1); % For reproducibility
[Response,Lower,Upper] = irf(Mdl,"E",E,"NumPaths",500,...
    "Confidence",0.9);

Постройте ортогонализируемый IRF с его доверительными границами курса облигаций, когда реальный доход будет потрясен во время 0.

irfshock2resp3 = Response(:,2,3);
IRFCIShock2Resp3 = [Lower(:,2,3) Upper(:,2,3)];

figure;
h1 = plot(0:19,irfshock2resp3);
hold on
h2 = plot(0:19,IRFCIShock2Resp3,'r--');
legend([h1 h2(1)],["IRF" "90% Confidence Interval"])
xlabel("Time Index");
ylabel("Response");
title("IRF of IB When Y Is Shocked");
grid on
hold off

Figure contains an axes object. The axes object with title IRF of IB When Y Is Shocked contains 3 objects of type line. These objects represent IRF, 90% Confidence Interval.

Когда реальный доход потрясен, курс облигаций реагирует и затем обосновывается между приблизительно 0 и 0.010 с 90%-м доверием.

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

свернуть все

Модель VEC в виде vecm объект модели создается vecm или estimate. Mdl должен быть полностью задан.

Аргументы name-value

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

Пример: 'NumObs',10,'Method',"generalized" задает оценку обобщенного IRF для 10 моментов времени, запускающихся во время 0, во время который irf применяет шок, и заканчивающийся в период 9.
Опции для всего IRFs

свернуть все

Количество периодов, для который irf вычисляет IRF в виде разделенной запятой пары, состоящей из 'NumObs' и положительное целое число. NumObs задает количество наблюдений, чтобы включать в IRF (количество строк в Response).

Пример: 'NumObs',10 задает включение 10 моментов времени в IRF, запускающемся во время 0, во время который irf применяет шок, и заканчивающийся в период 9.

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

Метод расчета IRF в виде разделенной запятой пары, состоящей из 'Method' и значение в этой таблице.

ЗначениеОписание
"orthogonalized"Вычислите импульсные характеристики с помощью ортогонализируемых, инновационных шоков с одним стандартным отклонением. irf использует факторизацию Холесского Mdl.Covariance для ортогонализации.
"generalized"Вычислите импульсные характеристики с помощью инновационных шоков с одним стандартным отклонением.

Пример: 'Method',"generalized"

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

Опции для оценки доверительной границы

свернуть все

Количество демонстрационных путей (испытания), чтобы сгенерировать в виде разделенной запятой пары, состоящей из 'NumPaths' и положительное целое число.

Пример: 'NumPaths',1000 генерирует 1000 демонстрационные пути, из которых программное обеспечение выводит доверительные границы.

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

Количество наблюдений для симуляции Монте-Карло или начальной загрузки на демонстрационный путь в виде разделенной запятой пары, состоящей из 'SampleSize' и положительное целое число.

  • Если Mdl предполагаемый vecm объект модели (объект, возвращенный estimate и немодифицированный после этого), затем значением по умолчанию является объем выборки данных, к которым модель является подходящей (см. summarize).

  • Если irf оценочные доверительные границы путем проведения симуляции Монте-Карло (для получения дополнительной информации смотрите E), необходимо задать SampleSize.

  • Если irf оценочные доверительные границы путем начальной загрузки остаточных значений, значением по умолчанию является длина заданной серии остаточных значений (size(E,1)).

Пример: Если вы задаете 'SampleSize',100 и не задавайте 'E' аргумент пары "имя-значение", программное обеспечение оценивает доверительные границы от NumPaths случайные пути длины 100 от Mdl.

Пример: Если вы задаете 'SampleSize',100,'E',E, программное обеспечение передискретизирует, с заменой, 100 наблюдения (строки) от E сформировать демонстрационный путь инноваций, чтобы проникнуть в Mdl. Программное обеспечение формирует NumPaths пути к случайной выборке, из которых это выводит доверительные границы.

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

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

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

numseries размерность модели Mdl.NumSeries входа VEC. Столбцы должны соответствовать переменным отклика в Mdl.SeriesNames.

  • Если Mdl предполагаемый vecm объект модели (объект, возвращенный estimate и немодифицированный после этого), irf наборы Y0 к преддемонстрационным данным об ответе, используемым для оценки по умолчанию (см. 'Y0').

  • В противном случае необходимо задать Y0.

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

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

numpreds количество переменных предикторов (size(Mdl.Beta,2)).

Строки соответствуют наблюдениям. X должен иметь, по крайней мере, SampleSize 'Строки' . Если вы предоставляете больше строк, чем необходимый, irf использование только последний SampleSize наблюдения. Последняя строка содержит последнее наблюдение.

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

Обеспечить непротиворечивость модели когда irf оценивает доверительные границы, хорошая практика должна задать X когда Mdl имеет компонент регрессии. Если Mdl предполагаемая модель, задайте данные о предикторе, используемые во время оценки модели (см. 'X').

По умолчанию, irf исключает компонент регрессии из оценки доверительной границы, независимо от ее присутствия в Mdl.

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

Серия остаточных значений, от которых можно чертить выборки начальной загрузки в виде разделенной запятой пары, состоящей из 'E' и числовая матрица, содержащая numseries столбцы. irf принимает тот E свободно от последовательной корреляции.

Столбцы содержат остаточный ряд, соответствующий серийным именам ответа в Mdl.SeriesNames.

Если Mdl предполагаемый vecm объект модели (объект, возвращенный estimate), можно задать E как выведенные остаточные значения оценки (см. E или infer).

По умолчанию, irf выводит доверительные границы путем проведения симуляции Монте-Карло.

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

Доверительный уровень для доверительных границ в виде числового скаляра в интервале [0,1].

В течение каждого периода случайным образом чертившие доверительные интервалы покрывают истинный ответ 100*Confidence% из времени.

Значением по умолчанию является 0.95, который подразумевает, что доверительные границы представляют 95% доверительных интервалов.

Типы данных: 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 отношениях, и детерминированные квадратичные тренды присутствуют на уровнях данных.

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

Для получения дополнительной информации о формах Йохансена смотрите estimate.

  • Если Mdl предполагаемый vecm объект модели (объект, возвращенный estimate и немодифицированный после этого), значением по умолчанию является форма Йохансена, используемая для оценки (см. 'Model').

  • В противном случае значением по умолчанию является "H1".

Совет

Лучшая практика должна обеспечить непротиворечивость модели во время симуляции, которая оценивает доверительные границы. Поэтому, если Mdl предполагаемый vecm объект модели (объект, возвращенный estimate и немодифицированный после этого), включите любые ограничения, наложенные во время оценки путем подчинения значению по умолчанию Model.

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

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

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

свернуть все

IRF, возвращенный как numobs- numseries- numseries числовой массив. numobs значение NumObs. Столбцы и страницы соответствуют переменным отклика в Mdl.SeriesNames.

Ответ (t + 1, jK) импульсная характеристика переменной k к инновациям с одним стандартным отклонением потрясают к переменной j во время 0, для t = 0, 1..., numObs – 1, j = 1,2..., numseries, и k = 1,2..., numseries.

Более низкие доверительные границы, возвращенные как numobs- numseries- numseries числовой массив. Элементы Lower соответствуйте элементам Response.

Ниже (t + 1, jK) нижняя граница 100*Confidence% интервал процентили на истинной импульсной характеристике переменной k к инновациям с одним стандартным отклонением потрясают к переменной j во время 0.

Верхние доверительные границы, возвращенные как numobs- numseries- numseries числовой массив. Элементы Upper соответствуйте элементам Response.

Верхний (t + 1, jK) верхняя граница 100*Confidence% интервал процентили на истинной импульсной характеристике переменной k к инновациям с одним стандартным отклонением потрясают к переменной j во время 0.

Больше о

свернуть все

Функция импульсной характеристики

impulse response function (IRF) модели временных рядов (или dynamic response of the system) измеряет изменения в будущих ответах всех переменных в системе, когда переменная потрясена импульсом. Другими словами, IRF во время, t является производной ответов во время t относительно инноваций во время t 0 (время, когда инновации были потрясены), t ≥ t 0.

Рассмотрите numseries- Модель D VEC (p - 1) для многомерной переменной отклика yt. В обозначении оператора задержки эквивалентный VAR (p) представление модели VEC (p - 1):

Γ(L)yt=c+dt+βxt+εt,

где Γ(L)=I−Γ1L−Γ2L2−...−ΓpLp и I является numseries- numseries единичная матрица.

В обозначении оператора задержки бесконечное представление MA задержки yt:

yt=Γ−1(L)(c+βxt+dt)+Γ−1(L)εt=Ω(L)(c+βxt+dt)+Ω(L)εt.

Общая форма IRF yt, потрясенного импульсом к переменной j одним стандартным отклонением его инноваций периоды m в будущее:

ψj(m)=Cmej.

  • ej является вектором выбора из длины numseries содержа 1 в элементе j и нули в другом месте.

  • Для ортогонализируемого IRF, Cm=ΩmP, где P является нижним треугольным фактором в факторизации Холесского Σ, и Ωm является задержкой коэффициент m Ω (L).

  • Для обобщенного IRF, Cm=σj−1ΩmΣ, где σj является стандартным отклонением инноваций j.

  • IRF свободен от постоянного, компонента регрессии модели и тренда времени.

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

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(B′yt−1+c0+d0t)+c1+d1t+Φ1Δyt−1+...+Φp−1Δyt−(p−1)+βxt+εt=c+dt+AB′yt−1+Φ1Δyt−1+...+Φp−1Δyt−(p−1)+βxt+εt.

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

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

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

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

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

  • Δyt = yt – y 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 ковариационная матрица Σ. Для t ≠ s, εt и εs независимы.

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

Φ(L)(1−L)yt=A(B′yt−1+c0+d0t)+c1+d1t+βxt+εt=c+dt+AB′yt−1+βxt+εt

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

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

Алгоритмы

  • NaN значения в Y0X, и E укажите на недостающие данные. irf удаляет недостающие данные из этих аргументов мудрым списком удалением. Каждый аргумент, если строка содержит по крайней мере один NaNто irf удаляет целую строку.

    Мудрое списком удаление уменьшает объем выборки, может создать неправильные временные ряды и может вызвать E и X не синхронизироваться.

  • Если Method "orthogonalized", затем получившийся IRF зависит от порядка переменных в модели временных рядов. Если Method "generalized", затем получившийся IRF является инвариантным к порядку переменных. Поэтому эти два метода обычно приводят к различным результатам.

  • Если Mdl.Covariance диагональная матрица, затем получившиеся обобщенные и ортогонализируемые IRFs идентичны. В противном случае получившиеся обобщенные и ортогонализируемые IRFs идентичны только, когда первая переменная потрясает все переменные (то есть, все остальное являющееся тем же самым, оба метода дают к тому же значению Response(:,1,:)).

  • Данные о предикторе X представляет один путь внешних многомерных временных рядов. Если вы задаете X и модель VAR Mdl имеет компонент регрессии (Mdl.Beta не пустой массив), irf применяет те же внешние данные ко всем путям, используемым для оценки доверительного интервала.

  • irf проводит симуляцию, чтобы оценить доверительные границы Lower и Upper.

    • Если вы не задаете остаточные значения Eто irf проводит симуляцию Монте-Карло путем выполнения этой процедуры:

      1. Симулируйте NumPaths пути к ответу длины SampleSize от Mdl.

      2. Подходящий NumPaths модели, которые имеют ту же структуру как Mdl к путям к симулированному отклику. Если Mdl содержит компонент регрессии, и вы задаете Xто irf соответствует NumPaths модели к путям к симулированному отклику и X (те же данные о предикторе для всех путей).

      3. Оцените NumPaths IRFs от NumPaths предполагаемые модели.

      4. Для каждого момента времени t = 0, …, NumObs, оцените доверительные интервалы путем вычисления 1 – Confidence и Confidence квантили (верхние и нижние границы, соответственно).

    • Если вы задаете остаточные значения Eто irf проводит непараметрическую начальную загрузку путем выполнения этой процедуры:

      1. Передискретизируйте, с заменой, SampleSize остаточные значения E. Выполните этот шаг NumPaths времена, чтобы получить NumPaths пути .

      2. Сосредоточьте каждый путь загруженных остаточных значений.

      3. Отфильтруйте каждый путь загруженных остаточных значений в центре через Mdl получить NumPaths загруженные пути к ответу длины SampleSize.

      4. Полные шаги 2 - 4 симуляции Монте-Карло, но замена пути к симулированному отклику с загруженными путями к ответу.

Ссылки

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

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

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

[4] Pesaran, H. H. и И. Шин. "Обобщенный Анализ Импульсной характеристики в Линейных Многомерных Моделях". Экономические Буквы. Издание 58, 1998, стр 17–29.

Смотрите также

Объекты

Функции

Введенный в R2019a

Документация Econometrics Toolbox

Поддержка