Этот пример показывает, как генерировать импульсные реакции на удар процентной ставки по реальному валовому внутреннему продукту (ВВП) США с использованием armairf.
Загрузить Data_USEconModel набор данных. В этом примере используются два квартальных временных ряда: логарифм реального ВВП и реальная трехмесячная ставка Т-векселя, которые различаются, чтобы быть приблизительно неподвижными. Предположим, что векторная модель авторегрессии, содержащая первые четыре последовательных лага (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).
Получение спущенного ряда импульсных откликов ВВП в результате шока к ряду процентных ставок путем отмены масштабирования реального ряда импульсных откликов ВВП.
defGDPImp = ret2price(YIRF(:,2,1));
Постройте график спущенных рядов импульсных откликов ВВП по первоначальной шкале.
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');
График показывает, что повышенная процентная ставка связана с небольшим падением дефлированной серии ВВП на короткое время. Затем реальный ВВП снова начинает подниматься, достигнув прежнего уровня примерно через три года.