Это основной вопрос о моделях Box-Jenkins MA. Как я понимаю, модель MA - это, в основном, линейная регрессия значений временного ряда относительно предыдущих слагаемых ошибок . То есть наблюдение сначала регрессирует к своим предыдущим значениям а затем одно или несколько значений используются в качестве терминов ошибки для MA модель.
Но как вычисляются члены ошибки в модели ARIMA (0, 0, 2)? Если модель MA используется без авторегрессионной детали и, следовательно, без расчетного значения, как я могу иметь ошибку?
regression
time-series
arima
box-jenkins
Роберт Кубрик
источник
источник
Ответы:
Оценка модели MA:
Давайте предположим ряд с 100 временными точками, и скажем, что это характеризуется моделью MA (1) без перехвата. Тогда модель задается
Термин ошибки здесь не соблюдается. Таким образом, чтобы получить это, Box et al. Анализ временных рядов: прогнозирование и контроль (3-е издание) , стр. 228 , позволяют предположить, что термин ошибки вычисляется рекурсивно
Таким образом, термин ошибки для : Теперь мы не можем вычислить это, не зная значения . Таким образом, чтобы получить это, нам нужно вычислить начальную или предварительную оценку модели, см. Box et al. упомянутой книги, раздел 6.3.2 стр. 202 утверждают, что,ε 1 = y 1 + θ ε 0 θt=1
Обратите внимание, что является предполагаемой автокорреляцией. Более подробное обсуждение приведено в разделе 6.3 «Начальные оценки параметров» . Теперь предположим, что мы получаем начальную оценку . Тогда Теперь другая проблема - у нас нет значения для потому что начинается с 1, и поэтому мы не можем вычислить . К счастью, есть два способа получить это,rk θ=0.5
Согласно Box и соавт. В разделе 7.1.3 стр. 227 значения могут быть заменены на ноль в качестве приближения, если является умеренным или большим, этот метод является условным правдоподобием. В противном случае используется безусловное правдоподобие, при котором значение получают путем обратного прогнозирования, Box et al. рекомендую этот метод. Подробнее об обратном прогнозировании читайте в разделе 7.1.4 на стр. 231 .ε0 n ε0
После получения начальных оценок и значения , наконец, мы можем приступить к рекурсивному вычислению члена ошибки. Затем последним этапом является оценка параметра модели , помните, что это уже не предварительная оценка.ε0 (1)
При оценке параметра я использую процедуру нелинейного оценивания, в частности алгоритм Левенберга-Марквардта, поскольку модели МА нелинейны по своему параметру.θ
В целом, я очень рекомендую вам прочитать Box et al. Анализ временных рядов: прогнозирование и контроль (3-е издание) .
источник
Гауссова модель MA (q) определяется (не только Боксом и Дженкинсом!) Как поэтому модель MA (q) является «чистой» моделью ошибок, степень определяющая, насколько далеко корреляция возвращается.
источник
источник
Смотрите мой пост здесь для объяснения того, как понимать термины нарушения в серии МА.
Вам нужны разные методы оценки, чтобы оценить их. Это связано с тем, что вы не можете сначала получить остатки линейной регрессии, а затем включить в качестве объясняющих переменных отстающие остаточные значения, поскольку процесс MA использует остатки текущей регрессии. В вашем примере вы создаете два уравнения регрессии и используете невязки из одного в другое. Это не то, чем является процесс МА. Это не может быть оценено с помощью OLS.
источник