Multi-Objective Optimization and Hyperparameter Tuning With Desirability Functions
Thomas Bartz-Beielstein
TL;DR
This work presents spotdesirability, a Python implementation of the desirability-function framework for multi-objective optimization, and demonstrates its use across classical optimization, surrogate-model-based optimization, and hyperparameter tuning. It translates multiple objectives into individual desirabilities and aggregates them with a geometric mean, enabling a unified objective $D = \left(\prod_{r=1}^R d_r\right)^{1/R}$ that guides optimization. The paper showcases both single- and multi-objective tuning strategies, compares direct Nelder-Mead optimization with surrogate- and desirability-based approaches, and introduces space-filling metrics via the Morris-Mitchell criterion to balance exploration and exploitation. Practical contributions include a Python ecosystem (spotdesirability, spotoptim) with concrete workflows for predictive modeling, neural-network hyperparameter tuning, and Pareto-front visualization, facilitating reproducible, multi-objective optimization in real-world settings. The results indicate that desirability-based methods can achieve competitive trade-offs, improve interpretability of objectives, and provide flexible tooling for hyperparameter tuning and design-space exploration, with formal handling of missing values, tolerances, and non-standard desirability shapes.
Abstract
The desirability-function approach is a widely adopted method for optimizing multiple-response processes. Kuhn (2016) implemented the packages desirability and desirability2 in the statistical programming language R, but no comparable packages exists for Python. The goal of this article is to provide an introduction to the desirability function approach using the Python package spotdesirability, which is available as part of the sequential parameter optimization framework. After a brief introduction to the desirability function approach, three examples are given that demonstrate how to use the desirability functions for (i) classical optimization, (ii) surrogate-model based optimization, and (iii) hyperparameter tuning. An extended Morris-Mitchell criterion, which allows the computation of the search-space coverage, is proposed and used in a fourth example to handle the exploration-exploitation trade-off in optimization. Finally, infill-diagnostic plots are introduced as a tool to visualize the locations of the infill points with respect to already existing points.
