Предскажите модель VAR Используя симуляцию Монте-Карло

Этот пример показывает, как использовать симуляцию Монте-Карло через simulate, чтобы предсказать модель VAR.

simulate позволяет вам сгенерировать симуляции временных рядов на основе вашей модели. Если у вас есть защищенный объект модели VAR, можно использовать эти симуляции в качестве демонстрационных прогнозов.

simulate требует:

  • Модель (EstMdl в дальнейшем)

  • Количество периодов для прогноза (numobs в дальнейшем)

simulate опционально берет:

  • Внешний ряд данных

  • Преддемонстрационные временные ряды (Y(end-3:end,:) в дальнейшем)

  • Будущие демонстрационные ответы для условной симуляции

  • Количество реализации или пути, чтобы моделировать (2000 в дальнейшем)

Загрузите набор данных Data_USEconModel. Этот пример использует два временных рядов: логарифм действительного GDP и действительного 3-месячного уровня Казначейского векселя, оба differenced, чтобы быть приблизительно стационарным. Для рисунка модель VAR (4) описывает временные ряды.

load Data_USEconModel
DEF = log(DataTable.CPIAUCSL);
GDP = log(DataTable.GDP);
rGDP = diff(GDP - DEF); % Real GDP is GDP - deflation
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);

Задайте горизонт прогноза.

numobs = 21;
FDates = dateshift(DataTable.Time(end),'end','quarter',1:numobs);

Моделируйте модель для шагов numobs вперед и сгенерируйте 2 000 путей. Задайте преддемонстрационные наблюдения от конца данных.

rng(1); %For reproducibility
Ysim = simulate(EstMdl,numobs,'Y0',Y(end-3:end,:),'NumPaths',2000);

Вычислите среднее и стандартное отклонение моделируемого ряда:

Ymean = mean(Ysim,3); % Calculate means
Ystd = std(Ysim,0,3); % Calculate std deviations

Постройте средние значения + стандартное отклонение/-1 для моделируемого ряда:

figure;
subplot(2,1,1)
plot(DataTable.Time(end-10:end),Y(end-10:end,1),'k')
hold('on')
plot([DataTable.Time(end) FDates],[Y(end,1);Ymean(:,1)],'r')
plot([DataTable.Time(end) FDates],[Y(end,1);Ymean(:,1)]+[0;Ystd(:,1)],'b')
plot([DataTable.Time(end) FDates],[Y(end,1);Ymean(:,1)]-[0;Ystd(:,1)],'b')
title('Transformed real GDP')
subplot(2,1,2)
plot(DataTable.Time(end-10:end),Y(end-10:end,2),'k')
hold('on')
plot([DataTable.Time(end) FDates],[Y(end,2);Ymean(:,2)],'r')
plot([DataTable.Time(end) FDates],[Y(end,2);Ymean(:,2)]+[0;Ystd(:,2)],'b')
plot([DataTable.Time(end) FDates],[Y(end,2);Ymean(:,2)]-[0;Ystd(:,2)],'b')
title('Transformed real 3-mo T-bill rate')

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

Объекты

Функции

Похожие темы

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