State-Level Youth Mental Health Need-vs-Access Measures Correlate with Within-Instrument Suicide Indicators and Show Modest Cross-Instrument Coupling: A Phase B Outcomes Correlation Analysis of the Trellison Need-vs-Access Framework
Authors: Rob Stillwell, DaedArch Corporation · Trellison Institute · DaedArch AI, DaedArch Platform
Corresponding author: [email protected]
Affiliations: DaedArch Corporation · Trellison Institute · LedgerWell Corporation
Version: Working paper · May 2026
Phase: B — outcomes correlation,
Methodology rating: Trellison Institute — passes 7 gates including external_validity_v1, location_authority_v1
Evidence-chain certificate: LedgerWell Corporation — verified against HRSA HPSA-MH (BCD_HPSA_FCT_DET_MH.csv) and AAMC-equivalent C&A psychiatrist rate computed independently from CMS NPPES May 2026 bulk
License: CC-BY-4.0
Working paper hub: https://trellison.com/research/youth-mental-health-outcomes-correlation-
Abstract
Background. The Trellison Need-vs-Access Framework identifies state-level workforce-capacity gaps in mental-health services relative to the prevailing insurance-coverage landscape. The May 2026 Youth Mental Health Access Gap V1 working paper applied the framework to under-18 populations across 35 U.S. states. The framework's analytical outputs have been treated by some readers as if they were predictors of poor mental-health outcomes. This Phase B analysis tests that interpretation directly, with two methodology corrections applied (state_abbr→practice_state denominator fix; NUCC-taxonomy-whitelist provider filter replacing a broken youth_specific flag).
Methods. We joined the 35-state gap dataset to three categories of state-level outcome data: (i) CDC Youth Risk Behavior Surveillance System (YRBSS) 2023 state-level self-reported suicide indicators among grades 9-12 (considered, planned, attempted in past 12 months); (ii) NCHS Underlying Cause of Death state-level age-adjusted suicide death rates (AADR) for 2017; (iii) NCHS Vital Statistics Rapid Release state-level drug overdose death counts (2023, all-age, 12-month-ending) converted to per-100K-population rates using ACS 1-year 2023 state denominators. We computed Pearson product-moment correlation coefficients between each of the five framework gap measures (need_value, covariate_value, access_value, log_gap_ratio, residual_z) and each of the five outcome measures, and we profiled the framework-identified outlier states alongside seven pre-specified comparator states.
Results. Within-YRBSS coupling is strong: the principal need metric (sad/hopeless prevalence) correlates with state-level YRBSS suicide indicators at r = +0.817 (considered, n=34), +0.800, +0.762 (attempted, n=33). Across-instrument coupling is moderate: need × all-age suicide AADR r = +0.338 (n=34); need × drug overdose r = +0.192 (n=35). The supply-side measure shows weak-but-protective correlation with YRBSS suicide indicators (r = -0.300 to -0.233) — a substantive direction-and-magnitude difference from the paper's near-zero supply-side correlations, which were artifacts of two methodology issues (Section 0). The framework's residual-classification signal correlates at r = -0.356 to -0.334 with YRBSS suicide indicators (weak-but-protective at state-level cross-section).
outlier classifications: positive outlier (best access) = VT (Vermont); negative outliers (worst access) = MS, PR, TN, WV (West Virginia, Puerto Rico, Tennessee, Mississippi). HRSA HPSA-MH cross-validation confirms the direction for all five outlier classifications. 's PR/NC/NJ-positive-outlier and VT/AK-negative-outlier claims were artifacts of the methodology issues; the classifications align with both HRSA HPSA-MH severity rankings and within-instrument suicide outcomes.
Conclusions. The framework's gap measures are valid cross-sectional proxies for state-level youth suicide ideation/attempt prevalence within the YRBSS instrument, and the supply-side measure shows weak-but-protective state-level correlation that was absent under 's broken methodology. The orthogonality finding in reflected the methodology issues, not a substantive framework property. The framework's residual classification, after, stratifies youth-suicide outcomes at state-level cross-section at r ≈ -0.3 (weak but meaningful). Drug overdose mortality has a state signature distinct from suicide outcomes and requires a separate framework binding. We publish the full joined dataset, per-state outcome profiles, methodology supplement, external-validity audit, sensitivity analysis, and replication recipe under CC-BY-4.0.
Keywords: state-level analysis, youth mental health, need-vs-access, suicide surveillance, drug overdose, provider supply, methodology audit, methodology, NUCC taxonomy, NPPES bulk file, Trellison Institute.
1. Introduction
This section locates the Phase B working paper within the Trellison Need-vs-Access Framework series, summarizes the methodology corrections, and lays out the questions the analysis answers.
1.1 The framework and its prior publication; the paper's methodology errors
The Need-vs-Access Framework was published in two prior Trellison Institute working papers in May 2026: an adult census-tract-level analysis [3] and a youth state-level analysis [2]. The release of Phase B: Outcomes Correlation (this paper's predecessor) tested whether the framework's residual-classification outputs predict adverse state-level mental-health outcomes. reported a non-protective supply signal and orthogonal residual signal at state-level cross-section.
In the weeks following, two methodology issues were identified by the Trellison Institute's own technical review stack (gate.location_authority_v1 and gate.external_validity_v1). Both are documented in Section 0 of this release. The corrections substantively change the outlier-classification conclusions and the central interpretive finding. The framework's contribution is preserved; the bounded interpretation is sharpened. The public-facing dartboard videos and supporting materials remain private until the supporting materials complete gate clearance.
1.2 The Phase B question
Multiple readers of asked: do the framework's gap-classification outputs predict adverse mental-health outcomes at state level? Specifically: do positive-outlier states (where the framework flags an under-built youth-serving workforce relative to the insurance landscape) show systematically worse state-level outcomes — higher YRBSS suicide ideation, higher attempted-suicide rates, higher state-level mortality — and do negative-outlier states show the inverse?
This question implies a specific interpretation of the framework's outputs: that the residual-classification signal operates as a state-level mental-health-outcome risk predictor. The analysis tests this interpretation rigorously, with both methodology issues fixed and the external-validity cross-checks now in place.
1.3 Why this matters
The Need-vs-Access Framework's residual analysis was designed for a specific question: which states have built or failed to build youth-serving workforce capacity relative to what their insurance-coverage landscape would predict? That question is policy-actionable for workforce-expansion programs (BHWET, CCBHC certification, Title V Maternal-Child Health Block Grants, state-level loan-repayment incentives for community youth mental-health workforce) [4]. Whether it is also policy-actionable for suicide-prevention deployment (988 Suicide and Crisis Lifeline capacity expansion, mobile rapid-response teams, firearm-access counseling, NHSC loan-repayment for rural rapid-response staff) [5] depends on whether the residual classification stratifies state-level outcomes. finds that it weakly does.
1.4 Contribution
This paper contributes (i) a state-level cross-sectional correlation analysis of the framework's gap measures against three categories of outcome data; (ii) the demonstration that the framework's need-metric is a strong cross-sectional proxy for within-instrument youth suicide indicators (r > +0.75); (iii) the demonstration that the supply-side metric shows weak-but-protective correlation with within-instrument youth suicide outcomes (r ≈ -0.2 to -0.3), reversing 's incorrect orthogonality claim; (iv) the demonstration that the residual-classification signal weakly stratifies state-level youth-suicide outcomes (r ≈ -0.3 for YRBSS attempted-suicide rate); (v) the demonstration that the same residual signal is essentially uncorrelated with all-age mortality outcomes (the cross-instrument signal-strength degradation is real); (vi) a clarifying policy interpretation that the framework's outputs are properly understood as a workforce-build-out priority list with weak cross-sectional suicide-outcome stratification; and (vii) an external-validity audit cross-checking outputs against HRSA HPSA-MH designations and AAMC-equivalent narrow C&A psychiatrist density (Section 8). The published dataset, methodology supplement, and reproduction recipe support independent verification.
1.5 Roadmap
Section 0 documents the → methodology corrections. Section 2 describes the data sources and their provenance. Section 3 details the joining, correlation, and methodology- methodology. Section 4 presents the results, organized into five subsections by outcome category. Section 5 discusses interpretation, including revised explanations for the supply signal and the policy implications. Section 6 lists ten specific limitations (eight from plus two new ones surfaced by cross-validation). Section 7 concludes. Section 8 (new in ) documents the external-validity audit.
2. Data sources
This section catalogs the input collections (framework gap measures, YRBSS suicide indicators, NCHS AADR, VSRR drug-OD) plus the external-validity benchmarks newly sourced for (HRSA HPSA-MH, AAMC-equivalent narrow C&A psychiatrist density). Each source is cited with retrieval date and access URL.
2.1 Gap measures (input)
Source collection: analysis_outputs.mh_gap_youth_v3_state_corrected_v1 (35 rows, 35 U.S. states with complete YRBSS 2023 Total-demographic + ACS + NPPES coverage). Each row contains:
need_value: state-level YRBSS sad/hopeless 2+ wks prevalence (grades 9-12, Total demographic, 2023)covariate_value: state-level ACS 1-year 2023 under-19 uninsured percentageaccess_value: state-level youth-serving NPI count (NUCC taxonomy whitelist) per 100K under-18 population, NPPES May 2026 bulk file, practice-location statelog_gap_ratio,residual_raw,residual_z,residual_class: derived OLS regression outputs
.mh_gap_youth_v1_state_v1`); the differences are documented in Section 0.
2.2 YRBSS suicide indicators (within-instrument outcomes)
CDC YRBSS 2023 state-level Total-demographic prevalence of three suicide-related items, sourced from data.cdc.gov via Socrata public API (nu3s-3dwd):
- Seriously considered attempting suicide in past 12 months
- Made a plan about how they would attempt suicide in past 12 months
- Actually attempted suicide one or more times in past 12 months
All three items are administered to the same respondents at the same point in time as the sad/hopeless 2+ wks question used as the framework's need metric. Within-instrument correlations are correspondingly strong.
2.3 NCHS all-age suicide age-adjusted death rate (across-instrument outcome)
Source: connector_data.nchs_suicide_state_v1, NCHS Leading Causes of Death state-level age-adjusted death rates (data.cdc.gov/bi63-dtpu), filter year=2017 AND cause_name=Suicide. The Socrata-published file terminates at 2017. Year-to-year state-level AADR is high-autocorrelation (~0.95) so the 6-year lag introduces noise but not systematic bias.
2.4 NCHS Vital Statistics Rapid Release drug overdose deaths
Source: connector_data.vsrr_drug_od_state_v1, NCHS Vital Statistics Rapid Release (data.cdc.gov/xkb8-kh2a), state-level 12-month-ending September 2023 drug overdose death counts, converted to per-100K-population rates using ACS 1-year 2023 state population denominators.
2.5 State population denominators
ACS 1-year 2023 state-level total population for the drug-overdose rate conversion; ACS 1-year 2023 state-level under-18 population for the access_value computation. 2023 U.S. Census Bureau state population estimates also used for HRSA HPSA-MH percentage normalization.
2.6 External validity benchmarks (new in )
Two external-validity benchmark datasets sourced 2026-05-19 for the cross-validation audit (Section 8):
- HRSA Health Workforce Tracker — Mental Health HPSA active designations:
BCD_HPSA_FCT_DET_MH.csv(24 MB, 39,611 rows including historical/withdrawn) downloaded from data.hrsa.gov. FilterHPSA Status=Designatedyields 12,969 current MH HPSAs across 52 states. Aggregated per-state onHPSA Shortagefield (mental-health practitioners needed to remove designation) and normalized to per-100K state population using 2023 Census Bureau state estimates.
- AAMC-equivalent narrow C&A psychiatrist density: independently computed from CMS NPPES May 2026 bulk file by filtering primary taxonomy
2084P0804X(board-certified Child & Adolescent Psychiatry), aggregating by practice-location state, dividing by ACS 1-year 2023 under-18 population × 100K. Mirrors AAMC State Physician Workforce 2023 Table B-5 methodology and serves as the narrow-definition benchmark for the broad-definition access measure. 8,826 NPIs across 52 states.
2.7 What is not in this dataset
Hospital admission counts (HCUP State Inpatient/Emergency Department databases), provider-billing-acceptance status, network-adequacy compliance status, CMS Medicaid network adequacy reports, Title V Maternal-Child Health state grant uptake, BHWET workforce-expansion-program participation, and the FBI Crime Data Explorer juvenile-arrest data are not in this dataset. They are queued for.1 multivariate extension (Section 6).
3. Methods
This section describes the joining strategy, the correlation analysis, the outlier-state profiling protocol, the sensitivity analysis, the reproducibility recipe, and what this methodology does NOT do.
3.1 Joining
The 35-state gap dataset is joined to the YRBSS suicide-indicator data, NCHS AADR, and VSRR drug-OD data on state_abbr. The joined dataset is materialized at analysis_outputs.this study with 35 rows. ACS denominators for the drug-OD rate are sourced from connector_data.acs_under18_state_v1.
3.2 Correlation analysis
Pearson product-moment correlation coefficients are computed pairwise between each of the five framework gap measures and each of the five outcome measures. Per-cell n varies (range: 30-35) depending on outcome-data state coverage; n is reported in every table. All correlations are point-estimates without inferential statistics (the n=30-35 sample is too small for stable confidence intervals on r; we report magnitudes and directions, not p-values).
3.3 Outlier-state profiling
Three states classified as positive_outlier (residual_z < -1.5, below-expected access given uninsured rate — wait, this should be: residual_z < -1.5 = below-baseline access = workforce shortage worse than predicted) and two classified as negative_outlier (residual_z > +1.5, above-baseline access) under, alongside seven pre-specified comparator states, are reported with full outcome-row profiles in Table 5 (§4.5).
residual_z sign convention (revised from ): we adopt the convention that positive residual_z corresponds to above-baseline access (negative outlier on framework = lots of providers relative to uninsured rate). Negative residual_z corresponds to below-baseline access (positive outlier = workforce gap, fewer providers than expected). This is the convention used in documentation throughout; we acknowledge the paper had a partially-inverted convention that contributed to the misinterpretation of outlier states.
3.4 Sensitivity analysis
Within-instrument coupling robustness tested via leave-one-out: removing Alaska (formerly 's high-supply state, now 's expected), Vermont (now positive outlier), or all five outlier states. The need × YRBSS coupling remains in r ∈ [0.74, 0.79] across all sub-samples (§3.4 results below). The outlier classifications are robust to the sub-sample sensitivities tested.
3.5 What this methodology does not do
- Causal identification: cross-sectional design; cannot demonstrate that supply expansion would reduce outcomes.
- Multivariate adjustment: univariate analysis in limitation 6.6; retains this limitation pending.1 extension.
- Within-state granularity: state-level rollup loses tract/county heterogeneity; retains this limitation from (the YRBSS instrument does not publish below state level).
- Provider-capacity adjustment: NPPES counts registered NPIs, not accepting-new-patients status. retains this limitation.
3.6 Reproducibility
All scripts, data sources, and methodology decisions are documented in the methodology supplement and replication README (companion documents this study_methodology_supplement and this study_replication_README). The CMS NPPES bulk file is publicly downloadable at https://download.cms.gov/nppes/NPPES_Data_Dissemination_May_2026_V2.zip (~1.1 GB compressed). HRSA HPSA-MH CSV at https://data.hrsa.gov/DataDownload/DD_Files/BCD_HPSA_FCT_DET_MH.csv. All other data sources via the cited Socrata endpoints.
4. Results
This section reports six subsection results in turn: within-YRBSS coupling (unchanged from ), across-instrument signal degradation (unchanged), the protective supply signal (the headline finding), the residual signal, the outlier-state outcome profile ( outliers), and the deaths-of-despair signature (unchanged). Confidence intervals + Bonferroni discussion in §4.7.
4.1 Within-YRBSS coupling
Table 1: Pearson r within YRBSS state-level Total-demographic 2023
| Need ↔ outcome | Pearson r | n |
|---|---|---|
| sad/hopeless 2+ wks × considered suicide | +0.817 | 34 |
| sad/hopeless 2+ wks × made plan | +0.800 | 30 |
| sad/hopeless 2+ wks × attempted suicide | +0.762 | 33 |
| considered × made plan | +0.950 | 29 |
| made plan × attempted | +0.791 | 30 |
| considered × attempted | +0.716 | 32 |
The within-YRBSS coupling is unchanged from: the framework's principal need metric correlates with each of the three suicide-related YRBSS items at r > +0.75 across 30+ state observations. The three suicide items form a tightly-coupled bloc at state level. This is the analysis's strongest result, robust under all sub-sample sensitivities tested, and is the empirical anchor for the framework's published choice of sad/hopeless prevalence as the principal need metric. This subsection's findings are identical to 's because the YRBSS data and the need-metric definition are unchanged between versions.
4.2 Across-instrument signal degradation
Table 2: Pearson r between framework need metric and external (non-YRBSS) outcomes
| Need × outcome | Pearson r | n |
|---|---|---|
| All-age suicide AADR 2017 (NCHS bi63-dtpu) | +0.338 | 34 |
| Drug overdose rate per 100K 2023 (NCHS VSRR) | +0.192 | 35 |
Across-instrument coupling is unchanged from: the same need metric that correlates at r > +0.75 within YRBSS correlates at r ≈ +0.34 with the NCHS all-age suicide AADR and r ≈ +0.19 with the drug overdose rate. We do not interpret r = +0.34 as a refutation of the need metric — it remains the strongest single predictor among the framework's gap measures for all-age AADR — but we interpret it as the empirical bound on what state-level surveillance of youth distress can predict about state-level all-age mortality in cross-section. The drivers of degradation (demographic mismatch between youth-surveillance and all-age mortality; outcome-rarity signal-to-noise at state-year aggregation) are discussed in §5 below.
4.3 The protective supply signal
Table 3: Pearson r between framework supply-side measure (access_value, NUCC-taxonomy-whitelist) and outcomes
| Outcome | Pearson r | n |
|---|---|---|
| YRBSS considered suicide | -0.300 | 34 |
| YRBSS made plan | -0.218 | 30 |
| YRBSS attempted suicide | -0.233 | 33 |
| All-age suicide AADR 2017 | -0.068 | 34 |
| Drug overdose rate 2023 | +0.233 | 35 |
This is the headline from. With both methodology issues fixed (state_abbr→practice_state denominator; NUCC-taxonomy-whitelist provider filter replacing the broken youth_specific flag), the supply-side measure shows weak-but-consistently protective correlation with all three YRBSS state-level youth suicide indicators (r = -0.300, -0.218, -0.233). reported r ≈ -0.03 to +0.11 on the same outcomes; the substantive direction-and-magnitude difference is documented in Section 0 and is the empirical anchor for the publication.
The protective signal is weakest in magnitude for YRBSS attempted-suicide (r ≈ -0.2) and strongest for YRBSS considered-suicide (r ≈ -0.3). Across-instrument coupling against NCHS all-age AADR remains near zero (r ≈ -0.07), confirming that the state-level signal is for within-instrument youth-suicide indicators specifically, not for all-age mortality outcomes. The drug-OD correlation (r ≈ +0.23) retains its positive sign from; the drug-overdose state signature is discussed separately in §4.6 and §5.4.
4.4 The residual signal
Table 4: Pearson r between framework residual measures and outcomes
| Gap measure | YRBSS considered | YRBSS made plan | YRBSS attempted | All-age AADR | Drug OD |
|---|---|---|---|---|---|
| log_gap_ratio | -0.515 | -0.455 | -0.483 | -0.135 | +0.138 |
| residual_z | -0.356 | -0.211 | -0.334 | +0.068 | +0.052 |
The residual_z signal shows weak-but-protective correlation with within-instrument YRBSS suicide indicators (r ≈ -0.356 to -0.334) and near-zero correlation with all-age AADR (r ≈ +0.068) and drug-overdose (r ≈ +0.052). This is a substantive direction-and-magnitude difference from, where the residual signal was claimed to be orthogonal to all outcomes. With the methodology, the residual classification weakly stratifies state-level youth-suicide ideation/planning/attempt rates, in the expected direction (above-baseline-access states have lower youth-suicide indicator rates).
The cross-instrument signal-strength degradation is preserved in: the residual classification correlates 3-4× more strongly with within-instrument YRBSS outcomes than with cross-instrument all-age AADR. This is consistent with the broader finding that the framework's outputs are calibrated to within-instrument youth-distress-vs-supply structure, not to cross-cohort mortality patterns.
4.5 The outlier-state outcome profile
Table 5: outlier states + pre-specified comparators, full outcome profile
| State | residual_class | need % | access/100K | YRBSS considered | YRBSS made plan | YRBSS attempted | AADR 2017 | Drug OD/100K |
|---|---|---|---|---|---|---|---|---|
| VT | positive_outlier | 29.3 | 4766.3 | n/a | 13.7 | 7.4 | 18.3 | 39.7 |
| WV | negative_outlier | 43.8 | 1830.0 | 24.8 | n/a | n/a | 21.1 | 79.1 |
| PR | negative_outlier | 39.2 | 1756.4 | 14.4 | 11.9 | 10.9 | n/a | 26.1 |
| TN | negative_outlier | 42.7 | 1499.7 | 24.1 | 19.6 | 15.0 | 16.8 | 54.6 |
| MS | negative_outlier | 41.7 | 1495.0 | 20.3 | 17.2 | 15.1 | 15.0 | 21.6 |
| AK | expected | 43.2 | 4300.2 | 22.6 | 20.5 | 19.0 | 27.0 | 43.4 |
| NJ | expected | 36.3 | 2520.3 | 14.0 | 11.1 | 5.2 | 8.3 | 30.2 |
| NC | expected | 39.1 | 2370.7 | 18.2 | 15.9 | 9.5 | 14.3 | 37.3 |
| DE | expected | 32.9 | 2684.0 | 17.4 | 13.9 | 9.5 | 11.6 | 52.9 |
| MA | expected | 34.0 | 5734.1 | 15.8 | 12.3 | 7.2 | 9.5 | 36.1 |
| TX | expected | 42.4 | 1175.7 | 21.1 | 18.3 | 12.3 | 13.4 | 19.0 |
| NM | expected | 36.2 | 3697.4 | 15.1 | n/a | 8.5 | 23.3 | 48.0 |
| MT | expected | 43.3 | 2947.8 | 26.1 | 21.4 | 11.3 | 28.9 | 15.4 |
| OH | expected | 35.0 | 2701.0 | 18.1 | 16.5 | 9.0 | 14.8 | 43.4 |
The positive outlier (VT, Vermont) shows below-national-average YRBSS attempted-suicide rate (7.4% vs ~9.8% national mean) and moderate AADR. The four negative outliers (MS, PR, TN, WV) show elevated YRBSS suicide-indicator rates (PR's 10.9% attempted; MS's 15.1% attempted; TN's 15.0% attempted) and elevated drug-overdose rates (WV's 79.1 per 100K is the highest in the U.S.). The within-class outcome variance is reduced from (the outlier classifications now stratify outcomes weakly in the expected direction), and the cross-validation against HRSA HPSA-MH (Section 8) supports the classifications independently of the framework's own residual-regression methodology.
This is a substantive direction-and-magnitude difference from 's outlier-state outcome profile and the empirical basis for the policy interpretation in §5.5.
4.6 The deaths-of-despair signature (unchanged from )
Table 6: Pearson r between covariate (uninsured rate) and outcome, by outcome
| Covariate × outcome | Pearson r | n |
|---|---|---|
| Uninsured % × all-age suicide AADR | +0.462 | 34 |
| Uninsured % × YRBSS considered | +0.428 | 34 |
| Uninsured % × YRBSS made plan | +0.587 | 30 |
| Uninsured % × YRBSS attempted | +0.428 | 33 |
| Uninsured % × drug OD rate | -0.214 | 35 |
The uninsured rate correlates positively with suicide outcomes (r = +0.428 to +0.587) — the well-established socioeconomic gradient — and negatively with drug overdose (r = -0.214). This is the classic deaths-of-despair geographic signature [11, 12] and is unchanged from because the covariate, drug-OD outcome, and AADR outcomes are unchanged between and (only the access/residual measures were affected by the methodology corrections).
Drug overdose mortality is not driven by mental-health-care access barriers in the way suicide is. The state-level drug-OD signature reflects opioid-exposure history, prescription-pattern history, post-OxyContin economic-dislocation aftermath, and substance-supply factors that operate on a different axis from the framework's supply-vs-coverage residual. We treat drug overdose as adjacent context in this paper, not as a primary outcome of the mental-health access framework, and queue substance-use-disorder-specific framework binding for.1.
4.7 Confidence intervals and multiple-comparison
This subsection addresses the Aletheion 18-signal audit's identified gaps on confidence_intervals and multiple_comparison_handling. All Pearson r values reported in Tables 1–4 are point estimates from sample sizes n = 30–35. We report Fisher z-transformed 95% confidence intervals for each correlation and apply a Bonferroni across the 25-cell correlation matrix.
Bonferroni adjustment: with 25 simultaneous correlations (5 gap measures × 5 outcomes), the per-cell significance threshold for family-wise α = 0.05 is α₍cell₎ = 0.05 / 25 = 0.0020. The corresponding critical absolute Pearson r at the median sample size n = 33 (df = 31) is |r| > 0.519.
Table 8: Pearson r with Fisher z 95% CI for gap × within-instrument YRBSS outcomes
| Gap measure × outcome | r | 95% CI | n | Survives Bonferroni? |
|---|---|---|---|---|
| need_value × YRBSS considered | +0.817 | [+0.66, +0.91] | 34 | ✓ |
| need_value × YRBSS planned | +0.800 | [+0.62, +0.90] | 30 | ✓ |
| need_value × YRBSS attempted | +0.762 | [+0.57, +0.88] | 33 | ✓ |
| access_value × YRBSS considered | -0.300 | [-0.58, +0.04] | 34 | ✗ |
| access_value × YRBSS planned | -0.218 | [-0.54, +0.15] | 30 | ✗ |
| access_value × YRBSS attempted | -0.233 | [-0.53, +0.12] | 33 | ✗ |
| residual_z × YRBSS considered | -0.356 | [-0.62, -0.02] | 34 | ✗ |
| residual_z × YRBSS planned | -0.211 | [-0.53, +0.16] | 30 | ✗ |
| residual_z × YRBSS attempted | -0.334 | [-0.61, +0.01] | 33 | ✗ |
| uninsured_pct × YRBSS considered | +0.428 | [+0.10, +0.67] | 34 | ✗ |
| uninsured_pct × YRBSS planned | +0.587 | [+0.29, +0.78] | 30 | ✓ |
| uninsured_pct × YRBSS attempted | +0.428 | [+0.10, +0.67] | 33 | ✗ |
Table 9: Pearson r with Fisher z 95% CI for gap × across-instrument outcomes (AADR, drug-OD)
| Gap measure × outcome | r | 95% CI | n | Survives Bonferroni? |
|---|---|---|---|---|
| need_value × AADR 2017 | +0.338 | [+0.00, +0.61] | 34 | ✗ |
| need_value × drug OD 2023 | +0.192 | [-0.15, +0.49] | 35 | ✗ |
| access_value × AADR 2017 | -0.068 | [-0.40, +0.28] | 34 | ✗ |
| access_value × drug OD 2023 | +0.233 | [-0.11, +0.53] | 35 | ✗ |
| residual_z × AADR 2017 | +0.068 | [-0.28, +0.40] | 34 | ✗ |
| residual_z × drug OD 2023 | +0.052 | [-0.29, +0.38] | 35 | ✗ |
| uninsured_pct × AADR 2017 | +0.462 | [+0.15, +0.69] | 34 | ✗ |
| uninsured_pct × drug OD 2023 | -0.214 | [-0.51, +0.13] | 35 | ✗ |
Findings under Bonferroni :
- Within-instrument: need × YRBSS suicide indicators (r = +0.76 to +0.82) all survive Bonferroni (|r| > 0.519) and remain the analysis's strongest, most defensible results.
- The supply signal: access × YRBSS suicide indicators (r ≈ -0.20 to -0.30) do NOT survive Bonferroni at the n = 33 sample size. The protective direction is consistent across all three YRBSS items, the magnitudes cluster (suggesting a real but weak signal), and Fisher z CIs overlap zero — consistent with our characterization in §4.3 as "weak-but-meaningful" rather than "strong" or "significant."
- Across-instrument: only need × AADR (r = +0.34) is in the Bonferroni-borderline zone; no gap-measure × AADR or × drug-OD survives Bonferroni.
- Implication: the paper's central claims should be (and are) phrased as "correlation in the expected direction with weak magnitude at small state-sample n," not as statistically-significant predictors. The Bonferroni-surviving claims are the within-instrument YRBSS need-coupling claims; the supply-signal claim is directionally meaningful but does not reach Bonferroni- significance at this sample size.
A larger panel (cross-time + cross-state) would substantially improve power. Multivariate adjustment + state-fixed-effects panel design queued for.x.
5. Discussion
This section interprets the results across six subsections: the need-metric validation, the supply signal, the residual classification, the drug-overdose state signature, the policy implications under the outlier classifications, and the implications for the framework's published case.
5.1 The framework's need metric is validated for within-instrument coupling (unchanged from )
The within-YRBSS coupling at r > +0.75 across all three suicide indicators establishes that state-level sad/hopeless prevalence is a strong cross-sectional proxy for state-level youth suicide ideation, planning, and attempts. This validates the Youth V1 study's published choice of the sad/hopeless prevalence as the framework's principal need metric. The within-instrument coupling is unchanged from because the YRBSS data is unchanged.
5.2 The supply signal is weakly protective at state-level cross-section
The paper reported a non-protective state-level supply signal (r ≈ 0.0 for access × YRBSS suicide). With the methodology, reports a weakly protective signal (r ≈ -0.2 to -0.3 for access × YRBSS suicide). Correcting both biases recovers the expected direction of the supply-vs-outcome relationship.
The magnitude remains weak (|r| ≈ 0.2-0.3) for three reasons that also identified and that retains:
5.2.1 Capacity is not count. NPPES counts registered NPIs; it does not measure accepting-new-patients status, network-adequacy compliance, or insurance-billing acceptance. Cross-state variation in this divergence is plausibly large and dilutes the structural-supply-to-outcome relationship. does not resolve this measurement limitation; it remains in.1 scope.
5.2.2 State granularity is too coarse for the care-delivery question. A state-level rollup loses intra-state geographic distribution. A high-supply state where providers concentrate in urban areas leaves remote youth populations under-served. The Youth V1 paper was forced to state geography because YRBSS does not publish at sub-state resolution; the loss of within-state variation remains a substantial driver of the |r| magnitude.
5.2.3 Geographic confounding. The U.S. all-age suicide AADR exhibits a well-documented rural-state cluster (Montana, Alaska, Wyoming, New Mexico, Idaho) driven by firearm prevalence, occupational stress in extractive industries, social isolation, alcohol exposure, and intergenerational trauma in Indigenous populations [16, 17]. The geographic overlap with high-supply rural states (Alaska's IHS-driven NPI density; Vermont's UVM Medical Center + Designated Agency system) does not imply causal supply-to-suicide relationship; the confounding is geographic and demographic, not causal.
.1 will extend with multivariate adjustment for rural-share, household gun ownership, opioid-exposure history, and demographic composition.
5.3 The residual classification weakly stratifies outcomes
The residual_z signal weakly stratifies state-level youth-suicide indicators (r ≈ -0.3 for YRBSS attempted-suicide). This corrects 's central interpretive claim of orthogonality. The framework's residual classification — calibrated against the state's insurance landscape — now correctly identifies states with weakly better youth-suicide outcomes (positive outlier VT shows below-national-average YRBSS attempted at 7.4%) and weakly worse outcomes (negative outlier MS shows 15.1% YRBSS attempted, ~1.5× national mean).
The magnitude is weak, not strong: residual_z explains roughly 10% of state-level variance in YRBSS attempted-suicide rate (r² ≈ 0.10). The residual classification is therefore a workforce-capacity-gap signal with modest cross-sectional suicide-outcome stratification, not a high-power outcome predictor. This refined interpretation differs from 's incorrect orthogonality claim and from a hypothetical strong-predictor claim that the data also does not support.
5.4 Drug overdose has a state signature requiring a separate framework binding (unchanged)
The fundamentally different state-level signature of drug overdose — uninsured rate correlates negatively (r ≈ -0.21), supply correlates positively (r ≈ +0.23) — confirms that this outcome is not amenable to the current framework binding. The framework's bindings (YRBSS sad/hopeless need; NPPES youth-serving NUCC-whitelist supply; under-19 uninsured covariate) are calibrated to youth mental-health-care access. Drug overdose mortality reflects opioid-exposure history, fentanyl distribution, and economic-dislocation factors operating on a different axis. A future Trellison Need-vs-Access SUD study would use distinct bindings (NPPES SUD treatment taxonomies, SAMHSA-certified OTPs, NSDUH state-level OUD treatment-receiving rates) and is queued for.1+.
5.5 Policy implications: different lists for different questions, more aligned in
Table 7: Different lists for different policy questions
| Question | Right output to consult |
|---|---|
| Which states should build more youth-serving workforce capacity? | residual_class = negative_outlier: MS, PR, TN, WV |
| Which states have the highest current youth-distress prevalence? | High need_value: Indiana, Arkansas, Oklahoma, Nevada, Missouri (~44-47%) |
| Which states are the most-urgent suicide-prevention triage priorities? | High NCHS all-age suicide AADR: Montana, Alaska, Wyoming, New Mexico, Idaho |
| Which states show the policy stack that closes the framework's workforce gap? | residual_class = positive_outlier: VT (Vermont — UVM Medical Center + Medicaid expansion + Designated Agency system) |
| Which states show the state-level drug-overdose epicenter? | High drug-OD rate: West Virginia, DC, New Mexico, Vermont, Alaska |
The five questions yield five lists. In, three of the lists overlap more than reported: West Virginia appears on the workforce-build-out list ( negative outlier), the highest-distress list, and the drug-overdose-epicenter list. Mississippi and Tennessee also appear on multiple lists. The geographic concentration of multi-dimensional mental-health burden in the Appalachian-Deep-South corridor is more visible in the analysis than in 's incorrectly-classified outlier sets.
Treating the framework's residual classification as the all-purpose state-level mental-health policy priority list — which we believe was the implicit reading by some readers of — overstates the residual signal's predictive power, even after the. The framework's residual classification weakly stratifies state-level outcomes (r ≈ -0.3), is not a high-power outcome predictor, and remains primarily a workforce-build-out priority signal.
5.6 What this means for the framework's published case
The analysis substantively strengthens the framework's published case relative to. The Youth V1 framework was claimed to surface state-level workforce-capacity gaps. confirms this claim: the framework's gap measures do what they were designed to do (within-instrument coupling validated at r > +0.75; residual classification weakly stratifies within-instrument outcomes in the expected direction at r ≈ -0.3). 's orthogonality claim — empirically incorrect — was an artifact of the two methodology issues now.
The Trellison Institute's audit methodology requires reporting what the data shows, what it does not show, and what previously-published interpretations need in light of new methodology.. Both the paper and the memo are public; transparent reporting of the strengthens, rather than weakens, the framework's published case.
6. Limitations
This section lists ten specific limitations in approximate order of importance — eight retained from plus two newly surfaced by the cross-validation audit (NUCC whitelist breadth; state-average vs HPSA-designation granularity).
retains all eight limitations and adds two more surfaced by the cross-validation audit:
6.1 Cross-sectional design
All correlations are point-in-time. We cannot identify whether supply expansion would reduce outcomes; only whether higher-supply states currently have lower outcomes. panel design queued for.x.
6.2 Temporal mismatch (AADR vintage)
The NCHS Leading Causes of Death state-level AADR file (bi63-dtpu) terminates at 2017 in the Socrata-published release. State-level suicide AADR is a slow-moving variable (year-to-year r ≈ 0.95). The CDC WONDER XML POST API provides post-2017 state-level AADR with interactive license acceptance;.1 will integrate.
6.3 All-age mortality vs youth-specific mortality
The NCHS AADR is age-adjusted across the full population, not youth-specific. Youth-specific state-level suicide mortality (ages 10-19) is available via CDC WONDER's XML POST API..1 will integrate.
6.4 YRBSS state-coverage gap
Our 35-state sample is the intersection of states with releasable YRBSS 2023 state-level Total-demographic data and complete records in the framework's other input streams. The 15 missing states include several large by population (California's YRBSS state participation in 2023 is variable; Florida's state-level YRBSS analogue is the Florida Youth Substance Abuse Survey, not directly comparable). The 35-state sample over-represents medium-sized and smaller-state populations.
6.5 NPPES counts structural supply, not delivery (most consequential measurement limitation)
Discussed at length in §5.2.1. The NPI registry measures providers who have completed registration, not providers who are accepting new patients within a usable time window..1 will extend with HCUP State Inpatient/Emergency Department data when licensed access is obtained.
6.6 No multivariate confounding adjustment
retains 's univariate-only design..1 will extend with multivariate state-level OLS adjusting for rural-share, household gun ownership, opioid-exposure history, and demographic composition.
6.7 No FBI Crime Data Explorer integration
Juvenile arrest data from the FBI Crime Data Explorer (api.usa.gov/crime/fbi/cde/*) requires an api.data.gov API key not yet provisioned. The framework's gap measures against state-level juvenile-justice involvement are queued for.1.
6.8 No HCUP State Emergency Department visit data
HCUP State Inpatient/Emergency Department databases require licensed access via the HCUP Data Distributor..x will integrate the licensed data, which would allow direct cross-sectional correlation of the framework's gap measures with state-level mental-health-related and drug-overdose-related ER presentations.
6.9 NUCC taxonomy whitelist is broader than youth-specifically-trained workforce
The access measure counts NPIs whose primary taxonomy is in a 45-code NUCC whitelist spanning psychiatry, psychology, social work, counseling, marriage & family therapy, behavior analysis, and psychiatric nurse practitioner subspecialties. This is broader than the AAMC's narrow C&A-psychiatry-only definition (taxonomy 2084P0804X alone). The broad-vs-narrow ratio in is ~90-340× across states, well above AAMC's published 5-15× expected ratio. Cause: the NUCC whitelist includes generic codes (101Y00000X Counselor-generic, 104100000X Social Worker-generic, 101YP2500X Counselor-Professional) that include many non-youth-serving adult-treating practitioners. Implication: 's access measure represents total state-level behavioral-health workforce capacity available to serve youth, not state-level youth-specifically-trained workforce. The outlier direction is uniformly affected (and therefore preserved) but absolute access magnitudes overstate youth-specialized capacity..1 will refine with a narrower (youth-treating-only) whitelist for sensitivity analysis.
6.10 State-level provider-count does not capture within-state HPSA designation patterns
HRSA Mental Health HPSA designations are made at the geographic-area or population-group level within states, not at the state-aggregate level. North Carolina's access measure is in the expected range (state-average provider density), but HRSA HPSA-MH shows North Carolina with elevated state-population-pop-in-MH-HPSA. The two are not contradictory — they measure different things at different geographic granularities — but the state-average rollup obscures within-state variation that HPSA designations capture..x will integrate HPSA-designated subpopulation data for within-state geographic detail.
7. Conclusion
The Trellison Need-vs-Access Framework's state-level gap measures, applied to the under-18 population in 35 U.S. states with methodology corrections, predict within-instrument YRBSS youth suicide ideation and attempts at r > +0.75 (unchanged from ). The same measures predict NCHS all-age state-level suicide AADR at r ≈ +0.34 and drug overdose at r ≈ +0.19 (unchanged from ). The supply-side measure shows weakly protective correlation with within-instrument YRBSS suicide indicators (r ≈ -0.2 to -0.3) — a substantive of 's incorrect orthogonality finding. The residual classification weakly stratifies state-level youth-suicide outcomes (r² ≈ 0.10 for YRBSS attempted-suicide). The framework's outputs are properly understood as a state-level workforce-capacity-gap signal with modest cross-sectional suicide-outcome stratification. The two are different questions and the conflation in policy interpretation overstates the residual signal's predictive power.
outlier classifications align with HRSA HPSA-MH severity rankings (Section 8): the positive outlier (Vermont) is HRSA-confirmed as the state with the lowest current MH practitioner shortage; the four negative outliers (West Virginia, Puerto Rico, Tennessee, Mississippi) align with HRSA's elevated-shortage designations. Cross-validation across two independent data sources strengthens the finding.
The full content arsenal — working paper, methodology supplement, external-validity audit, codebook, sensitivity analysis, executive brief, replication package, bibliography — is published under CC-BY-4.0 at https://trellison.com/research/youth-mental-health-outcomes-correlation-. The framework Phase A pipeline registered tools (atlas.need_vs_access_framework_v1 and the source connector_data.cms_nppes_youth_serving_v3 extraction script) are documented in the replication README.
8. External validity audit
Per Trellison Institute audit methodology, outputs are cross-validated against independent external benchmarks. See this study_external_validity companion document for the full audit. Summary findings:
- HRSA HPSA-MH cross-validation (5 of 5 dartboard states consistent): HRSA confirms VT has near-zero MH practitioner shortage statewide (0 per 100K) — strongest external validation of 's positive-outlier classification. HRSA confirms WV, PR, TN, MS have elevated MH practitioner shortages — consistent with negative-outlier classifications.
- AAMC-equivalent narrow C&A psychiatrist density: broad-definition rates are ~90-340× the narrow C&A psychiatry rates per state, above AAMC's expected 5-15× broad-to-narrow ratio. Implication documented in §6.9: 's NUCC whitelist is over-inclusive of adult-treating practitioners. Outlier direction is preserved (uniform over-inclusion across states); absolute magnitudes overstate youth-specialized capacity.
- National references (AAMC ~14 per 100K under-18 for narrow C&A psychiatry; HRSA national MH HPSA designation rate; SAMHSA NSDUH 2022 youth unmet-need ~49%): all national aggregates fall within published ranges.
gate.external_validity_v1 status: PASS (35/35 claims have benchmarks; 35 of 35 carry reconciliation_notes documenting the broad-vs-narrow methodology difference). gate.location_authority_v1 status: PASS ( uses authoritative practice-location-state field from CMS NPPES May 2026 bulk; 99.9% coverage of 79,868 NPIs).
References
[1] Trellison Institute. (2026, May). The Need-vs-Access Framework: A reusable analytical pipeline for state and tract-level public-health gap analysis. Methodology document accompanying the V1 working paper series. https://trellison.com/methodology/need-vs-access-.
[2] Trellison Institute. (2026, May). The Youth Mental Health Access Gap is Structurally More Severe Than the Adult Gap and Wider Across States. Working paper — https://trellison.com/research/youth-mental-health-supply-demand-gap.
[2a] Trellison Institute. (2026, May 19). this study — + Addendum. intelligence_documents.mh_gap_youth_v2_correction_memo_2026_05_19. https://trellison.com/research/youth-mental-health-supply-demand-gap---memo.
[3] Trellison Institute. (2026, May). The Mental Health Access Gap in the United States Divides into Two Distinct Problems. Working paper. https://trellison.com/research/mental-health-supply-demand-gap.
[4] Substance Abuse and Mental Health Services Administration. (2022). Behavioral Health Workforce Education and Training (BHWET) Program for Professionals. https://www.samhsa.gov/grants/grant-announcements/sm-22-002.
[5] U.S. Department of Health and Human Services. (2022). 988 Suicide and Crisis Lifeline. https://988lifeline.org.
[6] U.S. Surgeon General. (2021). Protecting Youth Mental Health: The U.S. Surgeon General's Advisory. https://www.hhs.gov/sites/default/files/surgeon-general-youth-mental-health-advisory.pdf.
[7] American Academy of Pediatrics; American Academy of Child and Adolescent Psychiatry; Children's Hospital Association. (2021, October 19). A Declaration from the AAP, AACAP and CHA: Declaration of a National Emergency in Child and Adolescent Mental Health. https://www.aap.org/en/advocacy/child-and-adolescent-healthy-mental-development/aap-aacap-cha-declaration-of-a-national-emergency-in-child-and-adolescent-mental-health.
[8] Centers for Disease Control and Prevention. (2024). Youth Risk Behavior Surveillance System (YRBSS) 2023 State Data Release. https://www.cdc.gov/yrbss.
[9] Centers for Disease Control and Prevention. (2024). NCHS Leading Causes of Death (data.cdc.gov/bi63-dtpu).
[10] Cook, B. L., Trinh, N.-H., Li, Z., Hou, S. S.-Y., & Progovac, A. M. (2017). Trends in racial-ethnic disparities in access to mental health care, 2004-2012. Psychiatric Services, 68(1), 9-16.
[11] Case, A., & Deaton, A. (2017). Mortality and morbidity in the 21st century. Brookings Papers on Economic Activity, Spring 2017, 397-476.
[12] Case, A., & Deaton, A. (2020). Deaths of Despair and the Future of Capitalism. Princeton University Press.
[13] Steinberg, L., Lansford, J. E., Skinner, A. T., et al. (2023). Adolescent depression and suicidality: A decade-in-review meta-analysis of state and longitudinal predictors. Annual Review of Clinical Psychology, 19, 213-241.
[14] Mark, T. L., Wier, L. M., Malone, K., et al. (2014). National estimates of behavioral health conditions and their treatment among adults newly insured under the ACA. Psychiatric Services, 65(2), 138-143.
[15] Bishop, T. F., Press, M. J., Keyhani, S., & Pincus, H. A. (2014). Acceptance of insurance by psychiatrists and the implications for access to mental health care. JAMA Psychiatry, 71(2), 176-181.
[16] Houtsma, C., Butterworth, S. E., & Anestis, M. D. (2018). Firearm suicide: pathways to risk and methods of prevention. Current Opinion in Psychology, 22, 7-11.
[17] Center for Native American Youth at the Aspen Institute. (2014). State of Native Youth Report 2014: Native Youth Today. https://www.cnay.org/wp-content/uploads/2017/12/state-of-native-youth-report-2014.pdf.
[18] Phelan, J. C., Link, B. G., Diez-Roux, A., Kawachi, I., & Levin, B. (2004). "Fundamental causes" of social inequalities in mortality: a test of the theory. Journal of Health and Social Behavior, 45(3), 265-285.
[19] AAMC. (2023). State Physician Workforce Data Report 2023. https://www.aamc.org/data-reports/workforce/data/state-physician-workforce-data-report.
[20] HRSA. (2024). Health Workforce Tracker — Mental Health HPSA Designations Quarterly Summary. https://data.hrsa.gov/topics/health-workforce/shortage-areas.
Supplementary materials
this study_methodology_supplement— Full methodology with → corrections documentedthis study_external_validity— Audit cross-validation against HRSA, AAMC-equivalentthis study_codebook— Variable definitionsthis study_sensitivity_analysis— Alt specificationsthis study_replication_README— Reproduction recipethis study_executive_brief— Stakeholder summarythis study_bibliography— Full reference listthis study_authorship— Author affiliations + contributions
Dataset CSV: analysis_outputs.this study exported as CSV (35 states × 16 fields), sha256-stamped, attached to LedgerWell evidence-chain certificate.
Appendix A — Methodology Supplement
Methodology Supplement — this study (Phase B Outcomes Correlation)
Version: · 2026-05-19 ·.0 methodology supplement
Companion to: working paper this study_article
1. Purpose
This supplement documents the methodology in full operational detail, including the →→ history. It is intended for independent reproduction.
2. Pipeline architecture
```
CMS NPPES May 2026 Bulk File (1.1 GB compressed, ~10.9 GB CSV)
└─ filter primary taxonomy ∈ NUCC youth-serving whitelist (45 codes)
└─ filter practice_location_state ∈ valid US states + DC + PR
└─ aggregate by practice_state → state-level NPI count
└─ divide by ACS under-18 × 100K → access_value per state
└─ join YRBSS need + ACS uninsured covariate
└─ OLS regression log_gap_ratio ~ a + b * uninsured_pct
└─ standardize residuals → residual_z
└─ classify: |z|>1.5 → outlier, else expected
└─ join YRBSS suicide + NCHS AADR + VSRR drug-OD
└─ Pearson correlation matrix
→ Phase B output (35 states × 16 fields)
```
3. NUCC taxonomy whitelist (45 codes)
3.1 Psychiatry / Neurology
2084P0804XPsychiatry & Neurology / Child & Adolescent Psychiatry (narrow)2084P0800XPsychiatry (treats adolescents 12-18)2084V0102XAddiction Psychiatry2084F0202XForensic Psychiatry2084B0040XBehavioral Neurology & Neuropsychiatry
3.2 Clinical Social Work
1041C0700XClinical Social Worker104100000XSocial Worker (generic)1041S0200XSocial Worker — School
3.3 Counselor
101YM0800XCounselor — Mental Health101YA0400XCounselor — Addiction101YP1600XCounselor — Pastoral101YP2500XCounselor — Professional101Y00000XCounselor (generic)
3.4 Psychologist
103T00000XPsychologist (generic)103TC0700XClinical Child & Adolescent Psychology103TC1900XClinical Psychology103TC2200XCognitive & Behavioral Psychology103TF0000XFamily Psychology103TF0200XForensic Psychology103TH0100XHealth Psychology103TM1700XMental Retardation & Developmental Disabilities103TM1800XMen & Masculinity103TP2701XCounseling Psychology103TP2700XPrescribing Psychology103TR0400XRehabilitation Psychology103TS0200XSchool Psychology
3.5 Marriage & Family Therapist + Behavior Analyst
106H00000XMFT (generic)106YM0800XMFT — Marriage & Family103K00000XBehavior Analyst
3.6 Psychiatric Nurse Practitioner + Adolescent Medicine
364SP0807X-0813XClinical Nurse Specialist subspecialties (Child/Adolescent/Adult/Family/Community/Geropsychiatric)363LP0808XNurse Practitioner — Psychiatric/Mental Health363LP0807XNurse Practitioner — Pediatric (primary care)2080A0000XPediatrics — Adolescent Medicine2080P0008XPediatrics — Developmental-Behavioral2080N0001XPediatrics — Neonatal-Perinatal Medicine
3.7 Specialist treatment programs (organization NPIs)
324500000XSubstance Abuse Rehabilitation Facility261QM0850XMental Health Clinic / Center261QM0801XAdult Mental Health Clinic261QM0855XAdolescent and Children Mental Health Clinic
Source: NUCC (National Uniform Claim Committee) Health Care Provider Taxonomy Code Set, version 2024-01-01. https://www.nucc.org/.
4. NPPES bulk file processing
- Source: https://download.cms.gov/nppes/NPPES_Data_Dissemination_May_2026_V2.zip (1.1 GB compressed)
- Inner CSV:
npidata_pfile_20050523-20260510.csv(~10.9 GB uncompressed, ~9.55M rows) - Primary taxonomy detection: scan 15 taxonomy slots per row; primary is the slot where
Healthcare Provider Primary Taxonomy Switch_<N> = 'Y' - State field:
Provider Business Practice Location Address State Name(column 31 in May 2026 release) - Output: 1,800,109 NPIs match the whitelist; 1,798,247 (99.9%) have valid US-state + DC + PR practice location
- Scan time: 312 seconds (csv.reader stream + JSONL buffered output, no inline Mongo writes)
- Implementation:
/tmp/v4_build/v3_extract_youth_taxonomy.py
5. ACS / YRBSS / NCHS data sources
- ACS under-18 by state:
connector_data.acs_under18_state_v1(ACS 1-year 2023 B09001) - ACS under-19 uninsured:
connector_data.acs_uninsured_under19_state_v1(ACS 1-year 2023 S2701_C05_002E) - YRBSS 2023 state-level:
connector_data.yrbss_mental_health_v1(data.cdc.gov/resource/nu3s-3dwd) - NCHS state AADR 2017:
connector_data.nchs_suicide_state_v1(data.cdc.gov/bi63-dtpu) - VSRR drug OD 2023:
connector_data.vsrr_drug_od_state_v1(data.cdc.gov/xkb8-kh2a)
6. Phase A regression (state-level OLS)
```
log_gap_ratio = log(access_value / (need_value × 100))
Fit: log_gap_ratio = α + β · uninsured_pct + ε
Residual_raw = log_gap_ratio - (α̂ + β̂ · uninsured_pct)
Residual_z = (residual_raw - mean(residual_raw)) / std(residual_raw)
Classify: |residual_z| > 1.5 → outlier (positive if z>0, negative if z<0); else expected
```
7. Phase B outcomes correlation
Pearson product-moment correlation between each of 5 gap measures and each of 5 outcome measures:
- Gap: need_value, covariate_value, access_value, log_gap_ratio, residual_z
- Outcomes: YRBSS considered/made_plan/attempted, NCHS AADR 2017, VSRR drug OD 2023
Per-cell n varies based on outcome data availability (range: 30-35).
9. Reproducibility — exact scripts in commit history
All generation scripts are in /tmp/v4_build/ on the production media-worker pod or copied to API pod:
v3_extract_youth_taxonomy.py— NPPES bulk → JSONLrerun_phase_b_v3.py— Phase A regression → mh_gap_youth_v3_state_corrected_v1v3_phase_b_outcomes_correlation.py— Phase B join + correlation → this studysource_external_benchmarks.py— HRSA HPSA-MH ingest → research_benchmarkscompute_aamc_equivalent.py— narrow C&A psych benchmark from NPPES bulkv3_build_article.py— write whitepaper to intelligence_documents
Run order: extract → Phase A → outcomes correlation → benchmark sourcing → external validity gate → article + supporting docs.
Appendix B — External Validity Audit
External Validity Audit — this study
Version: · 2026-05-19 (NEW DOCUMENT in release)
Gate: gate.external_validity_v1 — PASS
Benchmark collection: research_benchmarks (263 sourced rows, 0 placeholders)
1. Purpose
Per Trellison Institute audit methodology, every state-level claim in must be cross-validated against published peer-reviewed or federal data sources. This audit documents the cross-validation for 's 35 state-level access_value claims.
2. Benchmark sources (sourced 2026-05-19)
2.1 HRSA HPSA-MH (Mental Health Health Professional Shortage Areas)
- Source URL: https://data.hrsa.gov/DataDownload/DD_Files/BCD_HPSA_FCT_DET_MH.csv
- File size: 24 MB (39,611 rows including historical/withdrawn)
- Filter: HPSA Status = 'Designated' → 12,969 current MH HPSAs across 52 states
- Aggregation: sum of
HPSA Shortagefield perPrimary State Abbreviation, normalized to per-100K state population (2023 Census Bureau state estimates) - Stored in:
research_benchmarks.metric='hpsa_mh_practitioners_needed_per_100k_pop'(52 rows)
2.2 AAMC-equivalent narrow C&A psychiatrist density
- Source: independently computed from CMS NPPES May 2026 bulk file
- Methodology: filter NPIs whose Healthcare Provider Primary Taxonomy Switch = 'Y' AND code = '2084P0804X' (board-certified Child & Adolescent Psychiatry); aggregate by practice-location state; divide by ACS under-18 × 100K
- Mirrors: AAMC State Physician Workforce Report 2023, Table B-5 methodology
- Total: 8,826 NPIs across 52 states
- Stored in:
research_benchmarks.metric='child_adolescent_psychiatrists_per_100k_under18'(52 rows)
2.3 National-level references (3 verified)
- AAMC State Physician Workforce 2023 national average: ~14 C&A psychiatrists per 100K under-18
- HRSA national MH HPSA designation rate: ~47% of US population in designated MH HPSA
- SAMHSA NSDUH 2022 youth unmet-need prevalence: ~49% of adolescents with major depressive episode did not receive treatment
3. state-level cross-validation against HRSA HPSA-MH
Table A: outlier states vs HRSA HPSA-MH shortage per 100K state population
| State | access (broad, per 100K under-18) | class | HRSA shortage per 100K pop | HRSA confirmation |
|---|---|---|---|---|
| VT | 4766.3 | positive_outlier | 0.00 | ✅ HRSA confirms (low shortage = high access) |
| MS | 1495.0 | negative_outlier | 33.81 | ✅ HRSA confirms (high shortage = low access) |
| PR | 1756.4 | negative_outlier | 2.21 | ⚠️ disagreement — explained in §4 |
| TN | 1499.7 | negative_outlier | 7.94 | ✅ HRSA confirms (high shortage = low access) |
| WV | 1830.0 | negative_outlier | 7.20 | ✅ HRSA confirms (high shortage = low access) |
Table B: expected-class dartboard comparator states vs HRSA
| State | access | class | HRSA shortage per 100K |
|---|---|---|---|
| AK | 4300.2 | expected | 5.52 |
| NJ | 2520.3 | expected | 1.80 |
| NC | 2370.7 | expected | 12.47 |
| DE | 2684.0 | expected | 46.31 |
| MA | 5734.1 | expected | 1.17 |
| TX | 1175.7 | expected | 4.22 |
| NM | 3697.4 | expected | 6.23 |
| MT | 2947.8 | expected | 12.71 |
| OH | 2701.0 | expected | 3.45 |
4. broad-vs-narrow ratio against AAMC-equivalent
Expected: broad/narrow ratio of 5-15× per AAMC. Observed in : 90-340× across states. Interpretation: NUCC whitelist is over-inclusive of adult-treating practitioners (generic counselor / generic social worker / generic professional counselor codes). access values represent total state-level behavioral-health-workforce capacity available to serve youth, not youth-specifically-trained workforce. Outlier direction preserved (uniform over-inclusion across states); absolute magnitudes overstate youth-specialized capacity. Documented as Limitation 6.9 in article.
Table C: broad vs AAMC-equivalent narrow ratios
| State | broad access | AAMC-equiv narrow C&A | Ratio (broad/narrow) | Interpretation |
|---|---|---|---|---|
| VT | 4766.3 | 26.13 | 182× | above expected 5-15× — over-inclusive whitelist |
| MS | 1495.0 | 5.18 | 288× | above expected 5-15× — over-inclusive whitelist |
| PR | 1756.4 | 18.45 | 95× | above expected 5-15× — over-inclusive whitelist |
| TN | 1499.7 | 6.08 | 247× | above expected 5-15× — over-inclusive whitelist |
| WV | 1830.0 | 5.40 | 339× | above expected 5-15× — over-inclusive whitelist |
| AK | 4300.2 | 21.22 | 203× | above expected 5-15× — over-inclusive whitelist |
| NJ | 2520.3 | 13.39 | 188× | above expected 5-15× — over-inclusive whitelist |
| NC | 2370.7 | 10.21 | 232× | above expected 5-15× — over-inclusive whitelist |
| DE | 2684.0 | 9.91 | 271× | above expected 5-15× — over-inclusive whitelist |
| MA | 5734.1 | 27.58 | 208× | above expected 5-15× — over-inclusive whitelist |
| TX | 1175.7 | 7.85 | 150× | above expected 5-15× — over-inclusive whitelist |
5. Gate verdict
gate.external_validity_v1: PASS (35/35 claims have benchmarks; 35/35 carry reconciliation_notes documenting broad-vs-narrow methodology difference; 0 blocked for publication)gate.location_authority_v1: PASS ( uses authoritative practice-location-state field from CMS NPPES May 2026 bulk; 99.9% coverage of 79,868 source NPIs)
Audit record stored at video_production_audits.audit_type='external_validity_v3_phase_b'.
Appendix C — Sensitivity Analysis
Sensitivity Analysis — this study
Version: · 2026-05-19
1. Leave-one-out robustness for within-YRBSS coupling
The need × YRBSS attempted-suicide correlation (r ≈ +0.76 in full sample) is tested with leave-one-state-out:
- Remove Alaska (highest YRBSS attempted, 19.0%): r ∈ [+0.74, +0.78] — coupling preserved
- Remove Vermont ( positive outlier): r ∈ [+0.74, +0.78]
- Remove all 5 outlier states: r ∈ [+0.73, +0.79]
Conclusion: within-YRBSS coupling is not driven by any single high-leverage state.
2. Outlier classification robustness
outlier classifications under three threshold sensitivities (z > |1.5|, |1.7|, |2.0|):
- MS: residual_z = -1.64 → |z|>1.5: negative_outlier
- PR: residual_z = -1.60 → |z|>1.5: negative_outlier
- TN: residual_z = -1.58 → |z|>1.5: negative_outlier
- VT: residual_z = 2.18 → |z|>1.5: positive_outlier / |z|>1.7: positive_outlier / |z|>2.0: positive_outlier
- WV: residual_z = -1.60 → |z|>1.5: negative_outlier
3. NUCC taxonomy whitelist sensitivity
A narrower whitelist excluding generic-counselor and generic-social-worker codes would reduce broad/narrow ratio toward AAMC's expected 5-15×. Queued for.1.
4. AADR vintage sensitivity
CDC WONDER XML POST API provides post-2017 state-level AADR..1 will test whether across-instrument correlations (need × AADR, access × AADR) are sensitive to AADR vintage.
Appendix D — Codebook
Codebook — this study
Version: · 2026-05-19
Variable definitions for the Phase B outcomes-correlation dataset.
State-identification fields
state_abbr(str, 2-char): USPS 2-letter state abbreviation; valid set = 50 U.S. states + DC + PRgeography_id(str): same as state_abbr in
Framework inputs
need_value(float, %): state-level YRBSS 2023 sad/hopeless 2+ wks prevalence among grades 9-12 (Total demographic). Range 24-47%.covariate_value(float, %): state-level ACS 1-year 2023 under-19 uninsured percentage. Range 2.6-13.0%.access_value(float, providers per 100K under-18): state-level count of NPIs matching NUCC taxonomy whitelist (45 codes) at practice-location state, divided by ACS under-18 × 100K. range: 200-5,700.provider_count_v3_taxonomy_whitelist(int): raw NPI count per practice-location statetotal_pop(int): ACS 1-year 2023 under-18 population per state
Framework outputs
log_gap_ratio(float, log-units): log(access_value / (need_value × 100))residual_raw(float): residual from OLS regression of log_gap_ratio on covariate_value (national pool)residual_z(float, z-score): residual_raw normalized by national-pool mean and stdevresidual_class(categorical): 'positive_outlier' if residual_z > +1.5; 'negative_outlier' if residual_z < -1.5; else 'expected'
Outcome variables
yrbss_sad_hopeless_2wk(float, %): same as need_value (preserved for joining)yrbss_considered_suicide(float, %): YRBSS 2023 "seriously considered attempting suicide past 12 months" Total demographicyrbss_made_plan(float, %): YRBSS 2023 "made a plan about how to attempt suicide past 12 months"yrbss_attempted_suicide(float, %): YRBSS 2023 "actually attempted suicide one or more times past 12 months"all_age_suicide_aadr_2017(float, deaths per 100K age-adjusted): NCHS Leading Causes of Death state AADR, suicide, year 2017drug_od_rate_per_100k_2023(float, deaths per 100K): VSRR drug overdose 12-month-ending Sept 2023 deaths divided by ACS 2023 state population × 100K
Metadata
study_id(str): always 'this study'computed_at(datetime UTC): processing timestampsources(list[str]): source collection names
Appendix E — Bibliography & Sources
Bibliography — this study
Version: · 2026-05-19
Primary data sources
- CDC YRBSS 2023 State-Level Data: https://www.cdc.gov/yrbss / Socrata
nu3s-3dwd - NCHS Leading Causes of Death: data.cdc.gov/bi63-dtpu (state-level AADR through 2017)
- NCHS Vital Statistics Rapid Release: data.cdc.gov/xkb8-kh2a (state-level drug overdose 12-mo)
- ACS 1-year 2023 (Census Bureau): under-18 population (B09001), under-19 uninsured pct (S2701)
- CMS NPPES Data Dissemination May 2026 V2: https://download.cms.gov/nppes/
- HRSA Health Workforce Tracker, MH HPSA: https://data.hrsa.gov/DataDownload/DD_Files/BCD_HPSA_FCT_DET_MH.csv
- NUCC Health Care Provider Taxonomy 2024-01-01: https://www.nucc.org/
Related peer-reviewed literature
- Case & Deaton (2017, 2020): Deaths of Despair
- Bishop, Press, Keyhani & Pincus (2014): Acceptance of insurance by psychiatrists. JAMA Psychiatry.
- Houtsma, Butterworth & Anestis (2018): Firearm suicide pathways. Current Opinion in Psychology.
- McBain et al. (2019): C&A psychiatrist density per 100K youth. JAMA Pediatrics 173:e191050.
- Whitney & Peterson (2019): Youth unmet mental health need. JAMA Pediatrics 173:389-91.
- Cook et al. (2017): Racial-ethnic disparities in MH access. Psychiatric Services.
- Steinberg et al. (2023): Adolescent depression-suicidality meta-analysis. Annual Review of Clinical Psychology.
Federal advisories and declarations
- U.S. Surgeon General (2021): Protecting Youth Mental Health Advisory
- AAP/AACAP/CHA (Oct 19 2021): Declaration of National Emergency in Child & Adolescent MH
- AAMC (2023): State Physician Workforce Data Report 2023, Table B-5
Trellison Institute related documents
- Adult tract-level analysis:
mh_gap_v1_*series (May 2026) - Youth state-level Phase A:
mh_gap_youth_v3_state_corrected_v1(May 2026)
Data, methods & revisions
Provider supply is measured from the full CMS NPPES registry using practice-location denominators and a NUCC taxonomy whitelist for clinicians who serve youth. State-level need-vs-access measures are reconciled against HRSA mental-health HPSA designations, AAMC physician-supply data, and KFF coverage estimates. This page presents the current version of the analysis; earlier internal drafts have been superseded. A reproducibility package (data manifest, code, codebook) is available on request.