Como a pecuária contribui para a expansão do território brasileiro?

Este trabalho analisa a distribuição espacial do rebanho bovino nas microrregiões brasileiras entre 1995 e 2016, com enfoque especial na Amazônia Legal. A partir de uma análise espacial com dados em painel, a taxa de crescimento do rebanho se mostrou entre 5,6% a 15,68% mais elevada na Amazônia Legal, enquanto o crescimento da densidade bovina (rebanho por km²) foi superior em 11,11% a 21,47%, indicando uma possível intensificação da produção na região. Além disso, verificou-se uma desconcentração regional da pecuária brasileira, em especial devido à presença de externalidades espaciais positivas que auxiliaram no espraiamento do rebanho. O acesso ao crédito pecuário também se mostrou importante para explicar esse crescimento ao gerar forças de atração e propulsão para a produção. De maneira específica, a produtividade dos insumos utilizados na nutrição bovina (soja, milho e cana-de-açúcar) auxiliou no crescimento da densidade bovina, mas não do rebanho, corroborando sua relação com a intensificação da produção.

Palavras-chave:
pecuária; Amazônia Legal; análise espacial com dados em painel

Abstract

This paper analyses the spatial distribution of the bovine herd in the Brazilian micro-regions between 1995 and 2016, especially in the Legal Amazon. From a spatial analysis with panel data, the growth rate showed to be between 5.6% to 15.68% higher in the Legal Amazon, while the increase in bovine density (herd per km²) was higher in 11.11% to 21.47%, indicating a possible intensification of production in the region. In addition, there was a regional deconcentration of Brazilian livestock, especially due to positive spatial externalities that helped the spread of the herd. Credit for livestock also proved to be important to explain this growth by generating forces of attraction and propulsion for production. Specifically, the productivity of the inputs used in bovine nutrition (soybeans, corn and sugar cane) helped the growth of bovine density, but not in the herd, which corroborates its relationship with the intensification of production.

Keywords:
Cattle; Legal Amazon; Spatial Panel Data

1 Introdução

As atividades relativas à pecuária bovina no Brasil possuem grande importância, dado que o país tem o maior rebanho comercial do mundo, sendo o segundo maior produtor e o maior exportador mundial de carne bovina (Carvalho; Zen, 2017CARVALHO, T. B.; ZEN, Sérgio. A cadeia de pecuária de corte no Brasil: evolução e tendências. Revista iPecege, v. 3, n. 1, p. 85-99, 2017.). O país figura, no cenário internacional, como um dos principais atores na produção e comércio de carne bovina, reflexo de um estruturado processo de desenvolvimento que elevou não só a produtividade como também a qualidade do produto brasileiro e, consequentemente, sua competitividade e abrangência de mercado. Nas últimas quatro décadas, a pecuária brasileira sofreu uma modernização sustentada por avanços no nível tecnológico dos sistemas de produção e na organização da cadeia produtiva, com reflexo na qualidade da carne bovina (EMBRAPA, 2017). Segundo Mcmanaus et al. (2016), esse rápido crescimento e modernização ocorreu principalmente na região Norte e Centro-Oeste do Brasil, especialmente a partir do início desse século.

No que se refere ao avanço da produção pecuária no Brasil, é visível uma divisão na literatura especializada. Alguns autores, como Cezar et al. (2005CEZAR, I. M.; QUEIROZ, H. P.; THIAGO, L. R. L. S. Sistemas de produção de gado de corte no Brasil: uma descrição com ênfase no regime alimentar e no abate. Campo Grande, MS: Embrapa Gado de Corte, 2005.), Carvalho (2007CARVALHO, T. B. Estudo da elasticidade-renda da demanda de carne bovina, suína e de frango no Brasil. Tese (Doutorado) - Universidade de São Paulo, 2007.), McManaus et al. (2016) e Silva et al. (2018SILVA, M. M.; OLIVEIRA, F. A.; SANTANA, A. C. Mudanças socioambientais no uso da terra em Altamira, Amazônia Oriental. Novos Cadernos NAEA, v. 20, n. 3, 2018.) afirmam que o aumento na produção bovina no Brasil ocorreu principalmente devido à utilização do sistema de produção com base no regime de alimentação em pastagem. Desse modo, a expansão para áreas antes não ocupadas é uma importante ferramenta para a evolução da produção brasileira, sustentando a hipótese de um avanço da produção segundo um regime alimentar no formato extensivo.1 1 Regime alimentar exclusivo de pastagem. Por outro lado, autores como Fearnside (2005FEARNSIDE, P. M. Desmatamento na Amazônia brasileira: história, índices e consequências. Megadiversidade, v. 1, n. 1, p. 113-123, 2005.), Riveiro et al. (2009), e Assad (2016ASSAD, E. D. Amazônia Legal: propostas para uma exploração agrícola sustentável (sumário executivo). Relatório Técnico - FGVAgro, 2016), afirmam que a bovinocultura avançou especialmente de forma intensiva2 2 Regime alimentar de pastagem com suplementação e confinamento. a partir de um processo de desenvolvimento, integração e intensificação de sua produção, o que minimizou a abertura de novas áreas de pastagem. Em termos de rebanho, o efetivo mais que dobrou nas últimas quatro décadas, enquanto a área de pastagens não avançou na mesma proporção, indicando aumento da produtividade no período (Carvalho; Zen 2017).

Segundo Fearnside (2005FEARNSIDE, P. M. Desmatamento na Amazônia brasileira: história, índices e consequências. Megadiversidade, v. 1, n. 1, p. 113-123, 2005.), Riveiro et al. (2011) e Assad (2016ASSAD, E. D. Amazônia Legal: propostas para uma exploração agrícola sustentável (sumário executivo). Relatório Técnico - FGVAgro, 2016), o crescimento apresentado pelo rebanho bovino brasileiro, a partir da década de 1990, se deu principalmente na região da Amazônia Legal.3 3 Formam a Amazônia legal os estados da região Norte (Acre, Amapá, Amazonas, Pará, Rondônia, Roraima e Tocantins) além de Mato Grosso e Maranhão. Nessa região, em especial, o rebanho aumentou, em grande medida, em virtude do aumento da área de pastos, aliado a ganhos de produtividade da pecuária fruto do sucesso do controle da febre aftosa, além de subsídios públicos. Entretanto, o crescimento da produção, em especial daquela advinda do aumento das pastagens, também gera externalidades ambientais negativas, principalmente devido ao desmatamento para abertura de novas áreas de produção.

Nesse contexto, este artigo tem como finalidade verificar a evolução e o desempenho da produção bovina nas microrregiões brasileiras, assim como os determinantes do crescimento apresentado no período de 1996 a 2016, com especial enfoque para a Amazônia Legal. Ademais, o trabalho buscará verificar se o rebanho bovino, assim como sua densidade, apresentou taxas de crescimento superiores em regiões, como a Amazônia, que inicialmente não possuíam grandes níveis de produção. Dito de outra forma, há um processo de desconcentração da produção pecuária brasileira? Basicamente, será considerado o nível da produção inicial das microrregiões em relação à sua posterior taxa de crescimento no período de análise, além de diversos fatores que possam ter induzido ao crescimento da produção pecuária, como a produtividade dos insumos básicos, área destinada ao plantio e montante de crédito rural para pecuária. Por fim, o número de cabeças de gado criadas nas microrregiões será utilizado como proxy para a produção pecuária brasileira. Isso possibilitará, portanto, a verificação de vinte anos de evolução espacial e temporal na produção bovina no Brasil.

Para se atingir os objetivos propostos, serão utilizadas as técnicas da Análise Exploratória de Dados Espaciais (AEDE), que permitem analisar a distribuição espacial, além de possíveis processos espaciais adjacentes, a saber, a existência de dependência e da heterogeneidade espacial. O uso dessas técnicas de análise espacial se deve ao fato de a agropecuária ser frequentemente sujeita aos efeitos espaciais citados, pois a existência de diferentes técnicas de produção, clima, topografia e condições de solo podem causar diferenças regionais significativas.

Ademais, muitos estudos que visam analisar fatores econômicos e agropecuários, utilizam-se de dados regionais, sendo, portanto, espaciais. Entretanto, é comum a não consideração de possíveis efeitos espaciais contidos na amostra, o que pode invalidar as inferências estatísticas. Por esses motivos, além da AEDE e da estimação econométrica convencional de dados em painel, caso seja identificada a presença de dependência espacial, essa será devidamente considerada, adotando-se técnicas advindas da econometria espacial. Os procedimentos advindos das AEDE e da econometria espacial, portanto, possibilitarão o desenvolvimento de estimativas mais confiáveis, permitindo uma nova perspectiva da dinâmica espacial do crescimento do rebanho bovino brasileiro.

Sendo assim, o artigo se encontra estruturado em cinco seções, incluindo esta introdução. Na segunda será exposto o referencial teórico sobre a taxa de crescimento da produção bovina, especialmente sob a perspectiva brasileira, além da realização preliminar de um breve panorama sobre a produção no país. A terceira seção destina-se à descrição da metodologia, assim como da base de dados utilizada no trabalho. Os resultados encontrados e sua análise são realizados na quarta seção. Finalmente, na quinta seção têm-se as considerações finais.

A produção bovina brasileira, com a incorporação de novas tecnologias que visam ao aumento da produtividade, de acordo com Carvalho (2007CARVALHO, T. B. Estudo da elasticidade-renda da demanda de carne bovina, suína e de frango no Brasil. Tese (Doutorado) - Universidade de São Paulo, 2007.), cresceu especialmente devido à adoção de sistemas intensivos de produção. Na década de 1970, com auxílio do crédito rural, novas tecnologias foram incorporadas, como a importação de animais de raças europeias e indianas, da introdução de sal na alimentação; além da atenção especial dada ao controle de doenças como a febre aftosa. Desse modo, a partir de 1990, observou-se uma pecuária nacional com maiores índices de produtividade conjuntamente a uma indústria integrada nacionalmente.

O entendimento do sistema de produção bovino é importante para compreender a dinâmica de seu crescimento. Sendo assim, segundo a Cezar et al. (2005CEZAR, I. M.; QUEIROZ, H. P.; THIAGO, L. R. L. S. Sistemas de produção de gado de corte no Brasil: uma descrição com ênfase no regime alimentar e no abate. Campo Grande, MS: Embrapa Gado de Corte, 2005.) e Kamali et al. (2016KAMALI, F. P.; LINDEN, A. V. D.; MEUWISSEN, M. C. Environmental and economic performance of beef farming systems with different feeding strategies in southern Brazil. Agricultural Systems, v. 146, p. 70-79, 2016.), a pecuária de corte brasileira apresenta uma ampla gama de sistemas de produção e é desenvolvida em todos os estados e ecossistemas do país. Além disso, há uma heterogeneidade nos sistemas produtivos, nos mecanismos de gestão e de comercialização do gado, sendo características importantes no desenvolvimento dessa atividade no Brasil.

No que se refere aos regimes alimentares adotados pelos pecuaristas brasileiros, coexistem basicamente três. O primeiro é o de alta qualidade, caracterizado pela adoção de tecnologia avançada e padrões eficientes de gestão e de comercialização com o uso mais intensivo de pastagens cultivadas aliado à prática de confinamento na terminação. O segundo, de baixa qualidade, baseia-se na produção extensiva, pequena intensidade tecnológica e padrões precários de gestão e de comercialização do gado bovino, que são caracterizados pela utilização de pastagens nativas e cultivadas como únicas fontes de alimentos energéticos e proteicos. Por fim, o terceiro e menos utilizado na produção pecuária, é o sistema semi-intensivo, que também apresenta como base alimentar as pastagens (nativas e cultivadas) e os suplementos minerais, acrescidos de suplementos proteicos/energéticos (Cezar et al., 2005CEZAR, I. M.; QUEIROZ, H. P.; THIAGO, L. R. L. S. Sistemas de produção de gado de corte no Brasil: uma descrição com ênfase no regime alimentar e no abate. Campo Grande, MS: Embrapa Gado de Corte, 2005.)

A partir do regime alimentar adotado, é possível observar o aumento no consumo de subprodutos agroindustriais, segundo Filho (2001FILHO, R. A. M. Sistemas de terminação em confinamento: perspectivas para dinamização da cadeia produtiva da carne bovina em Minas Gerais. Revista Brasileira de Economia, v. 55, n. 1, 2001.), visando usufruir da sazonalidade do boi gordo. Por isso, um número considerável de pecuaristas passara a manter seus animais, no último período de engorda, em sistemas intensivos, principalmente no confinamento. Tal procedimento possibilita o aumento no giro do capital, dadas a redução no tempo para a obtenção do boi gordo e a possibilidade de utilização de subprodutos agroindustriais na alimentação desses animais, como o bagaço hidrolisado da cana, grãos de soja e milho.

De acordo com Martinelli e Filoso (2008MARTINELLI, L. A.; FILOSO, S. Expansion of sugarcane ethanol production in Brazil: environmental and social challenges. Ecological Applications, v. 18, n. 4, p. 885-898, 2008.), Mendonça et al. (2016MENDONÇA, F. S.; VAZ, R. Z.; COSTA, A. D. Fatores que afetam o bem-estar de bovinos durante o período pré-abate. Archivos de Zootecnia, v. 65, n. 250, p. 279-287, 2016.) e Carvalho e Zen (2017CARVALHO, T. B.; ZEN, Sérgio. A cadeia de pecuária de corte no Brasil: evolução e tendências. Revista iPecege, v. 3, n. 1, p. 85-99, 2017.), subprodutos agroindustriais são amplamente utilizados na nutrição bovina, especialmente no regime alimentar de alta qualidade. Por isso, ao avaliar a produção pecuária, os autores afirmam que é necessário analisar e considerar o comportamento da produção de tais culturas, assim como suas respectivas produtividades, pois estas podem impactar diretamente a evolução da densidade bovina regional. Dito de outra forma, a alta produtividade dessas culturas em regiões próximas à produção pecuária pode induzir indiretamente o aumento do rebanho bovino. Conforme destaca Cezar et al. (2005CEZAR, I. M.; QUEIROZ, H. P.; THIAGO, L. R. L. S. Sistemas de produção de gado de corte no Brasil: uma descrição com ênfase no regime alimentar e no abate. Campo Grande, MS: Embrapa Gado de Corte, 2005.), a produtividade da pecuária está associada basicamente ao sistema de produção intensivo, tornando, portanto, o regime alimentar um indicativo de seu aumento. Destaca-se, por exemplo, o crescimento do consumo de suplementos alimentares com o aumento da utilização subprodutos dos grãos de soja, milho e da cana de açúcar utilizada como volumoso.4 4 Suplementação para gado de corte em confinamento. Alimentos volumosos são basicamente aqueles que têm alto teor de fibra e baixo valor energético, geralmente utilizados para produção de silagem.

Outro fator importante para o avanço da produção bovina no Brasil é a disponibilidade de crédito. Avanços na política de expansão do crédito para a pecuária são fundamentais no processo de modernização da bovinocultura do país, sendo eles oriundos de financiamentos do PRONAF (Programa Nacional de Fortalecimento da Agricultura Familiar). Entretanto, o credito pecuário é distribuído de forma desigual ao longo do território brasileiro, além de algumas regiões receberem valores mais elevados para custeio5 5 O PRONAF Custeio: aquisição de sementes, fertilizantes, defensivos, vacinas, ração e outros itens necessários para o dia a dia da sua produção, seja ela agrícola ou pecuária (PRONAF, 2018). enquanto em outras há predominância de investimentos6 6 O PRONAF Investimento: financiar despesas de produção, investir em benfeitorias, tratores, máquinas e implementos agrícolas, além de comercializar a produção. As soluções de crédito rural contribuem para a diminuição de custos, aumento da produtividade e melhoria no fluxo de caixa, e capital de giro para atividade agropecuária (PRONAF, 2018). (Macedo, 2005MACEDO, L. O. B. Modernização da pecuária de corte bovina no Brasil e a importância do crédito rural. Agroanalysis, Rio de Janeiro, v. 25, n. 6, p. 35-36, 2005.).

Porém, vale ressaltar que o crescimento da agropecuária e especialmente da pecuária bovina no Brasil veio acompanhado do desmatamento dos três principais biomas brasileiros: a Mata Atlântica, o Cerrado e a Amazônia. A Mata Atlântica sofreu mais degradação devido à agricultura, à expansão da urbanização e, em menor grau, da pecuária, sendo que atualmente a maior parte desse bioma permanece em fragmentos, localizados no litoral, Serra do Mar e nas encostas da Serra Geral. O bioma Cerrado foi substituído principalmente por pasto e/ou outros cultivos diversos, como a soja. E, por fim, na Amazônia, a ocupação ocorreu principalmente pela incorporação de novas áreas de produção agropecuária, especialmente da pecuária, nas áreas da fronteira agrícola da região. Além disso, nas últimas décadas, esse bioma tem sofrido grande pressão para o aumento da produção pecuária (Alencar et al., 2015ALENCAR, A.; PEREIRA, C.; CASTRO, I.; CARDOSO, A. Desmatamento nos Assentamentos da Amazônia: histórico, tendências e oportunidades. Brasília, DF: IPAM, 2015.; Domingues, 2017DOMINGUES, M. S.; BERMANN, C.; SIDNEIDE, S.; A produção da soja no Brasil e sua relação com o desmatamento na Amazônia. Revista Presença Geográfica, v. 1, n. 1, 2017.).

Carvalho (2007CARVALHO, T. B. Estudo da elasticidade-renda da demanda de carne bovina, suína e de frango no Brasil. Tese (Doutorado) - Universidade de São Paulo, 2007.) e Martinelli et al. (2010MARTINELLI, L. A.; NAYLOR, R.; VITOUSEK, P. M.; MOUTINHO, P. Agriculture in Brazil: impacts, costs, and opportunities for a sustainable future. Current Opinion in Environmental Sustainability, v. 2, n. 5-6, p. 431-438, 2010.) também enfatizam que a pecuária bovina de corte se desenvolveu, nas últimas décadas, principalmente devido à expansão da fronteira agrícola na Amazônia, com a incorporação de novas terras. Segundo Barona et al. (2010BARONA, E; RAMANKUTTY, N; HYMAN, G; COMMES, O. The role of pasture and soybean in deforestation of the Brazilian. Environmental Research Letters, 5, 024002 (9pp), 2010.), o avanço do desmatamento na Amazônia pode estar ligado especialmente à produção da pecuária e, em menor grau, da soja; devido à pressão que essas exercem nas áreas de florestas próximas à fronteira agrícola. De fato, segundo Silva et al. (2018SILVA, M. M.; OLIVEIRA, F. A.; SANTANA, A. C. Mudanças socioambientais no uso da terra em Altamira, Amazônia Oriental. Novos Cadernos NAEA, v. 20, n. 3, 2018.), a expansão da pecuária na região amazônica tem como consequência direta uma maior pressão sobre o meio ambiente, refletida na diminuição das áreas de florestas primárias. Além disso, o crescimento da produção pecuária na região gerou diversas mudanças socioambientais no uso da terra, agravadas pela baixa capacidade de suporte das pastagens na região, assim como o baixo nível tecnológico dos sistemas de produção desenvolvidos. Por fim, segundo Brondizio (2009BRONDIZIO, E. S. Análise inter-regional de mudança do uso da terra na Amazônia. In: MORAN, E. F.; OSTROM, E. (Org.). Ecossistemas Florestais: interação homem-ambiente. São Paulo: EDUSP, p. 289-326, 2009.), o aumento do desmatamento na região se deve, em grande medida, aos incentivos dados pelo crédito rural destinado aos produtores.

O avanço do rebanho bovino na Amazônia Legal pode ser evidenciado, por exemplo, pelo número de cabeças de gado na região. Em 1995, por exemplo, havia 37 milhões de cabeças na região, valor equivalente a 23% do total nacional. Em 2016, por outro lado, esse valor saltou para 85 milhões - cerca de 40% da produção brasileira. Nascimento (2017NASCIMENTO, F. Desmatamento zero na Amazônia: como e por que chegar lá. Greenpeace, 2017.) destaca que o aumento da pastagem para a criação de gado é a atividade que mais contribui para o desmatamento na Amazônia, compreendendo 65% dessa área desmatada. Uma série de estudos vem apontando para os impactos negativos de práticas agropecuárias nos desmatamentos na Amazônia, em especial, no que se refere às atividades ligadas à criação de bovinos, que possivelmente geram maiores impactos ambientais do que outras atividades (Brondizio, 2009BRONDIZIO, E. S. Análise inter-regional de mudança do uso da terra na Amazônia. In: MORAN, E. F.; OSTROM, E. (Org.). Ecossistemas Florestais: interação homem-ambiente. São Paulo: EDUSP, p. 289-326, 2009.; Alencar et al. 2015ALENCAR, A.; PEREIRA, C.; CASTRO, I.; CARDOSO, A. Desmatamento nos Assentamentos da Amazônia: histórico, tendências e oportunidades. Brasília, DF: IPAM, 2015.).

Nesse sentido, o presente artigo buscar contribuir com a discussão sobre a evolução do rebanho bovino brasileiro, com especial enfoque para a Amazônia Legal. Isso será realizado com a estimação de um modelo econométrico visando verificar a consistência do avanço da produção bovina no Brasil. Para tanto, na seção seguinte será apresentada uma discussão sobre a metodologia adotada e a base de dados utilizada.

3.1 Base de dados e mecanismos metodológicos

A variável básica considerada é o número de cabeças de bovinos nas 588 microrregiões brasileiras no período entre 1995 e 2016, tendo como base de dados a Pesquisa Pecuária Municipal (PPM), extraída do Sistema de Recuperação Automática (SIDRA) do Instituto Brasileiro de Geografia e Estatística (IBGE). A utilização do número de cabeças de gado será utilizada para captar o desenvolvimento espaçotemporal da produção pecuária brasileira ao longo de suas regiões do país. Também serão empregadas variáveis referentes à produtividade do milho, soja e cana-de-açúcar.7 7 Rendimento médio da produção (quilogramas por hectare). Segundo Martinelli e Filoso (2008MARTINELLI, L. A.; FILOSO, S. Expansion of sugarcane ethanol production in Brazil: environmental and social challenges. Ecological Applications, v. 18, n. 4, p. 885-898, 2008.), Mendonça et al. (2016MENDONÇA, F. S.; VAZ, R. Z.; COSTA, A. D. Fatores que afetam o bem-estar de bovinos durante o período pré-abate. Archivos de Zootecnia, v. 65, n. 250, p. 279-287, 2016.), Carvalho e Zen (2017CARVALHO, T. B.; ZEN, Sérgio. A cadeia de pecuária de corte no Brasil: evolução e tendências. Revista iPecege, v. 3, n. 1, p. 85-99, 2017.), para caracterizar o sistema intensivo de produção, em caso da inexistência de dados específicos, é coerente utilizar como proxy a produtividade dos subprodutos agrícolas, insumos básicos para a produção em sistemas intensivos de produção. Ao reduzir o custo dessas matérias-primas, pode-se induzir a adoção desses insumos pela pecuária, intensificando a produção. Portanto, o presente trabalho busca contribuir com literatura no que se refere à discussão da produtividade das culturas utilizadas na alimentação bovina, principalmente na utilização desses insumos no sistema intensivo de engorda e seus efeitos na produção bovina.

Entretanto, a produtividade dos insumos tende a apresentar valores semelhantes para a mesma microrregião devido à semelhança de características físicas, naturais e humanas dessas localidades, resultando em colinearidade entre as variáveis. Para minimizar o problema, será construído um indicador de produtividade a partir da estatística multivariada, especialmente da análise fatorial. Essa técnica possibilitará evitar o viés de má especificação do modelo, ao mesmo tempo que trata o problema de multicolinearidade, pois o método permite sintetizar as variáveis em fatores ortogonais e não correlacionados entre si, melhorando o ajuste do modelo. A extração dos fatores será realizada por meio do método de componentes principais, relacionando linearmente as Zi variáveis e os k fatores comuns extraídos, conforme:

Z1=l11F1+…+l1kFk+ε1⋮⋮⋮Zp=lp1Fp+…+lpkFk+εp (1)

ou matricialmente por:

D(X−μ)=LF+ε (2)

onde F é um vetor aleatório contendok fatores, que resume as p variáveis incluídas; ε é um vetor de erros aleatório, que contém a parcela de Zi não explicada pelos Fj fatores, com j=1,2,…,k∈ℕ; L é uma matriz de parâmetros lij(loadings) a serem estimados, os quais representarão o grau de relação linear entre Zi e Fj. De posse do fator produtividade de insumos,8 8 A aplicação da análise fatorial para criar o fator de produtividade de insumos possibilitou sintetizar as três variáveis em apenas um fator com um valor de 0,5795 para o teste Kaiser-Meyer-Olkin (KMO), 1,71 de carga fatorial e 57,06% de variância acumulada. será testada a hipótese de que regiões com alta produtividade das culturas consideradas impactam positivamente no crescimento do rebanho bovino.

Outra variável importante, que visa captar mudanças estruturais das microrregiões, é a área total da microrregião destinada ao plantio. As variáveis de produtividade e área plantada provêm da Pesquisa Pecuária Municipal (PPM), extraídas do Sistema de Recuperação Automática (SIDRA) do IBGE. Também serão utilizados dados para o montante de crédito pecuário destinado à microrregião, variável extraída do Instituto de Pesquisa Econômica Aplicada (IPEA) e Banco Central do Brasil. Por fim, para investigar se a região Amazônica teve uma taxa de crescimento diferenciada em relação às demais localidades do país, uma variável dicotômica para a Amazônia Legal será incluída no modelo utilizado, a qual caberá maiores explicações nos próximos parágrafos. Em síntese, as variáveis utilizadas no presente artigo estão descritas no Quadro 1.

Thumbnail

Quadro 1:
Variáveis utilizadas

De posse desse banco de dados, efetuou-se uma média de dois anos para representar cada unidade de cross section que irá compor o painel de dados.9 9 Por exemplo, o primeiro período será composto pela média dos anos de 1995 e 1996. O procedimento é comumente adotado com a finalidade de se minimizar efeitos aleatórios, como o de variações de preços, climáticas, entre outros, que podem afetar o crescimento do rebanho. Entretanto, essas flutuações, muitas vezes, não representam a dinâmica da pecuária brasileira. Para verificar o comportamento do crescimento do rebanho e a densidade bovina na região Amazônia, será adotado um procedimento metodológico proposto por Wooldridge (2002WOOLDRIDGE, J. M. Econometric analysis of cross section and panel data. Cambridge: The MIT Press, 2002.), que se caracteriza como uma forma de se incluir variáveis dicotômicas no contexto de dados em painel. Primeiramente, criou-se uma variável temporal (t) para os quatro períodos analisados e uma dummy que assume 1 para a microrregião pertencente à Amazônia Legal e 0 para as demais. Por fim, estabeleceu-se uma dummy de interação entre elas, resultando, portanto, numa variável representativa da Amazônia Legal para cada período. De forma geral, o resultado pode ser formalizado como

yit=α+θit+β1Xit+βtAmazôniaLegal+ϵit (3)

em que X é o vetor de variáveis explicativas para a microrregião i no período t; θ é a heterogeneidade não observada. O procedimento busca isolar o efeito temporal das microrregiões pertencentes à Amazônia Legal com a finalidade de identificar se a taxa de crescimento do rebanho e da densidade bovino apresenta alguma dinâmica diferenciada nessa região em relação às demais localidades. Portanto, o modelo básico específico a ser estimado pelo presente artigo pode ser representado como

ln(Rebanho.Bovi,tRebanho.Bovi,t−1)=α+θit+β1ln(Rebanho.bovi,t−1)+δln(CEMi,t−1)+βtAmazôniaLegal+εi,t (4)

e

ln(Densidade.Bovi,tDensidade.Bovi,t−1)=α+θit+β1ln(Densidade.bovi,t−1)+δln(CEMi,t−1)+βtAmazôniaLegal+εi,t (5)

Onde ln(Rebanho.Bovi,tRebanho.Bovi,t−1)eln(Densidade.Bovi,tDensidade.Bovi,t−1) é o logaritmo natural da razão entre o rebanho e a densidade bovina, respectivamente, no período t pelo rebanho em t−1 e pode ser interpretada como a taxa de crescimento da produção e densidade pecuária na microrregião i; βln(Rebanho.bovi,t−1)eβln(Densidade.bovi,t−1) é o logaritmo natural do rebanho e densidade bovina no período inicial t−1; εi é o termo de erro. Sendo que a hipótese de desconcentração do rebanho no território nacional é determinada pelo sinal do parâmetro β1: se esse for negativo e significativo, o processo se confirma. Por fim, CEMi,t−1 se refere a um vetor de variáveis de controle, que busca representar as características estruturais da microrregião i no período t−1, sendo composto pelas seguintes variáveis: área total destinada ao plantio; crédito pecuário total e o indicador de produtividade dos insumos. A defasagem temporal das variáveis de controle busca evitar problemas de endogeneidade na estimação. As variáveis para a Amazônia Legal10 10 Para evitar a armadilha da variável dummy, excluiu-se a variável para o primeiro período, 1995. se referem às dummys para o período entre 1995 e 2016.

3.2 Análise Exploratória de Dados Espaciais (AEDE) e Análise Espacial com Dados em Painel

As técnicas advindas da Análise Exploratória de Dados Espaciais utilizadas neste trabalho para identificar efeitos espaciais, especificamente os de dependência e heterogeneidade espacial, foram o I de Moran global e o local. O primeiro capta o grau de correlação espacial de uma variável através das regiões,

It=(nS0)(zt'Wztzt'zt) (6)

onde W pode ser definida como uma matriz quadrada em que cada célula wij indica relação existente entre a região iej em um sistema de n regiões; o termo S0 é um escalar igual à soma de todos os elementos de W. O indicador acima é equivalente ao grau de associação linear entre o vetor de valores observados, z, e a média ponderada de valores vizinhos, Wz. A versão local do I de Moran, por sua vez, é

Ii=zi∑j−1Jwijzj (7)

onde zi representa a variável de interesse da regiãoi padronizada, wij é o elemento da matriz de ponderação espacial (W) e zj é o valor da variável de interesse na região j padronizada, podendo representar quatro tipos de clusters espaciais: Alto-Alto (AA), Baixo-Baixo (BB), Alto-Baixo (AB) e Baixo-Alto (BA). O cluster Alto-Alto, por exemplo, captura regiões com alto valor que estejam circundadas por outras com valores similares.

O método utilizado para investigar os determinantes do crescimento da produção bovina brasileira será o de painel de dados espacial. Tal escolha se deve ao fato de a produção bovina não estar, ao longo do período, relacionada simplesmente com os fatores explicitamente incluídos na regressão, mas também com características próprias de cada microrregião. O método de dados em painel possibilita captar esses efeitos fixos não observáveis que iriam, caso contrário, compor o termo de erro, enviesando as estimativas quando correlacionadas com as variáveis explicativas. A utilização da regressão com dados em painel tem como objetivo minimizar o viés e tornar o modelo consistente, ou seja, simétrico e com menor variância possível.

Entretanto, um dos pressupostos do modelo de dados em painel é que as unidades de corte transversal são independentes entre si. Segundo Almeida (2012ALMEIDA, E. Econometria espacial aplicada. Campinas: Alinea, 2012. 498 p.), ao se investigarem unidades espaciais, como as microrregiões, há grande probabilidade de a amostra apresentar autocorrelação espacial, especialmente ao se tratar de atividade economicas como a pecuária, que usualmente se concentra espacialmente, o que exige a incorporação desse efeito no modelo no sentido de evitar problemas advindos de omissão de variável relevante. Além disso, para efeitos de comparação, serão estimados os modelos espaciais por meio de efeito fixo (fixed effect - FE) e de efeito aleatório (random effect - RE). Posteriormente, o teste de Hausman espacial será realizado para identificar qual modelo é o mais adequado para conjunto de dados. Por fim, o critério de informação de Akaike será empregado para verificar a qualidade do ajuste dos modelos, considerando o melhor modelo aquele que apresentar o menor valor.

No presente artigo, serão estimados basicamente três modelos, considerando uma matriz de peso espacial n×n, W, para cada período t=1,…,T, um vetor n×1 da variável dependente, yt e Xt como uma matriz n×n de regressores. O primeiro modelo será o Spatial Autoregressive Model (SAR):

yt=ρWyt+Xtβ+μ+ϵtt=1,…,T (8)

onde μ~N(0,σμ2) se efeitos aleatórios, enquanto μ é um vetor de parâmetros se estimado por efeitos fixos. As hipóteses-padrão sobre o modelo são que ϵit~N(0,σϵ2) e E(∫it∫js)=0 para i≠j e/ou t≠s. O Spatial Error Model (SEM) foca na autocorrelação espacial no termo de erro,

yt=Xtβ+μ+vt (9)

vt=λMvt+ϵt (10)

em que M é a matriz de pesos espaciais que pode ser ou não igual a W. Por fim, o Spatial Durbin Model se caracteriza como uma generalização do modelo SAR, mas com a adição de variáveis independentes defasadas espacialmente como variáveis explicativas,

yt=ρWyt+Xtβ+WZtθ+μ+ϵt (11)

É possível que Zt≠Xt. A estimação dos referidos modelos adotará o método de Máximo Verossimilhança.

Uma análise descritiva prévia dos dados foi realizada com a finalidade de verificar a dinâmica apresentada pelo rebanho bovino brasileiro. Em 1990, o Brasil possuía aproximadamente 151 milhões de cabeças enquanto, em 2015, esse número aumentou para 215 milhões, um crescimento de 42,38% em 25 anos, uma taxa de 1,69% ao ano. Entretanto, autores como Carvalho e Zen (2017CARVALHO, T. B.; ZEN, Sérgio. A cadeia de pecuária de corte no Brasil: evolução e tendências. Revista iPecege, v. 3, n. 1, p. 85-99, 2017.) e Nascimento (2017NASCIMENTO, F. Desmatamento zero na Amazônia: como e por que chegar lá. Greenpeace, 2017.) sustentam que tal crescimento não aconteceu de forma homogênea ao longo do território nacional. Quando se consideram níveis regionais, por exemplo, de acordo com Rodrigues et al. (2018RODRIGUES, M.; SILVA, D. C.; JUNIOR, W. C. A redução do desmatamento na Amazônia: o papel dos mecanismos institucionais. Revista de Administração e Negócios da Amazônia, v. 10, n. 1, p. 108-123, 2018.), “o rebanho bovino nos estados da Amazônia Legal brasileira saltou de 26,2 milhões de cabeças em 1990 para 80,7 milhões de cabeças em 2013, um acréscimo de 207,38% no efetivo, sendo que os estados de Mato Grosso e Pará, juntos, corresponderam por quase 60% desta variação - 32,3 milhões de cabeças adicionais nos dois estados no mesmo período”. Fica evidente, portanto, um expressivo aumento do tamanho do rebanho bovino na Amazônia Legal. Para verificar esse processo em nível local, as Tabelas 1 e 2 trazem as dez microrregiões com maior número de cabeças no ano de 1990 e 2015, respectivamente. A finalidade é investigar a existência de alterações significativas entre as principais regiões produtoras do país, além de verificar como o crescimento do rebanho na Amazônia Legal alterou a disposição das microrregiões.

Thumbnail

Tabela 1:
Número de cabeças por microrregião em 1990

Por meio da Tabela 1, constata-se que, em 1990, as microrregiões com maior número de cabeças estavam, em sua maioria, localizadas na região Centro-Oeste do país, excetuando-se uma microrregião localizada no Rio Grande do Sul, com nenhuma representante da Amazônia Legal. Teixeira e Hespanhol (2006TEIXEIRA, C.;HESPANHOL, A. N. A região Centro-Oeste no contexto das mudanças agrícolas ocorridas no período após 1960. Revista da Associação de Geógrafos Brasileiros, p. 52, 2006.) apontam que a alta produção no Centro-Oeste se deve a uma acelerada modernização pela qual passou a região na década de 1970, especialmente incentivada por aumento na concessão de crédito para a realização de investimentos. Por fim, o somatório do estoque bovino das dez maiores microrregiões produtoras resulta em 14,66% do total nacional.

Nota-se, na Tabela 2, que apenas três das microrregiões que estavam no grupo das dez maiores produtoras em 1990, mantiveram-se entre as maiores em 2015. Desse modo, é evidente um deslocamento da produção bovina no Brasil, com uma mudança substancial na localização dos maiores rebanhos. Ademais, é evidente o avanço de microrregiões localizadas na Amazônia Legal, que passaram a contar com a maioria das microrregiões em 2015. Os números aqui apresentados corroboram, portanto, as hipóteses sustentadas por Carvalho (2007CARVALHO, T. B. Estudo da elasticidade-renda da demanda de carne bovina, suína e de frango no Brasil. Tese (Doutorado) - Universidade de São Paulo, 2007.), Martinelli et al. (2010MARTINELLI, L. A.; NAYLOR, R.; VITOUSEK, P. M.; MOUTINHO, P. Agriculture in Brazil: impacts, costs, and opportunities for a sustainable future. Current Opinion in Environmental Sustainability, v. 2, n. 5-6, p. 431-438, 2010.) e Alencar et al. (2015ALENCAR, A.; PEREIRA, C.; CASTRO, I.; CARDOSO, A. Desmatamento nos Assentamentos da Amazônia: histórico, tendências e oportunidades. Brasília, DF: IPAM, 2015.), que identificaram na região amazônica o principal vetor de atração da produção bovina no país nas últimas décadas. Dessa forma, confirma-se o deslocamento da produção pecuária do bioma Cerrado para o Amazônico, conforme destacado por Rodrigues et al. (2018RODRIGUES, M.; SILVA, D. C.; JUNIOR, W. C. A redução do desmatamento na Amazônia: o papel dos mecanismos institucionais. Revista de Administração e Negócios da Amazônia, v. 10, n. 1, p. 108-123, 2018.).

Thumbnail

Tabela 2:
Número de cabeças por microrregião em 2015

Outro ponto importante a ser destacado é a ocorrência de uma desconcentração da produção, com as principais regiões produtoras passando de 14,66% do total para 13,39%. Apesar de o número de cabeças pertencentes às principais microrregiões terem aumentado de 22.163.203 para 29.041.8653, sua representação no total diminuiu, indicando que as outras regiões apresentaram taxas de crescimento superiores. Para melhor observar a evolução e distribuição da produção bovina no Brasil ao longo do período de análise, foram construídos três mapas de quantile (Figura 1). O primeiro, da esquerda para a direita, refere-se ao estoque de gado nas microrregiões em 1990 enquanto o segundo é para 2015 e o terceiro mostra a taxa de crescimento entre os períodos. Fica evidente que, em 1990, os maiores produtores se localizavam especialmente no Centro-Oeste do país. Por outro lado, em 2015, verifica-se uma distribuição mais igualitária entre as microrregiões, indicando um possível processo de desconcentração do rebanho bovino ao longo do território brasileiro, em especial devido ao aumento da produção nas microrregiões pertencentes à Amazônia Legal.

Figura 1:
Rebanho bovino nas microrregiões brasileiras (1990 e 2015) e a taxa de crescimento no período

As regiões localizadas na Amazônia Legal apresentaram taxas de crescimento mais elevadas do que as demais, fato que evidencia o processo de avanço da produção bovina para essa região. Conforme salienta Dias Filho (2011), diversos elementos foram responsáveis para o avanço da pecuária na fronteira agrícola da Amazônia, destacando-se a abertura de estradas, o estabelecimento de incentivos fiscais e crédito rural subsidiado. Fatores que visavam ao desenvolvimento socioeconômico da região, assim como condições climáticas adequadas, com temperaturas mais uniformes e períodos de estiagem relativamente menos severos e extensos que em outras regiões do país. Ademais, é visível a concentração espacial do rebanho em ambos os anos, assim como da taxa de crescimento. Isso é comprovado por meio da estatística I de Moran da Tabela 3, cujos valores são positivos e estatisticamente significativos independentemente da matriz de convenção aplicada. Portanto, microrregiões com elevado rebanho e taxa de crescimento tenderam a estar rodeados também por regiões com altos valores (e vice-versa). Além disso, observa-se aumento na magnitude dos coeficientes, sinalizando uma intensificação da concentração espacial.

Thumbnail

Tabela 3:
I de Moran para o rebanho brasileiro em 1990, 2015 e para a taxa de crescimento no período

Ao utilizar os indicadores locais de associação espacial (Mapas Lisa) identificou-se a existência de clusters espaciais da produção e da taxa de crescimento da pecuária ao longo do Brasil. Pode-se notar, em 1990, um grande cluster com configuração espacial AA no Centro-Oeste do país, indicando que essa região concentrava espacialmente o maior rebanho do Brasil. Em 2015, por outro lado, verifica-se uma ampliação desse grande cluster para a fronteira agrícola na região Amazônica, em especial, englobando o estado do Mato Grosso, Pará e Rondônia. Em relação à taxa de crescimento do rebanho, destacam-se aglomerações com padrão espacial AA localizadas na região da Amazônia Legal, em especial, nos estados do Pará, Mato Grosso, Rondônia, Acre e Amazônia, o que indica uma expansão da pecuária na região. Tal dinâmica, entretanto, não está imune de externalidades negativas. De acordo com Rivero et al. (2009RIVERO, S.; ALMEIDA, O.; ÁVILA, S.; OLIVEIRA, W. Pecuária e desmatamento: uma análise das principais causas diretas do desmatamento na Amazônia. Nova Economia, Belo Horizonte, v. 19, n. 1, p. 41-66, 2009.), a produção agropecuária, especialmente da pecuária bovina, é a principal atividade causadora de desmatamento e danos ambientais no bioma amazônico. Sendo assim, é evidente a necessidade de maiores investigações dos determinantes do crescimento da produção pecuária no país, assim como de potenciais externalidades.

Figura 2:
Mapas LISA para produção bovina nas microrregiões brasileiras

Identificadas as características básicas da distribuição espacial e do avanço do rebanho no Brasil, o próximo passo será encontrar os determinantes desse avanço, com especial enfoque em suas características estruturais e na Amazônia Legal. Vale destacar a existência de dependência espacial para as variáveis dependentes, identificada por meio do I de Moran dos resíduos do painel de efeito fixo conforme destacado nas Tabelas 1A e 2A no Apêndice. Diante desse contexto, é importante a escolha da matriz de vizinhança, a qual deve respeitar certas propriedades desejáveis e certas particularidades do estudo em questão. Segundo Almeida (2012ALMEIDA, E. Econometria espacial aplicada. Campinas: Alinea, 2012. 498 p.), a escolha da matriz de pesos espaciais deve ser uma matriz W que capte a maior dependência espacial dos resíduos através do seguinte procedimento: a) Estima-se o modelo clássico de regressão linear, b) Testam-se os resíduos desse modelo para autocorrelação espacial, usando o I de Moran para um conjunto de matrizes W; c) seleciona-se a matriz de pesos espaciais que tenha gerado o mais alto valor estatisticamente significativo do I de Moran, conforme pode ser observado nas Tabelas A2 e A3 no Apêndice. O resultado do teste indicou que a matriz de peso espacial de três vizinhos é a mais indicada para a estimação dos modelos espaciais.

Nas Tabelas 4 e 5, têm-se os resultados para os modelos espaciais estimados11 11 Vale mencionar, por fim, que a correlação entre as variáveis utilizadas no modelo (Tabela A1 do apêndice) não apresentaram valores que pudessem comprometer as estimações. com a técnica de dados em painel. No que se refere aos problemas econométricos, tanto no modelo de rebanho bovino quanto no modelo de densidade bovina, foi identificada a presença de heterocedasticidade. Dessa forma, todos os modelos espaciais foram estimados usando o erro padrão robusto Disccoll-Kraay, com a finalidade de controlar o problema. Para a escolha entre o modelo de efeito fixo e o de efeito aleatório, utilizou-se o teste de Hausman espacial, que rejeitou a hipótese nula de que não há diferença sistemática entre os coeficientes estimados,12 12 Estatística Chi² para o modelo de rebanho bovino (Painel de dados: 936.02 SAR: 969.68; SDM: 960.52) e para o modelo de densidade bovina (Painel de dados: 1444.51; SAR: 1466.07; SDM: 47.51). conformando a estimação por efeito fixo como a mais indicada.

Posteriormente, definiu-se o modelo espacial que melhor representa e capta a dinâmica espaçotemporal do fenômeno em estudo. Tanto a partir do critério de informação de Akaike quanto da análise residual das regressões observadas, verificou-se que o modelo Durbin Espacial (SDM) é aquele com menor valor (Tabelas 4, 5 e A3 do Apêndice). Portanto, as análises posteriores serão realizadas tendo como referência esse modelo. Vale destacar que, apesar de ainda haver, em certa medida, autocorrelação espacial nos resíduos do modelo, essa se reduziu consideravelmente com a incorporação das defasagens espaciais. De qualquer maneira, é imperioso cautela ao considerar os coeficientes diante desse fato. Por outro lado, a grande quantidade de informações e anos considerados, conjuntamente ao alto grau de significância estatística de muitas das variáveis e a constância dos coeficientes em diferentes modelos, indica estabilidade estrutural para os resultados, minimizando o problema.

Os resultados das estimações apontam para a existência de um processo de desconcentração da pecuária de corte nas microrregiões brasileiras, ou seja, regiões com menor rebanho e menor densidade apresentaram taxas de crescimento superiores. O fato é captado pelo sinal negativo das variáveis Rebanho.Bov.i,t−1 e Densidade.Bov.i,t−1 que apresentaram significância estatística no nível de 1%. Desse modo, pode-se afirmar que no período de 1995 a 2016 houve avanço da pecuária bovina em microrregiões que anteriormente possuíam um rebanho relativamente menor, comprovando a hipótese levantada por diversos autores (Carvalho, 2007CARVALHO, T. B. Estudo da elasticidade-renda da demanda de carne bovina, suína e de frango no Brasil. Tese (Doutorado) - Universidade de São Paulo, 2007.; Brondizio, 2009BRONDIZIO, E. S. Análise inter-regional de mudança do uso da terra na Amazônia. In: MORAN, E. F.; OSTROM, E. (Org.). Ecossistemas Florestais: interação homem-ambiente. São Paulo: EDUSP, p. 289-326, 2009.; Barona et al., 2010BARONA, E; RAMANKUTTY, N; HYMAN, G; COMMES, O. The role of pasture and soybean in deforestation of the Brazilian. Environmental Research Letters, 5, 024002 (9pp), 2010.; Martinelli et al., 2010MARTINELLI, L. A.; NAYLOR, R.; VITOUSEK, P. M.; MOUTINHO, P. Agriculture in Brazil: impacts, costs, and opportunities for a sustainable future. Current Opinion in Environmental Sustainability, v. 2, n. 5-6, p. 431-438, 2010.; Alencar et al., 2015ALENCAR, A.; PEREIRA, C.; CASTRO, I.; CARDOSO, A. Desmatamento nos Assentamentos da Amazônia: histórico, tendências e oportunidades. Brasília, DF: IPAM, 2015.; Domingues, 2017DOMINGUES, M. S.; BERMANN, C.; SIDNEIDE, S.; A produção da soja no Brasil e sua relação com o desmatamento na Amazônia. Revista Presença Geográfica, v. 1, n. 1, 2017.; e Silva et al., 2018SILVA, M. M.; OLIVEIRA, F. A.; SANTANA, A. C. Mudanças socioambientais no uso da terra em Altamira, Amazônia Oriental. Novos Cadernos NAEA, v. 20, n. 3, 2018.).

Além disso, verificou-se que a taxa de crescimento do rebanho e densidade bovina é capaz de gerar transbordamentos espaciais positivos na dinâmica de crescimento de seus vizinhos próximos, fato captado pelo coeficiente espacial ρ, que se mostrou significativo no nível de 1% em ambos os casos. Conforme destaca Silva et al. (2012SILVA, J. C.; SOUZA, E. C.; MARTINS, H. E. P. Produção agropecuária em municípios de Minas Gerais (1996-2006): padrões de distribuição, especialização e associação espacial. Revista de Economia e Sociologia Rural, v. 50, n. 2, p. 333-349, 2012.), a existência de externalidades espaciais positivas torna o setor pecuário dinâmico, em especial, devido aos spillovers de produtividade entre as regiões, o que induz o espraiamento dessa produção ao longo do território, contribuindo para a formação de clusters espaciais de alta produção e/ou produtividade. Além disso, o autor destaca iniciativas governamentais,13 13 As principais contribuições governamentais passam por aumentos na cobertura de vacinas aplicadas contra a febre aftosa, intensificação de medidas sanitárias, apoio ao desenvolvimento e disseminação de melhoramentos genéticos, incentivo ao crédito e inspecionamento da produção. além do desenvolvimento da cadeia produtiva nas regiões, como fatores adicionais que podem induzir a concentração e transbordamento espacial da pecuária de corte. Por fim, de acordo com a Nova Geografia Econômica, tal fenômeno ocorre porque determinadas atividades são atraídas e aglomeradas em uma determinada localidade diante de forças atrativas (centrípetas) próprias, o que reforça a importância de se considerar o espaço e seus possíveis efeitos quando da análise da dinâmica de crescimento de determinada variável entre as regiões (Krugman, 1991KRUGMAN, P. Geography and trade. Cambridge, MA: MIT Press, 1991.), o que de fato foi confirmado pelo presente artigo para a pecuária bovina brasileira.

Para verificar a hipótese de que as regiões localizadas na Amazônia Legal apresentam uma taxa de crescimento estruturalmente maior do que as demais microrregiões brasileiras, mesmo após o controle de outros fatores relevantes, empregaram-se variáveis dicotômicas de interação espaçotemporal (descritas na subseção 3.1). Essas variáveis, referentes ao período de 1997 a 201514 14 Para não cair na armadilha da variável dummy, excluiu-se a variável dicotômica para o primeiro período (1995).

Como a pecuária contribuiu para a expansão territorial?

Como a pecuária influenciou a expansão do Brasil colonial?

Qual foi a importância da pecuária para o território brasileiro?

Que impacto a expansão da pecuária?