Прогнозирование модели VAR, симуляция и анализ

Прогнозирование модели VAR

Когда у вас есть модели параметрами (известный или предполагаемый), можно исследовать предсказания моделей. Для получения информации о создании моделей VAR смотрите Векторную Авторегрессию (VAR) Создание Модели. Для получения информации об оценке моделей смотрите Оценку Модели VAR.

Этот список описывает основные методы прогнозирования.

Эти функции основывают свои прогнозы на полностью заданном объекте модели и исходных данных. Функции отличаются по своим инновационным процессам:

  • forecast принимает инновации с нулевым знаком. Поэтому forecast дает к детерминированному прогнозу, условному выражению или в противном случае.

  • simulate принимает, что многомерные инновации совместно Гауссовы распределенный с ковариационной матрицей Σ. simulate псевдослучайные выражения, демонстрационные пути Монте-Карло.

  • filter требует инновационных путей к процессу. filter дает демонстрационный путь, который детерминировано основан на заданных инновационных путях к процессу.

forecast быстрее и требует меньшей памяти, чем генерация многого демонстрационного использования путей simulate или filter. Однако forecast не так гибко как simulate и filter. Например, предположите, что вы преобразовываете некоторые временные ряды прежде, чем сделать модель и хотите отменить преобразование при исследовании прогнозов. Ошибочные границы, данные преобразованиями forecast ошибочные границы не являются допустимыми границами. В отличие от этого ошибочные границы, данные статистикой преобразованных симуляций, допустимы.

Как Функции для прогнозирования работы

Для безусловного прогнозирования, forecast генерирует два количества:

  • Детерминированные временные ряды прогноза на основе 0 инноваций

  • Временные ряды матриц среднеквадратичной погрешности прогноза на основе Σ, инновационной ковариационной матрицы.

Для условного прогнозирования:

  • forecast требует массива будущих данных об ответе, которые содержат соединение пропавших без вести (NaN) и известные значения. forecast генерирует прогнозы для условного выражения отсутствующих значений на известных значениях.

  • Прогнозы сгенерированы forecast также детерминированы, но матрицы среднеквадратичной погрешности основаны на Σ и значениях имеющегося отклика в горизонте прогноза.

  • forecast использует Фильтр Калмана, чтобы сгенерировать прогнозы. В частности:

    1. forecast представляет модель VAR как модель в пространстве состояний (ssm объект модели) без ошибки наблюдения.

    2. forecast пропускает данные о прогнозе через модель в пространстве состояний. Таким образом, в период t в горизонте прогноза любой неизвестный ответ

      y^t=Φ^1y^t1+...+Φ^py^tp+c^+δ^t+β^xt,

      где y^s, s <t, отфильтрованная оценка y с периода s в горизонте прогноза. forecast преддемонстрационные значения использования в течение периодов перед горизонтом прогноза.

    Для получения дополнительной информации смотрите filter и [4], стр 612 и 615.

Для любого типа прогноза, Чтобы инициализировать модель VAR (p) в горизонте прогноза, forecast требует преддемонстрационных наблюдений p. Можно опционально задать больше чем один путь преддемонстрационных данных. Если вы действительно задаете разнообразные пути, forecast возвращает 3D массив предсказанных ответов, с каждой страницей, соответствующей пути преддемонстрационных значений.

Для безусловной симуляции, simulate:

  1. Генерирует случайные временные ряды на основе модели с помощью случайных путей многомерных Гауссовых инноваций, распределенных со средним значением нуля и ковариацией Σ

  2. Пропускает случайные пути инноваций через модель

Для условной симуляции:

  • simulate, как forecast, требует массива будущих данных об ответе, которые содержат соединение того, чтобы избегать и известных значений, и генерируют значения для недостающих ответов.

  • simulate выполняет условную симуляцию с помощью этого процесса. В каждый раз t в горизонте прогноза:

    1. simulate выводит (или, обратные фильтры) инновации (E (t,:)) от известных будущих ответов.

    2. Для недостающих будущих инноваций, simulate:

      1. Чертит Z1, который является случайным, стандартным условным выражением воздействий Распределения Гаусса на известных элементах E (t,:).

      2. Шкалы Z1 нижним треугольным Фактором Холецкого условной ковариационной матрицы. Таким образом, Z2 = L*Z1, где L = chol(Covariance,'lower') и Covariance ковариация условного Распределения Гаусса.

      3. Приписывает Z2 вместо соответствующих отсутствующих значений в E (t,:).

    3. Для отсутствующих значений в будущих данных об ответе, simulate пропускает соответствующие случайные инновации через модель VAR Mdl.

Для любого типа симуляции:

  • simulate не требует преддемонстрационных наблюдений. Для получения дополнительной информации на значениях по умолчанию преддемонстрационных данных, смотрите 'Y0'.

  • Чтобы выполнить вывод, сгенерируйте 1000-е путей к ответу, и затем оцените демонстрационную статистику от сгенерированных путей каждый раз в горизонте прогноза. Например, предположите Y 3D массив предсказанных путей. Точка Монте-Карло и оценки интервала прогноза во время t в прогнозе горизонт

    MCPointEst = mean(Y(t,:,:),3);
    MCPointInterval = quantile(Y(t,:,:),[0.025 0.975],3);

Таким образом, точечная оценка Монте-Карло является средним значением через страницы, и оценка интервала Монте-Карло состоит из 2.5th и 97.5th процентили, вычисленные через пути. Заметьте, что оценки Монте-Карло подвергаются ошибке Монте-Карло, и таким образом, оценки отличаются каждый раз, когда вы запускаете анализ при тех же условиях, но использование различного seed случайных чисел.

Масштабирование данных

Если вы масштабировали какие-либо временные ряды прежде, чем подобрать модель, можно не масштабировать получившиеся временные ряды, чтобы изучить его предсказания более легко.

  • Если вы масштабировали ряд с log, преобразуйте предсказания соответствующей модели с exp.

  • Если вы масштабировали ряд с diff(log) или, эквивалентно, price2ret, преобразуйте предсказания соответствующей модели с cumsum(exp), или, эквивалентно, ret2price. cumsum инверсия diff; это вычисляет совокупные суммы. Как в интегрировании, необходимо выбрать соответствующую аддитивную постоянную за совокупную сумму. Например, возьмите журнал итоговой записи в соответствующем ряду данных и используйте его в качестве первого термина в ряду перед применением cumsum.

Вычисление импульсных характеристик

Можно исследовать эффект impulse responses к моделям с armairf. Импульсная характеристика является детерминированным ответом модели временных рядов к инновационному процессу, который имеет значение одного стандартного отклонения в одном компоненте в начальное время и нули во всех других компонентах и времена. Основным компонентом функции импульсной характеристики является dynamic multipliers, то есть, коэффициенты представления VMA модели VAR.

Учитывая полностью заданный varm модель, необходимо предоставить коэффициенты авторегрессии к armairf. По умолчанию, armairf отправляет модульный шок через систему, которая приводит к forecast error impulse response. Можно опционально предоставить инновационную ковариационную матрицу и выбрать, сгенерировать ли импульсные характеристики orthogonalized или generalized. Обобщенные импульсные характеристики составляют фильтрацию шока одной стандартной погрешности каждых инноваций хотя модель VAR. Ортогонализируемые импульсные характеристики масштабируют динамические множители нижним треугольным Фактором Холецкого инновационной ковариации. Для получения дополнительной информации см. [2].

Для примера смотрите, Генерируют Импульсные характеристики Модели VAR.

Ссылки

[1] Lütkepohl, H. Новое введение в несколько анализ временных рядов. Берлин: Спрингер, 2005.

[2] Pesaran, H. H. и И. Шин. “Обобщенный анализ импульсной характеристики в линейных многомерных моделях”. Экономические буквы. Издание 58, 1998, 17-29.

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

Объекты

Функции

Похожие темы