A different view of ridge analysis from numerical optimization

This research falls into the gap between applied statistics and numerical optimization in a specific topic – Ridge Analysis (RA). This article proposes using the trust-region (TR) methods in numerical optimization to solve the RA problem, arising from the literature of response surface methodology (RSM) in applied statistics, where its goal is to help engineers for ‘process improvement’ to find the better response value of the predicted response function within the boundary of experimentation. In the field of numerical optimization, as the family of TR approaches always exhibits excellent mathematical properties during optimization steps, thus the algorithm presented in this study guarantees global optima for the RA problem. Two examples found in the RSM literature are included to illustrate the algorithm, demonstrating its capability of locating better operating conditions than existing computing methods and pointing out particular circumstances (termed the ‘hard case’) where the classical RA procedure fails. An important application to the response modeling problem arising from the philosophy of Taguchi's quality engineering illustrates the hard case. Finally, the utility of the presented TR algorithm is demonstrated through a sequential framework with iterative updates of the TR model under local approximation provided that the predicted response model is a high-order or even non-polynomial function.