Этот пример показывает, как сгенерировать импульсные ответы шока процентной ставки на США действительный валовой внутренний продукт (ВВП) с помощью 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 начинает подниматься снова, достигнув его прежнего уровня приблизительно через три года.