ARMA

Na análise estatística de séries temporais, modelos auto-regressivos de médias móveis (autoregressive-moving-average ou ARMA, na sigla em inglês) oferecem uma descrição parcimoniosa de um processo estocástico fracamente estacionário em termos de dois polinômios, um para a auto-regressão e outro para a média móvel. O modelo ARMA geral foi descrito pelo matemático neo-zelandês Peter Whittle em sua tese de 1951, Hypothesis testing in time series analysis, e popularizado pelos estatísticos britânicos George E. P. Box e Gwilym Jenkins em seu livro de 1970.

Dada uma série temporal de dados $X_t$, o modelo ARMA é uma ferramenta para entender e, talvez, prever valores futuros nesta série. O modelo consiste em duas partes, uma parte auto-regressiva (AR) ou uma parte de média móvel (MA). A parte AR envolve regressar a variável em seus próprios valores defasados, isto é, passados. A parte MA envolve modelar o termo de erro como uma combinação linear de termos de erro que ocorrem contemporaneamente e em vários momentos no passado.

O modelo é geralmente referido como o modelo ARMA($p,q$), em que $p$ é a ordem da parte auto-regressiva e $q$ é a ordem da parte de média móvel.

Modelos ARMA podem ser estimados seguindo a abordagem de Box–Jenkins.

Modelo auto-regressivo

A notação AR(

$p$

) se refere ao modelo auto-regressivo de ordem

$p$

. O modelo AR(

$p$

) é escrito como:

$X_t=c+\sum _{i=1}^p{\phi }_i{X}_{t-i}+{\epsilon }_t,$

em que

${\phi }_1...{\phi }_p$

são os parâmetros,

$c$

é uma constante e a variável aleatória

${\epsilon }_t$

é ruído branco.

Algumas restrições são necessárias sobre os valores dos parâmetros a fim de que o modelo permaneça estacionário. Por exemplo, processos no modelo AR(1) com $|{\phi }_1|\ge 1$ não são estacionários.

Modelo de média móvel

A notação MA(

$q$

) se refere ao modelo de média móvel de ordem

$q$

:

$X_t=\mu +{\epsilon }_t+\sum _{i=1}^q{\theta }_i{\epsilon }_{t-i},$

em que

${\theta }_1,...,{\theta }_q$

são os parâmetros do modelo,

$\mu$

é o valor esperado de

$X_t$

(frequentemente assumido como igual a zero) e

${ϵ}_t,{ϵ}_{t-1},...$

são novamente os termos de erro de ruído branco.

Modelo ARMA

A notação ARMA(

$p,q$

) se refere ao modelo com termos auto-regressivos

$p$

e termos de média móvel

$q$

. Este modelo contém os modelos AR(

$p$

) e MA(

$q$

):

$X_t=c+{\epsilon }_t+\sum _{i=1}^p{\phi }_i{X}_{t-i}+\sum _{i=1}^q{\theta }_i{\epsilon }_{t-i}.$

O modelo ARMA geral foi descrito por Peter Whittle em sua tese de 1951, em que usou análise matemática (séries de Laurent e análise de Fourier) e inferência estatística. Modelos ARMA foram popularizados por George E. P. Box e Gwilym Jenkins em um livro de 1970, em que expuseram um método iterativo (Box–Jenkins) para escolher e estimar estes modelos. Este método foi útil para polinômios de ordem baixa (grau três ou menos).

Termos de erro

Os termos de erro ${ϵ}_t$ são geralmente assumidos como variáveis aleatórias independentes e identicamente distribuídas amostradas a partir de uma distribuição normal com média zero, ${ϵ}_t\sim N(0,{\sigma }^2)$, em que ${\sigma }^2$ é a variância. Estes pressupostos podem ser enfraquecidos, mas com isto as propriedades do modelo são alteradas. Em particular, uma mudança no pressuposto da independência e da distribuição idêntica faz uma diferença bastante fundamental.

Especificação em termos de operador de defasagem

Em alguns textos, os modelos serão especificados em termos do operador de defasagem

$L$

. Nestes termos, o modelo AR(

$p$

) é dado por:

${\epsilon }_t=(1-\sum _{i=1}^p{\phi }_i{L}^i)X_t=\phi (L)X_t,$

em que

$\phi$

representa o polinômio:

$\phi (L)=1-\sum _{i=1}^p{\phi }_i{L}^i.$

O modelo MA(

$q$

) é dado por:

$X_t=(1+\sum _{i=1}^q{\theta }_i{L}^i){\epsilon }_t=\theta (L){\epsilon }_t,$

em que

$\theta$

representa o polinômio:

$\theta (L)=1+\sum _{i=1}^q{\theta }_i{L}^i.$

Finalmente, o modelo combinado ARMA(

$p,q$

) é dado por

$(1-\sum _{i=1}^p{\phi }_i{L}^i)X_t=(1+\sum _{i=1}^q{\theta }_i{L}^i){\epsilon }_t,$

ou, de forma mais concisa,

$\phi (L)X_t=\theta (L){\epsilon }_t,$

ou

$\frac{\phi (L)}{\theta (L)}{X}_t={\epsilon }_t.$

Notação alternativa

Alguns autores, incluindo George E. P. Box, Gwilym M. Jenkins e Gregory C. Reinsel, usam uma convenção diferente para os coeficientes de auto-regressão. Isto permite a todos os polinômios que envolvem o operador de defasagem aparecer de forma semelhante por toda a parte. Assim, o modelo ARMA é escrito como:

$(1-\sum _{i=1}^p{\varphi }_i{L}^i)X_t=(1+\sum _{i=1}^q{\theta }_i{L}^i){\epsilon }_t.$

Além disso, se configurarmos

${\varphi }_0={\theta }_0=1$

, então, obtemos uma formulação ainda mais sofisticada:

$\sum _{i=0}^p{\varphi }_i{L}^i{X}_t=\sum _{i=0}^q{\theta }_i{L}^i{\epsilon }_t.$

Modelos de ajuste

Escolha de $p$ e $q$

A descoberta dos valores apropriados de $p$ e $q$ no modelo ARMA($p,q$) pode ser facilitada ao mapear as funções de autocorrelação parcial para um valor estimado de $p$ e, da mesma forma, usar as funções de autocorrelação para um valor estimado de $q$. Mais informações podem ser obtidas considerando as mesmas funções para os resíduos de um modelo ajustado com uma seleção inicial de $p$ e $q$.

Peter J. Brockwell e Richard A. Davis recomendaram o uso do critério de informação de Akaike para encontrar $p$ e $q$.

Estimação dos coeficientes

Em geral, modelos ARMA não podem ser, depois de escolher $p$ e $q$, ajustados por regressão de mínimos quadrados para encontrar os valores dos parâmetros que minimizam o termo de erro. Geralmente, considera-se uma boa prática encontrar os menores valores de $p$ e $q$ que oferecem um ajuste aceitável para os dados. Para um modelo AR puro, as equações de Yule–Walker podem ser usadas para oferecer um ajuste.

Aplicações

O modelo ARMA é apropriado quando um sistema é uma função de uma série de choques não observados (a parte da média móvel ou MA), assim como seu comportamento. Por exemplo, preços de ações podem ser chocados por informação fundamental, assim como pela exibição de tendências técnicas e efeitos de reversão à média devido aos participantes do mercado.

Generalizações

A dependência de $X_t$ de valores passados e os termos de erro ${\epsilon }_t$ são assumidos como lineares, a não ser quando especificado o contrário. Se a dependência for não linear, o modelo é especificamente chamado de modelo de médias móveis não lineares (NMA), modelo auto-regressivo não linear (NAR) ou modelo auto-regressivo de médias móveis não lineares (NARMA).

Os modelos auto-regressivos de médias móveis podem ser generalizados de outras formas. Se séries temporais múltiplas devem ser ajustadas, um vetor ARIMA (ou VARIMA) pode ser ajustado. Se a série temporal em questão exibe memória longa, então, a modelagem ARIMA fracionária (FARIMA ou ARFIMA) pode ser apropriada. Se os dados contêm efeitos sazonais, pode ser modelada por um modelo ARIMA sazonal (SARIMA) ou um modelo ARIMA periódico.

Outra generalização é o modelo auto-regressivo multiescalar (MAR). Um modelo MAR é indexado pelos nós de uma árvore, enquanto um modelo auto-regressivo padrão (de tempo discreto) é indexado por números inteiros.

Perceba que o modelo ARMA é um modelo univariado. Extensões para o caso multivariado são a auto-regressão vetorial (VAR) e a auto-regressão vetorial de médias móveis (VARMA).

Modelo auto-regressivo de médias móveis com modelo de inputs exógenos (ARMAX)

A notação ARMAX(

$p,q,b$

) se refere ao modelo com

$p$

termos auto-regressivos,

$q$

termos de média móvel e

$b$

termos de input exógeno. Este modelo contém os modelos AR(

$p$

) e MA(

$q$

) e uma combinação linear dos últimos

$b$

termos de uma série temporal conhecida e externa

$d_t$

. É dado por:

$X_t={\epsilon }_t+\sum _{i=1}^p{\phi }_i{X}_{t-i}+\sum _{i=1}^q{\theta }_i{\epsilon }_{t-i}+\sum _{i=1}^b{\eta }_i{d}_{t-i}.$

em que

${\eta }_1,...,{\eta }_b$

são os parâmetros do input exógeno

$d_t$

.

Algumas variantes não lineares de modelos com variáveis exógenas foram definidas.

Pacotes estatísticos implementam o modelo ARMAX por meio do uso de variáveis "exógenas" ou "independentes". Deve-se tomar cuidado ao interpretar o output daqueles pacotes, porque os parâmetros estimados geralmente se referem à regressão:

$X_t-m_t={\epsilon }_t+\sum _{i=1}^p{\phi }_i(X_{t-i}-m_{t-i})+\sum _{i=1}^q{\theta }_i{\epsilon }_{t-i},$

em que

$m_t$

incorpora todas as variáveis exógenas (ou independentes):

$m_t=c+\sum _{i=0}^b{\eta }_i{d}_{t-i}.$

Ver também

• Modelos ARCH
• ARIMA
• Série temporal