Of 4,802 matched facility records, 3,045 hospitals (63.4%) appear in both source files; the remainder split between hospitals covered only by the broader Unplanned Visits measure set and a small number covered only by HRRP. The two files were reconciled on a common six-digit CCN, reshaped long-to-wide and back to long by condition, and modeled at the hospital-condition level.
Data Sources and File Construction
Two public CMS hospital-level files were used: the FY2026 HRRP Supplemental Data File and the Unplanned Hospital Visits — Hospital file. HRRP reports, for six conditions and procedures (AMI, CABG, COPD, heart failure, hip/knee replacement, and pneumonia), each hospital's excess readmission ratio — predicted-to-expected 30-day readmissions given that hospital's case mix — along with predicted and expected rates, discharge counts, and readmission counts. The Unplanned Visits file reports standardized scores, denominators, and national comparison flags across fourteen measures, including excess-days-in-acute-care (EDAC) for AMI, heart failure, and pneumonia; a hybrid hospital-wide readmission/mortality measure; and outpatient procedure follow-up measures for colonoscopy, chemotherapy, and surgery.
Facility identifiers were reconciled to a common six-digit, zero-padded CMS Certification Number; the HRRP file's numeric ID was converted to this format before merging. Non-numeric placeholder text ("Not Available," "Too Few to Report") was forced to missing during numeric conversion, reflecting CMS's low-volume suppression rules.
Each source file was reshaped long-to-wide so one row represented one hospital, with condition- or measure-specific variables distinguished by suffix. The two wide files were merged 1:1 on CCN, then the merged file was reshaped back to long format on the six HRRP conditions — producing one row per hospital per condition. Hospital-level covariates from the Unplanned Visits file (EDAC scores, the hybrid mortality measure, outpatient procedure measures) are invariant across a hospital's six condition-rows and were retained unchanged through this reshape.
Outcome and predictors
The outcome was a binary indicator of excess readmissions — coded 1 if a hospital's excess readmission ratio exceeded 1.0 for a given condition, 0 otherwise. Records with a missing ratio (typically low-volume suppression) were excluded from the corresponding model. Discharge volume was log-transformed for right skew; condition entered as a categorical predictor (AMI as reference); the eight hospital-level quality covariates were standardized (z-scored) so odds ratios are interpretable per one SD change, given their non-comparable native scales (days per 100 discharges, percentages, and standardized ratios).
Modeling Strategy
Three sets of models were estimated. First, condition-specific logistic regressions predicted excess readmission status from log discharge volume alone, fit separately for each of the six conditions. Second, a pooled mixed-effects logistic regression (melogit) was fit across all six conditions jointly — log discharge volume and condition as fixed effects, a random intercept for hospital — to test whether hospital performance correlates across conditions net of volume and condition. Third, the pooled model was re-estimated with the eight standardized quality covariates added as fixed effects, to see how much of the hospital-level random-effect variance from the second model is explained by these measured quality domains. Predicted probabilities and empirical Bayes posterior means of the random effects (BLUPs) were extracted from the final model for hospital-specific case profiling.
| Condition | N | OR | p | 95% CI |
|---|---|---|---|---|
| Heart attack (AMI) | 1,240 | 0.386 | <.001 | 0.313–0.475 |
| Bypass surgery (CABG) | 363 | 0.107 | <.001 | 0.060–0.190 |
| COPD | 1,545 | 0.549 | <.001 | 0.450–0.668 |
| Heart failure | 2,316 | 0.691 | <.001 | 0.624–0.764 |
| Hip/knee replacement | 253 | 0.049 | <.001 | 0.023–0.106 |
| Pneumonia | 2,320 | 0.818 | .001 | 0.729–0.917 |
Odds ratios below 1.0 indicate higher discharge volume associated with lower odds of excess readmissions — a volume effect strong enough to specifically flag hip/knee replacement and CABG as the most volume-sensitive programs.
| Predictor | OR | p | 95% CI |
|---|---|---|---|
| Discharge volume (log) | 0.564 | <.001 | 0.517–0.615 |
| Condition: CABG (ref. AMI) | 1.174 | .252 | 0.893–1.543 |
| Condition: COPD (ref. AMI) | 0.779 | .004 | 0.656–0.926 |
| Condition: Heart failure (ref. AMI) | 0.974 | .753 | 0.827–1.147 |
| Condition: Hip/knee (ref. AMI) | 1.099 | .563 | 0.798–1.513 |
| Condition: Pneumonia (ref. AMI) | 0.972 | .733 | 0.825–1.144 |
| Constant | 26.787 | <.001 | 16.883–42.504 |
N = 8,037 hospital-condition observations, 2,477 hospitals. Hospital-level variance = 0.962 (SE 0.096); LR test vs. standard logistic regression: χ²(1) = 246.02, p < .001 — hospital identity matters, net of volume and condition.
| Predictor | OR | p | 95% CI |
|---|---|---|---|
| Discharge volume (log) | 0.681 | <.001 | 0.596–0.778 |
| Condition: CABG (ref. AMI) | 1.256 | .151 | 0.920–1.714 |
| Condition: COPD (ref. AMI) | 0.723 | .004 | 0.578–0.904 |
| Condition: Heart failure (ref. AMI) | 0.815 | .092 | 0.642–1.034 |
| Condition: Hip/knee (ref. AMI) | 0.963 | .851 | 0.650–1.426 |
| Condition: Pneumonia (ref. AMI) | 0.856 | .190 | 0.679–1.080 |
| EDAC, heart attack (z) | 1.316 | <.001 | 1.201–1.441 |
| EDAC, heart failure (z) | 1.461 | <.001 | 1.283–1.665 |
| EDAC, pneumonia (z) | 1.555 | <.001 | 1.375–1.758 |
| Hospital-wide mortality (z) | 1.656 | <.001 | 1.528–1.795 |
| Outpatient colonoscopy follow-up (z) | 0.992 | .786 | 0.937–1.051 |
| Outpatient chemo admissions (z) | 0.992 | .806 | 0.928–1.060 |
| Outpatient chemo ED visits (z) | 0.912 | .016 | 0.846–0.983 |
| Outpatient surgery follow-up (z) | 0.997 | .925 | 0.929–1.070 |
| Constant | 9.039 | <.001 | 4.388–18.621 |
N = 3,879 hospital-condition observations, 929 hospitals. Hospital-level variance = 0.073 (SE 0.065); LR test vs. standard logistic: χ²(1) = 1.43, p = .116. "(z)" denotes a standardized (mean 0, SD 1) covariate.
| Model | N (obs) | N (hospitals) | Hospital variance | LR test vs. logistic |
|---|---|---|---|---|
| Model 1: volume + condition only | 8,037 | 2,477 | 0.962 (SE 0.096) | χ²=246.02, p<.001 |
| Model 2: adding quality covariates | 3,879 | 929 | 0.073 (SE 0.065) | χ²=1.43, p=.116 |
From Model to Hospital: Palomar Health Downtown Campus
A statistical model is only useful if it can answer a question about one real hospital, not just a population average. Palomar Health Downtown Campus in California was selected at random from hospitals with complete model data and a mixed performance profile — some conditions above the national benchmark, others below — specifically because it was representative, not because its results were unusually strong or weak.
Instead of asking whether the hospital's readmission rate is above or below benchmark, the model asks a sharper question: given this hospital's patient volume and quality-measure profile, how much readmission risk should we expect — and do the actual results line up?
Palomar Health Downtown Campus, California
Two things stand out. First, heart failure was expected to be the lowest-risk of the four modelable conditions, yet it's the one where the hospital exceeds benchmark — a clear gap between expectation and reality, and one that lines up with a specific driver identified below. Second, CMS suppresses case counts for bypass surgery and hip/knee replacement at this hospital due to low volume, so the model can't evaluate those two — their results are real, but rest on a thinner evidence base than the other four.
Is this hospital unusual?
The model separates readmission performance into measurable factors (volume, related quality scores) and whatever hospital-specific effect is left over once those are accounted for — in effect, asking whether something unique to this hospital is driving its results. For Palomar Downtown, that residual hospital-specific effect (its empirical Bayes random-effect estimate) is essentially zero and not statistically distinguishable from an average hospital. There's no sign of a hidden, unmeasured issue — its performance is fully explained by volume and the related quality measures already in the model.
What drives the hospital's risk?
Which measured factors move the needle
Patient volume is the strongest factor by a wide margin — it outweighs any single quality measure in the model. Overall hospital mortality performance is next, and works in the hospital's favor: hospitals with better mortality outcomes tend to have fewer excess readmissions, and that pattern holds here. The one factor working against the hospital is return visits after heart failure care — patients returning to the ED, observation, or inpatient care within 30 days of a heart failure discharge. That's the only measured factor raising this hospital's risk, and it's exactly what shows up in the heart-failure benchmark gap noted above. Return visits after heart attack care, by contrast, barely move this hospital's risk at all — even though they matter across hospitals generally, which is a useful reminder that an industry-level driver isn't automatically the driver for every individual hospital.
A plain-language version of this case profile, written for a non-technical hospital administrator audience, accompanies the full technical write-up in this portfolio's companion documents.