# Evolução temporal da relação entre os efetivos ovinos e caprinos nas grandes regiões brasileiras (2015–2025) **Conjunto de dados e scripts R para reprodução das análises estatísticas** ## 1. Descrição Este repositório disponibiliza o conjunto de dados e os scripts em linguagem R utilizados no estudo “Evolução temporal da relação entre os efetivos ovinos e caprinos nas grandes regiões brasileiras (2015–2025)”. O material permite reproduzir a preparação e auditoria dos dados, as análises descritivas, a estimação das tendências temporais regionais, a modelagem por mínimos quadrados generalizados (GLS) e a construção das figuras científicas. O banco contempla cinco grandes regiões brasileiras e onze anos consecutivos, totalizando 55 observações região-ano. ## 2. Fonte dos dados Os dados foram obtidos da Pesquisa da Pecuária Municipal (PPM), disponibilizada pelo Instituto Brasileiro de Geografia e Estatística (IBGE), por meio do Sistema IBGE de Recuperação Automática (SIDRA), Tabela 3939. Fonte: https://sidra.ibge.gov.br/tabela/3939 A planilha de entrada utilizada no processamento é `tabela3939_ovinos_caprinos_2015_2025.xlsx`. A variável central do estudo corresponde à razão entre os efetivos ovinos e caprinos: **Razão O/C = Efetivo de ovinos / Efetivo de caprinos** Valores superiores a 1 indicam predominância numérica de ovinos em relação aos caprinos na região e no ano considerados. ## 3. Arquivos disponibilizados O depósito contém sete arquivos: | Arquivo | Descrição | |---|---| | `tabela3939_ovinos_caprinos_2015_2025.xlsx` | Planilha de dados utilizada nas análises | | `01_preparacao_auditoria_analise_temporal_oc.R` | Preparação, auditoria, estatísticas descritivas e tendências temporais | | `02_figura_evolucao_temporal_oc.R` | Construção da Figura 1 | | `03_modelagem_gls_contrastes_oc.R` | Modelagem GLS e contrastes entre inclinações regionais | | `04_figura_inclinacoes_gls_oc.R` | Construção da Figura 2 | | `README.md` | Documentação e instruções de reprodução | | `dicionario_variaveis.csv` | Descrição das variáveis e respectivas codificações | As tabelas, os resultados intermediários e as figuras **não integram os arquivos depositados**, pois podem ser reproduzidos pela execução dos scripts. ## 4. Procedimentos estatísticos ### 4.1. Preparação e auditoria O banco foi organizado em estrutura longitudinal, com observações anuais para as cinco grandes regiões brasileiras entre 2015 e 2025. Foram realizadas verificações de integridade, duplicidade de combinações região-ano e presença de valores ausentes. ### 4.2. Análise das tendências temporais Foram ajustadas regressões lineares regionais utilizando o logaritmo da razão O/C como variável resposta. A inferência considerou erros-padrão robustos à autocorrelação temporal, estimados pelo procedimento de Newey–West com defasagem de uma unidade temporal. A variação percentual anual (VPA) foi calculada por: VPA = 100 × [exp(β) − 1] em que β representa o coeficiente angular estimado na escala logarítmica. ### 4.3. Modelagem GLS Foi utilizado um modelo de mínimos quadrados generalizados, incluindo região, tempo e interação região × tempo. O modelo incorporou estrutura de correlação residual autorregressiva de primeira ordem (AR1) dentro de cada região e variâncias residuais heterogêneas entre regiões. A estimação foi realizada por máxima verossimilhança (ML). A interação região × tempo foi avaliada por teste de razão de verossimilhanças entre modelos aninhados. As inclinações temporais regionais foram estimadas com o pacote `emmeans`, utilizando o método de Satterthwaite para aproximação dos graus de liberdade. Os contrastes entre inclinações foram avaliados com ajuste `mvt` para comparações múltiplas. O nível de significância adotado foi de 5%. ## 5. Ambiente computacional As análises foram desenvolvidas em linguagem R, com utilização do ambiente RStudio. Entre os pacotes empregados nos procedimentos estão `readxl`, `dplyr`, `tidyr`, `ggplot2`, `nlme`, `emmeans`, `sandwich` e `lmtest`. As dependências efetivas e suas versões devem ser verificadas no ambiente de execução e nos comandos de carregamento dos scripts. ## 6. Reprodução das análises Para reproduzir os procedimentos: 1. Instalar o R e os pacotes requeridos. 2. Manter a planilha de entrada e os quatro scripts R no mesmo diretório. 3. Definir esse diretório como diretório de trabalho. 4. Executar os scripts na ordem numérica. Comandos: ```r source("01_preparacao_auditoria_analise_temporal_oc.R") source("02_figura_evolucao_temporal_oc.R") source("03_modelagem_gls_contrastes_oc.R") source("04_figura_inclinacoes_gls_oc.R") ``` O Script 01 prepara o banco longitudinal e gera os arquivos necessários às etapas posteriores. Os Scripts 02 e 04 produzem as figuras científicas, enquanto o Script 03 executa a modelagem GLS e exporta as estimativas e os contrastes regionais. Os resultados são armazenados automaticamente na subpasta `resultados/`. Essa subpasta é produzida durante a execução e não precisa estar presente no depósito original. ## 7. Reprodutibilidade A sequência de scripts permite reproduzir os procedimentos estatísticos e os produtos analíticos do estudo a partir da planilha de entrada. O conjunto foi organizado para possibilitar a verificação independente dos cálculos, modelos, estimativas, intervalos de confiança, contrastes e representações gráficas. Pequenas diferenças numéricas podem ocorrer em procedimentos que utilizam aproximações ou integração numérica, dependendo das versões dos pacotes e do ambiente computacional. ## 8. Limitações Os dados representam efetivos pecuários agregados por grande região brasileira e ano. A razão O/C expressa a relação numérica entre os efetivos ovinos e caprinos, não constituindo uma medida direta de produtividade, eficiência econômica ou desempenho zootécnico. Os resultados descrevem tendências e associações temporais, sem estabelecer relações causais. ## 9. Autoria e instituição **Autores:** 1. Cyro Rego Cabral Junior 2. Tobyas Maia de Albuquerque Mariz **Afiliação institucional:** Universidade Federal de Alagoas (UFAL), Campus de Engenharias e Ciências Agrárias (CECA), Rio Largo, Alagoas, Brasil. ## 10. Identificador persistente e citação O conjunto de dados e scripts R associado a este estudo foi depositado no repositório Zenodo. **DOI:** https://doi.org/10.5281/zenodo.23901754 O depósito destina-se à disponibilização dos dados e códigos necessários à reprodução das análises estatísticas, sem incluir o manuscrito científico, as tabelas ou as figuras produzidas pelos scripts. A referência bibliográfica completa deverá incluir os autores, o ano de publicação, o título e a versão registrados no Zenodo. ## 11. Licença e reutilização A licença será definida antes da publicação do depósito, considerando as condições de reutilização dos dados originais do IBGE e dos scripts desenvolvidos pelos autores. O conjunto destina-se à transparência metodológica, à reprodução das análises e ao uso acadêmico e científico, respeitadas as condições da licença selecionada.