Background. County-level DALY disease burden estimates are rarely available for rural U.S. counties outside large academic centers. This analysis applies the DALY framework as a planning-grade educational model for seven major condition groups in Isabella County, Michigan — a rural county of 64,565 residents with a designated Medical Underserved Area status and multiple Health Professional Shortage Area (HPSA) designations. The DALY metric and disability weights are drawn from the IHME Global Burden of Disease Study — a body of work we deeply respect. This project does not replicate or claim to replicate a formal GBD estimate; it applies GBD's openly published disability weights and WHO reference standards for educational and local planning purposes only.
Methods. We constructed disability-adjusted life year (DALY) estimates from secondary data: CDC WONDER 2020–2022 (county mortality with rural adjustment factors), CDC PLACES 2024 (prevalence), IHME GBD 2021 (disability weights), and U.S. Census ACS 2022 (population denominators). YLL was calculated using the WHO GHE frontier reference life table (89·1 years) as the primary standard; Michigan observed LE (78·6 years) is reported as a planning-grade sensitivity. For mental health, a remission factor of 0·50 was applied to CDC PLACES lifetime-diagnosis prevalence to approximate active burden. Uncertainty was quantified via Monte Carlo simulation (n = 10,000 iterations). Economic burden was quantified by human capital (DALYs × Michigan GDP per capita) and value of statistical life (HHS 2026, income-adjusted).
Results. Primary burden (frontier LE, remission-adjusted MH): 15,430 DALYs (95% UI 14,459–16,481) · 9,915 YLL + 5,515 YLD. Cancer #1 (4,057), CVD #2 (3,302), SUD #3 (2,904), MH #4 (1,733). Planning-grade sensitivity (Michigan LE 78.6 yrs): 11,710 DALYs; MH remains #1 (2,713). Stroke mortality remains 49% above Michigan average. Total deaths: ~461/year. Economic burden: ~$957M/year (human capital, frontier) to ~$4.84B/year (VSL).
Limitations. Mixed-source modeled planning scenario - not a validated county burden estimate. Prevalence inputs are modeled by CDC PLACES (MLRP methodology). YLL uses simplified mean-age approach with an added age-banded sensitivity analysis; rural adjustment factors (×1.2–1.85) are applied to state-level mortality rates for suppressed county cells. Factors are derived from published rural/urban mortality differentials [8–12] and represent central estimates within reported ranges - county-specific validation data are not available. Results are directional planning approximations only; publication-grade extension requires age-sex-specific deaths, formal uncertainty propagation, and validation against claims or hospital discharge data.
Conclusion. Isabella County's burden is led by mental health access deficits (29.8% prevalence, 176,938:1 psychiatrist ratio), cancer, and substance use. Findings support telepsychiatry/collaborative care, AI-assisted cancer screening, MAT telehealth, and hypertension control as priorities. See ROI analysis → for investment return estimates.
This supplement uses the following abbreviations throughout: NCDs (non-communicable diseases, including cardiovascular disease, cancer, diabetes, and COPD), SUDs (substance use disorders), and MHCs (mental health conditions). Together these three groups account for the full disease burden estimated in this analysis.
Isabella County is a rural county in mid-Michigan (FIPS 26073), home to Central Michigan University (CMU) in Mount Pleasant. Population: 64,565 (Census ACS 2022). The presence of a large university creates an unusual demographic: substantial young adult population and elevated poverty rates under standard definitions that do not distinguish student poverty from structural community poverty.
Federal access designations: Primary Care HPSA (low-income), Dental HPSA (low-income), Mental Health HPSA (Central Michigan Service Area, multi-county), and Medically Underserved Area/Population (MUA/P) for low-income residents.
MDHHS 2020 PCNA ranks: Overall Rank 48 of 83 Michigan counties; Health Status Rank 50 of 83 (rank 1 = best). These are distinct indices - see Section 4.7. All MDHHS PCNA data are pre-pandemic baseline; post-COVID values may differ substantially.
Cause-specific mortality from CDC WONDER 2020–2022 (3-year pooled counts, Isabella County, FIPS 26073). Rural adjustment factors were applied to state-level rates for suppressed cells (<10 deaths/3-yr pool): SUD ×1.85, Mental Health ×1.40, CVD ×1.30, COPD ×1.35, Stroke ×1.25, Diabetes ×1.30, Cancer ×1.20. These factors are derived from published rural/urban mortality differentials in Michigan and nationally [8–12]; they represent central estimates within reported ranges, and county-specific validation data are not available. For all-cause comparison, MDHHS 2020 PCNA rate (819.9/100k) is used.
County mortality data for SUD and several other causes remain suppressed in CDC WONDER (cell count <10). State-level rates with rural adjustments were applied as proxies for suppressed cells. All substitutions are explicitly flagged in results.
CDC PLACES 2024 (county model-based estimates) provides primary prevalence inputs. CDC PLACES uses multilevel regression and poststratification (MLRP) applied to BRFSS data - all prevalence inputs are classified as modeled, not directly observed. Derived estimates (marked derived) used where CDC PLACES does not report the specific condition (SUD, CVD).
Disability weights: IHME GBD 2021. Moderate-severity weights applied as primary estimates.
Two reference life expectancy standards are used:
Primary (academic/international): WHO GHE frontier reference life table, sex-averaged 89·1 years (male 86·0 yr, female 92·0 yr). This is the standard used in GBD analyses globally and reflects how long a person should live given the best observed survival globally, rather than the burden-depressed observed LE of the study population. Using observed LE as the primary standard causes systematic underestimation of YLL for conditions with late mean age at death (CVD, cancer), because the reference is pulled down by the disease burden being measured.
Planning-grade sensitivity: Michigan observed life expectancy, sex-averaged 78·6 years (male 76·2 yr, female 80·9 yr, MDHHS 2024). lifeUSMI.asp → Used for direct comparison with prior county planning documents and the interactive dashboard. Under this standard, mental health ranks #1 due to its high YLD-driven burden, which is insensitive to the LE reference.
Cancer DW: Weighted composite DW ≈ 0.294 constructed from MDHHS 2020 site distribution: lung (0.279), colorectal (0.167), prostate (0.126), breast (0.049), other sites (0.400).
Total population 64,565; adult population 52,297 (≥18 yrs) (Census ACS 5-yr 2022). For conditions where CDC PLACES reports adult-only prevalence (COPD, stroke, mental health, CVD, SUD), adult population was used as the base. For population-wide conditions (diabetes, cancer), the full 64,565 was used.
Human capital: DALYs × $62,000 (Michigan GDP per capita 2024). VSL: HHS ASPE 2026 (central $13.4M), income-adjusted for Isabella County using elasticity 0.4 (Viscusi & Aldy 2003 [7]): $13.4M × (46,000/80,000)^0.4 = $10.5M per statistical life (county median income $46,000 per ACS 2022). See ROI analysis for intervention-level return estimates.
How to read each chart (Trends page). The line is age-adjusted prevalence (% of adults). The shaded band is the 95% confidence interval — the plausible range around each estimate. Dots mark years with a published county estimate; gaps mean no estimate was released that year (a data gap, not an improvement or worsening). The badge compares the earliest vs. latest reported year: ▲ Worsening · ▼ Improving · → Stable, and the line color encodes the same direction.
Abbreviations. COPD = chronic obstructive pulmonary disease · CHD = coronary heart disease · CI = confidence interval · PLACES = CDC small-area health estimates (model-based) · BRFSS = Behavioral Risk Factor Surveillance System · FIPS = federal county code · age-adjusted = standardized to a reference age distribution so counties compare fairly · DALY = disability-adjusted life year · YLL = years of life lost · YLD = years lived with disability.
Standard WHO DALY formula, no age-weighting, no discounting:
L = reference LE (89·1 yrs for primary frontier estimate; 78·6 yrs for Michigan planning sensitivity); Ā = mean age at death (condition-specific proxy). Deaths estimated as (WONDER rate × rural adj.) / 100,000 × population. Note: standard WHO/GBD YLL uses remaining life expectancy at the age of death - r(age) from a reference life table - summed over age-specific death counts, not LE at birth minus mean age at death. The formula above approximates remaining LE as L − Ā, a simplification applicable when county-level age-specific death distributions are unavailable.
P = prevalent cases (CDC PLACES 2024 prevalence × population base); DW = IHME GBD 2021 disability weight (moderate severity). Does not account for severity distributions, comorbidity adjustment, or uncertainty intervals. Incidence-based YLD would be preferable for conditions with well-characterized duration data but requires county-level incidence unavailable here.
Two methods with different conceptual frameworks: Human Capital captures productivity losses only; VSL reflects societal willingness to pay to prevent deaths. The ~6–7× difference between central estimates is methodologically expected, not an error. For grant applications use VSL (aligned with HHS/OMB); for academic GBD comparison use Human Capital.
CDC PLACES 2024 reports the prevalence of lifetime-diagnosed depressive disorder (BRFSS question ADDEPEV3: "Has a doctor, nurse, or other health professional ever told you that you have a depressive disorder, including depression, major depression, dysthymia, or minor depression?"). This is a lifetime diagnosis measure, not a current-episode measure.
For YLD calculation in the primary (frontier LE) estimate, a remission factor of 0·50 is applied to convert lifetime prevalence to active disease prevalence. Evidence basis: NSDUH 2021 reports national past-year MDE at 8·3%, while the PLACES national lifetime-diagnosed estimate is approximately 20·7% — implying approximately 60% of ever-diagnosed individuals are in remission at any given time. Published psychiatric epidemiology literature consistently reports 12-month remission rates of 40–60% for major depressive disorder. The central remission factor of 0·50 yields an active prevalence of 14·9% (range: 11·9%–17·9% for remission factors 0·40–0·60).
Active YLD (primary): 52,297 × 0·149 × 0·145 (DW) = 1,130 DALYs.
Upper-bound YLD (no remission adjustment): 52,297 × 0·298 × 0·145 = 2,260 DALYs.
The planning-grade Michigan LE sensitivity retains the raw 29·8% figure (YLD = 2,260) for continuity with v4 documents and the public dashboard.
A structured Monte Carlo simulation (n = 10,000 iterations) propagated uncertainty from three parameter sources. Prior to this version, uncertainty was reported informally as ±20% across all DALY estimates; this section replaces that approximation with simulation-based 95% uncertainty intervals.
Mortality uncertainty. For conditions with directly observed CDC WONDER death counts (cancer, CVD, COPD, stroke), annual deaths were sampled from a Poisson distribution: N ~ Poisson(λ = observed annual deaths). For proxy-estimated conditions (SUD, mental health, diabetes), where CDC WONDER cell counts are suppressed, rural adjustment factors were sampled from uniform distributions: SUD Uniform(1·65, 2·05); mental health Uniform(1·30, 1·50); diabetes Uniform(1·20, 1·40), reflecting the published ranges in Hedegaard et al. (2021) and Moy et al. (2017).
Prevalence uncertainty. CDC PLACES 2024 reports model-based 95% confidence intervals for each measure. Prevalence was sampled from Beta distributions parameterised using the method of moments to match the reported central estimate and 95% CI bounds. For mental health, an additional remission factor was sampled from Uniform(0·40, 0·60) to propagate active-prevalence uncertainty.
Disability weight uncertainty. IHME GBD 2021 disability weights carry uncertainty bounds of approximately ±15–25% at moderate severity. Disability weights were sampled from Beta distributions parameterised on the IHME-reported 95% CIs.
Reported 95% uncertainty intervals represent simulation-based percentile intervals (2·5th to 97·5th percentiles across 10,000 iterations). They do not constitute formal Bayesian posterior intervals; systematic biases in the underlying data models (e.g., MRP model error in CDC PLACES, rural adjustment factor misspecification) are not captured by this approach. Simulation implemented with NumPy seed 42.
Under the primary estimate (WHO GHE frontier LE 89·1 yrs with remission-adjusted MH prevalence), mental health ranks #4 (1,733 DALYs). Cancer is #1 (4,057), CVD #2 (3,302), and SUD #3 (2,904). Under the planning-grade Michigan LE sensitivity (raw PLACES 29·8%, no remission adjustment), mental health retains #1 (2,713 DALYs, 23·2% of total).
The rank shift from v4 to v5 under the primary estimate reflects two simultaneous changes: (1) the WHO frontier LE raises YLL more for cancer and CVD (which have late mean ages at death) than for mental health (which is YLD-dominant); (2) the remission adjustment reduces MH YLD by approximately 50%. Mental health retains the highest YLD:DALY ratio of any condition (65%) — it is the dominant source of avoidable disability regardless of LE standard. Analysts should note that CDC PLACES MH prevalence includes depression, anxiety, and other common conditions - a broader definition than severe mental illness.
Several cause-specific county mortality rates are suppressed in CDC WONDER (<10 deaths/3-yr pooled). Rural adjustment factors were applied to state-level rates: SUD ×1.85, Mental Health ×1.40, CVD ×1.30, COPD ×1.35, Stroke ×1.25, Diabetes ×1.30, Cancer ×1.20.
These factors are derived from published rural/urban mortality differentials. For drug overdose mortality, Hedegaard et al. (2021) [8] documented rural/urban rate differentials of 1.6–2.0× nationally; Mack et al. (2017) [9] reported nonmetropolitan excess for illicit drug use disorders of approximately 1.5–1.9× using MMWR WONDER data. For general rural excess mortality across leading causes, Moy et al. (2017) [10] and Garcia et al. (2017) [11] provide systematic documentation of condition-specific rural/urban differentials that informed the COPD, CVD, stroke, and diabetes adjustments. Cancer rural excess mortality (~1.2–1.4×) is documented in Zahnd et al. (2021) [12].
If MDHHS PCNA mortality rates are age-adjusted (common for geographic comparison dashboards), converting them to death counts by multiplying by county population is not mathematically valid - age-adjusted rates use a reference population structure, not the actual county structure. Death count estimates should be treated as indicative only. This should be resolved before publication-grade use.
Isabella vs. Michigan: (819.9 − 783.1) / 783.1 = +4.7%. Isabella vs. US: (819.9 − 723.6) / 723.6 = +13.3%. Earlier versions of this analysis used incorrect values (5% and 87% respectively), now corrected.
Three estimates exist for Isabella County: MDHHS 2020 PCNA 26.5% below 100% FPL; Census ACS 5-yr 2024: 21.9%; Census QuickFacts 2020–24: 19.0%. Sources differ in ACS vintage, reference year, and student-population treatment. CMU's ~14,000+ enrollment substantially inflates standard ACS poverty counts (students counted as residents in poverty without capturing parental support). All three estimates remain well above Michigan average; the range reflects measurement differences, not contradictions.
The 176,938:1 ratio is a HRSA HPSA designation metric - FTE-weighted for the Central Michigan Service Area (multi-county). It exceeds Isabella County's total population (64,565) because the denominator is an administrative service-area construction, not a county headcount. Interpret as a severe shortage designation indicator, not a literal one-psychiatrist-per-176,938 county ratio.
Health Status Rank 50/83: composite of health outcome indicators (mortality rates, disease prevalence, birth outcomes). Overall Rank 48/83: combines health status with social determinants (poverty, food insecurity, housing). These are distinct metrics from the same MDHHS 2020 PCNA - do not average or conflate them.
| Condition (ICD-10) | Prev. base | Active Prev* | Deaths/yr† | Mean age | YLL (frontier) | DW | YLD | DALYs | 95% UI | Rank | Quality |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Cancer, all sites (C00–C97) | 7.7% total | — | ~130 | 67 | 2,873 | 0.294‡ | 1,184 | 4,057 | 3,475–4,651 | #1 | WONDER 2020-22 DW comp. |
| Cardiovascular Disease (I20–I51) | 7.1% adults | — | ~178 | 72 | 3,044 | 0.070 | 258 | 3,302 | 2,857–3,781 | #2 | WONDER 2020-22 Derived |
| Substance Use Disorders (F10–F19) | 7.2% adults | — | ~37 | 44 | 1,669 | 0.329 | 1,235 | 2,904 | 2,511–3,349 | #3 | WONDER × rural adj. |
| Mental Health (F30–F48) | 29.8% adults | 14.9%* | ~14 | 46 | 603 | 0.145 | 1,130 | 1,733 | 1,396–2,133 | #4 | CDC PLACES DW: IHME |
| COPD & Respiratory (J44) | 8.2% adults | — | ~44 | 73 | 704 | 0.198 | 849 | 1,553 | 1,258–1,883 | #5 | WONDER 2020-22 CDC PLACES |
| Stroke (I60–I69) | 3.2% adults | — | ~34 | 73 | 544 | 0.316 | 529 | 1,073 | 833–1,347 | #6 | WONDER 2020-22 CDC PLACES |
| Type 2 Diabetes (E11) | 11.6% total | — | ~25 | 70 | 480 | 0.054 | 328 | 808 | 718–903 | #7 | CDC PLACES |
| Total | — | — | ~461 | — | 9,915 | — | 5,515 | 15,430 | 14,459–16,481 | — | — |
| Condition (ICD-10) | Prev. base | Deaths/yr | Mean age | YLL (MI) | DW | YLD | DALYs | Rank |
|---|---|---|---|---|---|---|---|---|
| Mental Health (F30–F48) | 29.8% adults† | ~14 | 46 | 453 | 0.145 | 2,260 | 2,713 | #1 |
| Cancer, all sites (C00–C97) | 7.7% total | ~130 | 67 | 1,509 | 0.294 | 1,184 | 2,693 | #2 |
| Substance Use Disorders (F10–F19) | 7.2% adults | ~37 | 44 | 1,280 | 0.329 | 1,235 | 2,516 | #3 |
| Cardiovascular Disease (I20–I51) | 7.1% adults | ~178 | 72 | 1,174 | 0.070 | 258 | 1,432 | #4 |
| COPD & Respiratory (J44) | 8.2% adults | ~44 | 73 | 244 | 0.198 | 849 | 1,093 | #5 |
| Stroke (I60–I69) | 3.2% adults | ~34 | 73 | 188 | 0.316 | 529 | 717 | #6 |
| Type 2 Diabetes (E11) | 11.6% total | ~25 | 70 | 218 | 0.054 | 328 | 545 | #7 |
| Total | — | ~461 | — | 5,067 | — | 6,643 | 11,710 | — |
| Condition | Mean age | YLL (MI 78.6) | YLL (Frontier 89.1) | DALYs (MI) | DALYs (Frontier) | Rank MI / Frontier |
|---|---|---|---|---|---|---|
| Mental Health | 46 | 453 | 603 | 2,713† | 1,733* | #1 / #4 |
| Cancer | 67 | 1,509 | 2,873 | 2,693 | 4,057 | #2 / #1 |
| SUD | 44 | 1,280 | 1,669 | 2,516 | 2,904 | #3 / #3 |
| CVD | 72 | 1,174 | 3,044 | 1,432 | 3,302 | #4 / #2 |
| COPD | 73 | 244 | 704 | 1,093 | 1,553 | #5 / #5 |
| Stroke | 73 | 188 | 544 | 717 | 1,073 | #6 / #6 |
| Diabetes | 70 | 218 | 480 | 545 | 808 | #7 / #7 |
| Total | — | 5,067 | 9,915 | 11,710 | 15,430 | — |
| Cause (ICD-10) | Isabella County | Michigan | United States | vs. Michigan | vs. US |
|---|---|---|---|---|---|
| All causes | 819.9 | 783.1 | 723.6 | +4.7% | +13.3% |
| IHD & HF (I20–I51) | 198.6 | 194.9 | 163.6 | +1.9% | +21.4% |
| Cancer (C00–C97) | 175.1 | 161.1 | 149.1 | +8.7% | +17.4% |
| COPD (J44) | 60.6 | 44.2 | 39.7 | +37.1% | +52.6% |
| Stroke (I60–I69) | 59.4 | 39.9 | 37.1 | +48.9% | +60.1% |
| Diabetes (E11) | 31.1 | 21.9 | 21.4 | +42.0% | +45.3% |
| Indicator | Isabella County | Michigan Avg | US Avg | Source & year |
|---|---|---|---|---|
| Below 100% FPL (MDHHS PCNA) | 26.5% | 15.0% | 14.1% | MDHHS 2020 |
| Below poverty line (Census ACS 5-yr) | 21.9% | - | 12.5% | Census ACS 2024 |
| Below poverty line (QuickFacts) | 19.0% | - | - | QuickFacts 2020–24 |
| Food insecurity | 16.4% | 13.7% | 12.5% | MDHHS 2020 |
| Adult obesity (BMI ≥30) | 43.1% | 36.7% | 35.3% | CDC PLACES 2024 |
| Current smoking (adults) | 18.6% | - | 14.0% | CDC PLACES 2024 |
| Hypertension (adults) | 32.1% | - | - | CDC PLACES 2024 |
| Maternal tobacco use | 22.0% | 14.3% | 6.5% | MDHHS 2020 |
| Uninsured (<65 yrs) | 7.9% | 6.1% | 10.6% | CDC PLACES 2024 |
| NAS rate (per 100k live births) | 583.1 | 835.8 | n/a | MDHHS 2020 |
| MH/SUD hospitalization (/100k) | 3,420 | 3,676 | 3,088 | MDHHS 2020 |
| Psychiatrist HPSA ratio‖ | 176,938:1 | 9,371:1 | - | HRSA 2020 |
| Method | Formula (frontier) | Central (frontier) | Formula (MI LE) | Central (MI planning) | Measures |
|---|---|---|---|---|---|
| Human Capital | 15,430 × $62,000 | ~$957M/yr | 11,710 × $62,000 | ~$726M/yr | Indirect productivity losses |
| VSL (HHS 2026) | 461 deaths × $10.5M | ~$4.84B/yr | 461 deaths × $10.5M | ~$4.84B/yr | Societal willingness to pay |
Under the WHO GHE frontier LE standard with remission-adjusted mental health prevalence, cancer emerges as the leading disease burden in Isabella County (4,057 DALYs), followed by CVD (3,302), SUD (2,904), and mental health (1,733). This ranking is consistent with GBD-standard methodology and enables cross-county comparison. Total primary burden is 15,430 DALYs (95% UI 14,459–16,481), with 9,915 YLL and 5,515 YLD.
Under the planning-grade Michigan LE sensitivity, mental health retains the #1 rank (2,713 DALYs, 23.2% of total) due to its high YLD-driven burden, which is insensitive to LE reference choice. This is the estimate reflected in prior county planning documents.
Mental health carries the highest YLD:DALY ratio of any condition (65% of MH burden is non-fatal disability), making it the dominant source of avoidable morbidity regardless of the LE standard applied. The CDC PLACES 29.8% estimate is a lifetime-diagnosis figure; the active disease approximation (14.9% after remission adjustment) is consistent with national NSDUH past-year MDE estimates (~15%), providing convergent validity for the adjustment.
Cancer has the highest absolute YLL under the frontier standard (2,873), reflecting 130 deaths per year at mean age 67 — giving a remaining LE of 22.1 years under frontier vs. 11.6 years under Michigan LE. This gap explains why cancer and CVD gain the most from adopting the frontier standard. Substance use disorders remain the third-highest burden and have the highest years lost per death under both standards due to their early mean age at death (44 years).
The poverty data discrepancy (19–26.5%) reflects real measurement complexity from CMU's student population, not data error. Analysts should consider a sensitivity analysis excluding the 18–24 age group to better characterize structural community poverty.
1. Age-specific YLL. Query CDC WONDER for Isabella County deaths by 5-year age group and ICD-10 cause (3–5 year pool, e.g., 2020–2024). Apply WHO reference life table for proper age-specific YLL.
2. Post-COVID mortality update. MDHHS 2020 PCNA predates pandemic and the recent opioid mortality surge. Incorporate 2023–2024 county data when available.
3. Severity-stratified YLD. Split prevalence into mild/moderate/severe using GBD severity distributions instead of single moderate-severity DW.
4. Mental health validation. The 29.8% MH prevalence (CDC PLACES 2024) warrants triangulation against BRFSS, SAMHSA NSDUH state estimates, and administrative claims data to confirm county-level estimate is not an artifact of model extrapolation.
5. Rural adjustment factor validation. Apply CDC WONDER 5-year pooled data (2019–2023) to reduce cell suppression and directly compare observed county mortality with adjusted state-level proxies. Where counts remain suppressed, consider Bayesian small-area estimation.
6. Uncertainty quantification — COMPLETED in v5.0. Monte Carlo propagation of input uncertainty — Monte Carlo simulation now propagates mortality, prevalence, and disability weight uncertainty to produce 95% simulation-based uncertainty intervals (Section 3.6). Remaining gap: systematic model biases in CDC PLACES MRP are not captured.
7. Student population sensitivity. Poverty and MH analysis with and without 18–24 CMU population to separate student vs. structural community burden.
8. GBD 2023 alignment. Published Oct 2025 in The Lancet. Verify if disability weights or reference LE changed vs. GBD 2021 and update accordingly.