An easy way to estimate a first-degree polynomial model is to use a factorial experiment or a fractional factorial design. This is sufficient to determine which explanatory variables affect the response variable(s) of interest. Once it is suspected that only significant explanatory variables are left, then a more complicated design, such as a central composite design can be implemented to estimate a second-degree polynomial model, which is still only an approximation at best. However, the second-degree model can be used to optimize (maximize, minimize, or attain a specific target for) the response variable(s) of interest.
Cubic designs are discussed by Kiefer, by Atkinson, Donev, and Tobias and by Hardin and Sloane.
Spherical designs are discussed by Kiefer and by Hardin and Sloane.
Mixture experiments are discussed in many books on the design of experiments, and in the response-surface methodology textbooks of Box and Draper and of Atkinson, Donev and Tobias. An extensive discussion and survey appears in the advanced textbook by John Cornell.
See also: Multiobjective optimization and Pareto efficiency
Some extensions of response surface methodology deal with the multiple response problem. Multiple response variables create difficulty because what is optimal for one response may not be optimal for other responses. Other extensions are used to reduce variability in a single response while targeting a specific value, or attaining a near maximum or minimum while preventing variability in that response from getting too large.
Response surface methodology uses statistical models, and therefore practitioners need to be aware that even the best statistical model is an approximation to reality. In practice, both the models and the parameter values are unknown, and subject to uncertainty on top of ignorance. Of course, an estimated optimum point need not be optimum in reality, because of the errors of the estimates and of the inadequacies of the model.
Nonetheless, response surface methodology has an effective track-record of helping researchers improve products and services: For example, Box's original response-surface modeling enabled chemical engineers to improve a process that had been stuck at a saddle-point for years. The engineers had not been able to afford to fit a cubic three-level design to estimate a quadratic model, and their biased linear-models estimated the gradient to be zero. Box's design reduced the costs of experimentation so that a quadratic model could be fit, which led to a (long-sought) ascent direction.34
Karmoker, J.R.; Hasan, I.; Ahmed, N.; Saifuddin, M.; Reza, M.S. (2019). "Development and Optimization of Acyclovir Loaded Mucoadhesive Microspheres by Box -Behnken Design". Dhaka University Journal of Pharmaceutical Sciences. 18 (1): 1–12. doi:10.3329/dujps.v18i1.41421. https://doi.org/10.3329%2Fdujps.v18i1.41421 ↩
Asadi, Nooshin; Zilouei, Hamid (March 2017). "Optimization of organosolv pretreatment of rice straw for enhanced biohydrogen production using Enterobacter aerogenes". Bioresource Technology. 227: 335–344. Bibcode:2017BiTec.227..335A. doi:10.1016/j.biortech.2016.12.073. PMID 28042989. https://www.researchgate.net/publication/311881656 ↩
Box & Wilson 1951 - Box, G.E.P.; Wilson, K.B. (1951). "On the Experimental Attainment of Optimum Conditions". Journal of the Royal Statistical Society, Series B. 13 (1): 1–45. doi:10.1111/j.2517-6161.1951.tb00067.x. https://doi.org/10.1111%2Fj.2517-6161.1951.tb00067.x ↩
Improving Almost Anything: Ideas and Essays, Revised Edition (Wiley Series in Probability and Statistics) George E. P. Box ↩