Acessibilidade / Reportar erro

Existe convergência espacial da produtividade agrícola no Brasil?

Resumos

Este trabalho analisa a convergência espacial da produtividade agrícola da terra para as microrregiões brasileiras entre o período de 1991 a 2003. Tal análise avalia se há evidências para a existência de convergência beta da produtividade agrícola, controlando-se explicitamente para efeitos espaciais. Para tanto, conduziu-se preliminarmente uma análise exploratória de dados espaciais, que detectou a presença de autocorrelação espacial para a produtividade da terra. Numa etapa posterior, desenvolveram-se modelos econométricos espaciais para a análise da convergência.

beta convergência; produtividade da terra; econometria espacial; autocorrelação espacial


This work analyzes the convergence of land productivity for Brazilian regions from 1991 to 2003. It evaluates if there is or not beta convergence of the agricultural productivity, controlling explicitly for spatial effects. To do so, an exploratory spatial data analysis was previously implemented and the presence of spatial autocorrelation for the land productivity was detected. In addition, spatial econometric models were developed for the convergence analysis.

beta convergence; land productivity; spatial econometrics; spatial autocorrelation


Existe convergência espacial da produtividade agrícola no Brasil?

Eduardo Simões de AlmeidaI; Fernando Salgueiro PerobelliII; Pedro Guilherme Costa FerreiraIII

IProfessor Adjunto da FEA/UFJF e pesquisador do CNPq, email: eduardo.almeida@ufjf.edu.br

IIProfessor Adjunto da FEA/UFJF e pesquisador do CNPq, email: fernando.perobelli@ufjf.edu.br

IIIMestrando em Economia do Curso de Mestrado em Economia Aplicada FEA/UFJF, email: pgcferreira@gmail.com

RESUMO

Este trabalho analisa a convergência espacial da produtividade agrícola da terra para as microrregiões brasileiras entre o período de 1991 a 2003. Tal análise avalia se há evidências para a existência de convergência b da produtividade agrícola, controlando-se explicitamente para efeitos espaciais. Para tanto, conduziu-se preliminarmente uma análise exploratória de dados espaciais, que detectou a presença de autocorrelação espacial para a produtividade da terra. Numa etapa posterior, desenvolveram-se modelos econométricos espaciais para a análise da convergência.

Palavras-chave: b convergência; produtividade da terra; econometria espacial; autocorrelação espacial.

Classificação JEL: R12, C31, Q10.

ABSTRACT

This work analyzes the convergence of land productivity for Brazilian regions from 1991 to 2003. It evaluates if there is or not b convergence of the agricultural productivity, controlling explicitly for spatial effects. To do so, an exploratory spatial data analysis was previously implemented and the presence of spatial autocorrelation for the land productivity was detected. In addition, spatial econometric models were developed for the convergence analysis.

Keywords: b convergence; land productivity, spatial econometrics, spatial autocorrelation.

JEL Classification: R12, C31, Q10.

1. Introdução

A importância da agricultura para o crescimento econômico tende a ser desigual para as diferentes regiões. Segundo Johnston & Mellor (1961), as funções da agricultura podem ser listadas como: a) fornecer alimentos; b) prover capital, especialmente para a expansão do setor não-agrícola; c) oferecer mão-de-obra para o crescimento e diversificação de atividades na economia; d) gerar ganhos cambiais; e, e) constituir mercado para os produtos do setor não-agrícola. Na visão de Johnston & Mellor (1961), as funções desempenhadas pela agricultura estão relacionadas aos primeiros estágios de desenvolvimento econômico. Isso ocorre, em parte, devido à diminuição da participação da agricultura no produto e no emprego regional. No entanto, a agricultura pode contribuir para o desenvolvimento da economia como um todo, mesmo nos estágios mais avançados.

Deste modo, fica evidente que analisar o comportamento da agricultura é de suma importância para a economia de um País, principalmente no Brasil, onde o setor participa de forma significativa tanto nas exportações quanto no mercado de trabalho. Uma das maneiras de avaliar tal desempenho é acompanhar o comportamento da produtividade agrícola ao longo do tempo e através das regiões, avaliando se a sua evolução mostra convergência. Conforme Lopes (2004), " a convergência é um processo em que uma mesma variável (por exemplo, renda per capita ou produtividade da terra) apresenta diferentes valores entre países, regiões ou estados, mas essa diferença se reduz ao longo do tempo, indicando que a desigualdade diminui" .

A discussão de convergência da renda atraiu muitos estudiosos e existe um grande número de trabalhos empíricos que discutem essa questão. Em outros termos, a questão do crescimento mais rápido das regiões mais pobres do que as mais ricas foi exaustivamente discutida na literatura.1 1 Embora a maioria dos trabalhos que analisou a questão da convergência de renda encontrou evidências de convergência com nivelamento para cima, é possível, teórica e empiricamente, admitir a existência de convergência com nivelamento para baixo. Todavia, como ficará claro adiante, não é razoável admitir convergência com nivelamento para baixo da produtividade agrícola no Brasil no período analisado. O trabalho seminal de Baumol (1986) examinou a convergência no período de 1870 a 1979 entre 16 países industrializados, conforme a seguinte expressão:

em que ln (Y / N) é o logaritmo da renda per capita, e é o termo de erro e i é o indexador para os diversos países.

Segundo Baumol (1986), se existir convergência absoluta, b será negativo, ou seja, os países com renda inicial maior terão menores taxas de crescimento. Com isso, Baumol procura mostrar que ao longo de um período de tempo, as rendas dos diversos países estariam convergindo para uma renda comum. Portanto, se os países de rendas menores crescem mais, a tendência é que, ceteris paribus, tais rendas se igualem no tempo. Cabe ressaltar que muito se avançou na discussão sobre a convergência absoluta desde o trabalho de Baumol (1986). Todavia, não é do escopo deste artigo discutir todas as variantes desta linha de pesquisa.2 2 Consultar Barro e Sala-i-Martin (1995) para uma discussão avançada sobre as convergências absoluta e condicional.

É importante elencar a existência de trabalhos, como o de Rey e Montouri (1999), que aplicam a análise da convergência da renda, adotando a abordagem da econometria espacial. Dall’erba (2003) estima a convergência da produtividade do trabalho com base em 48 regiões da Espanha, durante os anos de 1980-1996. Entre os artigos relacionados ao tema, no Brasil, destaca-se a tese de Lopes (2004), que analisa a evolução da produtividade da terra das onze principais culturas brasileiras, avaliando o processo de convergência dessa variável entre os estados produtores. Magalhães et al. (2000), por sua vez, discute a convergência da renda através dos estados brasileiros durante o período de 1970 a 1995.

A investigação de produtividade média agrícola sob o prisma do modelo de b convergência, controlando-se para os efeitos espaciais, é explorada no presente trabalho de forma inédita tanto para o Brasil quanto para outros países. Assim, o objetivo deste trabalho é adaptar a idéia apresentada por Baumol (1986) para entender o comportamento da produtividade média agrícola na terra. Nesse sentido, ao formalizar uma equação de convergência absoluta a la Baumol para o índice de produtividade média agrícola na terra, buscar-se-á verificar se as regiões com menor grau de produtividade agrícola estão convergindo para as regiões com maior grau de produtividade. Enfim, a idéia de utilizar a análise de convergência no trabalho é testar se os índices de produtividade agrícola estão ficando mais homogêneos no país.

De acordo com Lopes (2004), existem basicamente três razões para se verificar a convergência da produtividade agrícola, a saber, a ocorrência de mudanças estruturais no processo de produção (vinculadas, sobretudo a inovações tecnológicas e/ou incentivadas por políticas públicas, como, por exemplo, o processo de apreciação do câmbio na maior parte dos anos noventa e a abertura econômica, entre outros); o fenômeno da difusão tecnológica de novas sementes, novos sistemas de plantio, novos implementos agrícolas etc; e a eliminação de obstáculos ao crescimento da produção, como restrições ambientais, falta de adequada logística e carência de financiamento, dentre outros.

Outro aspecto relevante no artigo é o período selecionado para a análise (1991 a 2003). Esse período pode ser considerado promissor para a agricultura em termos de expansão da produção e de obtenção de ganhos de produtividade. Sem sombra de dúvida, neste período, consolida-se a posição do País como um dos mais importantes e competitivos produtores agrícolas do mundo.

Ademais, esse intervalo de tempo testemunhou várias mudanças estruturais que marcaram indelevelmente o desenvolvimento agrícola no Brasil, tais como a abertura econômica, a estabilização dos preços, a apreciação cambial (1994-1999) e a posterior depreciação cambial (a partir de fevereiro de 1999) que acarretou a expansão extraordinária das exportações, impulsionada também pela expansão do crescimento econômico do mundo. Deste modo, torna-se relevante dividir a estimativa da convergência em sub-períodos que abrangem tais medidas macroeconômicas.

A primeira estimativa testa se no período de 1991 a 1994, existe convergência da produtividade média agrícola. Isso porque esse é o período do início da abertura econômica brasileira até a adoção do Plano Real. Nessa época, de acordo com documento do DIEESE (2004), a moeda interna estava desvalorizada em relação ao dólar, constituindo um processo de superinflação crônica.

A segunda estimativa de convergência engloba o período de 1995 até 1999. Esse é o período após a implantação do Plano Real. Nesse período, houve um controle mais rigoroso da inflação, ao passo que o câmbio apreciou-se, transformando-se na âncora cambial da economia. É importante observar que esse é um período de altos e baixos na economia brasileira, considerando-se que, no ano de 1995, o país vivia um momento de estabilidade com crescimento econômico, fazendo com que o nível de atividade acelerasse e o desemprego diminuísse. No entanto, tal cenário econômico foi interrompido pelas sucessivas crises internacionais ocorridas no cenário mundial, que provocaram a fuga de capitais do Brasil. Isso levou o governo a promover sucessivos aumentos na taxa de juros, trazendo repercussões negativas para a economia brasileira.

A terceira estimativa de convergência aplica-se ao período de 2000 a 2003 e objetiva testar qual foi o impacto da moeda desvalorizada na produtividade média agrícola em relação ao insumo terra. Por último, a quarta estimativa engloba todo o período de 1991 a 2003 em que ocorreram profundas mudanças macroeconômicas na economia brasileira.

Outro ponto observado é a evolução da produtividade média no Brasil e nos estados ao longo dos sub-períodos. Isto porque poderia, eventualmente, ocorrer convergência da produtividade agrícola no país, caso houvesse queda na produtividade tanto no Brasil quanto nos estados, proporcionando uma conclusão enganosa dos dados. Todavia, como se verifica na Tabela 1, isso não acontece. Segundo Considera et al. (2002), são os ganhos de produtividade a principal alavanca do processo de crescimento do setor agrícola brasileiro na década de 1990.

Vale observar também a importância da utilização de econometria espacial nos estudos de convergência. Segundo Quah (1996), a maioria dos trabalhos sobre convergência utiliza dados regionais, enquanto grande parte dos trabalhos trata as regiões como unidades isoladas, ou seja, não considerando a interação com seu entorno. Ainda, de acordo com Rey e Montouri (1999), o desenvolvimento de ferramentas de análise exploratória de dados espaciais, bem como as técnicas econométricas espaciais, permitiu uma nova percepção da dinâmica geográfica dos padrões de crescimento da renda no tempo.3 3 De acordo com Anselin (1999, p. 2), " econometria espacial é um sub-campo da econometria que lida com as complicações causadas pela interação espacial (autocorrelação espacial) e pela estrutura espacial (heterogeneidade espacial), em modelos de regressão para dados na forma de seção cruzada e painel de dados" .

Rey e Montouri (1999) reconheceram que um modelo de convergência, composto por dados organizados em unidades espaciais, deveria levar em consideração os efeitos espaciais que poderiam resultar da interação espacial entre os agentes.

A diferença entre a econometria espacial e a econometria convencional concentra-se na preocupação de se incorporar na modelagem o padrão da interação sócio-econômica entre os agentes num sistema (a autocorrelação espacial), assim como as características da estrutura desse sistema no espaço (a heterogeneidade espacial). Essas interações e as características estruturais geram efeitos espaciais que influenciam vários processos econômicos (Anselin, 1988; Anselin e Bera, 1998).

Não é difícil de perceber que a agricultura é muito sensível aos efeitos espaciais. O desenvolvimento das culturas agrícolas através do espaço é heterogêneo, pois depende de diferentes técnicas de produção, condições de climas, solo e topografia. Por outro lado, os efeitos de interdependência entre diferentes regiões produtoras manifestam-se de diversas formas, tais como ordenação geográfica dos dados, difusão espacial de certos fenômenos que influenciam os vizinhos, processos de competição espacial no âmbito da expansão da fronteira agrícola ou na formação de cinturões agrícolas.

Em suma, este trabalho apresenta uma análise de convergência absoluta para a produtividade média agrícola, controlando-se os efeitos espaciais. Na presença dos efeitos espaciais, os métodos econométricos tradicionais levam à inferência espúria, com graves danos para a formulação de políticas públicas. Nesse sentido, busca-se verificar a existência de convergência da produtividade média para as microrregiões brasileiras, tendo dois objetivos centrais. O primeiro é fornecer uma nova percepção da dinâmica geográfica utilizando para isso métodos de análise exploratória de dados espaciais e de econometria espacial. O segundo, diz respeito à identificação e exposição do grau da convergência na produtividade agrícola que vem ocorrendo no País.

O artigo está organizado em cinco seções, incluindo esta primeira, de caráter introdutório. A próxima seção trata da descrição sucinta do banco de dados. A terceira seção traz a análise exploratória de dados espaciais. Na quarta seção, são apresentados a metodologia e os resultados das regressões de convergência, bem como sua discussão. Na última seção, são abordadas as conclusões de relevo e tecidas as considerações finais.

2. Dados

Este trabalho utiliza os dados da Pesquisa Agrícola Municipal do IBGE (Instituto Brasileiro de Geografia e Estatística), consultados por meio do Sistema IBGE de Recuperação Automática (SIDRA). Para a elaboração do artigo, foram compilados dados sobre a produção agrícola total e a área plantada total para os anos de 1991 a 2003 por microrregião.

Devido à incompatibilidade de unidade de medida, já que na consulta ao SIDRA não estavam disponíveis dados referentes à produção total, em toneladas, para os anos de 1991 a 2000, utilizou-se um conversor de unidade de medida, para transformar as diversas unidades de medidas de grandeza (i.e, mil cachos, mil frutos) para tonelada.

É importante, ainda, ressaltar que foram utilizadas na análise variáveis relativas, pois variáveis absolutas podem levar ao engano na interpretação dos resultados, pois costumam estar correlacionadas ao tamanho da população ou da área das regiões em estudo. Para tanto, foram construídos indicadores parciais de produtividade agrícola média da terra, formados pela razão entre a produção agrícola total (em tonelada) e a área plantada (em hectares) em nível de microrregião para todos os anos entre 1991 e 2003.

3. Análise exploratóra de dados espaciais

3.1. Idéias Gerais

A análise exploratória de dados espaciais (AEDE) trata diretamente de efeitos decorrentes da dependência espacial e da heterogeneidade espacial. Em outras palavras, o objetivo desse método é descrever a distribuição espacial, os padrões de associação espacial (clusters espaciais), verificar a existência de diferentes regimes espaciais ou outras formas de instabilidade espacial (não-estacionariedade) e identificar observações atípicas (i.e. outliers).

Para implementar a AEDE, assim como para aplicar as técnicas de econometria espacial, como será visto mais adiante, é preciso, primeiramente, definir uma matriz de pesos espaciais (W). Essa matriz é a forma de se expressar um determinado arranjo espacial das interações resultantes do fenômeno a ser estudado. Por exemplo, é razoável supor que, no estudo de vários fenômenos, regiões vizinhas tenham uma interação mais forte entre si do que regiões que não são contíguas. De modo similar, regiões distantes entre si teriam uma interação menor. Nesse caso, em que a distância entre as regiões importa na definição da força da interação, seria possível construir uma matriz W baseada na distância inversa entre as regiões a fim de capturar tal arranjo espacial da interação. A escolha da matriz de pesos espaciais é muito importante em uma AEDE, pois os resultados da análise são sensíveis a tal seleção.4 4 Qualquer matriz de pesos espaciais precisa atender às condições de regularidade impostas pela necessidade de invocar as propriedades assintóticas dos estimadores e dos testes. Segundo Anselin (1988), " isso significa que os pesos precisam ser não-negativos e finitos e que correspondam a uma determinada métrica" .

A matriz de peso espacial W, utilizada neste trabalho, está baseada na idéia dos k vizinhos mais próximos, e calculada utilizando a métrica baseada no grande círculo entre os centros das regiões.5 5 De acordo com Le Gallo e Ertur (2003), a escolha de um número fixo de vizinhos mais próximos ao invés do uso de uma matriz simples de contigüidade é preferível, pois é possível evitar alguns problemas metodológicos que podem ocorrer quando há variações no número de vizinhos. Formalmente, tal matriz de pesos espaciais é expressa como:

na qual dij é a distância, medida pelo grande círculo, entre os centros das regiões i e j. Di (k) denota um valor crítico que define o valor de corte para cada região, ou seja, a distância máxima para considerar regiões vizinhas à região i. Ou seja, distâncias acima deste ponto não serão tomadas como vizinhas da região em questão.6 6 Neste estudo, foi utilizado valor de k=15 (ie., quinze vizinhos mais próximos), todavia, a fim de verificar a robustez dos resultados foi testado k=10 e k=20. Vale destacar que os resultados k=10 e k=20 vizinhos mais próximos estão disponíveis para consulta.

Outro mecanismo muito útil para implementar a AEDE e a econometria espacial é o operador de defasagem espacial.7 7 O operador defasagem espacial é oriundo da idéia de defasagem temporal da econometria de séries de tempo. Contudo, o conceito de operador de defasagem espacial não é tão simples e direto como o operador de defasagem em séries de tempo devido à natureza multidirecional do processo de interação no espaço. A solução apresentada na literatura para se construir esse operador consiste na utilização de pesos espaciais para se obter uma medida de variável defasada para uma dada região. O operador de defasagem espacial para uma variável qualquer y, formalmente Wy, pode ser interpretado como sendo a média do valor dessa variável nas regiões vizinhas.8 8 Note-se que a matriz de pesos espaciais W tem dimensão n por n, sendo n o número de unidades espaciais (por ex., microrregiões), ao passo que y é um vetor-coluna com os valores da variável de interesse (por exemplo, produtividade agrícola média).

3.2. Autocorrelação Espacial Global

A autocorrelação espacial pode ser calculada por meio da estatística I de Moran. Essa estatística fornece a indicação formal do grau de associação linear entre os vetores de valores observados de uma variável de interesse no tempo t (zt) e a média ponderada dos valores da vizinhança, ou as defasagens espaciais (Wzt). Valores de I maiores (ou menores) do que o valor esperado E(I)= –1/(n – 1) significa que há autocorrelação positiva (ou negativa). Seguindo Cliff e Ord (1981), em termos formais, a estatística I de Moran pode ser expressa como:

em que zt é o vetor de n observações para o ano t na forma de desvio em relação à média. Mais uma vez, W é a matriz de pesos espaciais, sendo que os elementos wii na diagonal são iguais a zero, enquanto que os elementos wij indicam a forma como a região i está espacialmente conectada com a região j. O termo So é um escalar igual à soma de todos os elementos de W.

Quando a matriz de pesos espaciais é normalizada na linha, isto é, quando os elementos de cada linha somam um, a expressão (3) assume a seguinte forma:

A autocorrelação espacial positiva indica que microrregiões que apresentam elevada produtividade média agrícola são vizinhas de outras microrregiões que também apresentam elevada produtividade média agrícola. Ou, alternativamente, que microrregiões com baixa produtividade média agrícola são circundadas por outras microrregiões também ostentando baixa produtividade média.

A Tabela 2 mostra os resultados do Índice de Moran, o desvio padrão e o z-value para a produtividade média agrícola para os anos selecionados de 1991, 1997 e 2003. Como podem ser observados, são positivos e crescentes todos os valores de uma variável de interesse qualquer z quanto do I de Moran, o que identifica uma autocorrelação espacial global positiva em todos os períodos.

Cabe ressaltar que a estatística I de Moran, utilizando o conceito de k vizinhos mais próximos para os valores de vizinhos utilizados neste trabalho (isto é, quinze e vinte), leva ao mesmo resultado em relação ao sinal e à significância da autocorrelação espacial, ou seja, os resultados são robustos em relação à escolha da matriz de pesos.

3.3. Diagramas de Dispersão de Moran

Segundo Anselin (1996), o diagrama de dispersão de Moran é uma das formas de interpretar a estatística I de Moran. Trata-se de uma representação do coeficiente de regressão e permite visualizar a correlação linear entre z e Wz através do gráfico de duas variáveis. A regressão é especificada de forma que o coeficiente I de Moran é a inclinação da curva de regressão de Wz contra z , sendo que essa inclinação indica o grau de ajustamento.

O diagrama de dispersão de Moran é dividido em quatro quadrantes (AA, BB, AB e BA). Esses quadrantes correspondem a quatro padrões de associação local espacial entre as regiões e seus vizinhos.

As regiões localizadas nos quadrantes AA e BB apresentam autocorrelação espacial positiva, ou seja, essas regiões formam clusters de valores similares. Por outro lado, os quadrantes BA e AB apresentam autocorrelação espacial negativa, ou seja, estas regiões formam clusters com valores diferentes.

O diagrama de dispersão de Moran tem uma limitação devido à ausência da avaliação da incerteza estatística dos clusters AA, BA, BB e AB.

O diagrama de dispersão de Moran9 9 Todos os diagramas de dispersão de Moran estão disponíveis para consulta. (Figura 1) apresenta no eixo horizontal a produtividade média para 1991 (PME91),respectivamente,e, na vertical, a defasagem espacial da variável de interesse para o mesmo período (W_PME91).


Os resultados de autocorrelação espacial (I de Moran), mostrados na seção 3.2, podem ser corroborados pelo fato de que a maioria das regiões está localizada nos quadrantes AA e BB para todos os anos avaliados no presente trabalho. Para o ano de 1991, por exemplo, observa-se que 76% das microrregiões são caracterizadas por valores similares (23% no quadrante AA e 52% no quadrante BB). Desse modo, a autocorrelação espacial positiva manifesta-se mais na forma de clusters com microrregiões de baixa produtividade do que clusters de alta produtividade agrícola.

4. Análise de convergência espacial

A definição de um modelo para captar, por exemplo, questões inerentes a efeitos de transbordamento entre as regiões, devem considerar explicitamente componentes espaciais em sua forma funcional. Portanto, os modelos tradicionais de regressão linear, que não levam em conta os efeitos espaciais (autocorrelação e heterogeneidade espaciais), não devem ser estimados por intermédio do método dos mínimos quadrados ordinários, pois as estimativas serão inconsistentes e/ou ineficientes.10 10 A esse respeito, ver Anselin (1988) e Anselin e Bera (1998).

Assim sendo, ao estimar por mínimos quadrados ordinários (MQO), o modelo de b convergência, busca-se apenas identificar qual é a melhor maneira de se estimar a Equação 5, dada como:

na qual é o logaritmo natural da razão entre a produtividade média entre dois anos em análise (e.g., produtividade média em 2003 e produtividade média em 1991); b ln(Pme t – n)é o logaritmo natural da produtividade média no período inicial e uié o termo de erro.

A fim de identificar a melhor especificação do modelo de b convergência, levando em conta autocorrelação espacial, segue-se a proposta feita por Florax, Folmer e Rey (2003), ou seja, como primeiro estágio, os autores recomendam que seja seguido o roteiro:

a) Estimar o modelo clássico de análise de regressão linear por MQO;

b) Testar a hipótese de ausência de autocorrelação espacial devido a uma defasagem ou a um erro por meio do valor do multiplicador de Lagrange para defasagem espacial (MLr) e Multiplicador de Lagrange para o erro espacial (MLl);

c) Caso ambos os testes não sejam significativos, a utilização do modelo clássico é mais apropriada. Caso contrário, é necessário seguir o próximo passo;

d) Caso ambos sejam significativos, estima-se o modelo apontado como o mais significante de acordo com as versões robustas desses testes, ou seja, o multiplicador de Lagrange robusto para a defasagem espacial (MLRr) e o multiplicador de Lagrange robusto para o erro espacial (MLRl). Caso MLRr > MLRl identifica-se o modelo com defasagem espacial como o mais apropriado. Caso contrário, MLRr < MLRl, adota-se o modelo de erro espacial como o mais apropriado.

4.1. Modelos

4.1.1. Modelo de Erro Espacial

A primeira modificação pode ser o caso em que o termo de erro ui na Equação 5 segue um processo espacial autoregressivo, como mostrado na Equação 6:

na qual l representa o coeficiente escalar do erro espacial, enquanto o termo de erro ei é normalmente distribuído com média zero e variância constante. Substituindo a Equação (6) na Equação (5) obtém-se a forma funcional do modelo de regressão do erro espacial:

Note-se que a matriz W é a mesma matriz de contigüidade utilizada no cálculo da estatística I de Moran. Quando l assumir o valor nulo, não existe autocorrelação espacial do erro. De acordo com Rey e Montouri (1999), quando l ¹ 0, um choque ocorrido em uma unidade geográfica se espalha não só para os seus vizinhos imediatos, mas por todas as outras unidades. Esse tipo de dependência espacial poderia ser resultante de efeitos não-modelados que não fossem aleatoriamente distribuídos através do espaço.

4.1.2. Modelo de Defasagem Espacial

Neste modelo, a autocorrelação espacial é considerada como sendo gerada pela interação atual entre as regiões. Nesse caso, é introduzida uma defasagem espacial como variável independente na equação de convergência b. O modelo é especificado da seguinte forma:

Note-se que r é o coeficiente de defasagem espacial (um escalar). O elemento novo nesta forma funcional pode ser entendido como uma média dos valores da taxa de crescimento das regiões vizinhas. Espera-se que r > 0, sugerindo a existência de autocorrelação espacial positiva.

4.1.3. Modelo Regressivo Cruzado Espacial

Trata-se de um modelo que inclui efeitos de transbordamento espacial. No contexto da convergência, o efeito de transbordamento é representado pela defasagem espacial da produtividade média agrícola na terra do período inicial (1991). Assim, formalmente, o modelo assume a seguinte expressão:

na qual p é o coeficiente de transbordamento espacial, W ln (Pmet – n) denota a defasagem espacial da produtividade no período inicial (e.g. produtividade média em 1995) e ei representa o termo de erro bem comportado.

Nesse modelo, conforme apontado por Rey e Montouri (1999, p. 151), a dependência espacial remanescente toma a forma da média da produtividade média agrícola do começo do período nas microrregiões vizinhas, que seria o termo de transbordamento cruzado.11 11 Considerando ln (P me91 i) como uma variável exógena, assim o é a variável defasada espacialmente W ln (P me91 i). Portanto, é possível estimar esse modelo por intermédio dos mínimos quadrados ordinários (Rey e Montouri, 1999, p. 151). Note que p, nesse caso específico, um escalar.12 12 Se houvesse efeitos de transbordamentos de outras variáveis explicativas, p seria um vetor e não um escalar.

4.2. Resultados e Discussão13 13 É importante observar que o presente trabalho estima apenas a convergência para produtividade agrícola, e não quais são os fatores responsáveis pelo aumento dessa produtividade.

A estimação da Equação (5), apresentada na seção 4, permite obter a estimativa do parâmetro da convergência b, que estuda a hipótese de que microrregiões com produtividade média agrícola baixa tendem a crescer mais rapidamente que microrregiões com produtividade média agrícola maior. Convém relembrar que a Equação 5 faz referência ao modelo de convergência absoluta.

Ao seguir os passos do procedimento de especificação, citado na seção 4, estima-se a Equação 5 por MQO. Os resultados estão reportados na Tabela 3.

O conjunto de testes para averiguar a presença de autocorrelação espacial é útil tanto para auxiliar a identificação do modelo econométrico espacial quanto para a tarefa de validação ou diagnóstico desse modelo. Para tal, fez-se uso dos testes de multiplicador de Lagrange e multiplicador de Lagrange robusto.

É importante salientar dois problemas que podem surgir nos modelos estimados. Primeiro, a existência (da) heterocedasticidade, conforme se verifica aplicando o teste Breusch-Pagan. Segundo, os erros não seguem uma distribuição normal, de acordo com o teste Jarque-Bera. Ambos os problemas sugerem cautela com as inferências.

Para os sub-períodos 1991 a 1994 e 1991 a 2003, o melhor modelo é o do erro espacial. Tal resultado pode ser corroborado pela observação dos resultados quando se aplicam critérios de informação Akaike Information Criterion (AIC), Schwarz Criterion (SC) e do valor da função de verossimilhança (LIK), apresentados na Tabela 3.

Contudo, para o período de 1991-03 foram encontrados ambos os problemas citados acima (não-normalidade dos erros e heterocedasticidade) e para o sub-período 1991-94, apenas os erros não são normais.

Para corrigir a ausência de homocedasticidade, a especificação do modelo de erro espacial foi modificada para acomodar a heterocedasticidade na forma de grupos (groupwise heteroskedasticity).

Quanto ao problema da ausência de normalidade nos erros, adotou-se o método de momentos generalizados (MG), proposto por Kelejian e Prucha (1999), para estimar o modelo de erro espacial, pois prescinde da hipótese da normalidade dos erros.

Já para os sub-períodos 1995-1999 e 2000-2003, o melhor modelo é o de defasagem espacial como podem ser confirmados por intermédio dos critérios de informação Akaike Information Criterion (AIC), Schwarz Criterion (SC) e o valor da função de verossimilhança (LIK). Tais resultados podem ser observados na Tabela 3.

Existem, no entanto, dois problemas recorrentes nos modelos estimados para os dois sub-períodos. Em primeiro lugar, existe heterocedasticidade, e em segundo, os erros não são normais.

Para corrigir a ausência de homocedasticidade, a especificação do modelo de defasagem espacial foi modificada para acomodar a heterocedasticidade na forma de grupos (groupwise heteroskedasticity), como feito para os sub-períodos 1991 a1994 e 1991 a 2003.

Quanto ao problema da ausência de normalidade nos erros, estimou-se o modelo por Variáveis Instrumentais (VI), uma vez que essa estimação prescinde da normalidade dos erros.

Após seguir todos os critérios brevemente explicitados, foi possível estimar os modelos de convergência b, tal como pode ser observado na Tabela 4.

Assim, ao analisar a estimativa do parâmetro b para ambos os modelos, torna-se possível observar a existência ou não de convergência em relação à produtividade média agrícola brasileira. 14 14 A taxa de convergência q foi computada segundo a fórmula: q = ln(b + 1)/(-k). A respeito disso, consulte Rey e Montouri (1999, p. 152).

Primeiramente, ao analisar os sub-períodos 1995-1999 e 2000-2003 observa-se que os parâmetros estimados pelo modelo não são significativos estatisticamente. Portanto, pode-se afirmar que em ambos os sub-períodos, a taxa de convergência é nula.

Em segundo lugar, analisa-se a taxa de convergência para os períodos de 1991 a 1994 e 1991 a 2003, já que, como podem ser observados na tabela 4, os resultados de convergência para tais períodos são significativos estatisticamente e por isso, cabe fazer a análise sobre a velocidade de convergência.

No período entre os anos de 1991 a 2003, a taxa de convergência é de 0,00635, portanto, ocorre a um passo lento. Todavia, é importante observar que o sub-período entre os anos de 1991 a 1994 foi o que mais contribuiu para a convergência da produtividade agrícola no Brasil, com uma taxa de convergência de 0,02193. Ou seja, o período pós-abertura econômica foi o que mais contribuiu para diminuir as desigualdades de produtividade agrícola no Brasil.

5. Considerações finais

O artigo proporciona evidências sobre a evolução da produtividade agrícola brasileira, aplicando-se uma análise do grau de convergência b. No período entre 1991 e 2003, o estudo apresenta evidências em favor da existência de convergência absoluta da produtividade agrícola brasileira através das microrregiões. Isso sugere que a produtividade agrícola brasileira está convergindo para a média, porém a uma taxa muito lenta.

Portanto, cabe indagar quais são as implicações desse resultado para a economia brasileira e, principalmente, em qual proporção tal constatação pode influenciar o setor agrícola no longo prazo.

Ao observar uma convergência lenta na produtividade média agrícola, evidencia-se que a disparidade relacionada à produtividade média agrícola no Brasil, observada no item 2, está diminuindo de forma extremamente lenta durante os anos, colaborando para o aumento das disparidades regionais.

Portanto, a mensuração da convergência b no presente artigo pode ser tomada como inovadora, com relação ao Brasil, sob dois aspectos: a) a incorporação dos efeitos espaciais explícitos no modelo; ou seja, mensura a convergência do setor agrícola, levando em conta a interação das regiões com o seu entorno, e b) em contraste com a grande maioria dos trabalhos, a utilização da produtividade como variável de convergência em detrimento da variável renda per capita.

Os resultados sugerem que pode ser interessante a realização de uma análise de convergência b condicional, incorporando variáveis explicativas como estoques de capital físico, de capital humano e de fatores edafo-climáticos.

6. Referências bibliográficas

  • ANSELIN, L. Interactive techniques and exploratory spatial data analisys. Longley P. A, Goodchild M.F, Maguire D.J e Wind D. W (eds). Geographical information system: principles, techniques, management and applications Wiley, New York. p 253-365. 1998
  • ANSELIN, L. Spatial econometrics: methods and models Kluwer Academic, Boston, 1988.
  • ANSELIN, L. e BERA, A. Spatial dependence in linear regression models with an introduction to spatial econometrics. In: Ullah A. and Giles D. E. (eds.) Handbook of Applied Economic Statistics Marcel Dekker, New York, p. 237-289, 1998.
  • ANSELIN, L The Moran scatterplot as an ESDA tool to assess local instability in spatial association. Fisher, M, Scholten, H.J e Unwin, D W (eds). Spatial analytical perspectives in GIS Taylor&Francis. London. p 111-125. 1996.
  • ANSELIN, L. Local indicators of spatial association LISA. Geographical Analysis, v. 27, n. 2, abril. p. 93-115, 1995.
  • BARRO, R. e SALA-I-MARTIN, X. Economic growth. MacGraw Hill. New York, 1995.
  • BAUMOL, W. J. Productivity growth, convergency, and welfare: What the long-run show. American Economic Review, v. 76, n. 5, p. 1072-85, 1986.
  • CLIFF, A. D. e ORD, J.K. Spatial processes: models and applications Pion, London. 1981.
  • CONSIDERA, C. M., SOUSA, E. L. L. e BRACALE, G. Âncora verde: o papel da agricultura no ajuste econômico. SEAE, Brasília, Distrito Federal, 2002.
  • DALLERBA, S. Productivity convergence and spatial dependence among Spanish regions. Texto para discussão REAL/UIUC, Illinois, 2003.
  • DEPARTAMENTO INTERSINDICAL DE ESTATÍTICAS E ESTUDOS SÓCIO-ECONÔMICOS DIEESE. Dez anos do Real São Paulo, 2004.
  • FIGUEIREDO, A. M. R. Resposta da produção agrícola aos preços no Centro-Oeste brasileiro: uma análise de econometria espacial para o período 1975-1995/1996 Tese de Doutorado, Universidade Federal de Viçosa, 2002. (mimeo).
  • FLORAX, R. J. G. M., FOLMER, H., REY, S. J. Specification searches in spatial econometrics: The relevance of Hendrys methodology. Regional Science and Urban Economics, vol. 33, n. 5, p. 557-79, 2003.
  • JOHNSTON, B.F.; MELLOR, J.W. The role of agriculture in economic development. American Economic Review, vol. 51, n.4, p. 566-93, 1961.
  • KELEJIAN, H. H. e PRUCHA, I. R. (1999). A generalized moments estimator for the autoregressive parameter in a spatial model. International Economic Review, vol. 40, n. 2.
  • LE GALLO, J and ERTHUR, C. Exploratory spatial data analysis of the distribution of regional per capita. GDP in Europe, 1980-1995. Papers in Regional Science, vol. 82, n. 2, p. 175-201. 2003.
  • LOPES, J. L. Avaliação do processo de convergência da produtividade da terra na agricultura brasileira no período de 1960 a 2001 Tese de Doutorado, ESALQ/USP, 2004. (mimeo).
  • MAGALHÃES, A. HEWINGS, G. J. D. AZZONI, C. R. Spatial dependence and regional convergence in Brazil. Texto para discussão, REAL/UIUC, Illinois, 2000.
  • Ministério da Agricultura, Pecuária e Abastecimento (MAPA). Brasília, 2005. Disponível em: < http://www.agricultura.gov.br/> Acesso em: 10 fev 2005.
    » link
  • QUAH, D. T. Regional convergence clusters across Europe. European Economic Review, vol. 40, p. 951-958, 1996.
  • REY, J. S., e MONTOURI, B. D. US Regional income convergence: a spatial sconometric perspective. Regional Studies, vol. 33, n. 2, p. 143-156, 1999.
  • Sistema IBGE de Recuperação Automática (SIDRA). Brasília, 2005. Disponível em: < http://www.sidra.ibge.gov.br/bda/ > Acesso em: 5 fevereiro de 2005.
    » link
  • 1
    Embora a maioria dos trabalhos que analisou a questão da convergência de renda encontrou evidências de convergência com nivelamento para cima, é possível, teórica e empiricamente, admitir a existência de convergência com nivelamento para baixo. Todavia, como ficará claro adiante, não é razoável admitir convergência com nivelamento para baixo da produtividade agrícola no Brasil no período analisado.
  • 2
    Consultar Barro e Sala-i-Martin (1995) para uma discussão avançada sobre as convergências absoluta e condicional.
  • 3
    De acordo com Anselin (1999, p. 2),
    " econometria espacial é um sub-campo da econometria que lida com as complicações causadas pela interação espacial (autocorrelação espacial) e pela estrutura espacial (heterogeneidade espacial), em modelos de regressão para dados na forma de seção cruzada e painel de dados" .
  • 4
    Qualquer matriz de pesos espaciais precisa atender às condições de regularidade impostas pela necessidade de invocar as propriedades assintóticas dos estimadores e dos testes. Segundo Anselin (1988), "
    isso significa que os pesos precisam ser não-negativos e finitos e que correspondam a uma determinada métrica" .
  • 5
    De acordo com Le Gallo e Ertur (2003), a escolha de um número fixo de vizinhos mais próximos ao invés do uso de uma matriz simples de contigüidade é preferível, pois é possível evitar alguns problemas metodológicos que podem ocorrer quando há variações no número de vizinhos.
  • 6
    Neste estudo, foi utilizado valor de
    k=15 (ie., quinze vizinhos mais próximos), todavia, a fim de verificar a robustez dos resultados foi testado k=10 e k=20. Vale destacar que os resultados k=10 e k=20 vizinhos mais próximos estão disponíveis para consulta.
  • 7
    O operador defasagem espacial é oriundo da idéia de defasagem temporal da econometria de séries de tempo. Contudo, o conceito de operador de defasagem espacial não é tão simples e direto como o operador de defasagem em séries de tempo devido à natureza multidirecional do processo de interação no espaço.
  • 8
    Note-se que a matriz de pesos espaciais
    W tem dimensão
    n por
    n, sendo
    n o número de unidades espaciais (por ex., microrregiões), ao passo que
    y é um vetor-coluna com os valores da variável de interesse (por exemplo, produtividade agrícola média).
  • 9
    Todos os diagramas de dispersão de Moran estão disponíveis para consulta.
  • 10
    A esse respeito, ver Anselin (1988) e Anselin e Bera (1998).
  • 11
    Considerando ln (P
    me91
    i) como uma variável exógena, assim o é a variável defasada espacialmente
    W ln (P
    me91
    i). Portanto, é possível estimar esse modelo por intermédio dos mínimos quadrados ordinários (Rey e Montouri, 1999, p. 151).
  • 12
    Se houvesse efeitos de transbordamentos de outras variáveis explicativas, p seria um vetor e não um escalar.
  • 13
    É importante observar que o presente trabalho estima apenas a convergência para produtividade agrícola, e não quais são os fatores responsáveis pelo aumento dessa produtividade.
  • 14
    A taxa de convergência q foi computada segundo a fórmula: q = ln(b + 1)/(-k). A respeito disso, consulte Rey e Montouri (1999, p. 152).
  • Datas de Publicação

    • Publicação nesta coleção
      10 Jun 2008
    • Data do Fascículo
      Mar 2008
    Sociedade Brasileira de Economia e Sociologia Rural Av. W/3 Norte, Quadra 702 Ed. Brasília Rádio Center Salas 1049-1050, 70719 900 Brasília DF Brasil, - Brasília - DF - Brazil
    E-mail: sober@sober.org.br