Mostrar código-fonte de ARMA

{{Minidesambig|a série de jogos eletrônicos|ARMA (série)|}}

Na análise [[estatística]] de [[Série temporal|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.<ref>{{Citar livro|url=https://books.google.com.br/books?id=FubniGJ0ECQC&printsec=frontcover&dq=Spectral+Analysis+for+Physical+Applications&hl=pt-BR&sa=X&ved=0ahUKEwiZ6Zv0w7nVAhWFFZAKHfCYBJwQ6AEIJzAA#v=onepage&q=Spectral%20Analysis%20for%20Physical%20Applications&f=false|título=Spectral Analysis for Physical Applications|ultimo=Percival|primeiro=Donald B.|ultimo2=Walden|primeiro2=Andrew T.|data=1993-06-03|editora=Cambridge University Press|lingua=en|isbn=9780521435413}}</ref> 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 <math>X_t</math>, 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.<ref>{{Citar livro|url=https://books.google.com.br/books?id=YxwW4gMGv5AC&printsec=frontcover&dq=gujarati+econometria+b%C3%A1sica&hl=pt-BR&sa=X&ved=0ahUKEwjF2LrPxLnVAhWJE5AKHfU7A2oQ6AEIJzAA#v=onepage&q=gujarati%20econometria%20b%C3%A1sica&f=false|título=Econometria Basica|ultimo=Gujarati|primeiro=Dadomar N.|data=2006|editora=Elsevier|lingua=pt|isbn=9788535216646}}</ref>

O modelo é geralmente referido como o modelo ARMA(<math>p,q</math>), em que <math>p</math> é a ordem da parte auto-regressiva e <math>q</math> é 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(<math>p</math>) se refere ao modelo auto-regressivo de ordem <math>p</math>. O modelo AR(<math>p</math>) é escrito como:<blockquote><math> X_t=c+\sum_{i=1}^p\varphi_i X_{t-i}+\varepsilon_t,</math></blockquote>em que <math>\varphi_1...\varphi_p</math> são os [[Parâmetro (estatística)|parâmetros]], <math>c</math> é uma constante e a variável aleatória <math>\varepsilon_t</math> é [[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 <math>|\varphi_1|\geq1</math> não são estacionários.<ref name=":0">{{Citar livro|url=https://books.google.com.br/books?id=OKUGARAXKMwC&printsec=frontcover&dq=Ruey+Tsay+Analysis+of+Financial+Time+Series&hl=pt-BR&sa=X&redir_esc=y#v=onepage&q=Ruey%20Tsay%20Analysis%20of%20Financial%20Time%20Series&f=false|título=Analysis of Financial Time Series|ultimo=Tsay|primeiro=Ruey S.|data=2010-10-26|editora=John Wiley & Sons|lingua=en|isbn=9781118017098}}</ref>

==Modelo de média móvel==
A notação MA(<math>q</math>) se refere ao modelo de média móvel de ordem <math>q</math>:<blockquote><math> X_t=\mu+\varepsilon_t+\sum_{i=1}^q\theta_i \varepsilon_{t-i},</math></blockquote>em que <math>\theta_1,...,\theta_q</math> são os parâmetros do modelo, <math>\mu</math> é o [[valor esperado]] de <math>X_t</math> (frequentemente assumido como igual a zero) e <math>\epsilon_t,\epsilon_{t-1},...</math> são novamente os termos de erro de ruído branco.<ref name=":0" />

==Modelo ARMA==
A notação ARMA(<math>p,q</math>) se refere ao modelo com termos auto-regressivos <math>p</math> e termos de média móvel <math>q</math>. Este modelo contém os modelos AR(<math>p</math>) e MA(<math>q</math>):<blockquote><math> X_t=c+\varepsilon_t+\sum_{i=1}^p\varphi_i X_{t-i}+\sum_{i=1}^q\theta_i\varepsilon_{t-i}.</math></blockquote>O modelo ARMA geral foi descrito por Peter Whittle em sua tese de 1951, em que usou análise matemática ([[Série de Laurent|séries de Laurent]] e [[análise de Fourier]]) e [[inferência estatística]].<ref>{{Citar livro|url=https://books.google.com.br/books?id=INhGjn_4f-QC&printsec=frontcover&dq=Hannan+Multiple+time+series&hl=pt-BR&sa=X&redir_esc=y#v=onepage&q=Hannan%20Multiple%20time%20series&f=false|título=Multiple Time Series|ultimo=Hannan|primeiro=Edward James|data=2009-09-25|editora=John Wiley & Sons|lingua=en|isbn=9780470317136}}</ref><ref>{{Citar livro|url=https://books.google.com.br/books?id=Ghc4ngAACAAJ&dq=Whittle+Prediction+and+Regulation+by+Linear+Least-Square+Methods&hl=pt-BR&sa=X&ved=0ahUKEwjA_tHCvbnVAhVChpAKHZBLDz8Q6AEIJzAA|título=Prediction and Regulation by Linear Least-square Methods|ultimo=Whittle|primeiro=Peter|data=1983|editora=University of Minnesota Press|lingua=en|isbn=9781452938554}}</ref> 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).<ref>{{Citar livro|url=https://books.google.com.br/books?id=n_5JCmPCQ7QC&printsec=frontcover&dq=Hannan+Deistler+Statistical+theory+of+linear+systems&hl=pt-BR&sa=X&ved=0ahUKEwjZnMjVvbnVAhVOlpAKHbJhCuQQ6AEIJzAA#v=onepage&q=Hannan%20Deistler%20Statistical%20theory%20of%20linear%20systems&f=false|título=The Statistical Theory of Linear Systems|ultimo=Hannan|primeiro=E. J.|ultimo2=Deistler|primeiro2=Manfred|data=2012-05-31|editora=SIAM|lingua=en|isbn=9781611972184}}</ref>

==Termos de erro==
Os termos de erro <math>\epsilon_t</math> 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, <math>\epsilon_t\sim N(0,\sigma^2)</math>, em que <math>\sigma^2</math> é 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.<ref>{{Citar livro|url=https://books.google.com.br/books?id=D6-CAAAAIAAJ&q=McCleary+Hay+Applied+Time+Series+Analysis+for+the+Social+Sciences&dq=McCleary+Hay+Applied+Time+Series+Analysis+for+the+Social+Sciences&hl=pt-BR&sa=X&ved=0ahUKEwiKv__-xLnVAhXKDpAKHf6ZBWkQ6AEIJzAA|título=Applied time series analysis for the social sciences|ultimo=McCleary|primeiro=Richard|ultimo2=Hay|primeiro2=Richard|data=1980-06-01|editora=Sage Publications|lingua=en|isbn=9780803912052}}</ref>

==Especificação em termos de operador de defasagem==
Em alguns textos, os modelos serão especificados em termos do [[operador de defasagem]] <math>L</math>. Nestes termos, o modelo AR(<math>p</math>) é dado por:<blockquote><math> \varepsilon_t=\left(1-\sum_{i=1}^p\varphi_iL^i\right)X_t=\varphi(L)X_t,</math></blockquote>em que <math>\varphi</math> representa o polinômio:<blockquote><math> \varphi(L)=1-\sum_{i=1}^p\varphi_iL^i.</math></blockquote>O modelo MA(<math>q</math>) é dado por:<blockquote><math> X_t=\left(1+\sum_{i=1}^q\theta_iL^i\right)\varepsilon_t=\theta(L)\varepsilon_t,</math></blockquote>em que <math>\theta</math> representa o polinômio:<blockquote><math> \theta(L)=1+\sum_{i=1}^q\theta_iL^i.</math></blockquote>Finalmente, o modelo combinado ARMA(<math>p,q</math>) é dado por<blockquote><math> \left(1-\sum_{i=1}^p\varphi_iL^i\right)X_t=\left(1+\sum_{i=1}^q\theta_iL^i\right)\varepsilon_t,</math></blockquote>ou, de forma mais concisa,<blockquote><math> \varphi(L)X_t=\theta(L)\varepsilon_t, </math></blockquote>ou<blockquote><math> \frac{\varphi(L)}{\theta(L)}X_t=\varepsilon_t.</math><ref>{{Citar livro|url=https://books.google.com.br/books?id=B8_1UBmqVUoC&dq=Hamilton+Time+Series+Analysis&hl=pt-BR&sa=X&redir_esc=y|título=Time Series Analysis|ultimo=Hamilton|primeiro=James Douglas|data=1994|editora=Princeton University Press|lingua=en|isbn=9780691042893}}</ref></blockquote>

===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.<ref>{{Citar livro|url=https://books.google.com.br/books?id=rNt5CgAAQBAJ&printsec=frontcover&dq=Box+Jenkins+Reinsel+Time+Series+Analysis&hl=pt-BR&sa=X&ved=0ahUKEwi35c-UvrnVAhWKhpAKHff1BPwQ6AEIJzAA#v=onepage&q=Box%20Jenkins%20Reinsel%20Time%20Series%20Analysis&f=false|título=Time Series Analysis: Forecasting and Control|ultimo=Box|primeiro=George E. P.|ultimo2=Jenkins|primeiro2=Gwilym M.|ultimo3=Reinsel|primeiro3=Gregory C.|ultimo4=Ljung|primeiro4=Greta M.|data=2015-05-29|editora=John Wiley & Sons|lingua=en|isbn=9781118674925}}</ref> 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:<blockquote><math> \left(1-\sum_{i=1}^p\phi_i L^i\right)X_t=\left(1+\sum_{i=1}^q\theta_iL^i\right)\varepsilon_t.</math></blockquote>Além disso, se configurarmos <math>\phi_0=\theta_0=1</math>, então, obtemos uma formulação ainda mais sofisticada:<blockquote><math> \sum_{i=0}^p\phi_iL^iX_t=\sum_{i=0}^q\theta_iL^i\varepsilon_t.</math></blockquote>

==Modelos de ajuste==

===Escolha de <math>p</math> e <math>q</math>===
A descoberta dos valores apropriados de <math>p</math> e <math>q</math> no modelo ARMA(<math>p,q</math>) pode ser facilitada ao mapear as funções de [[autocorrelação]] parcial para um valor estimado de <math>p</math> e, da mesma forma, usar as funções de autocorrelação para um valor estimado de <math>q</math>. 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 <math>p</math> e <math>q</math>.

Peter J. Brockwell e Richard A. Davis recomendaram o uso do critério de informação de Akaike para encontrar <math>p</math> e <math>q</math>.<ref>{{Citar livro|url=https://books.google.com.br/books?id=DJ_lBwAAQBAJ&printsec=frontcover&dq=Brockwell+Davis+Time+Series&hl=pt-BR&sa=X&ved=0ahUKEwiZk5bNvrnVAhXHE5AKHTm4AgYQ6AEIJzAA#v=onepage&q=Brockwell%20Davis%20Time%20Series&f=false|título=Time Series: Theory and Methods|ultimo=Brockwell|primeiro=Peter J.|ultimo2=Davis|primeiro2=Richard A.|data=2013-11-11|editora=Springer Science & Business Media|lingua=en|isbn=9781489900043}}</ref>

===Estimação dos coeficientes===
Em geral, modelos ARMA não podem ser, depois de escolher <math>p</math> e <math>q</math>, ajustados por regressão de [[Método dos mínimos quadrados|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 <math>p</math> e <math>q</math> 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.<ref>{{Citar livro|url=https://books.google.com.br/books?id=cNe3xrFg3PcC&printsec=frontcover&dq=Time+Series+Techniques+for+Economists&hl=pt-BR&sa=X&redir_esc=y#v=onepage&q=Time%20Series%20Techniques%20for%20Economists&f=false|título=Time Series Techniques for Economists|ultimo=Mills|primeiro=Terence C.|data=1991-06-13|editora=Cambridge University Press|lingua=en|isbn=9780521405744}}</ref>

==Generalizações==
A dependência de <math>X_t</math> de valores passados e os termos de erro <math>\varepsilon_t</math> 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(<math>p,q,b</math>) se refere ao modelo com <math>p</math> termos auto-regressivos, <math>q</math> termos de média móvel e <math>b</math> termos de ''input'' exógeno. Este modelo contém os modelos AR(<math>p</math>) e MA(<math>q</math>) e uma combinação linear dos últimos <math>b</math> termos de uma série temporal conhecida e externa <math>d_t</math>.<ref>{{Citar periódico|ultimo=Yang|primeiro=Hong-Tzer|ultimo2=Huang|primeiro2=Chao-Ming|ultimo3=Huang|primeiro3=Ching-Lien|ano=1996|mes=02|titulo=Identification of ARMAX model for short term load forecasting: an evolutionary programming approach|jornal=IEEE Transactions on Power Systems|volume=11|numero=1|paginas=403–408|issn=0885-8950|doi=10.1109/59.486125|url=http://ieeexplore.ieee.org/document/486125/|acessadoem=}}</ref> É dado por:<blockquote><math> X_t=\varepsilon_t+\sum_{i=1}^p\varphi_i X_{t-i}+\sum_{i=1}^q\theta_i \varepsilon_{t-i}+\sum_{i=1}^b\eta_id_{t-i}.</math></blockquote>em que <math>\eta_1,...,\eta_b</math> são os parâmetros do ''input'' exógeno <math>d_t</math>.

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:<blockquote><math> X_t-m_t=\varepsilon_t+\sum_{i=1}^p \varphi_i(X_{t-i}-m_{t-i})+\sum_{i=1}^q \theta_i\varepsilon_{t-i},</math></blockquote>em que <math>m_t</math> incorpora todas as variáveis exógenas (ou independentes):<blockquote><math>m_t=c+\sum_{i=0}^b\eta_id_{t-i}.</math></blockquote>

==Ver também==
*[[Modelos ARCH]]
*[[ARIMA]]
*[[Série temporal]]

== Referências ==
{{Reflist}}

{{Econometria}}
{{Processos estocásticos}}

[[Categoria:Econometria]]
[[Categoria:Processos estocásticos]]