Симуляция Монте-Карло модели векторного исправления ошибок (VEC)
Y = simulate(Mdl,numobs)
Y = simulate(Mdl,numobs,Name,Value)
[Y,E] =
simulate(___)
дополнительные опции использования заданы одним или несколькими аргументами пары "имя-значение". Например, Y
= simulate(Mdl
,numobs
,Name,Value
)'NumPaths',1000,'X',X
задает симуляцию 1 000 путей и X
как внешние данные о предикторе для компонента регрессии.
simulate
выполняет условную симуляцию с помощью этого процесса для всех страниц k
= 1..., numpaths
и в течение каждого раза t
= 1..., numobs
.
simulate
выводит (или обратные фильтры) инновации
от известных будущих ответов E(t,:,k)
. Для YF(t,:,k)
E(t,:,k)
simulate
подражает шаблону значений NaN
, который появляется в
.YF(t,:,k)
Для недостающих элементов E(t,:,k)
simulate
выполняет эти шаги.
Чертите Z1
, случайное, стандартное условное выражение воздействий Распределения Гаусса на известных элементах
.E(t,:,k)
Масштабируйте Z1
нижним треугольным Фактором Холесского условной ковариационной матрицы. Таким образом, Z2
= L*Z1
, где L
= chol(C,'lower')
и C
являются ковариацией условного Распределения Гаусса.
Припишите Z2
вместо соответствующих отсутствующих значений в
.E(t,:,k)
Для отсутствующих значений в YF(t,:,k)
simulate
пропускает соответствующие случайные инновации через модель Mdl
.
simulate
использует этот процесс, чтобы определить источник времени t 0 из моделей, которые включают линейные тренды времени.
Если вы не задаете Y0
, то t 0 = 0.
В противном случае simulate
устанавливает t 0 на size(Y0,1)
– Mdl.P
. Поэтому временами в компоненте тренда является t = t 0 + 1, t 0 + 2..., t 0 + numobs
. Это соглашение сопоставимо с поведением по умолчанию образцовой оценки, по которой estimate
удаляет первые ответы Mdl.P
, уменьшая эффективный объем выборки. Несмотря на то, что simulate
явным образом использует первые преддемонстрационные ответы Mdl.P
в Y0
, чтобы инициализировать модель, общее количество наблюдений в Y0
(исключая любые отсутствующие значения) определяет t 0.
[1] Гамильтон, J. D. Анализ timeseries. Принстон, NJ: Издательство Принстонского университета, 1994.
[2] Йохансен, S. Основанный на вероятности вывод в векторных авторегрессивных моделях Cointegrated. Оксфорд: Издательство Оксфордского университета, 1995.
[3] Juselius, K. Модель VAR Cointegrated. Оксфорд: Издательство Оксфордского университета, 2006.
[4] Lütkepohl, H. Новое введение в несколько анализ временных рядов. Берлин: Спрингер, 2005.