Relatório · Fase 2

Câncer de esôfago: o que as áreas com excesso têm em comum

Geração de hipóteses com dados públicos de produção agrícola, água, ruralidade, PIB e geografia, para orientar um estudo com pacientes. English version.

Autor: Clériston Almeida Capistrano (ORCID 0009-0008-5280-1869) · pesquisador independente, sem formação na área da saúde · cledon.org

Resumo

As Fases 1 e 1b mostraram um excesso persistente de mortalidade por câncer de esôfago em áreas do interior de Minas Gerais, do Sul e da Bahia. A Fase 2 descreve que características municipais, em dados públicos, acompanham esse excesso, com plano de análise registrado no git antes dos dados (commit 0d88889). Esta versão incorpora uma revisão técnica assistida por IA (não é revisão independente por pares): dados ausentes deixaram de ser tratados como zero, e as análises foram refeitas.

• Comparando cada município com o esperado pelas taxas do seu estado, a mortalidade por câncer de esôfago é maior onde há mais água de fonte, nascente ou poço, extração de erva-mate, população rural, produção de lenha e de fumo, e menor onde o PIB municipal per capita é maior. As direções se mantêm com efeitos fixos de estado, com a regra antiga de dados ausentes e com medidas de 2000 e 2010.
• Para pulmão e para os outros cânceres, as variáveis rurais e de água têm associação em direção oposta, e a diferença em relação ao esôfago é estatisticamente clara em 8 das 10 variáveis. Isso indica que o padrão não é o mesmo para todos os cânceres; não afasta vieses de registro ou de diagnóstico que afetem cânceres de formas diferentes.
• Por região (análise não pré-especificada): no Sul, o excesso acompanha a erva-mate extraída e o fumo; em Minas Gerais, onde quase não há mate nem fumo, acompanha ruralidade, menor PIB per capita e menor altitude, e tem associação inversa com a área de café.
• A comparação dos 24 agrupamentos com seus vizinhos é mais fraca: após correção para as 10 comparações, nenhuma variável fica com q < 0,05.

1. Método

Desfecho. Óbitos por câncer de esôfago de 2015 a 2024 por município, com esperados calculados pelas taxas do próprio estado por sexo e faixa etária, usados como offset. Isso compara cada município com o padrão do seu estado, mas não equivale a incluir efeitos de estado na regressão: por isso há uma sensibilidade com efeitos fixos de estado.

Exposições e unidades. Fumo, erva-mate cultivada e café: área colhida da cultura como percentual da área territorial do município (PAM/IBGE). Erva-mate extraída: toneladas por mil habitantes; lenha: m³ por habitante, somando extração e silvicultura (PEVS/IBGE). Produção em média de 2005 a 2010. Água de poço e de fonte ou nascente: % dos domicílios, Censo 2022. População rural: %, Censo 2022. PIB municipal per capita de 2010 (não é renda domiciliar). Altitude do centroide (Copernicus DEM). As variáveis de produção entram como log(1 + x) e todas são padronizadas: a razão de taxas (RT) é por desvio-padrão da variável transformada, não por unidade original.

Dados ausentes. No IBGE, “-” é zero absoluto; “...” é dado não disponível. Na primeira versão, “...” e municípios ausentes viravam zero, o que confundia falta de informação com falta de produção. Agora: “-” = 0; “...”, “X” e ausência na resposta = ausente; a média usa só os anos com dado; sem nenhum ano, o município fica fora da análise daquela variável. A extração vegetal tem “...” no mesmo conjunto de municípios e anos para erva-mate e lenha, concentrado em SP, RJ, PR, SC e MG; a silvicultura omite muitos municípios da resposta, por isso a lenha (extração + silvicultura, como no plano) tem dado em poucos municípios e há uma sensibilidade só com a extração.

Tabela 1. Auditoria dos símbolos nas tabelas de produção (6 anos × 5.570 municípios).

Modelo. Regressão de Poisson com offset, erros-padrão sanduíche agrupados por região imediata do IBGE (para dependência espacial e sobredispersão), Benjamini–Hochberg entre as 10 variáveis.

2. Associações por variável

RT de mortalidade por desvio-padrão de cada variável, com IC 95%: esôfago (azul) e os desfechos de comparação pulmão (laranja) e outros cânceres, exceto esôfago (verde).

Tabela 2. Esôfago: estimativa principal e sensibilidades.

3. Desfechos de comparação

Pulmão e todos os cânceres exceto esôfago. Não são controles negativos no sentido formal: o pulmão compartilha exposições (tabaco), e vieses de registro ou diagnóstico podem afetar cânceres de formas diferentes. Servem para ver se o padrão é próprio do esôfago. A diferença foi testada formalmente num modelo conjunto, nos mesmos municípios.

A razão das RT acima de 1 indica associação mais forte (ou mais positiva) para o esôfago que para o desfecho de comparação. Não houve evidência estatística de diferença para café e altitude.

4. Modelos conjuntos

Com as 10 variáveis, só entram os municípios com dado em todas (a lenha limita muito a amostra); por isso há também um modelo sem lenha e sem erva-mate extraída. Várias variáveis são correlacionadas (erva-mate cultivada e extraída, ρ = 0,81; ruralidade e PIB, ρ = −0,47).

5. Agrupamentos versus vizinhos

Média dentro de cada um dos 24 agrupamentos menos a média dos vizinhos do mesmo estado (até 100 km, fora dos agrupamentos), em desvios-padrão. O número de agrupamentos varia com a disponibilidade de dados (por exemplo, 21 para lenha e 23 para erva-mate extraída; 24 nas demais). Teste do sinal sem empates e correção de Benjamini–Hochberg entre as 10 variáveis. Os agrupamentos compartilham vizinhos, então as comparações não são independentes.

Os agrupamentos tendem a ter PIB per capita menor que os vizinhos (19 de 24), mas, após a correção, nenhuma variável fica com q < 0,05 (PIB: q = 0,07). Esta análise acrescenta pouca evidência às anteriores.

6. Sul e Minas Gerais (não pré-especificado)

Modelos separados por região; variáveis com menos de 20 municípios com valor acima de zero na região foram omitidas (—), e a lenha ficou fora dos modelos conjuntos.

Sul: erva-mate extraída e fumo seguem associados no modelo conjunto, com água de poço e menor PIB; a área de café tem associação inversa. Minas Gerais: ruralidade, menor PIB e menor altitude; a área de café também tem associação inversa; a água de poço, forte na análise univariada, perde força no modelo conjunto. Nenhuma das variáveis de produção medidas tem associação positiva estatisticamente clara com o excesso em Minas Gerais.

7. Cronologia das medidas

Água e ruralidade vêm do Censo 2022, depois do início dos óbitos (2015), e não documentam as condições anteriores ao adoecimento. Como comparação histórica (não pré-especificada), a população rural de 2000 e de 2010 e os domicílios com poço ou nascente em 2010 (o Censo 2010 não separa poço de nascente). A ruralidade é muito estável (correlação com 2022: 0,93 para 2000 e 0,95 para 2010); poço em 2022 e poço ou nascente em 2010 têm correlação 0,73.

Maior presença de poços não demonstra contaminação da água. Nitrato e outros contaminantes são hipóteses a investigar, não achados desta análise.

8. Os 16 municípios

Valor no município e, entre parênteses, a mediana do estado; “s/d” = sem dado. 14 dos 16 têm proporção de domicílios com água de poço acima da mediana do estado em 2022 (14 com poço ou nascente em 2010).

9. Hipóteses para um estudo com pacientes

Apoiadas, em nível municipal, por estes dados:

• Sul: erva-mate (produção extrativa como indicador regional) e fumo (área de cultivo).
• Minas Gerais e Sul: ruralidade e menor PIB municipal per capita.
• Abastecimento por poço: no Sul, a associação se mantém no modelo conjunto; em Minas Gerais, aparece na análise individual, mas perde precisão no modelo conjunto (o intervalo de confiança inclui 1). Fonte ou nascente: associação na análise individual e no modelo conjunto nacional sem lenha, mas não nos modelos conjuntos do Sul e de Minas Gerais.

Motivadas pela literatura, não testadas aqui (sem dado municipal adequado): temperatura do chimarrão e de outras bebidas (a IARC classificou bebidas muito quentes, acima de 65 °C, como provavelmente carcinogênicas, grupo 2A); fumo de rolo, de corda ou de palha; cachaça e destilados caseiros; fogão a lenha dentro de casa; nitrato e outros contaminantes da água; dieta pobre em frutas e verduras. A Fase 2b (https://cledon.org/fase2b) aprofunda Minas Gerais.

10. Limitações

• Estudo ecológico: associações entre municípios não indicam o que acontece com as pessoas (falácia ecológica).
• Produção local é aproximação grosseira do consumo; lenha produzida não mede uso de fogão a lenha.
• Dados ausentes não aleatórios: a extração vegetal falta sobretudo em alguns estados, o que pode enviesar as estimativas dessas variáveis.
• Confundimento residual por pobreza e ruralidade é provável; PIB per capita não é renda domiciliar.
• Análises por região e sensibilidades foram feitas depois de ver os resultados nacionais.
• Análises por pesquisador independente, sem formação na área da saúde, com assistência de IA (Claude, da Anthropic) para código e texto, sob revisão humana. Nenhum modelo de aprendizado de máquina foi treinado ou avaliado.

11. Reprodução

Código em https://cledon.org/codigo. Plano: PLANO-fase2.md (commit 0d88889). Depois das Fases 1 e 1b: exposicoes.py, hipoteses.py, gerar_fase2.py. Resultados em CSV na página do mapa e no depósito do Zenodo (fase2_*.csv).

Em memória de Francisco Cledio Capistrano (F.C.C.).