armairf

Сгенерируйте или постройте импульсные характеристики модели ARMA

Описание

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

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

Объект моделиФункция IRF
arimaimpulse
regARIMAimpulse
varmirf
vecmirf

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

пример

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

armairf функция:

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

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

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

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

пример

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

пример

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

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

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

Примеры

свернуть все

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

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

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

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

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

armairf(AR0,MA0);

Импульсная характеристика исчезает после четырех периодов.

В качестве альтернативы создайте модель ARMA, которая представляет yt. Задайте 1 для отклонения инноваций и никакой константы модели.

Mdl = arima('AR',AR0,'MA',MA0,'Variance',1,'Constant',0);

Mdl arima объект модели.

Постройте IRF использование Mdl.

impulse(Mdl);

impulse использует диаграмму стебель-листья, тогда как armairf использует линейный график. Однако IRFs в этих двух реализациях равны, потому что отклонение модели ARMA равняется 1.

Постройте целый обобщенный IRF одномерной модели ARMA(2,1)

(1-0.3L+0.1L2)yt=(1+0.05L)εt.

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

AR0Lag = LagOp([1 -0.3 0.1])
AR0Lag = 
    1-D Lag Operator Polynomial:
    -----------------------------
        Coefficients: [1 -0.3 0.1]
                Lags: [0 1 2]
              Degree: 2
           Dimension: 1
MA0Lag = LagOp([1 0.05])
MA0Lag = 
    1-D Lag Operator Polynomial:
    -----------------------------
        Coefficients: [1 0.05]
                Lags: [0 1]
              Degree: 1
           Dimension: 1

AR0Lag и MA0Lag LagOp изолируйте полиномы оператора, представляющие полиномы оператора задержки авторегрессивного и скользящего среднего значения, соответственно.

Постройте обобщенный IRF путем передачи в полиномах оператора задержки.

armairf(AR0Lag,MA0Lag,'Method','generalized');

IRF эквивалентен IRF в Графике Ортогонализируемый IRF Одномерной Модели ARMA.

Постройте целый IRF структурной векторной модели скользящего среднего значения авторегрессии (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);

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

figure;
armairf(VARLag,VMALag,'Method','generalized');

armairf возвращает три фигуры. Рисунок k содержит обобщенный IRF переменной k к шоку, применился ко всем другим переменным во время 0. Поскольку все IRFs исчезают после конечного числа периодов модель VARMA устойчива.

Вычислите целый ортогонализируемый IRF одномерной модели ARMA(2,1)

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

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

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

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

y = armairf(AR0,MA0)
y = 5×1

    1.0000
    0.3500
    0.0050
   -0.0335
   -0.0105

y вектор 5 на 1 из импульсных характеристик. y(1) импульсная характеристика в течение времени t=0, y(2) импульсная характеристика в течение времени t=1, и так далее. IRF исчезает после периода t=4.

В качестве альтернативы создайте модель ARMA, которая представляет yt. Задайте 1 для отклонения инноваций и никакой константы модели.

Mdl = arima('AR',AR0,'MA',MA0,'Variance',1,'Constant',0);

Mdl arima объект модели.

Постройте IRF модели ARIMA Mdl.

y = impulse(Mdl)
y = 5×1

    1.0000
    0.3500
    0.0050
   -0.0335
   -0.0105

IRFs в этих двух реализациях эквивалентны.

Вычислите обобщенный IRF 2D модели 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];

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

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

    31     2     2

Y(10,1,2)
ans = -0.0116

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

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

Вычислите и отобразите обобщенные импульсные характеристики в течение первых 10 периодов.

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

    0.7071   -0.2000
    0.7354   -0.3000
    0.2135   -0.1340
    0.0526   -0.0112
    0.2929   -0.0772
    0.3717   -0.1435
    0.1872   -0.0936
    0.0730   -0.0301
    0.1360   -0.0388
    0.1841   -0.0674


Y10(:,:,2) =

   -0.1414    0.5000
   -0.1131    0.1700
   -0.0509   -0.0040
    0.0058   -0.0113
    0.0040   -0.0003
   -0.0300    0.0100
   -0.0325    0.0133
   -0.0082    0.0054
   -0.0001   -0.0003
   -0.0116    0.0028

Y10 10 массивом 2 на 2 импульсных характеристик. Строки соответствуют временам 0 до 9 в горизонте прогноза.

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

Copyright 2018 The MathWorks, Inc.

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

свернуть все

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

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

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

  • Для 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)).

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

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

ar0 = [0.5 -0.8];
ma0 = [-0.6 0.08];
y = armairf(ar0,ma0);

Модель 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 = armairf(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.

  • Для 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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

свернуть все

Импульсные характеристики, возвращенные как числовой вектор-столбец или числовой массив.

Yt + 1, jK) импульсная характеристика переменной k к инновациям с одним стандартным отклонением потрясают к переменной j во время 0, для t = 0, 1..., numObs – 1, j = 1,2..., numVars, и k = 1,2..., numVars. Столбцы и страницы Y соответствуйте переменному порядку в ar0 и ma0.

Указатели на нанесенные на график графические объекты, возвращенные как numVars- numVars матрица графических объектов. h (jK) соответствует IRF переменной 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- размерная единичная матрица, для моделей в уменьшаемой форме.

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

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

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

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

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

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

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

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

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

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

ψj(m)=Cmej.

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

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

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

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

Модель временных рядов, написанная в 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 кажутся отрицаемыми относительно соответствующих условий в обозначении разностного уравнения. Знаки коэффициентов скользящего среднего значения являются тем же самым и появляются на той же стороне.

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

Советы

  • Чтобы вычислить forecast error impulse responses, используйте значение по умолчанию InnovCov, который является numVars- numVars единичная матрица. В этом случае, все доступные методы расчета (см. Method) приведите к эквивалентному IRFs.

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

  • Поскольку многомерный ортогонализировал IRFs, расположите переменные согласно Wold causal ordering [2]:

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

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

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

Алгоритмы

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

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

Вопросы совместимости

развернуть все

Поведение изменяется в R2018b

Поведение изменяется в R2018b

Ссылки

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

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

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

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