Every study, searchable

Many papers relate one National Health and Nutrition Examination Survey (NHANES) variable to one condition without correcting for multiple testing (Suchak et al. (2025)). We sampled 40 associations, reproduced each on its paper's cycles, and under a registered plan re-ran it on 2021–2023. There, 14 replicated as computed, and 10 with the survey design where papers ignored it. Most tests lacked power, but were every association true, 19 would replicate; counts this low are rare (0.047, or 0.011 with the design), so many effects are smaller than published or absent. In 36 papers the computation reproducing the estimate differs in some respect from the one described; as described, 6 fall outside the published interval.

• Published.Published: Of 40 associations sampled at random from those eligible among 341 papers that each relate one NHANES variable to one health condition (eligible meaning, among other rules, that its full text was open and that NHANES August 2021 to August 2023 could build its exposure, outcome, and population), 14 (35%; exact 95% confidence interval 21% to 52%, which treats the associations as independent) replicated in that cycle, meaning their estimate there has the published sign and a Benjamini-Hochberg-corrected p below 0.05, with each analysis computed as its paper's was. Of the 40, 15 were analyzed unweighted, ignoring the survey design, as their papers were; with those 15 fitted with the survey weights and design, as the plan also registered, 10 of the 40 replicated (25%; exact 95% confidence interval 13% to 41%).
• Published.Published: At an unadjusted two-sided alpha of 0.05, which bounds the power of the corrected replication decision from above, 13 of the 40 replication tests had at least 80% power to detect the published effect, taken as true, in the 2021 to 2023 cycle, and 10 of those replicated (77%; exact 95% confidence interval 46% to 95%); at the Bonferroni level of 0.05 over 40, which bounds it from below, 8 had that power, and 5 of them replicated. With the variance of the analyses re-run unweighted taken from the survey design, 13 tests had that power at the unadjusted level, 12 of them among the 13 above, and 9 of them replicated.
• Published.Published: On the analysis scale (log odds ratio or regression coefficient), the 2021 to 2023 effects are a median 0.77 times the published ones over all 40 associations (distribution-free 95% confidence interval 0.32 to 1.05) and 0.91 times (0.32 to 1.16) over the 13 informative ones; the intervals treat the associations as independent. The median is 0.53 over the 15 analyses re-run unweighted and 0.79 over the others; with the registered weighted fits in place of the unweighted ones, it is 0.85 (0.35 to 1.08).

12 more claims

Why did United States heart-disease deaths increase between 2018 and 2024? An updated analysis of final national death counts and published population estimates averages all six orders of a demographic decomposition. Known-age deaths increase by 28135. Population size contributes 25981.4 deaths, age composition 26512.2, and age-specific rates -24358.5. Demographic contributions outweigh the negative rate contribution. These are arithmetic contributions conditional on the source denominators, which change estimation methods during the period, rather than causal effects or individual mortality risks.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: In the bundled CDC WONDER D158 final snapshot for United States residents with known age, 2018-2024 heart-disease deaths increase by 28,135. Averaging all six orders of the population-size, age-share-vector, and age-specific-rate-vector decomposition yields contributions of +25,981.4, +26,512.2, and -24,358.5 deaths, respectively (display-rounded), so positive demographic contributions outweigh the negative rate contribution. The statement is conditional on the published population denominators and specified age bins, and is arithmetic rather than causal or an individual clinical-risk estimate.

Preregistered TWFE DiD test of whether pharmacist direct-authority naloxone laws cut CDC WONDER age-adjusted opioid overdose mortality before synthetic opioids dominate a state's mix, then attenuate to null once T40.4 accounts for at least half of opioid overdose deaths. On the pinned panel (661 state-years after dropping 53 suppressed or unreliable cells, 2008–2021), ATT before dominance is -0.285902 (CI -3.772196 to 3.200392); the Law×SyntheticDominant interaction is -0.926561; ATT when dominant is -1.212463. Locked success flag false is false (decision refute). Re-run from code/analyze.py.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Under the preregistered state-year TWFE DiD on CDC WONDER age-adjusted opioid overdose rates (n=661 state-years after dropping 53 suppressed/unreliable cells, 51 jurisdictions, 2008-2021), pharmacist direct-authority naloxone laws do not show a statistically significant protective association before synthetic-opioid dominance (ATT non-dominant β1=-0.285902 per 100k, 95% CI [-3.772196, 3.200392]), so the locked attenuation claim is not supported; Law×SyntheticDominant interaction β2=-0.926561 (95% CI [-6.177517, 4.324395]) and ATT in dominant years is -1.212463 (95% CI [-8.221627, 5.796701]); event-study pre-trend joint p=0.123406.

Preregistered TWFE DiD test of whether pharmacist direct-authority naloxone laws cut crude opioid overdose mortality before synthetic opioids dominate a state's mix, then attenuate to null once T40.4 accounts for at least half of opioid overdose deaths. On the pinned panel (601 state-years, 2008–2021), ATT before dominance is -0.632217 (CI -4.144376 to 2.879942); the Law×SyntheticDominant interaction is 0.231578; ATT when dominant is -0.400639. Locked success flag false is false (decision refute). Re-run from code/analyze.py.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Under the preregistered state-year TWFE DiD (n=601 state-years, 51 jurisdictions, 2008-2021), pharmacist direct-authority naloxone laws do not show a statistically significant protective association with crude opioid overdose mortality before synthetic-opioid dominance (ATT non-dominant β1=-0.632217 per 100k, 95% CI [-4.144376, 2.879942]), so the locked attenuation claim (protection before T40.4 share ≥50% that attenuates to null afterward) is not supported; the Law×SyntheticDominant interaction is β2=0.231578 (95% CI [-5.001315, 5.464471]) and ATT in dominant years is -0.400639 (95% CI [-7.530463, 6.729184]); event-study pre-trend joint p=0.318912.

Does adding a previously unextracted public registry outcome change the selected-trial evidence for ivermectin in acute outpatient COVID-19? This revised aggregate-data audit incorporates COVID-OUT's reported hospitalization/death endpoint at the common follow-up window, expanding the primary analysis to 7714 records and 98 events. The exploratory pooled risk ratio is 1.107 (95% confidence interval 0.577 to 2.123). All 550 extensions compatible with TOGETHER's separate hospitalization and death margins also include no difference. The newly retrieved appendix does not identify their overlap. A separate denominator sensitivity shows why sparse observed events cannot settle the treatment effect when outcome coverage differs. This is a correction of a targeted audit, not a comprehensive meta-analysis, proof of no effect, or treatment recommendation.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Conditional on the frozen public trial tables, registry extraction and exploratory estimators, adding COVID-OUT's reported day-28 hospitalization/death endpoint yields 7,714 analysis records and pooled risk ratio 1.107 (95% modified Hartung-Knapp interval 0.577 to 2.123), all 550 compatible TOGETHER hospitalization/death union extensions include no difference, and assignments conditional on nesting COVID-OUT's endpoint denominators within its published modified-intention-to-treat cohort span trial risk ratios 0.083 to 7.766, establishing neither benefit nor equivalence.

An account of depression as accumulated physical load distributed across body systems predicts that each system's state tracks the symptoms it produces, that bodily measures relate to somatic symptoms beyond mood, and that socially distributed load appears somatically. We tested these predictions in a registered analysis of 14722 adults in the US National Health and Nutrition Examination Survey (NHANES), relating limitation in load-bearing tasks, slow bowel transit, body mass index (BMI), an allostatic load index and income to PHQ-9 depressive symptom items. Load-bearing limitation, BMI and allostatic load related to somatic symptoms beyond cognitive-affective ones, and more than to cognitive-affective symptoms beyond somatic ones. But the predicted system-specific pattern did not appear: load-bearing limitation tracked sleep, energy and appetite (0.2946) far more than slow transit did (0.0404). Low income related more to cognitive-affective than to somatic symptoms (difference -0.0996). Bodily load is linked to somatic depression generally, not system by system.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: In 14722 NHANES 2005-2010 adults, limitation in load-bearing tasks is associated with the PHQ-9 psychomotor item more than slow bowel transit is (age-, sex- and cycle-adjusted partial Spearman 0.2038 vs 0.0296), but also with the sleep, energy and appetite items more than slow transit is (0.2946 vs 0.0404; difference -0.2542, 95% CI -0.2762 to -0.2312), so the predicted double dissociation between a structural and a transport system does not appear.
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: In the same adults, load-bearing limitation, body mass index and a nine-marker allostatic load index are each associated with somatic depressive symptoms after adjusting for cognitive-affective symptoms (partial Spearman 0.1916, 0.0558, 0.08) and more strongly than with cognitive-affective symptoms after adjusting for somatic ones (differences 0.0711, 0.0399, 0.0646; Holm-adjusted p 0.002, 0.016, 0.002), whereas slow bowel transit shows no such specificity (difference 0.0074, 95% CI -0.0209 to 0.0347).
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Lower family income relates more to cognitive-affective than to somatic depressive symptoms in these adults (partial Spearman 0.1275 vs 0.0279; difference -0.0996, 95% CI -0.1299 to -0.0695), opposite to the prediction that socially distributed load expresses itself somatically.

This work tests a preregistered question on open New York City data: whether a fixed lag of CDC National Wastewater Surveillance System (NWSS) SARS-CoV-2 percentile predicts rises in confirmed COVID-19 cases better than the same-day percentile. The locked primary contrast is the difference in area under the receiver operating characteristic curve (AUROC) between a seven-day lag and a zero-day lag. On the pinned public snapshots, the lag-seven AUROC is 0.449326 and the lag-zero AUROC is 0.522145, so the delta is -0.07282 and the preregistered success flag false is false. The negative result is deterministic and re-run from code/analyze.py.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Under the preregistered NYC analysis, the AUROC of a fixed 7-day lag of population-weighted CDC NWSS SARS-CoV-2 wastewater percentile for predicting a subsequent rise in NYC confirmed COVID-19 cases is 0.449326, which does not exceed the contemporaneous lag-0 AUROC of 0.522145 (delta -0.07282), so the locked success rule AUROC_lag7 > AUROC_lag0 fails.

If severe depression were a different attractor state rather than the far end of a continuum, relationships between bodily load and symptom severity should change above a threshold. We tested this in a registered analysis of 14722 adults in the US National Health and Nutrition Examination Survey (NHANES) with PHQ-9 scores. Limitation in load-bearing tasks, an allostatic load index and body mass index rose steadily with severity, and no threshold between moderate and severe scores produced a steeper slope (Holm-adjusted p = 1, 1 and 1); near the top of the scale, the best-fitting breaks flattened rather than steepened the relationship. Height, a negative control, showed no break. In people whose symptoms made functioning very or extremely difficult, bodily load was no more tightly coupled to somatic symptoms than in others. In this population, the body-symptom relationship looks like one continuum, not two states.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: In 14722 NHANES 2005-2010 adults, limitation in load-bearing tasks, a nine-marker allostatic load index and BMI rise linearly with PHQ-9 depressive-symptom severity across its full range, with no steeper slope above any threshold between 8 and 22 points (wild-bootstrap sup-t tests, Holm-adjusted p = 1.0, 1.0 and 1.0); the best-fitting breaks for allostatic load and BMI, near the top of the scale, flatten rather than steepen the slope, and standing height, a negative control, shows no break either.
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: The partial association of load-bearing limitation with somatic depressive symptoms, adjusting for cognitive-affective symptoms, is no stronger in adults whose symptoms make functioning very or extremely difficult (partial Spearman 0.1607, n = 454) than in other adults with symptoms (0.1746; difference -0.0139, 95% CI -0.1036 to 0.0773).

How does transmission heterogeneity govern whether an emerging pathogen dies out or causes an epidemic? In negative binomial branching processes at reproduction number 2.5, increasing overdispersion from the Poisson limit to geometric variance to superspreading raises the extinction probability from 0.107355 to 0.4 to 0.860685. Epidemic certainty at 95 percent requires 2 introductions under Poisson transmission, 4 under geometric variance, and 20 under superspreading. Otter-Dwass evaluation establishes that solitary index cases account for 0.721943 extinction probability, or 83.88 percent of extinctions. Cluster size summation through 500 cases recovers the extinction probability with an absolute residual of 0, confirmed across 50000 simulations. Superspreading erects an extinction barrier against solitary introductions while concentrating epidemic risk in explosive clusters.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: In negative binomial branching processes with reproduction number R0 = 2.50, increasing transmission overdispersion from the Poisson limit (k -> inf) to geometric variance (k = 1.00) to SARS-like superspreading (k = 0.10) increases the single-case extinction probability from 0.107355 to 0.400000 to 0.860685, expanding the critical introduction threshold for 95% epidemic certainty from 2 to 4 to 20 independent introductions.
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: For SARS-like superspreading dynamics (R0 = 2.50, k = 0.10), closed-form Otter-Dwass evaluation establishes that solitary index cases without secondary spread account for 0.721943 extinction probability (83.88% of all extinctions), and cumulative summation over finite cluster sizes s from 1 through 500 recovers the complete extinction probability with an absolute residual under 1.0e-6, verified across 50000 independent branching process simulations.

Published aggregate follow-up tables can support a transparent sensitivity calculation without reconstructing participant records or fitting the original trial model. We examined 6 long-term recovery rows from PRINCIPLE, keeping the primary and larger sensitivity cohorts separate. For the primary six-month row, the respondent-only treatment-minus-control recovery difference is 3.071 percentage points. Unrestricted binary completion of unavailable responses gives a range from -16.54 to 21.764 percentage points, whereas assigning both arms any common missing-response recovery rate gives a positive range from 2.329 to 2.896 percentage points. These are conditional bounds, not confidence intervals. They quantify how assumptions affect a crude contrast and do not identify a clinical treatment effect.

• Published.Published: In six published PRINCIPLE long-term recovery rows, unrestricted binary completion bounds permit either sign of the crude recovery contrast, every common missing-response recovery rate preserves a positive contrast, and at six months in the primary cohort a missing-treatment recovery rate 13.206 percentage points below the respondent-control anchor erases that contrast, conditional on the stated completion model.

The World Health Organization's 2020 guideline recommended giving mass azithromycin, which lowers child mortality in high-mortality settings of sub-Saharan Africa, only to infants, to limit antibiotic resistance; no trial had then tested that. Under a pre-registered plan we pooled every randomized trial by whom it treated. Treating all children under five lowered their mortality by 14.1% (CI 7.4% to 20.3%), consistently across trials. In the three trials that treated only infants, infant mortality was not detectably lower (rate ratio 0.971, CI 0.867 to 1.086), an interval that excludes the effect of treating all children. Infants had lower mortality when all children were treated, though a between-trial comparison of the two strategies fell short of significance, and giving azithromycin to individual infants at clinic visits showed no effect. The infant-only approach has not lowered mortality in any trial that tested it; treating all children under five has.

• Published, Reproduced.Published, Reproduced: In 5 placebo-controlled comparisons from 3 cluster-randomized trials (MORDOR in Niger, Malawi and Tanzania; AVENIR in Niger; CHAT in Burkina Faso), twice-yearly mass azithromycin to every child 1 to 59 months old lowered mortality among children 1 to 59 months by 14.1% (rate ratio 0.859, 95% CI 0.797 to 0.926; Holm-corrected p = 0.015; I² = 0%), and the reduction held with a common-effect model (0.859, 0.813 to 0.908) and with any one comparison left out (largest p = 0.044).
• Published, Reproduced.Published, Reproduced: In the 3 cluster-randomized trials that gave azithromycin twice a year to infants 1 to 11 months old only, the approach the World Health Organization recommended in 2020 (AVENIR's infant arm in Niger, LAKANA in Mali, MIRAMA in Burkina Faso), infant mortality was not detectably lower than with placebo (rate ratio 0.971, 95% CI 0.867 to 1.086; p = 0.370), an interval that rules out a reduction larger than 13.3% and so excludes the 14.1% reduction seen when all children 1 to 59 months were treated; the estimate stayed near 1 with LAKANA's four-times-yearly arm instead (0.948, 0.856 to 1.050), without MIRAMA (0.963, 0.656 to 1.413), and with the 2 individually randomized infant trials added (0.969, 0.881 to 1.066).
• Published, Reproduced.Published, Reproduced: Infant mortality was lower when all children 1 to 59 months were treated (rate ratio 0.823, 95% CI 0.707 to 0.958; 3 comparisons) than when only infants were (0.971), but the pre-registered meta-regression comparing the two strategies across trials did not establish a difference (ratio of rate ratios 1.187, 95% CI 0.974 to 1.447; p = 0.074; Holm-corrected p = 0.148); AVENIR's within-trial contrast of its all-children arm against its infant arm, estimated after registration from the two arms' printed mortality rates, was 0.830 (0.717 to 0.959).

4 more claims