Сгенерируйте импульсные характеристики модели VAR
В этом примере показано, как сгенерировать импульсные характеристики шока процентной ставки на США действительный валовой внутренний продукт (ВВП) с помощью armairf.
Загрузите Data_USEconModel набор данных. Этот пример использует два ежеквартальных временных рядов: логарифм действительного GDP и действительного трехмесячного уровня Казначейского векселя, оба differenced, чтобы быть приблизительно стационарным. Предположим, что векторная модель авторегрессии, содержащая первые четыре последовательных задержки (VAR (4)), является соответствующей, чтобы описать временные ряды.
load Data_USEconModel DEF = log(DataTable.CPIAUCSL); GDP = log(DataTable.GDP); rGDP = diff(GDP - DEF); % Real GDP is GDP - DEF TB3 = 0.01*DataTable.TB3MS; dDEF = 4*diff(DEF); % Scaling rTB3 = TB3(2:end) - dDEF; % Real interest is deflated Y = [rGDP,rTB3];
Подбирайте модель VAR (4) к данным.
Mdl = varm(2,4);
Mdl.SeriesNames = {'Transformed real GDP','Transformed real 3-mo T-bill rate'};
EstMdl = estimate(Mdl,Y);Вычислите и постройте ортогонализируемые импульсные характеристики в одном armairf вызов путем дополнительного возврата указателя на нанесенные на график графические объекты (второй выход). Задайте предполагаемую инновационную ковариационную матрицу.
[YIRF,h] = armairf(EstMdl.AR,{},'InnovCov',EstMdl.Covariance);
YIRF 15 массивом 2 на 2 импульсных характеристик. Каждая строка соответствует времени в горизонте прогноза (0..., 14), столбец j соответствует шоку для переменной Y(:,j) во время 0 и страницу k содержит импульсные характеристики переменной Y(:,k).
Получите выкачанный ряд импульсной характеристики GDP, следующий из шока для ряда процентной ставки путем отмены масштабирования действительного ряда импульсной характеристики GDP.
defGDPImp = ret2price(YIRF(:,2,1));
Постройте выкачанный ряд импульсной характеристики GDP по исходной шкале.
figure; numobs = size(defGDPImp,1); plot(0:numobs - 1,defGDPImp); title('Impact of Shock to Interest Rate on Deflated GDP') ylabel('Deflated GDP'); xlabel('Time');
График показывает, что увеличенная процентная ставка сопоставлена с маленьким падением в выкачанной серии GDP в течение короткого времени. Затем действительный GDP начинает подниматься снова, достигнув его прежнего уровня приблизительно через три года.