Altas habilidades · neurodivergência · saúde
Auditoria de perdas, pesos e variância do ciclo 001
2026-09-09 · NHANES 2021–2023: 6.337 adultos após junções por hipótese, 4.194–4.522 casos completos de 20–80 anos; HA/autismo/TDAH não avaliados.
Reconstituição das amostras e comparação Taylor com svy: erros-padrão reproduzidos e perdas após junções. Os ciclos 003–004 registram atualizações posteriores de mensuração e inferência.
Escopo e registro
Auditoria pós-resultados do ciclo 001 corrigido, registrada localmente em 2026-09-09 antes da execução desta rodada. Usa os mesmos seis XPT, participantes e especificações. Não é replicação, nova evidência de eficácia ou validação em outra população.
Foram reproduzidas as chaves dos participantes analíticos e todas as estimativas do ciclo anterior. Os hashes dos insumos permanecem iguais. O pacote svy 0.28.0 foi instalado como dependência de desenvolvimento pelo Pipenv; o lock registra suas dependências.
Cobertura e seleção
DEMO contém 11.933 pessoas, sendo 8.153 com idade ≥18. As ligações necessárias a cada pergunta deixam 6.337 adultos. A tabela usa este último denominador: as retenções não representam participação no inquérito ou na população inteira.
Os 273 participantes de 18–19 anos nas junções têm DMDEDUC2 ausente por elegibilidade da variável (20+). Portanto, a amostra final 20–80 não inclui os mais jovens do protocolo 18+. Ausência estrutural não é recusa. A idade máxima também é truncada na divulgação.
A renda/PIR falta em 831 dos 6.337 adultos. H1/H3 têm 1.351 recordatórios indisponíveis e pesos dietéticos zero. Essas perdas se sobrepõem: contagens marginais não podem ser somadas. Atribuições sequenciais dependem da ordem dos filtros e não revelam causas clínicas de ausência.
| Modelo | Após junções | Analisados | Retenção | Excluídos após junções |
|---|---|---|---|---|
| H1 | 6.337 | 4.211 | 66,45% | 2.126 |
| H2 | 6.337 | 4.522 | 71,36% | 1.815 |
| H3 | 6.337 | 4.194 | 66,18% | 2.143 |
Desigualdade dos pesos
Kish = (soma dos pesos)²/soma dos pesos² descreve apenas a desigualdade de ponderação. Não incorpora clusters, estratificação, desfecho ou covariáveis e não deve ser interpretado como tamanho efetivo completo da regressão.
Não foram truncados, recalibrados ou substituídos pesos. As somas após seleção não recuperam automaticamente a população-alvo: o ajuste original de não resposta do inquérito não elimina seleção adicional por casos completos.
| Modelo | Peso | N | Kish | CV | Fração dos pesos no 1% maior |
|---|---|---|---|---|---|
| H1 | WTDRD1 | 4.211 | 2460,2 | 0,844 | 5,45% |
| H2 | WTMEC2YR | 4.522 | 2957,1 | 0,727 | 4,39% |
| H3 | WTDRD1 | 4.194 | 2450,7 | 0,843 | 5,36% |
Variância: comparação independente executada
A matriz de covariáveis foi compartilhada para isolar o cálculo estatístico. Em svy, o desenho foi definido com todos os registros de peso positivo, estrato e PSU antes de selecionar o domínio de casos completos com where. Foi usada aproximação com reposição, sem correção de população finita.
Todos os 15 estratos e 30 PSUs do desenho têm participantes nas três amostras analíticas. Nesse caso concreto, o cálculo Taylor anterior coincide numericamente com o cálculo que preserva o desenho completo. Isso não autoriza reutilizar a rotina anterior em domínios que esvaziem PSUs: ela ignora estratos com menos de duas PSUs observadas.
A comparação verifica cada coeficiente e seu erro-padrão, não a matriz completa de covariância ou testes conjuntos. Não valida mensuração, escolha de confundidores, seleção de casos completos, causalidade ou transporte para HA/autismo/TDAH.
| Modelo | N no desenho completo | Colunas/rank da matriz | Maior diferença absoluta: coeficiente | Maior diferença relativa: erro-padrão |
|---|---|---|---|---|
| H1 | 6.754 | 17/17 | 5e-13 | 4.64e-14 |
| H2 | 8.860 | 17/17 | 1.13e-13 | 5.48e-14 |
| H3 | 6.754 | 18/18 | 2.16e-13 | 1.73e-13 |
Sensibilidade aos graus de liberdade
O ciclo 001 usa t com 15 graus (30 PSUs − 15 estratos). O padrão de svy 0.28.0 subtrai k−1 e aplica piso de 1: max(1, 15−(k−1)). Com k=17 em H1/H2 e k=18 em H3, todos chegam a 1 grau. Esse piso é comportamento do pacote, não demonstração de que 1 seja a escolha correta para estes modelos.
Os coeficientes e erros-padrão são os mesmos, mas p-valores e intervalos dependem da distribuição de referência. A documentação de R survey descreve tanto a convenção residual quanto a opção de usar os graus do desenho; alerta que a primeira pode ser muito conservadora para covariáveis individuais. Não foi executado R nesta rodada.
Recalculando p a partir de svy com 15 graus, reproduzimos o ciclo 001 e a mesma família FDR de oito termos. Com o padrão de 1 grau, nenhum teste tem q<0,05. Essa comparação é sensibilidade à convenção, não oito novas descobertas. Nenhuma das duas saídas substitui silenciosamente a outra; futuras especificações devem justificar complexidade e graus antes dos resultados.
| Modelo/termo | p: 15 graus | p: svy padrão (1) | q: 15 graus | q: svy padrão (1) |
|---|---|---|---|---|
| H1 / CAFF100 | 0.00362143 | 0.179938 | 0.00724286 | 0.359875 |
| H1 / CAFF_X_AGE | 0.0352676 | 0.259702 | 0.0564281 | 0.415523 |
| H2 / C(SLEEP_CAT)[T.short] | 2.72277e-05 | 0.106218 | 0.000108911 | 0.359875 |
| H2 / C(SLEEP_CAT)[T.long] | 0.00126081 | 0.157518 | 0.00336216 | 0.359875 |
| H2 / SED_HOURS | 6.50067e-07 | 0.0774049 | 5.20054e-06 | 0.359875 |
| H3 / DIET_MAG100 | 0.463415 | 0.589329 | 0.463415 | 0.589329 |
| H3 / SUPP_MAG100 | 0.394445 | 0.541747 | 0.463415 | 0.589329 |
| H3 / SUPP_MAG_UNQUANT | 0.426915 | 0.56404 | 0.463415 | 0.589329 |
Magnésio: limites da ausência registrada
Na amostra H3, 1.401 pessoas declaram uso de algum suplemento, mas DSQTMAGN está ausente. Isso não comprova que tomaram magnésio com dose desconhecida: pode envolver produto sem magnésio ou incompletude da composição. O nome operacional SUPP_MAG_UNQUANT do ciclo anterior não resolve essa ambiguidade.
Há 28 pessoas com DSD010=2 e magnésio quantificado; suplementos e antiácidos têm indicadores distintos. Não concluir inconsistência nem zerar automaticamente esses valores. A recodificação e a distinção entre não exposição e quantidade desconhecida exigem detalhamento de produtos/antiácidos em uma próxima auditoria.
O ciclo 001 permanece numericamente preservado. H3 é estimável com limites de mensuração; essa checagem não confirma forma química, uso de magnésio por todos os usuários de suplemento, ou indicação clínica.
| Categoria observada em H3 | N |
|---|---|
| use_yes_mag_missing | 1401 |
| use_no_mag_missing | 1445 |
| use_no_mag_reported | 28 |
| mag_reported | 1348 |
Conferência de leitura e códigos
Foram conferidos totais do codebook DSQTOT_L: 3.762 usuários, 2.969 não usuários, 1.613 valores presentes de magnésio e 2.106 pesos WTDRD1 zero. Os 2.106 zeros de peso estavam representados como floats SAS diminutos e foram normalizados pela regra existente. Essa checagem cobre esses campos, não todos os zeros de todos os componentes.
Decisão e próximos testes
H1, H2 e H3: estimable_with_limits para associações na amostra de casos completos, com estimadores e erros-padrão conferidos. A robustez da inferência depende da escolha de graus, da seleção e da mensuração. Os grupos HA, autismo e TDAH continuam not_assessed.
Estado na data desta auditoria: definir complexidade e graus de liberdade antes de novas especificações; auditar composição de suplementos/antiácidos e códigos de ausência; avaliar seleção/dados ausentes; só então registrar formas não lineares e validação por PSU/onda. Não executar uma busca de especificações que produzam significância.
Atualização posterior: o ciclo 003 auditou produtos e antiácidos; o ciclo 004 registrou 15 graus do desenho como convenção principal segundo o NCHS e testou reponderação pela seleção observável. As pendências atuais são a nova especificação de H3, seleção não observada, formas não lineares e validação independente.
Resultados
| Pergunta | Estimativa | Incerteza | Amostra | Interpretação |
|---|---|---|---|---|
| Quanto se perde após as junções por hipótese? | H1: 4.211/6.337 (66,45%); H2: 4.522/6.337 (71,36%); H3: 4.194/6.337 (66,18%) | Contagens exatas dos arquivos locais; não são estimativas de prevalência e não requerem IC. | 6.337 adultos após junções; antes delas, DEMO contém 8.153 adultos. | 273 pessoas de 18–19 anos não têm a escolaridade 20+ exigida. Ausência de renda, dieta e sintomas explica outras perdas observáveis, com sobreposição. |
| Qual é a concentração dos pesos analíticos? | Kish: H1=2.460,2; H2=2.957,1; H3=2.450,7 | Indicadores descritivos dos pesos; não são N efetivos completos do desenho. | WTDRD1 em H1/H3; WTMEC2YR em H2; sem truncamento de pesos. | O 1% de maiores pesos concentra aproximadamente 4,4–5,4% da soma; seleção por casos completos permanece uma limitação. |
| A implementação anterior reproduz a variância de um pacote independente? | Diferença relativa máxima nos erros-padrão <2×10⁻¹³; coeficientes diferem <6×10⁻¹³ | Concordância numérica por coeficiente; não é IC de efeito nem validação causal. | svy 0.28.0, domínio de casos completos no desenho completo; 15 estratos/30 PSUs representados em cada modelo. | Coeficientes e erros-padrão conferidos neste domínio; a rotina anterior não deve ser generalizada a domínios com PSUs vazias. |
| Os testes dependem da convenção de graus de liberdade? | 15 graus no ciclo 001; 1 no padrão svy. Com o padrão svy, nenhum dos oito q-valores é <0,05. | Mesmos coeficientes e erros-padrão; distribuições t diferentes. Com 15 graus, svy reproduz os p/q originais. | 17 coeficientes em H1/H2 e 18 em H3; svy aplica max(1, 15−(k−1)). | Relatar as duas convenções. O piso de 1 não prova superioridade; definir complexidade e inferência antes de novas análises. |
Gráficos


Limitações
- Auditoria na mesma onda e nos mesmos participantes, sem replicação independente.
- Inferência sensível à convenção de graus; modelos têm 17–18 coeficientes para 15 graus do desenho.
- Casos completos podem introduzir seleção não resolvida pelos pesos originais.
- DSQTMAGN ausente não identifica por si só consumo de magnésio com dose desconhecida.
- Não identifica altas habilidades, autismo ou TDAH; não orienta dose individual.
Fontes
- CDC: variância e subpopulações NHANES
- CDC: ponderação NHANES
- DEMO_L: escolaridade de adultos 20+
- DSQTOT_L: suplementos e antiácidos
- svy: regressão e análise de domínio
- R survey: convenções de graus de liberdade
Reprodução
{
"command": "MPLCONFIGDIR=/tmp/science-matplotlib PIPENV_VENV_IN_PROJECT=1 pipenv run python research/discoveries/cycle-002/run_audit.py",
"code_revision": "working tree; run_audit.py sha256=f190c0e3a52d60fe23e7e398bcf3dd8e81e06ee540b4425fe0ee5a566f9d2aef; svy 0.28.0; additional hashes in results.json",
"input_sha256": {
"DEMO_L.xpt": "ca4374a158b493b8b0163e1388da21d57a18d1b9cecff2aa4e2fa2bec494fe23",
"DPQ_L.xpt": "25605a02685035fe997477b31a0991cca2776b0f72f24a8e61ac5b018456304b",
"DR1TOT_L.xpt": "73a3097aab64000f9c1d138367363ca605968b691e349981132bd6032ca6c9ef",
"DSQTOT_L.xpt": "8ee3e55578f56918a190100bb321a0b580f824ca3cc439901f2fd707dbe695f8",
"PAQ_L.xpt": "de34acfed4523f10f52bea06e64ed33c8db973cb38d3539eeceb2b68039248c8",
"SLQ_L.xpt": "918c3258b2f466ac1cae55ea45330c736dbdc573ee4dc72306e3d2b222764741"
}
}