Glauber Machado Rodrigues

gmrodrigues / pesquisa aberta

Altas habilidades, neurodivergência e saúde

Perguntas, dados e descobertas em construção.

← Site principal · Forte e Leve

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.

ModeloApós junçõesAnalisadosRetençãoExcluídos após junções
H16.3374.21166,45%2.126
H26.3374.52271,36%1.815
H36.3374.19466,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.

ModeloPesoNKishCVFração dos pesos no 1% maior
H1WTDRD14.2112460,20,8445,45%
H2WTMEC2YR4.5222957,10,7274,39%
H3WTDRD14.1942450,70,8435,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.

ModeloN no desenho completoColunas/rank da matrizMaior diferença absoluta: coeficienteMaior diferença relativa: erro-padrão
H16.75417/175e-134.64e-14
H28.86017/171.13e-135.48e-14
H36.75418/182.16e-131.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/termop: 15 grausp: svy padrão (1)q: 15 grausq: svy padrão (1)
H1 / CAFF1000.003621430.1799380.007242860.359875
H1 / CAFF_X_AGE0.03526760.2597020.05642810.415523
H2 / C(SLEEP_CAT)[T.short]2.72277e-050.1062180.0001089110.359875
H2 / C(SLEEP_CAT)[T.long]0.001260810.1575180.003362160.359875
H2 / SED_HOURS6.50067e-070.07740495.20054e-060.359875
H3 / DIET_MAG1000.4634150.5893290.4634150.589329
H3 / SUPP_MAG1000.3944450.5417470.4634150.589329
H3 / SUPP_MAG_UNQUANT0.4269150.564040.4634150.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 H3N
use_yes_mag_missing1401
use_no_mag_missing1445
use_no_mag_reported28
mag_reported1348

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

PerguntaEstimativaIncertezaAmostraInterpretaçã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,7Indicadores 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

Participantes retidos e excluídos nas três hipóteses
NHANES 2021–2023: contagens não ponderadas após junções de adultos; não representa taxa populacional de resposta.
Comparação entre N observado e indicador de Kish
WTDRD1 para H1/H3, WTMEC2YR para H2; Kish reflete somente desigualdade dos pesos, sem IC ou correção por clusters.

Limitações

Fontes

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"
  }
}

Dados agregados do ciclo (JSON)