# Plano de análise — Fase 1b (registrado antes de examinar os resultados)

Registrado no git antes de rodar as análises abaixo. Até este registro, da base do RHC foram vistos apenas
os nomes dos campos e três registros de exemplo (de Salvador, BA), para conhecer o formato. Nenhum
resultado dos 16 municípios, de agrupamentos espaciais ou de tendências foi calculado.

Isto é um plano registrado localmente, não um pré-registro formal em plataforma pública.
Tudo o que estiver fora deste plano será rotulado como "não pré-especificado".

## A. Cruzamento com o Registro Hospitalar de Câncer (RHC/INCA) — `rhc.py`

**Pergunta.** Os 16 municípios com sinal de mortalidade por câncer de esôfago (após os critérios geográficos)
também têm, entre os casos de câncer de seus residentes atendidos em hospitais com RHC, uma fração de câncer de
esôfago maior que a esperada?

**Dados.** Bases nacionais do Integrador RHC (INCA), anos de primeira consulta 2015–2022 (análise principal) e
2015–2023 (sensibilidade, porque o ano mais recente tende a estar incompleto). Município de residência = campo
`PROCEDEN` (código IBGE de 7 dígitos, reduzido a 6). Incluídos: casos analíticos (`TPCASO` = 1), com sexo e idade
informados, código CID-O de localização `LOCTUDET` entre C00 e C97, excluindo C44 (pele não melanoma, de registro
muito variável). Câncer de esôfago: `LOCTUDET` = C15.

**Limitação de partida.** O RHC não é um registro de base populacional: cobre só os hospitais que enviam dados, e o
mesmo paciente pode aparecer em mais de um hospital. Por isso a medida principal não usa a população.

**Medida principal — razão de incidência proporcional (RIP).** Para cada município, esôfago observado ÷ esperado,
em que esperado = Σ (casos de câncer de residentes no estrato sexo × faixa etária de 5 anos, 80+) × (fração de
esôfago entre os casos do mesmo estrato no próprio estado). IC 95% exato (Garwood); teste de Poisson unilateral.
A RIP é robusta a diferenças de cobertura do RHC entre municípios, desde que a cobertura não dependa da localização
do tumor.

**Critério principal.** O RHC "corrobora" o sinal de um município se RIP > 1 e p < 0,05. Resumo: quantos dos 16
corroboram; RIP agrupada por região (MG, Sul = PR/SC/RS, BA), como na validação.

**Contraste.** A mesma RIP para todos os municípios dos estados envolvidos. Correlação de Spearman, entre os
municípios com ≥ 5 óbitos esperados por câncer de esôfago no SIM 2020–2024, entre a RMP de mortalidade (sem
suavização, em relação ao estado) e a RIP do RHC.

**Descritivos (sem limiar).** Fração de carcinoma epidermoide (`TIPOHIST` 8050–8089) entre os casos de esôfago
nos 16 municípios e no estado; razão casos no RHC ÷ óbitos no SIM (como indicador de cobertura); casos com
informação de álcool e tabaco, se houver (esperamos muitos "sem informação").

## B. Agrupamentos espaciais — `varredura.py`

**Método.** Estatística de varredura espacial de Kulldorff, modelo de Poisson, sobre os centroides dos municípios.
Janelas circulares centradas em cada município, com os vizinhos mais próximos, até 100 km de raio e até
1 milhão de habitantes (Censo 2022). Óbitos esperados em relação ao próprio estado (como na análise principal),
de modo que os agrupamentos são excessos dentro do estado, não diferenças entre estados. Só janelas com excesso
(observado > esperado). Valor de p por Monte Carlo com 999 réplicas (distribuição multinomial dos óbitos
condicionada ao total). Agrupamento mais provável e secundários sem sobreposição com p < 0,05.

**Análise principal.** Câncer de esôfago, 2020–2024. **Secundárias.** Os outros 12 grupos (13 varreduras, sem
correção entre elas; resultados tratados como exploratórios). **Validação.** Os agrupamentos de esôfago
encontrados em 2020–2024 são testados em 2015–2019 (RMP agrupada da área, em relação ao estado; replicação se
RMP > 1 e p < 0,05).

## C. Tendências anuais — `tendencias.py`

**Dados.** Óbitos por ano (2015–2024), município, sexo e faixa etária, lidos de novo dos mesmos arquivos do SIM.
População por ano: interpolação linear entre os Censos 2010 e 2022 para 2015–2021; Censo 2022 para 2022–2024.
Esperados em relação ao estado no mesmo ano.

**Análise.** Para os 16 municípios (agrupados por região e no total) e para os agrupamentos de esôfago da parte B:
RMP por ano e regressão de Poisson dos óbitos observados no ano, com log dos esperados como offset. A inclinação
estima a variação anual da RMP. Sem limiar de decisão: é descritivo. Como os municípios foram escolhidos por
terem excesso em 2020–2024, a inclinação tende a ser positiva por construção (regressão à média no sentido
inverso); a comparação de interesse é se o excesso já existia em 2015–2019 (ver validação) e se é estável.

## Ressalvas que valem para tudo

Estudos ecológicos e exploratórios. Nenhum resultado permite estimar risco individual nem estabelecer causa.
Concordância entre SIM e RHC pode refletir vieses comuns (diagnóstico, codificação, endereço).