Exponenta Banner

fevd

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

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

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

FEVD предоставляет информацию об относительной важности каждых инноваций во влиянии на отклонение ошибки прогноза всех переменных отклика в системе. Напротив, импульсная функция отклика (IRF) прослеживает эффекты инновационного шока для одной переменной на ответе всех переменных в системе. Чтобы оценить IRF модели VEC, охарактеризованной объектом модели vecm, смотрите irf.

Синтаксис

Decomposition = fevd(Mdl)
Decomposition = fevd(Mdl,Name,Value)
[Decomposition,Lower,Upper] = fevd(___)

Описание

пример

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

пример

Decomposition = fevd(Mdl,Name,Value) дополнительные опции использования заданы одним или несколькими аргументами пары "имя-значение". Например, 'NumObs',10,'Method',"generalized" задает оценку обобщенного FEVD в течение многих времен 1 - 10.

пример

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

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

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

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

Примеры

свернуть все

Скрипт Open Live Script

Соответствуйте 4-D модели VEC(2) двумя cointegrating отношениями к датским деньгам и поступите ряд уровня. Затем оценка и график ортогонализируемый FEVD из предполагаемой модели.

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

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).

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

Decomposition = fevd(Mdl);

Decomposition является 20 массивом 4 на 4, представляющим FEVD Mdl. Строки соответствуют последовательным моментам времени со времени 1 - 20, столбцы соответствуют переменным, получающим инновационный шок с одним стандартным отклонением во время 0, и страницы соответствуют переменным, отклонение ошибки прогноза которых fevd анализирует. Mdl.SeriesNames задает переменный порядок.

Поскольку Decomposition представляет ортогонализируемый FEVD, строки должны суммировать к 1. Эта характеристика иллюстрирует, что ортогонализируемые FEVDs представляют пропорции вкладов отклонения. Подтвердите, что все строки Decomposition суммируют к 1.

rowsums = sum(Decomposition,2);
sum((rowsums - 1).^2 > eps)
ans = 
ans(:,:,1) =

     0


ans(:,:,2) =

     0


ans(:,:,3) =

     0


ans(:,:,4) =

     0

Суммы строки среди страниц близко к 1.

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

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

    0.0694
    0.1744
    0.1981
    0.2182
    0.2329
    0.2434
    0.2490
    0.2522
    0.2541
    0.2559
      ⋮

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

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

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

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

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

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

Скрипт Open Live Script

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

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

load Data_JDanish

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

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

Decomposition = fevd(Mdl,"Method","generalized","NumObs",100);

Decomposition является 100 массивом 4 на 4, представляющим обобщенный FEVD Mdl.

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

figure;
plot(1:100,Decomposition(:,2,3))
title("FEVD of IB When Y Is Shocked")
xlabel("Forecast Horizon")
ylabel("Variance Contribution")
grid on

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

Скрипт Open Live Script

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

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

load Data_JDanish

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

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

rng(1); % For reproducibility
[Decomposition,Lower,Upper] = fevd(Mdl);

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

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

fevdshock2resp3 = Decomposition(:,2,3);
FEVDCIShock2Resp3 = [Lower(:,2,3) Upper(:,2,3)];

figure;
h1 = plot(1:20,fevdshock2resp3);
hold on
h2 = plot(1:20,FEVDCIShock2Resp3,'r--');
legend([h1 h2(1)],["FEVD" "95% Confidence Interval"],...
    'Location',"best")
xlabel("Forecast Horizon");
ylabel("Variance Contribution");
title("FEVD of IB When Y Is Shocked");
grid on
hold off

В конечном счете, и когда реальный доход потрясен, пропорция отклонения ошибки прогноза курса облигаций обосновывается между приблизительно 0 и 0.7 с 95%-й уверенностью.

Скрипт Open Live Script

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

Загрузите датские деньги и поступите набор данных, затем оцените модель 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.

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

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

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

fevdshock2resp3 = Decomposition(:,2,3);
FEVDCIShock2Resp3 = [Lower(:,2,3) Upper(:,2,3)];

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

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

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

свернуть все

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

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

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

Пример: 'NumObs',10,'Method',"generalized" задает оценку обобщенного FEVD в течение периодов 1 - 10.

Опции для всего FEVDs

свернуть все

Количество периодов, в течение которых fevd вычисляет FEVD (горизонт прогноза), заданный как пара, разделенная запятой, состоящая из 'NumObs' и положительного целого числа. NumObs задает количество наблюдений, чтобы включать в FEVD (количество строк в Decomposition).

Пример: 'NumObs',10 задает оценку FEVD в течение многих времен 1 - 10.

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

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

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

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

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

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

свернуть все

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Предположим Confidence = c. Затем 100 (1 – c)/2 процент импульсных ответов лежат вне доверительных границ.

Значением по умолчанию является 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

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

свернуть все

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

Decomposition(t,j,k) является вкладом в разложение отклонения переменной k, относящейся к инновационному шоку с одним стандартным отклонением для переменной j во время t, для t = 1,2, …, numobs, j = 1,2..., numseries и k = 1,2..., numseries.

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

Lower(t,j,k) является нижней границей интервала процентили % 100*Confidence на истинном вкладе в разложение отклонения переменной k, относящейся к инновационному шоку с одним стандартным отклонением для переменной j во время 0.

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

Upper(t,j,k) является верхней границей интервала процентили % 100*Confidence на истинном вкладе в разложение отклонения переменной k, относящейся к инновационному шоку с одним стандартным отклонением для переменной j во время 0.

Больше о

свернуть все

Разложение отклонения ошибки прогноза

forecast error variance decomposition (FEVD) многомерной, динамической системы показывает относительную важность шока для каждых инноваций во влиянии на отклонение ошибки прогноза всех переменных в системе.

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

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

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

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

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

Общая форма FEVD ykt (переменная k) периоды m в будущее, относящееся к инновационному шоку с одним стандартным отклонением для yjt,

γmjk=∑t=0m−1(ek′Ctej)2∑t=0m−1ek′ΩtΣΩt′ek.

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

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

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

  • Числитель является вкладом инновационного шока для переменной j к отклонению ошибки прогноза m - неродной вперед прогноз переменной k. Знаменатель является среднеквадратичной погрешностью (MSE) m - неродной вперед прогноз переменной k [4].

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

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) в первых различиях ответов.

Алгоритмы

  • Если Method является "orthogonalized", то fevd ортогонализирует инновационные шоки путем применения факторизации Холесского образцовой ковариационной матрицы Mdl.Covariance. Ковариация ортогонализируемых инновационных шоков является единичной матрицей и FEVD каждой переменной суммы одной, то есть, сумма вдоль любой строки Decomposition является той. Поэтому ортогонализируемый FEVD представляет пропорцию отклонения ошибки прогноза, относящегося к различным шокам в системе. Однако ортогонализируемый FEVD обычно зависит от порядка переменных.

    Если Method является "generalized", то получившийся FEVD, то получившийся FEVD является инвариантным к порядку переменных, и не основан на ортогональном преобразовании. Кроме того, получившийся FEVD суммирует одному для конкретной переменной только, когда Mdl.Covariance диагональный [5]. Поэтому обобщенный FEVD представляет вклад в отклонение ошибки прогноза мудрых уравнением шоков для переменных отклика в модели.

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

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

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

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

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

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

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

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

      3. Оцените NumPaths, FEVDs от NumPaths оценил модели.

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

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

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

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

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

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

Ссылки

[1] Гамильтон, J. D. Анализ timeseries. Принстон, NJ: Издательство Принстонского университета, 1994.

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

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

[4] Lütkepohl, H. Новое введение в несколько анализ временных рядов. Нью-Йорк, Нью-Йорк: Springer-Verlag, 2007.

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

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

Объекты

Функции

Введенный в R2019a

Документация Econometrics Toolbox
Поддержка