# Plano de análise — Fase 2b: Minas Gerais (registrado antes dos dados)

Registrado no git antes de baixar as variáveis abaixo e antes de cruzá-las com a mortalidade. Até aqui foram consultados
apenas catálogos (nomes de tabelas, códigos, anos disponíveis) do IBGE e do CNES. Plano registrado localmente, não um
pré-registro formal. O que não estiver aqui será rotulado como "não pré-especificado".

## Pergunta

Na Fase 2, o excesso de mortalidade por câncer de esôfago em Minas Gerais acompanhou ruralidade, menor PIB municipal per
capita e menor altitude, sem nenhuma variável de produção associada de forma positiva. A Fase 2b pergunta: (1) exposições
ainda não medidas, mais antigas ou específicas, acompanham esse excesso? (2) o acesso a diagnóstico e tratamento ajuda a
explicá-lo como artefato? Objetivo: gerar hipóteses para um estudo com pacientes, não testar causas.

## Segurança dos dados

Somente dados públicos e agregados: por município (IBGE) ou por estabelecimento de saúde (CNES, sem pacientes). Nenhum
microdado com identificação de pessoas é baixado ou publicado. Arquivos brutos ficam em `dados/brutos/fase2b/`, fora do
git. Nenhuma credencial é usada.

## Variáveis (definidas antes de ver os dados)

| Código | Variável | Fonte |
|---|---|---|
| `aguardente` | aguardente de cana produzida na agroindústria rural, por mil habitantes de 2010 | Censo Agropecuário 2006, tabela 1724 (produto 111701); sensibilidade com 2017, tabela 6960 |
| `milho` | área colhida de milho em grão como % da área municipal, média 2005–2010 (micotoxinas como hipótese) | PAM, tabela 1612, produto 2711 |
| `agua_pn_1991` | % de domicílios abastecidos por poço ou nascente, 1991 | Censo 1991, tabela 179 |
| `agua_pn_2000` | % de domicílios abastecidos por poço ou nascente, 2000 | Censo 2000, tabela 1436 |
| `rural_1991` | % da população rural, 1991 | tabela 202 |
| `endoscopia` | endoscópios digestivos (CNES, equipamento 0333) por 100 mil habitantes na região imediata do IBGE, dezembro de 2015 | CNES, arquivos EQ; sensibilidade com dezembro de 2020 |
| `dist_oncologia` | distância (km, entre centroides) até o município mais próximo com hospital habilitado em oncologia (CACON ou UNACON: códigos 1701–1709, 1712, 1713), dezembro de 2015 | CNES, arquivos HB |

Mesmas regras de dados ausentes da Fase 2 ("-" = 0; "...", "X" e ausência = ausente). Variáveis de produção e
`endoscopia` entram como log(1 + x); todas são padronizadas (RT por desvio-padrão).

## Análises

**A. Minas Gerais (principal).** Desfecho e modelo da Fase 2: óbitos por câncer de esôfago 2015–2024, esperados do
estado como offset, Poisson com erros-padrão agrupados por região imediata. RT por desvio-padrão de cada variável;
Benjamini–Hochberg entre as 7 variáveis.

**B. Acréscimo ao modelo da Fase 2 em MG.** Cada variável nova acrescentada, uma de cada vez, ao modelo com ruralidade
(2022), PIB per capita (2010) e altitude: a variável nova continua associada após esse ajuste?

**C. Desfechos de comparação em MG.** As mesmas RT para pulmão e para outros cânceres (exceto esôfago), com teste formal da
diferença (modelo empilhado), como na Fase 2.

**D. Brasil (secundária).** As mesmas variáveis no Brasil, sem e com efeitos fixos de estado.

**E. Leitura prevista das variáveis de acesso.** Se o excesso fosse artefato de mais diagnóstico, esperaríamos mais
mortalidade registrada onde há mais endoscópios e hospitais de oncologia por perto. Associação na direção oposta (mais
mortalidade com menos acesso) torna essa explicação menos provável, mas é compatível com diagnóstico tardio; como o câncer
de esôfago tem letalidade muito alta em qualquer cenário, o efeito do acesso sobre a mortalidade tende a ser pequeno. O
acesso também pode mudar a codificação da causa de morte (mais causas mal definidas onde há menos acesso), por isso o
modelo B é repetido com a % de causas mal definidas.

## Ressalvas

Estudo ecológico; produção local não mede consumo; distância entre centroides não é tempo real de viagem; o CNES registra
equipamentos existentes, não exames feitos; vários testes (corrigidos dentro de cada análise). O nitrato na água
(SISAGUA) ficou fora porque o portal de dados abertos do Ministério da Saúde não está acessível deste ambiente; pode entrar
depois, se os arquivos públicos forem obtidos por outro meio, como análise não pré-especificada.