1 Objetivo e principais resultados

Este relatório avalia a associação entre renda domiciliar per capita e posse de plano de saúde médico na PNS 2019, considerando o desenho amostral complexo. Os resultados são recalculados a cada compilação. O estudo é transversal e não permite interpretar as associações como efeitos causais.

A amostra analítica inclui 90.814 pessoas. A cobertura ponderada de plano médico foi de 26,55% (IC95%: 25,79–27,33%).

No modelo ajustado, a razão de chances por duplicação da renda positiva foi de 2,27 (IC95%: 2,19–2,36). O intervalo está inteiramente acima de 1, indicando associação positiva entre renda e posse de plano, condicionada às covariáveis.

A OR compara chances, não probabilidades; seu valor não representa a multiplicação da probabilidade de possuir plano.

2 Cobertura e renda

2.0.1 Formação da amostra analítica

Etapa N
Registros importados 293.726
Elegiveis com peso valido 90.846
Renda ausente 22
Cor ou raca ausente ou ignorada 10
Idade ausente ou fora de 15 a 130 0
Casos completos 90.814

Fonte: elaboração própria com os microdados da PNS 2019, conforme o script original.

2.0.2 Renda média segundo posse de plano

Plano N Media Inferior Superior
0 70.256 1.005,811 989,7731 1.021,850
1 20.558 3.191,318 3.054,7398 3.327,897

Fonte: elaboração própria com os microdados da PNS 2019, conforme o script original.

2.0.3 Cobertura por decis de renda e grupo zero

Grupo N Minimo Maximo Percentual_ponderado Cobertura Inferior Superior Tipo
Zero 1.191 0 0 0,8800 9,5007 7,0598 12,6704 Faixa exclusiva
D1 11.105 1 300 10,0160 2,0142 1,6325 2,4829 Faixa exclusiva
D2 8.925 301 480 9,9027 5,7326 4,9502 6,6300 Faixa exclusiva
D3 9.230 481 625 10,1224 9,7606 8,5984 11,0608 Faixa exclusiva
D4 7.675 626 792 9,6091 13,6593 12,2326 15,2236 Faixa exclusiva
D5 11.982 793 998 11,8446 16,2712 14,9352 17,7017 Faixa exclusiva
D6 6.617 999 1.165 8,0366 24,3818 22,3357 26,5513 Faixa exclusiva
D7 8.057 1.166 1.499 9,9486 28,7001 26,9407 30,5265 Faixa exclusiva
D8 8.093 1.500 1.999 9,8697 37,2468 35,2663 39,2712 Faixa exclusiva
D9 8.494 2.000 3.127 9,8585 52,6983 50,6540 54,7337 Faixa exclusiva
D10 9.445 3.129 203.196 9,9116 78,5420 76,9622 80,0415 Faixa exclusiva
Positiva (agregado) 89.623 1 203.196 99,1200 26,7031 25,9335 27,4871 Agregado: nao somar as faixas

Fonte: elaboração própria com os microdados da PNS 2019, conforme o script original.

A cobertura por grupos de renda descreve a desigualdade observada antes do ajuste. Diferenças entre grupos podem refletir tanto renda quanto diferenças de composição por idade, escolaridade e outras características. As contagens são não ponderadas; as medidas de cobertura e renda são ponderadas conforme o script.

3 Modelo ajustado

3.0.1 Razões de chances do modelo principal

Termo Beta OR Inferior Superior p
L 1,1821 2,2691 2,1859 2,3556 < 0,001
renda_zero -0,8527 0,4263 0,3080 0,5900 < 0,001
sexoMulher 0,1361 1,1458 1,0676 1,2298 < 0,001
cor_racaPreta -0,2436 0,7838 0,6996 0,8782 < 0,001
cor_racaAmarela 0,0165 1,0166 0,6633 1,5583 0,940
cor_racaParda -0,2206 0,8021 0,7411 0,8681 < 0,001
cor_racaIndigena -0,6433 0,5256 0,3538 0,7807 0,001
escolaridadeSem instrucao -1,6237 0,1972 0,1460 0,2663 < 0,001
escolaridadeFundamental incompleto -0,9804 0,3752 0,3388 0,4154 < 0,001
escolaridadeFundamental completo -0,4330 0,6486 0,5633 0,7467 < 0,001
escolaridadeMedio incompleto -0,3539 0,7019 0,6026 0,8177 < 0,001
escolaridadeSuperior incompleto 0,3463 1,4139 1,2156 1,6445 < 0,001
escolaridadeSuperior completo 0,7502 2,1175 1,9210 2,3341 < 0,001
regiaoNorte -0,6644 0,5146 0,4511 0,5870 < 0,001
regiaoNordeste -0,5193 0,5949 0,5413 0,6539 < 0,001
regiaoSul -0,3115 0,7323 0,6599 0,8127 < 0,001
regiaoCentro-Oeste -0,3611 0,6969 0,6234 0,7792 < 0,001

Fonte: elaboração própria com os microdados da PNS 2019, conforme o script original.

O modelo controla idade, sexo, cor ou raça, escolaridade e região. Para renda positiva, a OR corresponde à duplicação da renda. O indicador de renda zero compara esse grupo à renda positiva de R$ 998, mantendo constantes os demais controles. Nas variáveis categóricas, cada OR deve ser lida em relação à respectiva referência. Os resultados não permitem concluir que uma mudança individual na renda causaria a alteração estimada na cobertura.

4 Probabilidades padronizadas

Renda (R$) Probabilidade (%) IC95% inferior (%) IC95% superior (%)
0 11,35 8,62 14,81
250 5,71 5,16 6,32
500 11,69 10,95 12,48
1.000 22,02 21,20 22,86
2.000 37,03 35,98 38,09
5.000 60,44 58,62 62,23

Para renda de R$ 500, a probabilidade média padronizada foi de 11,69% (IC95%: 10,95–12,48%).

Para renda de R$ 1.000, a probabilidade média padronizada foi de 22,02% (IC95%: 21,20–22,86%).

Para renda de R$ 2.000, a probabilidade média padronizada foi de 37,03% (IC95%: 35,98–38,09%).

O contraste 1000 menos 500 corresponde a 10,33 pontos percentuais (IC95%: 9,95–10,70).

O contraste 2000 menos 1000 corresponde a 15,01 pontos percentuais (IC95%: 14,27–15,76).

As probabilidades mantêm a distribuição das demais covariáveis da amostra analítica, fixando apenas o cenário de renda. São mais diretas para avaliar a magnitude da associação na escala de cobertura. Os intervalos são pontuais; o grupo de renda zero é tratado separadamente da curva de renda positiva.

5 Gráficos da análise

5.1 Cobertura por escolaridade

As barras representam a cobertura ponderada por escolaridade e as hastes indicam os intervalos de confiança. A comparação é descritiva e não isola o efeito da escolaridade.

5.2 Logito empírico e ajuste dos modelos

A proximidade entre os valores observados e os previstos ajuda a avaliar a forma funcional da renda. Diferenças sistemáticas entre os pontos merecem atenção; o gráfico complementa os testes formais.

5.3 Probabilidades médias padronizadas por renda

A curva apresenta as probabilidades médias padronizadas e os intervalos pontuais de 95%. A renda positiva está em escala logarítmica, enquanto a renda zero aparece em painel separado.

5.4 Resíduos de resposta por grupos de renda

Os pontos representam a média de observado menos previsto por grupo de renda. Valores positivos indicam subestimação da cobertura; valores negativos indicam superestimação. Os intervalos são descritivos e não constituem um teste simultâneo de ajuste.

6 Robustez e ajuste

6.0.1 Sensibilidade aos pesos e controles

Modelo N OR Inferior Superior
Principal 90.814 2,2691 2,1859 2,3556
Pesos unitarios 90.814 2,2371 2,1892 2,2861
Pesos aparados P99 90.814 2,2706 2,1931 2,3507
Pos-estratificacao explicita 90.814 2,2691 2,1880 2,3533
  • ocupacao
90.814 2,3204 2,2316 2,4127
  • situacao urbana/rural
90.814 2,2350 2,1524 2,3207

Fonte: elaboração própria com os microdados da PNS 2019, conforme o script original.

Nas especificações apresentadas, a OR da duplicação variou de 2,23 a 2,32.

6.0.2 archer lemeshow

N Grupos GL_desenho Wald F GL1 GL2 p
90.814 10 7.453 15,8462 1,5829 9 7.445 0,114

Fonte: elaboração própria com os microdados da PNS 2019, conforme o script original.

O teste de Archer–Lemeshow não rejeitou a hipótese de médias residuais nulas nos grupos de risco ao nível de 5%. Isso não comprova ajuste perfeito.

A comparação de especificações ajuda a avaliar a estabilidade da associação diante de mudanças nos pesos e nos controles. Estabilidade não elimina confundimento não observado. O diagnóstico de ajuste deve ser interpretado junto da forma funcional e dos resíduos, sem depender de um único valor-p.

7 Interações e ajuste de Holm

7.0.1 Testes de interação e ajuste de Holm

Familia Componente Termo F GL1 GL2 p p_Holm Rejeita_Holm_5pct
regiao Conjunto L e Z regiao 5,4119 8 7.425 < 0,001 < 0,001 TRUE
escolaridade Conjunto L e Z escolaridade 1,4934 12 7.421 0,118 0,355 FALSE
regiao Renda positiva (L) Inclinacao da renda positiva 7,9867 4 7.425 < 0,001 < 0,001 TRUE
regiao Renda zero (Z) Indicador de renda zero 2,6977 4 7.425 0,029 0,116 FALSE
escolaridade Renda positiva (L) Inclinacao da renda positiva 1,5160 6 7.421 0,168 0,355 FALSE
escolaridade Renda zero (Z) Indicador de renda zero 1,4345 6 7.421 0,197 0,355 FALSE

Fonte: elaboração própria com os microdados da PNS 2019, conforme o script original.

Após Holm, 2 dos 6 testes rejeitaram suas hipóteses nulas a 5%.

  • regiao: Conjunto L e Z.
  • regiao: Renda positiva (L).

Holm (1979) controla a probabilidade de pelo menos uma rejeição falsa na família dos seis testes exploratórios. A decisão considera o valor-p ajustado. A interação com região avalia diferenças na associação na escala do logito; ela não equivale, automaticamente, a diferenças de efeitos na escala das probabilidades. A não rejeição de uma interação não demonstra igualdade exata entre grupos.

8 Limitações e conclusão

Os resultados indicam associação positiva entre renda domiciliar per capita e posse de plano médico, mesmo após os controles observados.

Atenção à execução: foram registrados 3 avisos. Examine avisos_analise no ambiente de execução antes de finalizar a interpretação.

A interpretação se restringe à população e ao desfecho definidos no estudo. O delineamento transversal, a análise de casos completos e as covariáveis disponíveis limitam a generalização e a interpretação causal. Os cenários padronizados são comparações estatísticas, não previsões do impacto de uma política de renda.

9 Referências

HOLM, Sture. A simple sequentially rejective multiple test procedure. Scandinavian Journal of Statistics, v. 6, n. 2, p. 65–70, 1979.

LUMLEY, Thomas; SCOTT, Alastair. Fitting regression models to survey data. Statistical Science, v. 32, n. 2, p. 265–278, 2017.