Untitled - UFRJ

Untitled - UFRJ Untitled - UFRJ

11.07.2015 Views

Processo ARMA: Um Estudo BayesianoAdriana Strieder PhilippsenInstituto de Ciências Matemáticas e de Computação - ICMC-USPMarcos H. CasconeUniversidade Federal de São Carlos - UFSCarMarinho G. AndradeInstituto de Ciências Matemáticas e de Computação - ICMC-USPRicardo Sandes EhlersInstituto de Ciências Matemáticas e de Computação - ICMC-USPNa análise de uma série temporal, primeiramente deseja-se modelar o fenômeno de interesse paraque possamos descrever o seu comportamento, fazendo estimativas, avaliando quais fatores influenciamno comportamento da série, buscando definir relações de causa e efeito entre duas ou mais séries, e porúltimo fazer previsões de valores futuros da série com base em valores passados fornecendo informaçõese subsídios para uma consequente tomada de decisão.Na análise do comportamento de uma série livre de tendência e de sazonalidade podem ser utilizadosmodelos auto-regressivos (AR) ou que incorporem médias móveis (MA), sendo que o modelo ARMAcombina as características dos dois processos anteriores. Este processo relaciona os valores futuros comas observações passadas, assim como também com os erros passados apurados entre os valores reais e osprevistos.Neste trabalho, utilizando dados simulados e um conjunto de dados reais, é conduzida uma análiseclássica e bayesiana para o processo ARMA(p,q). Para o estudo clássico, as estimativas dos parâmetrosdo modelo adotado são obtidas fazendo uso do procedimento de estimadores de máxima verossimilhança,que equivale ao método de mínimos quadrados ponderados reiterados, conhecido por IRLS ou Scoring deFisher. Na estimação bayesiana, será utilizada uma priori não informativa que juntamente com a funçãode verossimilhança do modelo, fornecerá uma distribuição a posteriori, que é empregada no cálculo damédia da distribuição a posteriori como estimativa de um parâmetro. Isto tem justificativa teórica, nocontexto na Teoria de Decisão, quando considera-se a função de perda quadrática pois, neste caso, é ovalor que minimiza o risco de Bayes. Os modelos ajustados serão comparados usando o erro quadráticomédio (MSE) e o erro percentil absoluto da média (MAPE). Como resultado final das análises, seráapresentado um resumo descritivo das distribuições condicionais a posteriori dos parâmetros de interesse,bem como os valores das médias e seus gráficos. Também serão apresentadas as previsões k passos àfrente. As amostras da distribuição a posteriori serão geradas fazendo uso das técnicas de simulaçãoMCMC, em particular, o algoritmo de Gibbs sampling e para verificar a convergência será utilizado ocritério proposto por Geweke (1992).117

Estudo das Matrizes de Covariância de Modelos Bayesianoscom Interação Espaço-TemporalLetícia Cavalari PinheiroDepartamento de Estatística – Instituto de Ciências Exatas –Universidade Federal de Minas Gerais (UFMG)Renato Martins AssunçãoDepartamento de Estatística – Instituto de Ciências Exatas –Universidade Federal de Minas Gerais (UFMG)Nesse trabalho estudamos as matrizes de covariância de modelos bayesianos com efeitos de interaçãoespaço-temporal. Para isso, apresentamos os possíveis tipos de efeitos aleatórios espaciais e temporaisem modelos bayesianos e atribuimos a eles distribuições a priori comumente utilizadas. Construimospossíveis efeitos aleatórios espaço-temporais a partir da interação entre um efeito temporal e um espacial.Calculamos as matrizes de covariância a priori para os modelos com interação espaço-tempo eas escrevemos na forma de produto de Kronecker entre as matrizes de covariância a priori dos efeitostemporal e espacial. Conseguimos visualizar mais claramente o efeito de cada tipo de interação possível,relacionando as matrizes de covariância a priori com as estruturas de dependência espacial e/ou temporalenvolvidas nos modelos estudados. Como exemplo, apresentamos o estudo das matrizes de covariânciaa priori de dois modelos específicos existentes na literatura e fazemos suas interpretações.118

Processo ARMA: Um Estudo BayesianoAdriana Strieder PhilippsenInstituto de Ciências Matemáticas e de Computação - ICMC-USPMarcos H. CasconeUniversidade Federal de São Carlos - UFSCarMarinho G. AndradeInstituto de Ciências Matemáticas e de Computação - ICMC-USPRicardo Sandes EhlersInstituto de Ciências Matemáticas e de Computação - ICMC-USPNa análise de uma série temporal, primeiramente deseja-se modelar o fenômeno de interesse paraque possamos descrever o seu comportamento, fazendo estimativas, avaliando quais fatores influenciamno comportamento da série, buscando definir relações de causa e efeito entre duas ou mais séries, e porúltimo fazer previsões de valores futuros da série com base em valores passados fornecendo informaçõese subsídios para uma consequente tomada de decisão.Na análise do comportamento de uma série livre de tendência e de sazonalidade podem ser utilizadosmodelos auto-regressivos (AR) ou que incorporem médias móveis (MA), sendo que o modelo ARMAcombina as características dos dois processos anteriores. Este processo relaciona os valores futuros comas observações passadas, assim como também com os erros passados apurados entre os valores reais e osprevistos.Neste trabalho, utilizando dados simulados e um conjunto de dados reais, é conduzida uma análiseclássica e bayesiana para o processo ARMA(p,q). Para o estudo clássico, as estimativas dos parâmetrosdo modelo adotado são obtidas fazendo uso do procedimento de estimadores de máxima verossimilhança,que equivale ao método de mínimos quadrados ponderados reiterados, conhecido por IRLS ou Scoring deFisher. Na estimação bayesiana, será utilizada uma priori não informativa que juntamente com a funçãode verossimilhança do modelo, fornecerá uma distribuição a posteriori, que é empregada no cálculo damédia da distribuição a posteriori como estimativa de um parâmetro. Isto tem justificativa teórica, nocontexto na Teoria de Decisão, quando considera-se a função de perda quadrática pois, neste caso, é ovalor que minimiza o risco de Bayes. Os modelos ajustados serão comparados usando o erro quadráticomédio (MSE) e o erro percentil absoluto da média (MAPE). Como resultado final das análises, seráapresentado um resumo descritivo das distribuições condicionais a posteriori dos parâmetros de interesse,bem como os valores das médias e seus gráficos. Também serão apresentadas as previsões k passos àfrente. As amostras da distribuição a posteriori serão geradas fazendo uso das técnicas de simulaçãoMCMC, em particular, o algoritmo de Gibbs sampling e para verificar a convergência será utilizado ocritério proposto por Geweke (1992).117

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!