Previsão da volatilidade no mercado interbancário de câmbio

ARTIGOS

Previsão da volatilidade no mercado interbancário de câmbio

Clayton Peixoto Goulart^I; Hudson Fernandes Amaral^II; Luiz Alberto Bertucci^III; Aureliano Angel Bressan^IV

^IUFMG-CEPEAD

^IIUFMG-CEPEAD

^IIIUFMG-CEPEAD

^IVUFMG-CEPEAD

RESUMO

O artigo apresenta um estudo comparativo da capacidade preditiva dos modelos EWMA, GARCH (1,1), EGARCH (1,1) e TARCH (1,1), quando utilizados para prever a volatilidade das taxas de câmbio praticadas no mercado interbancário brasileiro. A amostra é composta pelas cotações diárias de fechamento da taxa de câmbio real/dólar estadunidense observadas no período de 20 de agosto de 2001 a 30 de setembro de 2003. Os resultados demonstraram que o modelo TARCH (1,1) apresentou o melhor desempenho preditivo para o período, acompanhado de perto pelo modelo EGARCH (1,1), seguindo-se o modelo GARCH (1,1) e, por último, o modelo EWMA. Constatou-se também que todos os modelos revelaram uma propensão a superestimar a volatilidade futura, e que a Clearing de Câmbio da BM&F atua de forma excessivamente conservadora e subjetiva na definição dos índices de variação da taxa de câmbio contratualmente garantidos e, conseqüentemente, na exigência de garantias.

Palavras-chave: Volatilidade, câmbio, risco, modelos econométricos.

ABSTRACT

This paper presents a comparative study of the predictive capacity of the models EWMA, GARCH (1,1) EGARCH (1,1) and TARCH (1,1) when applied to forecast the volatility of the exchange rates in the Brazilian inter-bank market. The sample consists of the daily closing quotations of the exchange rate real/US dollar obtained in the period from August 20, 2001 to September 30, 2003. The results showed that the TARCH (1,1) model, provided the most accurate forecast performance for this period, followed closely by the EGARCH (1,1) model, then the GARCH (1,1) model, and finally the EWMA model. There was also evidence that all the models revealed a tendency to overestimate the future volatility. It was confirmed as well that the Brazilian Exchange Clearinghouse operates in an extremely traditional and subjective way concerning the definition of the exchange rate variation indexes guaranteed under contract and, consequently, in collateral requirements.

Key words: Volatility, exchange, risk, econometrics models.

INTRODUÇÃO

Desde o trabalho inicial de Markowitz (1952), a volatilidade é uma variável presente em grande parte da moderna teoria de finanças e dos processos de tomada de decisão financeira. Nesse sentido, Bollerslev, Chou e Kroner (1992, p. 46) afirmam que a "volatilidade é uma variável-chave que permeia a maioria dos instrumentos financeiros e que exerce um papel central em diversas áreas de finanças". Entretanto, passados mais de 50 anos de pesquisas e discussões acadêmicas, observa-se que ainda não há um consenso sobre o modelo mais adequado para o cálculo da volatilidade, gerando ainda calorosas discussões e abrindo um amplo campo de pesquisas.

Durante muitos anos a proposição dominante nas publicações científicas financeiras e nas aplicações de mercado era baseada na pressuposição de que se poderia considerar a volatilidade dos retornos dos ativos financeiros como um fator constante no tempo. Embora Mandelbrot (1963) e Fama (1965) já tivessem alertado para as mudanças na volatilidade dos ativos ao longo do tempo, destaca-se a ênfase dada por Merton (1980, p 354) no que se refere à necessidade do "desenvolvimento de modelos mais precisos de estimação de variâncias que levassem em conta os erros de estimativas passadas".

Como preconizado por Merton (1980), Engle (1982) propôs uma nova classe de processos estocásticos denominada Autoregressive Conditional Heteroskedasticity [heteroscedasticidade condicional auto-regressiva] (ARCH), na qual a variância de uma série temporal modifica-se, com o passar do tempo, de forma condicional aos erros de previsão observados no passado. Posteriormente, Bollerslev (1986) generalizou o tratamento da variância condicional por meio do modelo intitulado Generalized Autoregressive Conditional Heteroskedasticity [heteroscedasticidade condicional auto-regressiva generalizada] (GARCH), no qual a variância condicional em um determinado tempo t depende não somente dos erros de previsão observados no passado, como também das variâncias condicionadas observadas no passado (Alexander 2002; Brooks, 2002).

Durante os anos que se seguiram, diversas variantes dos modelos originais ARCH e GARCH foram desenvolvidas no sentido de melhor representar o comportamento dos diversos ativos nos mercados financeiros, por vezes apresentando maior complexidade matemática e, conseqüentemente, de implementação; por vezes adotando simplificações que, por pouco, não descaracterizam os princípios básicos dos modelos originais. Dentre os modelos simplificados, destacase o adotado pelo RiskMetrics^TM, denominado Exponentially Weighted Moving Averages [modelo de médias móveis ponderadas exponencialmente] (EWMA), que pode ser interpretado como um caso especial do modelo Integrated GARCH (I-GARCH), e que obteve ampla aceitação pelo mercado financeiro (Alexander, 1998).

No campo acadêmico, foram e são efetuados diversos estudos empíricos testando a capacidade preditiva dos modelos da família ARCH, principalmente nos mercados de câmbio e de ações, originando um grande número de artigos, teses e dissertações sobre o assunto. No âmbito comercial, embora no início dos anos 1990 fosse possível considerar que os modelos ARCH - como ferramentas de decisão - estivessem em um estágio ainda experimental, atualmente muitas instituições já detêm a competência necessária para utilizálos (Gouriéroux, 1997). No Brasil, infelizmente, tanto no âmbito acadêmico como no comercial, ainda não se observa o mesmo volume de publicações, pesquisas e aplicativos que se constata nos mercados financeiros dos países desenvolvidos, principalmente no que se refere ao mercado de câmbio. Vários motivos podem ser levantados para esse fato, sendo o mais relevante as diversas mudanças de regime cambial ocorridas nos últimos 30 anos, inseridas ou não no contexto de um plano de estabilização financeira.

Entretanto, recentemente diversos fatores vieram a contribuir para a criação de um ambiente propício à consecução de novos estudos sobre o tema no mercado cambial brasileiro. Podem ser destacados: a) adoção da política de câmbio livre, a partir de janeiro de 1999, com a redução da participação do Estado - em especial do Banco Central do Brasil (Bacen) - na definição do patamar da taxa de câmbio corrente; b) a implementação do novo Sistema de Pagamentos Brasileiro (SPB), a partir de 22 de abril de 2002, com a formalização e institucionalização das Câmaras de Registro, Compensação e Liquidação - Clearings -, em particular da Clearing de Câmbio da Bolsa de Mercadorias e Futuros (BM&F); e, c) a necessidade de incremento do comércio exterior, visando a obtenção de superávits comerciais crescentes.

Dessa forma, este artigo apresenta um estudo comparativo da capacidade preditiva de modelos da família ARCH (GARCH, EGARCH e TARCH) em relação ao modelo EWMA - adotado pelo RiskMetrics^TM e pela Clearing de Câmbio da BM&F -, no sentido de contribuir não apenas para preencher parte da lacuna observada nos estudos relativos aos modelos de predição de volatilidade no mercado de câmbio interbancário brasileiro, como também para apresentar um estudo científico da capacidade preditiva dos modelos aplicados a processos de gerenciamento de riscos e controle de garantias no referido mercado. Ademais, o estudo pode servir como referência para discussões futuras, por parte da autoridade monetária, da Clearing de Câmbio da BM&F e dos operadores de mercado, sobre a pertinência da utilização de modelos econométricos mais elaborados no procedimento de definição do índice de variação da taxa de câmbio contratualmente garantido pela câmara.

REFERENCIAL TEÓRICO

Ao contrário dos preços, a volatilidade não é diretamente observável no mercado, podendo apenas ser estimada no contexto de um modelo. Dessa forma, é plenamente possível a obtenção de diferentes valores de volatilidade utilizando-se um mesmo conjunto de dados, bastando utilizar modelos distintos, que são, em sua grande maioria, teoricamente bem fundamentados (Alexander, 2002). Conforme mostram Poon e Granger (2003), em pesquisa realizada em 93 artigos e working papers publicados no período de 1976 a 2002, a diversidade de resultados ainda não permite que se firme um consenso sobre a capacidade preditiva dos diversos modelos de predição de volatilidade, configurando-se um amplo campo para pesquisas e debates. Nesse sentido, e no caso específico desta pesquisa, optou-se por comparar a capacidade preditiva dos modelos EWMA, GARCH (1,1), EGARCH (1,1) e TARCH (1,1), sobre os quais será efetuada uma breve discussão a seguir.

Segundo Brooks (2002, p. 442), "o modelo de média móvel com ponderação exponencial (EWMA) é, essencialmente, uma simples extensão do modelo de volatilidade histórica ponderada, no qual as observações mais recentes têm um impacto maior na previsão da volatilidade que os dados mais antigos". De outra sorte, Alexander (1998) postula que o modelo EWMA pode ser entendido como uma simplificação do modelo IGARCH, no qual considera-se a volatilidade de longo prazo igual a zero. Formalmente, o modelo de médias móveis com ponderação exponencial, para n períodos de uma série temporal constituída dos retornos de um ativo, pode ser definido como (Alexander, 2002):

onde:

é o retorno médio;

r _t-i é o retorno da série com t-i defasagens;

n é o número de observações utilizadas;

λ é o fator de decaimento exponencial, tal que 0 < λ < 1; e

σ _t²é a variância prevista para o tempo t.

O modelo EWMA recebeu especial impulso junto ao mercado financeiro a partir de 1994, com o lançamento do software de administração de riscos RiskMetrics^TM pelo J. P. Morgan. A grande aceitação desse software e a divulgação livre do respectivo documento técnico (RiskMetrics: Technical Document, Morgan Guaranty Trust Company & Reuters Ltd., 1996) rapidamente tornaram a metodologia um padrão de mercado.

Entretanto, é importante observar que o referido documento técnico tece algumas considerações e estabelece premissas para a utilização do modelo no âmbito do software RiskMetrics^TM, dentre as quais se destacam: a) os retornos das séries financeiras são calculados de forma contínua, tomando-se o logaritmo da relação entre os preços observados nos períodos t e t-1, considerando o comportamento log-normal dos preços; b) o retorno médio é definido como sendo nulo, para qualquer série financeira e em qualquer intervalo de tempo utilizado; c) o valor do fator de decaimento (λ), quando se utilizam retornos diários, obedece a relação predeterminada, definindo-se um total de 149 observações (n) e λ = 0,94, para qualquer mercado ou ativo analisado; e, d) a volatilidade para múltiplos dias não pode ser estimada diretamente por meio do modelo EWMA (tipicamente um modelo do tipo "um passo à frente"), assumindo-se a regra da "raiz quadrada do tempo" para se calcular a volatilidade em períodos maiores.

Em uma outra perspectiva, encontram-se os modelos da classe ARCH. A proposição original, elaborada por Engle (1982), mereceu extensos debates e diversos aperfeiçoamentos ao longo dos anos. A primeira, e mais significativa, foi introduzida por Bollerslev (1986), ao propor que a volatilidade condicionada fosse função não apenas dos quadrados dos erros passados (ε²_t-i), como também dos seus próprios valores passados (σ²_t-j), passando os modelos assim construídos a ser denominados Generalized ARCH (GARCH). Em termos matemáticos, um modelo GARCH(p,q) pode ser expresso como:

onde:

σ²_t é a é variância condicionada em t

ω é a constante;

ε_t-i é o erro observado no tempo t-i, em que i denota a defasagem;

α_i é a constante multiplicadora do termo de defasagem i;

β_j é a constante multiplicadora do termo de defasagem j;

σ² _t-j é a variância condicionada observada em t-j; e,

ν _t é um ruído branco [N~(0,1)].

Saliente-se que, para se garantir que a variância condicionada não seja negativa, e também a estacionariedade do processo, as seguintes restrições devem ser impostas à Equação 1: ω > 0; α_i> 0, para i = 1, ... , p; β_j> 0, para j = 1, ... , q ; e

Outro ponto a se salientar é que, apesar de o modelo GARCH(p,q) captar corretamente diversas características observadas nas séries históricas de finanças, tais como a leptocurtose e o agrupamento de volatilidades, não capta o efeito de alavancagem, uma vez que a variância condicional é função apenas da magnitude das inovações e não dos seus sinais (Brooks, 2002). Assim, surgiram algumas extensões com a finalidade de modelar o problema da assimetria, tais como os modelos EGARCH e TARCH.

O modelo EGARCH (Exponential GARCH) foi proposto por Nelson (1991), apresentando duas vantagens em relação ao modelo GARCH: a incorporação dos efeitos assimétricos de mercado aos modelos auto-regressivos de volatilidade condicionada, e a não imposição artificial de restrições aos parâmetros da equação, dada a sua formulação em termos logarítmicos. Por outro lado, o referido modelo necessita de adaptações para sua implementação computacional e, por vezes, apresenta um alto grau de persistência após choques significativos de volatilidade. Segundo Brooks (2002), há vários modos de se expressar a equação da variância condicional, dentre os quais:

onde:

σ²_t é a variância condicionada em t;

ω é a constante;

ε_t-1 é o erro observado no tempo t-1;

α, β, γ são coeficientes da equação;

σ_t-1 é o desvio padrão observado em t-1; e

ν _t é um ruído branco [N~(0,1)].

Um modelo mais simples e intuitivo para a captação do comportamento assimétrico das séries financeiras foi apresentado por Glosten, Jagannathan e Runkle (1993), sendo denominado modelo GJR ou TARCH (Threshold ARCH). Nesse novo modelo, é acrescentado um termo à equação-padrão GARCH, cuja função é o reconhecimento dos efeitos assimétricos à volatilidade condicional, ou seja:

onde:

σ² t é a variância condicionada em t;

ω é a constante;

ε_t-i é o erro observado no tempo t-i, onde i denota a defasagem;

α_i é a constante multiplicadora do termo de defasagem i;

β_j é a constante multiplicadora do termo de defasagem j;

σ²_t-j é a variância condicionada observada em t-j;

γ é o coeficiente do termo de assimetria;

d_t-1 é a variável dummy que assume o valor igual a 1 se ε_t-1 < 0 (más notícias), e valor igual a 0 se ε_t-1 > 0 (boas notícias); e

ν_t é um ruído branco [N~(0,1)].

Destacando que, como no modelo GARCH, as condições para não negatividade e estacionariedade são: ω> 0; α_i> 0, para i = 1, ... ,q; β_j> 0, para j = 1, ... , p; e Σα_i+γ> 0 (Brooks, 2002).

Evidentemente, a importância dos modelos acima citados está diretamente relacionada ao uso da volatilidade prevista nos processos de tomada de decisão financeira, precificação de ativos e nas metodologias de gerenciamento de risco. No caso específico da Clearing de Câmbio da BM&F, a utilização da volatilidade prevista está descrita no Anexo I do seu Manual de Gerenciamento de Risco, qual seja: os índices de variação da taxa de câmbio contratualmente garantidos - base para a exigência de garantias a serem depositadas pelos participantes - são determinados pela conjunção de duas análises. A primeira, a estimação da volatilidade, por meio de modelos econométricos (avaliação objetiva); e a segunda, a avaliação da conjuntura macroeconômica pelo Comitê da Câmara (avaliação subjetiva). A parcela estatística do índice (avaliação objetiva) segue a abordagem paramétrica sugerida por Jorion (1998) para a determinação do "valor em risco" com um nível de confiança de 99%. Assim, para operações negociadas pelo prazo de um dia, o desvio padrão calculado pelo modelo EWMA é multiplicado por 2,33; já para as operações negociadas pelo prazo de dois dias, o desvio padrão é multiplicado por 2,33.

. Ressaltando-se que a decisão final do Comitê da Clearing baseia-se na análise econômica dos fatos (BM&F, 2002, p. 83).

Observa-se, portanto, que duas questões estão envolvidas nesta pesquisa: a comparação da capacidade preditiva dos modelos anteriormente citados e os possíveis reflexos que uma substituição do modelo EWMA poderia acarretar na definição do índice de taxa de câmbio contratualmente garantido pela Clearing. A próxima seção discorrerá sobre a metodologia empregada para se avaliar essas questões.

PROCEDIMENTOS METODOLÓGICOS

Neste estudo foram utilizadas duas bases de dados. A primeira refere-se às cotações diárias de fechamento da taxa de câmbio real/dólar norte-americano (R$/US$) observadas no período de 20 de agosto de 2001 a 30 de setembro de 2003, num total de 536 observações. Os dados foram fornecidos pela Clearing de Câmbio da BM&F e são os mesmos por ela utilizados nas rotinas internas de cálculo de volatilidade. É importante ressaltar que a utilização de dados no período de 20 de agosto de 2001 a 19 de abril de 2002 foi devida a uma necessidade prática. Conforme observado anteriormente, a Clearing adotou o modelo EWMA para o cálculo, qual seja, a janela móvel com 149 observações e fator de decaimento (1) de 0,94. Assim, para viabilizar o início das atividades da Clearing em 22 de abril de 2002, a direção da Câmara decidiu utilizar as taxas de câmbio do período de 20/8/2001 a 19/3/2002 (que perfazem um total de 149 observações) como valores iniciais do modelo. Já o intervalo compreendido entre 20 de março de 2002 e 19 de abril de 2002 foi utilizado para o acompanhamento do comportamento da volatilidade no período imediatamente anterior ao início das atividades da Câmara. A escolha da data final para a amostra - 30 de setembro de 2003 - foi arbitrada pelos pesquisadores, não tendo, portanto, nenhum significado especial.

A segunda base de dados refere-se aos índices diários de variação da taxa de câmbio (percentual garantido), estipulados no período de 22 de abril de 2002 a 30 de setembro de 2003, para operações a serem liquidadas em D+1 e D+2. Também nesse caso, as informações foram disponibilizadas pela Câmara. Para efeito deste estudo, assume-se que as diferenças entre os valores calculados e os valores estipulados são devidas à parcela subjetiva referente à avaliação macroeconômica.

Com base na premissa de log-normalidade dos preços dos ativos, a série de taxas de câmbio (R$/US$) foi transformada por intermédio da equação r_t= lnP_t - lnP_t-1, resultando em uma série de retornos de taxas de câmbio (um total de 535 observações). Dessa forma, a série de retornos substituiu a série original de cotações de taxas de câmbio como dados de entrada para os testes estatísticos e para os modelos de predição de volatilidade.

O primeiro passo da investigação empírica foi efetuado por meio de uma análise estatística descritiva da série de retornos da taxa de câmbio. Com o intuito de analisar possíveis influências de fatores externos no comportamento dos retornos da taxa de câmbio, os testes foram realizados na série completa (de 21 de agosto de 2001 a 30 de setembro de 2003, perfazendo 535 observações) e em três subperíodos, a saber: a) de 21 de agosto de 2001 a 19 de abril de 2002, perfazendo 170 observações, cuja motivação foi verificar o comportamento da série no período de antecede o início de operação da Câmara e do novo SPB; b) de 22 de abril de 2002 a 30 de dezembro de 2002, perfazendo 177 observações, buscando observar o comportamento da série no período pré-eleitoral brasileiro; e c) de 2 de janeiro de 2003 a 30 de setembro 2003, perfazendo 188 observações, no intuito de observar o comportamento da série em um período de estabilização de um novo governo, marcado por uma política fiscal e monetária restritiva.

O segundo passo visou constatar a presença de heteroscedasticidade condicional auto-regressiva na série de retornos, pois tal característica é condição fundamental para justificar a aplicação de modelos da família ARCH. Sobre os mesmos períodos acima especificados, foi aplicado o teste ARCH-LM (ARCH - Lagrange Multipliers), conforme sugerido por Engle (1982). Os testes foram realizados com o auxílio do software Eviews^TM 4.0.¹

O modelo EWMA foi implementado na sua forma não-recursiva, conforme apresentada na Equação 1, com n = 149 e λ = 0,94. Os cálculos relativos a esse modelo foram efetuados com o uso da planilha eletrônica Microsoft^® Excel 2000, por meio de rotina desenvolvida pelos autores deste artigo.

Antes de descrever os modelos da família ARCH, é importante ressaltar que sua utilização exigiu a definição em conjunto de uma equação para a média condicional, uma vez que os erros apurados na equação da média são os insumos necessários à operacionalização dos cálculos da volatilidade condicionada. Alexander (1998) postula que, quanto mais parcimoniosa for a definição da equação da média condicional, melhor será possível avaliar o desempenho do modelo. Dessa forma, neste estudo, foi adotada a equação r_t = c + ε _t (regressão do retorno apenas sobre uma constante) como padrão para o cálculo da média condicional em todas as oportunidades em que se fez necessário.

Definida a equação da média condicional, os erros gerados serviram de insumo para as seguintes parametrizações dos modelos de volatilidade condicionada da família ARCH:

Por fim, todas as estimações dos parâmetros e cálculos dos modelos de volatilidade condicionada foram efetuadas com o auxílio do software econométrico Eviews^TM (versão 4.0), ressaltando-se que: a) os modelos foram estimados pelo método de log-máxima verossimilhança, sob a assunção de que os resíduos condicionais são normalmente distribuídos; b) para minimizar as possibilidades de erros oriundos da detecção de máximos locais em vez de máximos globais da função de log-máxima verossimilhança, ajustou-se o critério de convergência do processo iterativo para 10^-5 e o número máximo de iterações para 8.000; c) o algoritmo de otimização escolhido foi BHHH (Berndt-Hall-Hall-Hausman), com o critério selecionado para privilegiar a precisão e não a velocidade; d) foi utilizada a opção de backcasting, sugerida pelo Manual do Usuário de Eviews^TM; e, e) as previsões foram efetuadas com a rotina forecast, ajustando-se para o método dinâmico.

Como a volatilidade não pode ser diretamente observada, faz-se necessária a definição de um processo para se determinar esse parâmetro ou, em outras palavras, para se calcular o que será considerado como a "volatilidade real". Para este trabalho, como o foco de observação será a previsão da volatilidade diária para um horizonte máximo de dois dias úteis, optou-se pelo procedimento utilizado por Pagan e Schwert (1990), Day e Lewis (1992) e Bertucci (1999), que consiste em calcular a volatilidade instantânea por intermédio da seguinte equação:

onde:

σ²_it é a volatilidade instantânea para o dia t;

r_t é o retorno observado no dia t; e

é o retorno médio observado em cada período das subamostras, consideradas aqui como janelas de 149 observações.

No que tange à capacidade preditiva dos modelos de volatilidade, este estudo empírico restringiu-se a comparações de desempenho "fora da amostra", com ênfase nas previsões para D+1 e D+2, seguindo o atual contexto de negociação realizado no âmbito da Clearing de Câmbio da BM&F. Os cálculos com o modelo EWMA foram realizados sempre com a utilização de uma janela móvel de 149 observações. Para os modelos da família ARCH, foram empregados dois procedimentos: a utilização de uma janela móvel de 149 observações, visando comparar a capacidade preditiva em condições iguais às fixadas para o modelo EWMA, e a utilização de todas as observações disponíveis até a data imediatamente anterior ao processo de cálculo, buscando-se obter todo o potencial dos modelos auto-regressivos de heteroscedasticidade condicional.

A operacionalização dos cálculos de predição de volatilidade, para o caso em que se utilizou uma janela móvel para todos os modelos, seguiu os seguintes passos: a) o período de 21 de agosto de 2001 a 20 de março de 2002, correspondendo aos primeiros 149 retornos da série, foi utilizado para o cálculo inicial dos modelos; b) como o modelo EWMA não necessita de estimação de parâmetros, as 149 observações iniciais foram diretamente utilizadas para se calcular a previsão para D+1 (21 de março de 2002). Entretanto, como o referido modelo possibilita apenas a previsão do tipo "um passo a frente", seguindo o procedimento adotado pela Câmara, o valor da volatilidade atribuído para D+2 (22 de março de 2002) foi idêntico ao calculado para D+1 (σ²_D+2 = σ²_D+1); c) para os modelos da família ARCH, a estimação dos parâmetros foi efetuada com os dados da subamostra, gerando-se, a partir daí, as equações que foram utilizadas na predição da volatilidade. Como os modelos dessa família permitem a projeção para períodos superiores a um dia, foi utilizada a rotina forecast do Eviews^TM nas previsões para D+1 e D+2; d) na segunda iteração, deslocou-se a amostra em um dia útil, retirando-se a observação mais antiga e acrescentando-se uma nova, estabelecendo-se um novo conjunto com 149 observações. Com base nessa "nova amostra", foram refeitos os passos descritos nos itens b) e c); e) esse processo iterativo de retirada da observação mais antiga e inclusão de uma nova repetiu-se até 26 de setembro de 2003, dois dias úteis antes do final da amostra. O resultado final desse procedimento foi uma matriz de oito colunas (previsões para D+1 e D+2 realizadas por cada um dos modelos) e 385 linhas (número total de iterações).

Já os cálculos de predição de volatilidade utilizando-se todas as informações anteriormente disponíveis, efetuados apenas com os modelos ARCH, seguiram um procedimento similar ao anteriormente descrito. A diferença reside no fato de que, a cada iteração, uma observação era acrescentada à subamostra, não sendo retirada nenhuma. O resultado final desse procedimento foi uma matriz de seis colunas (previsões para D+1 e D+2 realizadas por cada um dos modelos da família ARCH) e 385 linhas (número total de iterações).

Dos dois procedimentos descritos anteriormente, realizados por meio de rotina desenvolvida pelos autores deste artigo no software Eviews^TM, obteve-se uma matriz de previsões de volatilidade de 14 colunas e 385 linhas, contendo os dados necessários para se avaliar o desempenho preditivo dos modelos selecionados. Passou-se, então, para o cálculo das medidas de avaliação de desempenho preditivo.

Conforme descrito anteriormente, adotou-se a "volatilidade instantânea" (Equação 8) como medida representativa (proxy) da volatilidade diária. Assim, para o período compreendido entre 21 de março de 2002 e 30 de setembro de 2003, foram calculadas as volatilidades diárias por meio do seguinte procedimento: a) definiu-se uma subamostra composta pelos 149 retornos compreendidos entre os dias 22 de agosto de 2001 e 21 de março de 2002; b) calculou-se a média aritmética da subamostra, decrementando-a do retorno observado no último dia. Elevou-se ao quadrado o resultado dessa operação para se obter a volatilidade instantânea para o último dia da amostra; c) deslocou-se a subamostra retirando-se a observação mais antiga e acrescentando-se uma nova, mantendo-se um total de 149 observações; d) repetiram-se os itens b) e c) até o dia 30 de setembro de 2003, obtendo-se um vetor composto pelas 386 volatilidades instantâneas do período.

O terceiro passo consistiu em estabelecer as medidas de avaliação de desempenho da capacidade preditiva dos modelos. Com base na matriz de previsões de volatilidade e do vetor de volatilidades instantâneas, foi possível efetuar a avaliação acerca do desempenho dos modelos na predição por intermédio das seguintes medidas:

onde:

EAM_D+n é o erro absoluto médio nas previsões para D+n;

REQM_D+n é a raiz do erro quadrático médio nas previsões para D+n;

Theil-U_D+n é a estatística Theil-U referente às previsões para D+n; σ²in_D+k é a volatilidade instantânea para a D+k (k=[1,2]) em relação à observação n;

σ2 pn_D+k é a volatilidade prevista para a D+k (k=[1,2]) em relação à observação n;e

σ²RWn_D+k é a volatilidade segundo o modelo escolhido como benchmark - random walk - para a D+k (k=[1,2]), em relação à observação n.

Foram ainda efetuados dois testes adicionais: a verificação de viés de predição, e presença de previsões de volatilidades negativas por parte dos modelos. O viés de predição foi verificado obtendo-se, para cada modelo, o percentual de vezes que a volatilidade prevista foi superior à volatilidade real. O teste de volatilidades negativas foi necessário uma vez que os testes de desempenho utilizados não capturam esta anomalia, o que poderia ocorrer nos modelos GARCH e TARCH.

Finalmente, procurou-se comparar os valores apontados pelos índices econométricos (I_eco) com os valores estipulados para o índice da taxa de câmbio contratualmente garantido (I_est), estipulado pela Câmara. Os I_eco (avaliação objetiva) são calculados segundo a abordagem paramétrica sugerida por Jorion, e adotada pela Câmara, para a determinação do valor em risco com um nível de confiança de 99%, considerando-se a volatilidade prevista por um modelo. Assim, foram calculados os valores dos I_eco, tanto para D+1 quanto para D+2, utilizando-se os valores previstos para a volatilidade por cada um dos modelos estudados. Os valores estipulados (I_est) foram obtidos diretamente junto à Câmara, sendo definidos pelo seu Comitê após uma análise concomitante do índice econométrico calculado e da conjuntura macroeconômica (avaliação subjetiva). De posse dos referidos valores, procedeuse a uma avaliação do comportamento dos índices durante o período da amostra, cujos resultados serão discutidos na próxima seção.

ANÁLISE DOS RESULTADOS

A Tabela 1 apresenta os resultados da primeira fase da investigação empírica, na qual se procurou observar as principais características estatísticas da série de retornos da taxa de câmbio. Para tanto, foram calculadas as estatísticas descritivas e efetuados testes ARCH-LM para os três subperíodos anteriormente definidos e para o conjunto de observações da amostra.

A análise da estatística Jarque-Bera demonstra que, excetuando-se o subperíodo de 2 de janeiro de 2003 a 30 de setembro de 2003, no qual a economia brasileira esteve sob um rigoroso processo de ajuste monetário e fiscal, a hipótese de normalidade dos retornos é fortemente rejeitada. Mesmo no referido subperíodo, a assimetria positiva (indicando uma cauda direita longa) e a diferença entre a média e a mediana sugerem a rejeição da hipótese de normalidade, ainda que ao nível de significância de 34,48%. Ademais, pelos valores de curtose apresentados, a função densidade de probabilidade da série é sempre leptocúrtica.

Percebe-se, também, que o período de maior instabilidade é aquele que antecede as eleições majoritárias brasileiras de 2002 e se estende até a posse do novo governo (período de 22 de abril de 2002 a 30 de dezembro de 2002). Nesse período, observa-se o maior desvio padrão e os valores extremos (máximo e mínimo) de toda a amostra. A razão desse comportamento pode ser atribuída às incertezas econômicas e políticas vivenciadas em virtude da troca de comando no governo e a ascensão, pela primeira vez, de um partido de oposição à presidência da República até então no poder.

Em relação à heteroscedasticidade da série, os p-values observados nas estatísticas ARCH-LM indicam uma forte presença de autocorrelação dos resíduos quadráticos no conjunto completo da amostra e em dois dos subperíodos. Apenas no subperíodo que antecede a entrada em operação da Câmara (21 de agosto de 2001 a 19 de abril de 2002) é que a hipótese de heteroscedasticidade teria que ser aceita em um nível de significância um pouco menor, ainda que representativo (grau máximo de certeza de 77,02%). Dessa forma, os testes sugerem a rejeição da hipótese de homoscedasticidade na série de retornos da taxa de câmbio no contexto deste trabalho.

Além disso, a partir da inspeção visual da Figura 1, que expressa o comportamento dos retornos diários da taxa de câmbio para todo o período da amostra, pode-se perceber que a série temporal demonstra, ainda, a presença de "clusters de volatilidade" e de um certo grau de assimetria na intensidade das respostas a eventos positivos e negativos. Verifica-se, portanto, que as características estatísticas da série de retornos da taxa de câmbio avalizam plenamente a utilização de modelos de volatilidade condicionada como os da classe ARCH.

Os resultados do desempenho preditivo dos diversos modelos estudados podem ser observados na Tabela 2. Uma primeira análise dessa tabela revela a superioridade dos modelos da família ARCH em relação ao modelo EWMA. Apenas o modelo GARCH (1,1), quando utilizado de forma restrita (ou seja, com a parametrização sugerida para o modelo EWMA pelo Manual do RiskMetrics^TM - janela móvel com 149 observações), apresenta uma capacidade preditiva inferior em todos os critérios. Além da restrição imposta, uma possível hipótese para o fraco desempenho do modelo GARCH é a assimetria dos retornos, a qual não é capturada pelo modelo.

Corrobora essa hipótese a supremacia dos modelos assimétricos (EGARCH e TARCH), segundo todos os critérios e implementações, restritas ou não. Observa-se que para D+1, todos os critérios apontam o modelo TARCH (1,1) como o de melhor desempenho preditivo. Já para D+2, pode-se considerar que há um "empate técnico" entre os modelos TARCH (1,1) e EGARCH (1,1). Segundo o erro absoluto médio (EAM), o melhor desempenho é do modelo EGARCH. Por sua vez, segundo a raiz do erro quadrático médio (REQM), o melhor desempenho é do modelo TARCH. Finalmente, tomando-se por base a estatística Theil-U, a diferença de desempenho não é suficientemente representativa para que se possa concluir pela superioridade de um modelo sobre o outro.

É importante observar que a estatística Theil-U proporciona, adicionalmente, uma comparação entre os modelos pesquisados e o modelo escolhido como benchmark, nesse caso, o Random Walk (RW). Por construção, o valor da estatística tem o seu limite inferior igual a zero - caso em que o modelo em teste apresenta todas as predições corretas -, sendo sempre igual a 1 para o modelo escolhido como benchmark. Dessa forma, pode-se verificar que todos os demais modelos apresentam um desempenho superior ao modelo RW, uma vez que as respectivas estatísticas Theil-U são sempre menores que a unidade. Entretanto, deve-se destacar que a diferença de desempenho para D+2 decresce substancialmente em comparação aos resultados observados para D+1.

A Tabela 2 permite ainda uma análise específica em relação aos três modelos da classe ARCH. Comparando-se o desempenho preditivo quando se emprega um número restrito de observações (149 neste caso) e quando são utilizadas todas as informações anteriormente disponíveis, verifica-se que todos os modelos se beneficiam com o aumento do número de observações utilizadas para a determinação de seus estimadores. Assim, as evidências empíricas apontam para a impropriedade de implementar restrições ao número de observações para o caso em que se optar pela utilização dos referidos modelos.

Conforme discutido anteriormente, devido ao fato de todos os critérios utilizados de avaliação de desempenho serem simétricos, fazem-se necessárias duas verificações adicionais - em relação às predições de volatilidade ponto a ponto - de forma a melhor validar os resultados obtidos: a ocorrência de volatilidades negativas e a presença de viés de estimação. No que tange à presença de volatilidades negativas, testou-se a presença de ocorrências da espécie nas predições efetuadas pelos modelos GARCH e TARCH, que, dentre os modelos concorrentes, eram os únicos suscetíveis a esse tipo de erro. Como não foi constatada nenhuma falha nesse sentido, não há sob esse aspecto nenhuma objeção em relação aos resultados obtidos pelos critérios de avaliação.

Em relação à presença de viés de predição, todos os modelos revelaram uma tendência a superestimar a volatilidade futura. Como os valores observados não apresentaram diferenças significativas entre os modelos, também sob esse aspecto não se podem oferecer objeções em relação aos resultados obtidos pelos critérios de avaliação.

Em resumo, os resultados até aqui evidenciados indicam que todos os modelos pesquisados tendem a superestimar a volatilidade futura, sendo que os da família ARCH apresentam, no contexto desta pesquisa, um desempenho preditivo superior ao do modelo EWMA. Ademais, analisando-se em conjunto os resultados para D+1 e D+2, conclui-se também pela superioridade dos modelos assimétricos - em especial do modelo TARCH (1,1) - na tarefa de predizer a volatilidade dos retornos da taxa de câmbio.

No que se refere à aplicação, no contexto da Clearing de Câmbio da BM&F, das volatilidades previstas, a comparação do comportamento diário entre os índices de taxa de câmbio contratualmente garantidos estipulados pela Câmara (I_est) e os índices calculados com base no modelo de valor em risco (I_eco) permitiu a identificação de três fases distintas, independentemente do modelo analisado: o subperíodo compreendido entre 22 de abril de 2002 e 26 de julho de 2002, representando o início das atividades da Câmara; o subperíodo compreendido entre 29 de julho de 2002 e 22 de outubro de 2002, que concentra os maiores valores calculados para o I_eco; e o subperíodo compreendido entre 23 de outubro de 2002 a 26 de setembro de 2003, correspondendo ao final da amostra.

Na primeira fase, para todos os modelos, nota-se certo equilíbrio entre as avaliações subjetivas e objetivas na determinação do índice final. Percebe-se uma tendência, tanto para D+1 quanto para D+2, que as elevações no índice econométrico sejam acompanhadas por elevações no índice estipulado pela Câmara. Em relação aos impactos esperados caso a Câmara optasse pela substituição do modelo EWMA, verificase que para D+1 os valores calculados para I_eco pelos demais modelos não ensejariam mudanças no I_est. Já para D+2, os modelos da família ARCH apontam para uma elevação no valor do I_est.

A segunda fase apresenta a maior dissonância entre as avaliações objetivas e subjetivas. Logo no início, enquanto todos os modelos econométricos sinalizam para um aumento substancial nos valores do I_est, a avaliação subjetiva parece considerar tal aumento por demais exagerado. Assim, verifica-se, principalmente para D+2, que os valores estipulados ficam aquém, em maior ou menor grau, dependendo do modelo, dos valores indicados pela avaliação objetiva. No que tange aos impactos esperados em razão da substituição de modelos, os da família ARCH apontam, com mais intensidade que o modelo EWMA, para a necessidade de elevação do valor I_est, o que, de alguma forma, poderia sensibilizar o Comitê da Câmara.

A última fase é marcada pelo predomínio crescente da avaliação subjetiva em relação à avaliação objetiva, na medida em que se aproxima do final da amostra. Enquanto todos os modelos econométricos apresentam uma tendência decrescente para os valores do I_eco, o valor do I_est se mantém constante, seja para D+1 ou D+2. Nessa fase, dada a diferença entre os valores calculados para o I_eco e os valores definidos pela Câmara, não há razão para se supor que uma mudança no modelo econométrico utilizado pudesse causar algum impacto no I_est.

Com base nas análises anteriores e considerando que o índice econométrico é calculado segundo procedimentos estatísticos que, em princípio, predizem a variação da taxa de câmbio em um nível de confiança de 99%, a tendência observada de distanciamento progressivo entre o índice de variação da taxa de câmbio estipulado pela Câmara e o apontado pelo método econométrico sugere um excesso de conservadorismo ou um relativo descrédito nos procedimentos estatísticos adotados. Nota-se, também, que no período de maior volatilidade, a avaliação subjetiva subestimou, por vezes, a volatilidade real, contrariando as indicações dos modelos econométricos. Dessa forma, podese concluir que, mantida essa tendência, nenhum efeito prático deve ser esperado de uma eventual substituição do modelo EWMA por qualquer dos modelos concorrentes. A Figura 2 apresenta os gráficos comparativos entre a volatilidade real, I_est e I_eco, para D+1 e D+2, segundo cada um dos modelos estudados.

CONCLUSÕES

No contexto desta pesquisa, conclui-se que o modelo TARCH (1,1) é o que apresenta melhor desempenho preditivo, sendo acompanhado de perto pelo modelo EGARCH (1,1). Em seguida, tem-se o modelo GARCH (1,1) e, por último, o modelo EWMA, ressaltando-se que todos os modelos apresentam uma tendência a superestimar a volatilidade futura. Conclui-se, ainda, pela inadequação, nos modelos da classe ARCH, de implementar restrições em relação ao número de observações utilizadas no processo de estimação dos parâmetros da equação da volatilidade condicional, haja vista a queda acentuada no desempenho preditivo, em especial do modelo GARCH. Ademais, verificouse que o modelo Random Walk não deve ser considerado como uma alternativa viável aos modelos econométricos utilizados neste estudo, uma vez que, conforme demonstrado pela estatística Theil-U, estes últimos apresentam um desempenho significativamente superior.

Os resultados acima confirmam o fato, preconizado pela teoria, de que a abordagem econométrica representada pelos modelos da classe ARCH é a mais adequada para o tratamento de séries que apresentam heteroscedasticidade, agrupamentos de volatilidade e leptocurtose, características verificadas na série de retornos das taxas de câmbio utilizada neste trabalho. Constatou-se também a característica de reação assimétrica a boas e más notícias, o que guarda relação direta com o melhor desempenho preditivo observado para os modelos TARCH e EGARCH, de resposta assimétrica, comparativamente aos modelos GARCH e EWMA, de resposta simétrica.

Em relação ao contexto operacional da Clearing, verificou-se um excesso de garantias em boa parte da amostra, sendo que, no período de maior volatilidade, a avaliação subjetiva adotada pela Câmara subestimou, por vezes, a volatilidade real, contrariando as indicações dos modelos econométricos. Essa supremacia da avaliação subjetiva sobre a objetiva pode estar fundamentada em um relativo descrédito nos procedimentos estatísticos ora adotados pela Câmara. Nesse caso, à semelhança do ocorrido com a Clearing de Derivativos da BM&F, sugere-se a reavaliação dos modelos utilizados de forma a minimizar o volume total de garantias exigidas sem comprometer a segurança do sistema.

Para finalizar, devem ser ressaltadas as limitações deste trabalho, como o restrito número de observações decorrente do curto prazo de operação da Câmara; a presença de pontos controversos e ainda inconclusos no meio acadêmico, inclusive a definição da volatilidade real e a determinação da medida de desempenho preditivo mais apropriada para modelos de volatilidade e o caráter probabilístico dos resultados. Dessa forma, os resultados aqui obtidos não têm a pretensão de sinalizar o fim das discussões acerca do melhor modelo de predição de volatilidade, mas sim recolocar em discussão um tema de crucial importância para o gerenciamento de riscos e a precificação de operações no mercado interbancário brasileiro. Nesse sentido, sugere-se a replicação do presente trabalho, ampliando-se o intervalo da amostra e utilizando-se dados de alta freqüência, bem como o desenvolvimento de estudos visando caracterizar a tendência a superestimar a volatilidade dos modelos avaliados.

NOTAS

1 Eviews^TM é um pacote estatístico desenvolvido pela empresa Quantitative Micro Software, que proporciona um conjunto de funções predefinidas para suporte à análise econométrica de dados. Mais informações podem ser obtidas no site www.eviews.com.

Artigo recebido em 23.08.2004.

Aprovado em 18.07.2005.

Clayton Peixoto Goulart

Pesquisador do NUFI/UFMG-CEPEAD, do CEGE/UFMG-CEPEAD e do IBMEC-MG. Mestre em Administração pela UFMG. Interesses de pesquisa nas áreas de mercado financeiro, de capitais e de derivativos e métodos e modelos matemáticos, econométricos e estatísticos.

E-mail: [e-mail]

Endereço: Rua Francisco Bicalho, 2385, ap. 301, Caiçaras, Belo Horizonte - MG, 30720-340.

Hudson Fernandes Amaral

Professor da UFMG-CEPEAD. Doutor em Administração pela Université Pierre Mendes - France. Interesses de pesquisa nas áreas de finanças corporativas, mercado financeiro, de capitais e de derivativos, finanças públicas.

E-mail: [e-mail]

Endereço: Rua Curitiba 832, 10º andar, Centro, Belo Horizonte - MG, 30170-120.

Luiz Alberto Bertucci

Professor da UFMG-CEPEAD. Doutor em Administração pela FGV-EAESP. Interesses de pesquisa nas áreas de mercado financeiro, de capitais e de derivativos, finanças corporativas, estratégia.

E-mail: [e-mail]

Endereço: Rua Curitiba 832, 10º andar, Centro, Belo Horizonte - MG, 30170-120.

Aureliano Angel Bressan

Professor da UFMG-CEPEAD. Doutor em Economia Aplicada pela UFV. Interesses de pesquisa nas áreas de técnicas de avaliação de investimentos em renda variável, análise de séries temporais em finanças, métodos e modelos matemáticos, econométricos e estatísticos.

E-mail: [e-mail] Endereço: Rua Curitiba 832, 10º andar, Centro, Belo Horizonte - MG, 30170-120.

• ALEXANDER, C. (Ed.). Risk Management and Analysis. v. 1: Measuring and Modeling Financial Risk. West Sussex: John Wiley, 1998.
• ALEXANDER, C. Market Models: A Guide to Financial Data Analysis West Sussex: John Wiley, 2002.
• BERTUCCI, L. A. Avaliação de opções sob consideração de volatilidades históricas, implícitas e condicionadas - o caso Telebrás na Bovespa 1999. 150 p. Tese (Doutorado em Administração) - Escola de Administração de Empresas de São Paulo: Fundação Getúlio Vargas. São Paulo.
• BM&F: Bolsa de Mercadorias & Futuros. Câmara de Registro, Compensação e Liquidação de Operações de Câmbio BM&F: regulamento, manual de operações, manual de gerenciamento de risco e órgãos auxiliares de administração. São Paulo: BM&F, 2002.
• BOLLERSLEV, T. Generalized autoregressive conditional heteroskedasticity. Journal of Econometrics, v. 31, p. 307-327, 1986.
• BOLLERSLEV, T.; CHOU, R. Y.; KRONER, K. F. ARCH modeling in finance: a review of theory and empirical evidence. Journal of Econometrics, v. 52, p. 5-59, 1992.
• BROOKS, C. Introductory Econometrics for Finance Cambridge, UK: Cambridge University Press, 2002.
• DAY, T.; LEWIS, C. Stock market volatility and the information content of stock index options. Journal of Econometrics, v. 52, p. 267-287, 1992.
• ENGLE, R. F. Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation. Econometrica, v. 50, n. 4, p. 987-1007, 1982.
• FAMA, E. The behavior of stock market prices. Journal of Business, n. 38, p. 34-105, 1965.
• GLOSTEN, L. R.; JAGANATHAN, R.; RUNKLE, D. E. On the relation between the expected value and the volatility of the nominal excess returns on stocks. Journal of Finance, v. 48, p. 1779-1801, 1993.
• GOURIÉROUX, C. Arch Models and Financial Applications New York: Springer-Verlag, 1997.
• JORION, P. Value at risk: a nova fonte de referência para o controle do risco de mercado São Paulo: Cultura, 1998.
• MANDELBROT, B. The variation of certain speculative prices. Journal of Business, v. 36, p. 394-419, 1963.
• MERTON, R. C. On estimating the expected return on market: an exploratory investigation. Journal of Financial Economics, v. 8, p. 323-361, 1980.
• MARKOWITZ, H. Portfolio selection. Journal of Finance, v. 7, n. 1, p. 77-91, 1952.
• MORGAN GUARANTY TRUST COMPANY & REUTERS LTD. RiskMetrics: Technical Document 4^th. ed., New York: RiskMetrics, 1996.
• NELSON, D. B. Conditional heteroskedasticity in asset return: a new approach. Econometrica, v. 59, n. 2, p. 347-370, 1991.
• PAGAN, A. R.; SCHWERT, G. W. Alternative models for conditional stock volatility. Journal of Econometrics, v. 45, p. 267-290, 1990.
• POON, S.; GRANGER, C. W. J. Forecasting volatility in financial markets: a review. Journal of Economic Literature, v. 41, p. 478-539, 2003.

Datas de Publicação

• Publicação nesta coleção
  10 Dez 2014
• Data do Fascículo
  Dez 2005

Histórico

• Aceito
  18 Jul 2005
• Recebido
  23 Ago 2004