exponenta event banner

Монте-Карло моделирование регрессионных моделей с ошибками ARIMA

Что такое моделирование Монте-Карло?

Моделирование Монте-Карло - процесс генерации независимых случайных розыгрышей из заданной вероятностной модели. При моделировании моделей временных рядов один розыгрыш (или реализация) представляет собой полный путь выборки заданной длины N, y1, y2,...,yN. При создании большого числа розыгрышей, например М, генерируются M путей образца, каждый длиной N.

Примечание

Некоторые расширения моделирования Монте-Карло основаны на генерации зависимых случайных розыгрышей, таких как Markov Chain Monte Carlo (MCMC). simulate функция в Econometrics Toolbox™ генерирует независимые реализации.

Некоторые применения моделирования Монте-Карло:

  • Демонстрация теоретических результатов

  • Прогнозирование будущих событий

  • Оценка вероятности будущих событий

Создание путей образца Монте-Карло

Часть временных рядов модели определяет динамическую эволюцию процесса безусловного возмущения во времени через условную среднюю структуру. Чтобы выполнить моделирование моделей регрессии с ошибками ARIMA по Монте-Карло:

  1. Укажите инновации предварительной выборки или безусловные нарушения (или используйте данные предварительной выборки по умолчанию).

  2. Создание некоррелированного инновационного ряда из распределения вероятностей.

  3. Фильтрация нововведений с помощью модели ошибок ARIMA для получения смоделированных безусловных возмущений.

  4. Используйте регрессионную модель, данные предиктора и смоделированные безусловные возмущения для получения ответов.

Например, рассмотрим моделирование N ответов из регрессионной модели с ошибками ARMA (2,1):

yt = Xtβ + utut = δ 1ut 1 + δ 2ut 2 + αt + start1αt − 1,

где αt - гауссов со средним 0 и дисперсией start2. Учитывая предварительные безусловные нарушения (u0 и u-1) и нововведения (α0), следуя следующим шагам:

  1. Генерировать N независимых инноваций из гауссова дистрибутива:

    {ε^1,ε^2,...,ε^N}.

  2. Фильтрация нововведений рекурсивно для получения безусловных возмущений:

    1. u^1 =ϕ1u0 +ϕ2u−1 +ε^1 +ε0

    2. u ^ 2 = 1 1u ^ 2 + 1 2u0 +

    3. u ^ 3 = 2 1u ^ 1 + 3 2u ^ 2 +

    4. ...

    5. u^N =ϕ1u^N−1 +ϕ2u^N−2 +ε^N +ε^N−1.

  3. Получить смоделированные ответы с использованием безусловных возмущений, регрессионной модели и предикторов:

    y ^ t = Xtβ + u ^ t.

Econometrics Toolbox автоматизирует этот процесс с помощью simulate. Передача в полностью заданной регрессионной модели с ошибками ARIMA (regARIMA), количество ответов на моделирование и, необязательно, количество путей и данных предварительной выборки, и simulate моделирует ответы.

Примечание

Econometrics Toolbox рассматривает предикторы в регрессионной модели как фиксированные нестохастические ряды. Поэтому, чтобы генерировать пути выборки Монте-Карло ответа, нужно знать значения предикторов.

Ошибка Монте-Карло

Используя множество моделируемых путей, можно оценить различные особенности модели. Однако оценка Монте-Карло основана на конечном количестве симуляций. Поэтому оценки Монте-Карло подвержены некоторой погрешности. Можно уменьшить количество ошибок Монте-Карло в расчетном исследовании, увеличив число путей образца M, создаваемых на основе модели.

Например, чтобы оценить вероятность будущего события:

  1. Создайте M путей образцов из модели.

  2. Оцените вероятность будущего события, используя выборочную долю возникновения события в M-моделировании,

    p^ =# событие времен происходит в M drawsM.

  3. Рассчитайте стандартную ошибку Монте-Карло для оценки,

    se = p ^ (1 p ^) М.

Можно уменьшить ошибку Монте-Карло оценки вероятности, увеличив число реализаций. Если вы знаете требуемую точность оценки, вы можете решить количество реализаций, необходимых для достижения этого уровня точности.

См. также

|

Связанные примеры

Подробнее