Benchmarking uncertainty quantification for engineering regression under graded distribution shift
Organizations: School of Computer Science and Informatics, De Montfort University, Leicester, UK
Abstract
Machine learning models in engineering are often used beyond their training conditions, where reliable uncertainty quantification (UQ) matters most. Existing comparisons of regression UQ methods rely on random splits, which assess calibration only in distribution. We introduce a three-rung shift ladder (in distribution, soft shift and hard shift) built along physically meaningful shift variables, and evaluate 12 methods from five families on seven public engineering datasets: five extrapolation settings, an unseen turbofan fault mode and a winter air-quality sensor deployment. All tuning uses source data only. In distribution, conformal methods cover at least 0.946 of test points on every dataset. Under hard shift, mean coverage falls to 0.40-0.75 across the 11 methods that always return an interval, and no systematic difference in method ranks is detected (Friedman p = 0.63). Weighted conformalised quantile regression stays near nominal on the points it answers, but abstains on 85% of hard-shift points on average and covers only 0.67 of the few Concrete points it answers. Source-validation gains predict improvement in distribution and under soft shift, but not under hard shift (Spearman's rho = 0.04). Engineering evaluations should split calibration data by deployment unit, report subgroup coverage and sharpness, and treat abstention as an explicit operational outcome.
Figures & tables
| Dataset | Rows | Feat. | Target | Shift variable | Hard shift: source target | Target rows | Groups for | Deployment scenario |
| Concrete | 1,030 | 8 | Compressive strength (MPa) | Curing age (log days) | age 28 d age 28 d | 281 | 4 age quantile bands | Predicting long-term strength from early-age tests |
| Energy | 768 | 8 | Heating load (kWh/m 2 ) | Glazing area ratio | glazing 25% 40% | 240 | 4 relative-compactness bands | Assessing a more glazed design than any simulated |
| CCPP | 9,568 | 4 | Net electrical output (MW) | Ambient temperature (°C) | coolest 75% of hours hottest 25% | 2,388 | 5 relative-humidity bands | Plant output during a heatwave |
| Yacht | 308 | 6 | Residuary resistance per unit weight | Froude number | Fr 0.35 Fr 0.35 | 88 | 3 Froude bands | Hull resistance at higher speed than tank-tested |
| Naval | 11,934 | 14 | Gas-turbine (GT) compressor decay coefficient | Ship speed (knots) | 21 kn 24–27 kn | 2,652 | 3 speed bands | Condition monitoring at top-speed regimes |
| C-MAPSS | 45,351 | 15 | RUL, capped at 125 cycles | Fault mode | FD001 (HPC fault) FD003 (HPC + fan faults) | 24,720 | 5 RUL bands (0–25 … 100–125) | Prognostics under an unseen fault mode |
| Method | Family | Uncertainty represented | Interval | Fixed settings; search space (default first) |
| GPR | Bayesian | Epistemic + homoscedastic aleatoric | Type-II maximum likelihood (GPy); exact up to 2,000 rows, else sparse with 200 -means inducing points. Automatic relevance determination (ARD) kernel {radial basis function (RBF), Matérn-5/2, Matérn-3/2} | |
| BNN (VI) | Bayesian | Epistemic (weights) + heteroscedastic aleatoric | , moment-matched over 100 samples | Mean-field Gaussian posterior, Bayes-by-backprop, up to 600 epochs, patience 40. Width {64, 32, 128}; prior sd {1, 0.1, 0.5}; learning rate { , } |
| MC Dropout | Bayesian (approx.) | Epistemic + heteroscedastic aleatoric | As BNN | 100 stochastic passes, patience 30. Dropout {0.1, 0.05, 0.2, 0.3}; width {64, 128}; weight decay { , }; learning rate { , } |
| Deep Ensemble | Ensemble | Member disagreement + heteroscedastic aleatoric | Gaussian mixture, moment-matched | 5 members with Gaussian negative log-likelihood heads, patience 30. Width {64, 32, 128}; learning rate { , }; weight decay {0, } |
| NGBoost | Distributional | Heteroscedastic aleatoric | Normal distribution, up to 1,000 trees, early stopping after 50 rounds. Learning rate {0.03, 0.01, 0.1}; tree depth {3, 2, 4, 5}; minibatch fraction {1, 0.5} | |
| Quantile RF | Quantile | Conditional quantiles (aleatoric) | Empirical 2.5% and 97.5% quantiles | 300 trees. Minimum leaf {5, 1, 10, 20, 40}; features per split {all, half, } |
| PICP | PINAW | NRMSE | Winkler | ||||||
| Method | Family | IID | Soft | Min, soft | Hard | Hard | IID | Hard | rank, hard |
| GPR | Bayesian | 0.93 | 0.87 | 0.79 | 0.70 | 0.56 | 0.044 | 0.209 | 5.7 |
| BNN (VI) | Bayesian | 0.96 | 0.90 | 0.81 | 0.54 | 0.29 | 0.061 | 0.303 | 7.0 |
| MC Dropout | Bayesian | 0.96 | 0.93 | 0.85 | 0.65 | 0.31 | 0.052 | 0.229 | 4.9 |
| Deep Ensemble | Ensemble | 0.96 | 0.91 | 0.81 | 0.60 | 0.41 | 0.051 | 0.362 | 6.0 |
| NGBoost | Distributional | 0.88 | 0.81 | 0.72 | 0.40 | 0.15 | 0.050 | 0.185 | 7.9 |
| IID | Soft shift | Hard shift | |||||||
| Dataset | ESS | Abstain | PICP | ESS | Abstain | PICP | ESS | Abstain | PICP |
| Concrete | 73 | 1% | 0.95 | 38 | 26% | 0.95 | 60 | 97% | 0.67 |
| Energy | 64 | 0% | 0.98 | 12 | 54% | 0.94 | 47 | 100% | – |
| CCPP | 593 | 0% | 0.95 | 172 | 0% | 0.97 | 46 | 95% | 0.93 |
| Yacht | 19 | 31% | 0.99 | 7 | 76% | 0.90 | 10 | 100% | – |
| Naval | 594 | 0% | 0.95 | 130 | 0% | 0.94 | 86 | 100% | – |
| PICP | PINAW | ||||||
| Method | Engine | Row | Engine | Row | Engine | Row | |
| Split CP | 0.96 | 0.90 | 0.90 | 0.78 | 0.66 | 0.50 | 0.002 |
| CQR | 0.95 | 0.93 | 0.85 | 0.81 | 0.54 | 0.48 | 0.016 |
| Conformal GPR | 0.96 | 0.91 | 0.91 | 0.80 | 0.74 | 0.52 | 0.004 |
| Weighted CQR | 0.95 | 0.94 | 0.86 | 0.81 | 0.54 | 0.49 | 0.027 |
Uncertainty Quantification
Uncertainty quantification estimates how reliable predictions and scientific inferences are, and how uncertainty propagates into decisions.
See all 663 papers