armafevd

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

Описание

armafevd функция возвращает или строит разложение отклонения ошибки прогноза переменных в одномерном или векторном (многомерном) авторегрессивном скользящем среднем значении (ARMA или VARMA) модель, заданная массивами коэффициентов или полиномов оператора задержки.

В качестве альтернативы можно возвратить FEVD из полностью заданный (например, оцененный) объект модели при помощи функции в этой таблице.

Объект моделиФункция IRF
varmfevd
vecmfevd

FEVD предоставляет информацию об относительной важности каждых инноваций во влиянии на отклонение ошибки прогноза всех переменных в системе. В отличие от этого функция импульсной характеристики (IRF) прослеживает эффекты инновационного шока для одной переменной на ответе всех переменных в системе. Чтобы оценить IRFs одномерных или многомерных моделей ARMA, смотрите armairf.

пример

armafevd(ar0,ma0) графики, на отдельных рисунках, FEVD numVars переменные временных рядов, которые составляют ARMA (p, q) модель, с авторегрессивным (AR) и коэффициенты скользящего среднего значения (MA) ar0 и ma0, соответственно. Каждая фигура соответствует переменной и содержит numVars линейные графики. Линейные графики являются FEVDs той переменной по горизонту прогноза, следование из инновационного шока с одним стандартным отклонением применилось ко всем переменным в системе во время 0.

armafevd функция:

  • Принимает векторы или векторы ячейки из матриц в обозначении разностного уравнения

  • Принимает LagOp изолируйте полиномы оператора, соответствующие AR и полиномам MA в обозначении оператора задержки

  • Вмещает модели временных рядов, которые являются одномерными или многомерными, стационарными или интегрированными, структурными или в уменьшаемой форме, и обратимыми или необратимыми

  • Принимает, что постоянный c модели 0

пример

armafevd(ar0,ma0,Name,Value) строит numVars FEVDs с дополнительными опциями, заданными одним или несколькими аргументами пары "имя-значение". Например, 'NumObs',10,'Method',"generalized" задает горизонт прогноза с 10 периодами и оценку обобщенного FEVD.

пример

Y = armafevd(___) возвращает numVars FEVDs, использующий любую из комбинаций входных аргументов в предыдущих синтаксисах.

armafevd(ax,___) графики к осям заданы в ax вместо осей в последних данных. Опция ax может предшествовать любой из комбинаций входных аргументов в предыдущих синтаксисах.

[Y,h] = armafevd(___) дополнительно возвращает указатели на нанесенные на график графические объекты. Используйте элементы h изменить свойства возвращенных графиков.

Примеры

свернуть все

Постройте FEVD одномерной модели ARMA(2,1)

yt=0.3yt-1-0.1yt-2+εt+0.05εt-1.

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

AR0 = [0.3 -0.1];
MA0 = 0.05;

Постройте ортогонализируемый FEVD yt.

armafevd(AR0,MA0);

Figure contains an axes. The axes with title Orthogonalized FEVD of Variable 1 contains an object of type line. This object represents Shock to Variable 1.

Поскольку yt является одномерным, FEVD тривиален.

Постройте FEVD модели VARMA(3,1)

yt=[-0.50.20.10.30.1-0.1-0.40.20.05]yt-1+[-0.050.020.010.10.010.001-0.040.020.005]yt-3+εt+[-0.020.030.30.0030.0010.010.30.010.01]εt-1

где yt=[y1ty2ty3t] и εt=[ε1tε2tε3t].

Модель VARMA находится в обозначении разностного уравнения, потому что текущий ответ изолируется от всех других условий в уравнении.

Создайте вектор ячейки, содержащий коэффициенты матрицы VAR. Положение матрицы коэффициентов в векторе ячейки определяет свою задержку. Поэтому задайте 3х3 матрицу нулей как второй элемент вектора.

var0 = {[-0.5 0.2 0.1; 0.3 0.1 -0.1; -0.4 0.2 0.05],...
    zeros(3),...
    [-0.05 0.02 0.01; 0.1 0.01 0.001; -0.04 0.02 0.005]};

Создайте вектор ячейки, содержащий коэффициенты матрицы VMA.

vma0 = {[-0.02 0.03 0.3; 0.003 0.001 0.01; 0.3 0.01 0.01]};

Постройте ортогонализируемый FEVDs модели VARMA.

armafevd(var0,vma0);

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

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

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

armafevd возвращает три фигуры. Рисунок k содержит обобщенный FEVD переменной k к шоку, применился ко всем другим переменным во время 0.

  • Можно приписать большую часть отклонения ошибки прогноза переменного 1 к шоку для переменного 1. Шок для переменных 2 не способствует очень отклонению ошибки прогноза переменного 1.

  • Можно приписать большую часть отклонения ошибки прогноза переменных 2 к шоку для переменных 2. Шок для переменных 3 не способствует очень отклонению ошибки прогноза переменных 2.

  • Можно приписать большую часть отклонения ошибки прогноза переменных 3 к шоку для переменных 3. Шок для переменных 2 не способствует очень отклонению ошибки прогноза переменных 3.

Постройте целый FEVD структурной модели VARMA(8,4)

{[10.2-0.10.031-0.150.9-0.251]-[-0.50.20.10.30.1-0.1-0.40.20.05]L4-[-0.050.020.010.10.010.001-0.040.020.005]L8}yt={[100010001]+[-0.020.030.30.0030.0010.010.30.010.01]L4}εt

где yt=[y1ty2ty3t] и εt=[ε1tε2tε3t].

Модель VARMA находится в обозначении оператора задержки, потому что ответ и инновационные векторы находятся на противоположных сторонах уравнения.

Создайте вектор ячейки, содержащий коэффициенты матрицы VAR. Поскольку эта модель является структурной моделью в обозначении оператора задержки, начните с коэффициента yt и введите остальных в порядок задержкой. Создайте вектор, который указывает на степень термина задержки для соответствующих коэффициентов (задержка структурного коэффициента 0).

var0 = {[1 0.2 -0.1; 0.03 1 -0.15; 0.9 -0.25 1],...
    -[-0.5 0.2 0.1; 0.3 0.1 -0.1; -0.4 0.2 0.05],...
    -[-0.05 0.02 0.01; 0.1 0.01 0.001; -0.04 0.02 0.005]};
var0Lags = [0 4 8];

Создайте вектор ячейки, содержащий коэффициенты матрицы VMA. Поскольку эта модель находится в обозначении оператора задержки, начните с коэффициента εt и введите остальных в порядок задержкой. Создайте вектор, который указывает на степень термина задержки для соответствующих коэффициентов.

vma0 = {eye(3),...
    [-0.02 0.03 0.3; 0.003 0.001 0.01; 0.3 0.01 0.01]};
vma0Lags = [0 4];

Создайте отдельные полиномы оператора задержки, которые описывают VAR и компоненты VMA модели VARMA.

VARLag = LagOp(var0,'Lags',var0Lags);
VMALag = LagOp(vma0,'Lags',vma0Lags);

Постройте обобщенный FEVDs модели VARMA.

armafevd(VARLag,VMALag,'Method','generalized');

Figure contains an axes. The axes with title Generalized FEVD of Variable 1 contains 3 objects of type line. These objects represent Shock to Variable 1, Shock to Variable 2, Shock to Variable 3.

Figure contains an axes. The axes with title Generalized FEVD of Variable 2 contains 3 objects of type line. These objects represent Shock to Variable 1, Shock to Variable 2, Shock to Variable 3.

Figure contains an axes. The axes with title Generalized FEVD of Variable 3 contains 3 objects of type line. These objects represent Shock to Variable 1, Shock to Variable 2, Shock to Variable 3.

armafevd возвращает три фигуры. Рисунок k содержит обобщенный FEVD переменной k к шоку, применился ко всем другим переменным во время 0.

  • Можно приписать большую часть отклонения ошибки прогноза переменного 1 к шоку для переменного 1. Шоки для переменных 2 и 3 способствуют так же отклонению ошибки прогноза переменного 1.

  • Можно приписать большую часть отклонения ошибки прогноза переменных 2 к шоку для переменных 2. Шок для переменных 3 не способствует очень отклонению ошибки прогноза переменных 2.

  • Можно приписать большую часть отклонения ошибки прогноза переменных 3 к шокам для переменных 1 и 3, каждый способствующие подобные суммы. Шок для переменных 2 не способствует очень отклонению ошибки прогноза переменных 3.

Вычислите обобщенный FEVDs двумерной модели VAR (3)

yt=[1-0.2-0.10.3]yt-1-[0.75-0.1-0.050.15]yt-2+[0.55-0.02-0.010.03]yt-3+εt.

В уравнении, yt=[y1,ty2,t], εt=[ε1,tε2,t], и, для всего t, εt является Гауссовым со средней нулевой и ковариационной матрицей

Σ=[0.5-0.1-0.10.25].

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

AR1 = [1 -0.2; -0.1 0.3];
AR2 = -[0.75 -0.1; -0.05 0.15];
AR3 = [0.55 -0.02; -0.01 0.03];
ar0 = {AR1 AR2 AR3};

InnovCov = [0.5 -0.1; -0.1 0.25];

Вычислите обобщенный FEVDs yt. Поскольку никакие условия MA не существуют, задают пустой массив ([]) для второго входного параметра.

Y = armafevd(ar0,[],'Method','generalized','InnovCov',InnovCov);
size(Y)
ans = 1×3

    31     2     2

Y(10,1,2)
ans = 0.1302

Y 31 массивом 2 на 2 FEVDs. Строки соответствуют временам 1 - 31 в горизонте прогноза, столбцы соответствуют переменным что armafevd шоки во время 0 и страницы соответствуют FEVD переменных в системе. Например, вкладом в отклонение ошибки прогноза переменных 2 во время 10 в горизонте прогноза, относящемся к шоку для переменного 1, является Y(10,1,2) = 0.1302.

armafevd удовлетворяет останавливающемуся критерию после 31 периода. Можно задать, чтобы прекратить раньше использовать 'NumObs' аргумент пары "имя-значение". Эта практика выгодна, когда система имеет много переменных.

Вычислите и отобразите обобщенный FEVDs в течение первых 10 периодов.

Y10 = armafevd(ar0,[],'Method','generalized','InnovCov',InnovCov,...
    'NumObs',10)
Y10 = 
Y10(:,:,1) =

    1.0000    0.0800
    0.9912    0.1238
    0.9863    0.1343
    0.9863    0.1341
    0.9873    0.1294
    0.9874    0.1313
    0.9864    0.1342
    0.9864    0.1343
    0.9866    0.1336
    0.9867    0.1336


Y10(:,:,2) =

    0.0800    1.0000
    0.1157    0.9838
    0.1235    0.9737
    0.1236    0.9737
    0.1237    0.9736
    0.1264    0.9709
    0.1296    0.9679
    0.1298    0.9677
    0.1298    0.9677
    0.1302    0.9673

Y10 10 массивом 2 на 2 FEVDs. Строки соответствуют временам 1 - 10 в горизонте прогноза. Во всем FEVDs вклады, кажется, стабилизировались, прежде чем 10 периодов протекают.

Для каждой переменной (страница) вычислите суммы строки.

sum(Y10,2)
ans = 
ans(:,:,1) =

    1.0800
    1.1150
    1.1206
    1.1204
    1.1167
    1.1187
    1.1206
    1.1207
    1.1202
    1.1203


ans(:,:,2) =

    1.0800
    1.0995
    1.0972
    1.0973
    1.0973
    1.0973
    1.0975
    1.0975
    1.0975
    1.0975

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

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

свернуть все

Авторегрессивные коэффициенты ARMA (p, q) модель в виде числового вектора, вектора ячейки из квадратных числовых матриц или LagOp изолируйте объект полинома оператора. Если ar0 вектор (числовой или ячейка), затем коэффициент yt является идентичностью (eye(numVars)).

Для модели MA задайте пустой массив или ячейку ([] или {}).

  • Для одномерных моделей временных рядов, ar0 числовой вектор, вектор ячейки из скаляров или одномерный LagOp изолируйте полином оператора. Для векторов, ar0 имеет длину p, и элементы соответствуют изолированным ответам, которые составляют полином AR в обозначении разностного уравнения. Другими словами, ar0(j) или ar0{j} коэффициент yt-j, j = 1, …, p. Разложения отклонения одномерных моделей тривиальны; смотрите Y.

  • Для numVars- размерные модели временных рядов, ar0 вектор ячейки из numVars- numVars числовые матрицы или numVars- размерный LagOp изолируйте полином оператора. Для векторов ячейки:

    • ar0 имеет длину p.

    • ar0 и ma0 каждый должен содержать numVars- numVars матрицы. Для каждой матрицы, строка k и столбец k соответствуют переменной k в системе k = 1, …, numVars.

    • Элементы ar0 соответствуйте изолированным ответам, которые составляют полином AR в обозначении разностного уравнения. Другими словами, ar0{j} матрица коэффициентов векторного yt-j, j = 1, …, p. Для всех содействующих матриц AR строка k содержит коэффициенты AR в уравнении переменной ykt и столбец, k содержит коэффициенты переменной ykt в рамках уравнений. Порядок строки и столбца всех коэффициентов авторегрессивного и скользящего среднего значения должен быть сопоставимым.

  • Для LagOp полиномы оператора задержки:

    • Коэффициенты в Coefficients свойство соответствует задержкам yt в Lags свойство.

    • Задайте модель в уменьшаемой форме путем предоставления идентичности для первого коэффициента (eye(numVars)).

    • armafevd составляет модель с помощью обозначения оператора задержки. Другими словами, когда вы работаете из модели в обозначении разностного уравнения, инвертируете коэффициенты AR изолированных ответов, чтобы создать эквивалентный полином оператора задержки.

      Например, рассмотреть yt=0.5yt10.8yt2+εt0.6εt1+0.08εt2. Модель находится в форме разностного уравнения. Чтобы вычислить FEVD, введите следующее в командной строке.

      y = armafevd([0.5 -0.8], [-0.6 0.08]);

      Модель ARMA, написанная в обозначении оператора задержки, (10.5L+0.8L2)yt=(10.6L+0.08L2)εt. Коэффициенты AR изолированных ответов отрицаются по сравнению с соответствующими коэффициентами в формате разностного уравнения. Чтобы получить тот же результат с помощью обозначения оператора задержки, введите следующее в командной строке.

      ar0 = LagOp({1 -0.5 0.8});
      ma0 = LagOp({1 -0.6 0.08});
      y = armafevd(ar0, ma0);

Коэффициенты скользящего среднего значения ARMA (p, q) модель в виде числового вектора, вектора ячейки из квадратных числовых матриц или LagOp изолируйте объект полинома оператора. Если ma0 вектор (числовой или ячейка), затем коэффициент εt является идентичностью (eye(numVars)).

Для модели AR задайте пустой массив или ячейку ([] или {}).

  • Для одномерных моделей временных рядов, ma0 числовой вектор, вектор ячейки из скаляров или одномерный LagOp изолируйте полином оператора. Для векторов, ma0 имеет длину q, и элементы соответствуют изолированным инновациям, которые составляют полином AR в обозначении разностного уравнения. Другими словами, ma0(j) или ma0{j} коэффициент εt-j, j = 1, …, q. Разложения отклонения одномерных моделей тривиальны; смотрите Y.

  • Для numVars- размерные модели временных рядов, ma0 вектор ячейки из числового numVars- numVars числовые матрицы или numVars- размерный LagOp изолируйте полином оператора. Для векторов ячейки:

    • ma0 имеет длину q.

    • ar0 и ma0 каждый должен содержать numVars- numVars матрицы. Для каждой матрицы, строка k и столбец k соответствуют переменной k в системе k = 1, …, numVars.

    • Элементы ma0 соответствуйте изолированным ответам, которые составляют полином MA в обозначении разностного уравнения. Другими словами, ma0{j} матрица коэффициентов εt-j, j = 1, …, q. Для всех содействующих матриц MA строка k содержит коэффициенты MA в уравнении переменной εkt и столбец, k содержит коэффициенты εkt в рамках уравнений. Порядок строки и столбца всех содействующих матриц авторегрессивного и скользящего среднего значения должен быть сопоставимым.

  • Для LagOp изолируйте полиномы оператора, коэффициенты в Coefficients свойство соответствует задержкам εt в Lags свойство.

    Чтобы задать модель в уменьшаемой форме, предоставьте идентичность (eye(numVars)) для коэффициента, который соответствует задержке 0.

Оси, на которых можно построить FEVD каждой переменной в виде вектора из Axes объекты с длиной равняются numVars.

По умолчанию, armafevd разложения отклонения графиков на осях на отдельных рисунках.

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

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

Пример: 'Method',"generalized",'NumObs',10 задает, чтобы вычислить обобщенный FEVD каждой переменной в течение 10 периодов.

Ковариационная матрица ARMA (p, q) инновации модели εt в виде разделенной запятой пары, состоящей из 'InnovCov' и числовой скаляр или numVars- numVars числовая матрица. InnovCov должна быть положительная скалярная величина или положительная определенная матрица.

Значением по умолчанию является eye(numVars).

Пример: 'InnovCov',0.2

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

Предскажите горизонт или количество периодов для который armafevd вычисляет FEVD в виде разделенной запятой пары, состоящей из 'NumObs' и положительное целое число. Другими словами, NumObs задает количество наблюдений, чтобы включать в FEVD (количество строк в Y).

По умолчанию, armafevd определяет NumObs критерием остановки mldivide.

Пример: 'NumObs',10

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

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

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

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

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

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

свернуть все

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

Yt JK) вклад в разложение отклонения переменной k относящийся к инновациям потрясают к переменной j во время t, для t = 1,2, …, numObsJ = 1,2..., numVars, и k = 1,2..., numVars. Столбцы и страницы Y соответствуйте переменному порядку в ar0 и ma0.

Для одномерных моделей, Y ones(numObs,1) потому что разложение отклонения один в течение каждого периода в горизонте прогноза.

Указатели на нанесенные на график графические объекты, возвращенные как numVars- numVars матрица графических объектов. h (jK) соответствует FEVD k относящийся к инновациям потрясают к переменной j во время 0.

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

Больше о

свернуть все

Обозначение разностного уравнения

Линейная модель временных рядов, написанная в difference-equation notation, располагает приведенную стоимость ответа и его структурного коэффициента на левой стороне уравнения. Правая сторона уравнения содержит сумму изолированных ответов, существующих инноваций, и изолировала инновации с соответствующими коэффициентами.

Другими словами, линейные временные ряды, написанные в обозначении разностного уравнения,

Φ0yt=c+Φ1yt1+...+Φpytp+Θ0εt+Θ1εt1+...+Θqεtq,

где

  • yt является numVars- размерный вектор, представляющий ответы numVars переменные во время t, для всего t и для numVars ≥ 1.

  • εt является numVars- размерный вектор, представляющий инновации во время t.

  • Φj является numVars- numVars матрица коэффициентов AR ответа yt-j, для j = 0..., p.

  • Θk является numVars- numVars матрица коэффициентов MA инноваций εt-k., k = 0..., q.

  • c является n - размерная константа модели.

  • Φ 0 = Θ 0 = I numVars, который является numVars- размерная единичная матрица, для моделей в уменьшаемой форме.

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

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

Предположим, что yt является ARMA (p, q) модель, содержащая numVars переменные отклика

Φ(L)yt=Θ(L)εt.

  • Φ (L) является полиномом оператора задержки авторегрессивных коэффициентов, другими словами, Φ(L)=Φ0Φ1LΦ2L2...ΦpLp.

  • Θ (L) является полиномом оператора задержки коэффициентов скользящего среднего значения, другими словами, Θ(L)=Θ0+Θ1L+Θ2L2+...+ΘqLq.

  • εt является вектором из numVars- D серия инноваций. Примите, что инновации имеют нулевое среднее значение и постоянную, положительно-определенную ковариационную матрицу Σ для всего t.

Представление MA бесконечной задержки yt

yt=Φ1(L)Θ(L)εt=Ω(L)εt.

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

γmjk=t=0m1(ekCtej)2t=0m1ekΩtΣΩtek.

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

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

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

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

Изолируйте обозначение оператора

Модель временных рядов, написанная в lag operator notation, располагает p - полином оператора задержки степени на существующем ответе на левой стороне уравнения. Правая сторона уравнения содержит константу модели и q - полином оператора задержки степени на существующих инновациях.

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

Φ(L)yt=c+Θ(L)εt,

где

  • yt является numVars- размерный вектор, представляющий ответы numVars переменные во время t, для всего t и для numVars ≥ 1.

  • Φ(L)=Φ0Φ1LΦ2L2...ΦpLp, который является авторегрессивным, полиномом оператора задержки.

  • L является оператором подготовительной смены, другими словами, Ljyt=ytj.

  • Φj является numVars- numVars матрица коэффициентов AR ответа yt-j, для j = 0..., p.

  • εt является numVars- размерный вектор, представляющий инновации во время t.

  • Θ(L)=Θ0+Θ1L+Θ2L2+...+ΘqLq, который является скользящим средним значением, полиномом оператора задержки.

  • Θk является numVars- numVars матрица коэффициентов MA инноваций εt-k., k = 0..., q.

  • c является numVars- размерная константа модели.

  • Φ 0 = Θ 0 = I numVars, который является numVars- размерная единичная матрица, для моделей в уменьшаемой форме.

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

Для получения дополнительной информации об обозначении оператора задержки смотрите Обозначение Оператора Задержки.

Советы

  • Вмещать структурный ARMA (p, q) модели, LagOp предоставления изолируйте полиномы оператора для входных параметров ar0 и ma0. Задавать структурный коэффициент, когда вы вызываете LagOp, установите соответствующую задержку на 0 при помощи 'Lags' аргумент пары "имя-значение".

  • Для ортогонализируемого многомерного FEVDs расположите переменные согласно Wold causal ordering [3]:

    • Первая переменная (соответствие первой строке и столбцу обоих ar0 и ma0) скорее всего, окажет мгновенное влияние (t = 0) на всех других переменных.

    • Вторая переменная (соответствие второй строке и столбцу обоих ar0 и ma0) скорее всего, окажет мгновенное влияние на остающиеся переменные, но не первую переменную.

    • В общем случае переменная j (соответствующий строке j и столбец j обоих ar0 и ma0) наиболее вероятное должно оказать мгновенное влияние на последний numVars – переменные j, но не предыдущие переменные j - 1.

Алгоритмы

  • armafevd FEVDs графиков только, когда это не возвращает выходных аргументов или h.

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

    Если Method "generalized"то:

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

    • Получившийся FEVD не основан на ортогональном преобразовании.

    • Получившийся FEVD переменной суммирует к одному единственному когда InnovCov диагональный [4].

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

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

Ссылки

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

[2] Lütkepohl, H. "Асимптотические Распределения Функций Импульсной характеристики и Разложения Отклонения Ошибки прогноза Векторных Авторегрессивных Моделей". Анализ Экономики и Статистики. Издание 72, 1990, стр 116–125.

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

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

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

| | |

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