Краткий аудит программного кода квартальных прогнозных моделей от ЦБ

Всем привет, данная работа по воспроизведению квартальной прогнозных моделей ЦБ на основе программного кода долго лежала в бэклоге, так как является трудоемкой, а с появлением LLM она стала возможной. 


Начнем с  kpm_dip.rar - Программный код квартальной прогнозной модели Департамента исследований и прогнозирования

Содержание:

Очень старая и жёсткая привязка среды

В инструкции прямо требуются:

• MATLAB 2014;
• IRIS 2014;
• MiKTeX;
• размещение IRIS и модели на диске

Модель представляет собой осмысленную линейную semi-structural QPM с:

• гибридной кривой Филлипса;
• эндогенными инфляционными ожиданиями;
• разрывом выпуска;
• процентным каналом;
• валютным каналом;
• бюджетным блоком;
• кредитным и финансовым блоком;
• внешним сектором.

Но архив:

• не содержит базовую инфляцию как отдельный показатель;
• не рассчитывает инфляцию из корзины Росстата;
• не содержит фактических данных;
• не воспроизводит официальный прогноз ЦБ;
• содержит как минимум две вероятные технические ошибки в российской ветке запуска.

Далее смотрим kpm_ddkp.rar - Программный код квартальной прогнозной модели Департамента денежно-кредитной политики

Архив КПМ ДДКП действительно содержит модель с отдельной базовой инфляцией. В этом смысле он гораздо ближе к вашему первоначальному воспоминанию, чем kpm_dip. 

Как базовая инфляция используется в решении по ставке

ЦБ не подставляет в правило только headline CPI. Сначала формируется специальная ожидаемая инфляция:

<code>E3_PIE4_MP =
    w_mp_pie_c*E3_PIE_C4
    + (1-w_mp_pie_c)*E3_PIE4;</code>

Здесь:

<code>w_mp_pie_c = 0,75</code>

То есть в индикаторе для денежно-кредитного правила:

• 75% — ожидаемая базовая инфляция;
• 25% — ожидаемая общая инфляция.

Правило ставки:

<code>RS =
    gamma1*RS{-1}
    + (1-gamma1)*(
        RS_NEUTRAL
        + gamma2*E3_PIE4_DEV
        + gamma3*Y_GAP
      )
    + RES_RS;</code>

Базовые параметры:

<code>сглаживание ставки gamma1 = 0,75
реакция на инфляцию gamma2 = 1,50
реакция на разрыв выпуска gamma3 = 0,50</code>

В более жёстком сценарии:

<code>gamma2 = 1,90</code>

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

Архив:

• моделирует динамику уже подготовленного индекса базовой инфляции;
• разделяет headline CPI на базовую и несколько небазовых компонент;
• строит импульсные отклики;
• не содержит исторических данных, перечня товаров Росстата, процедуры сезонной корректировки, оценки параметров или воспроизведения официального прогноза. 

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

Полностью запустить модель в текущем окружении не удалось: необходимы MATLAB 2019a+ и IRIS Toolbox 2015+, а они отсутствуют. При этом я распаковал архив, раскрыл макросы IRIS, проверил структуру уравнений, параметры и приложенные результаты MATLAB.