Optimal Designs for Rational Function Regression
针对含异方差噪声的多项式和有理函数回归模型，提出统一方法求解D、E、A和Φp最优设计，通过优化模型确定支撑点，算法高效且数值稳定。
We consider the problem of finding optimal nonsequential designs for a large class of regression models involving polynomials and rational functions with heteroscedastic noise also given by a polynomial or rational weight function. Since the design weights can be found easily by existing methods once the support is known, we concentrate on determining the support of the optimal design. The proposed method treats D-, E-, A-, and Φ p -optimal designs in a unified manner, and generates a polynomial whose zeros are the support points of the optimal approximate design, generalizing a number of previously known results of the same flavor. The method is based on a mathematical optimization model that can incorporate various criteria of optimality and can be solved efficiently by well-established numerical optimization methods. In contrast to optimization-based methods previously proposed for the solution of similar design problems, our method also has theoretical guarantee of its algorithmic efficiency; in concordance with the theory, the actual running times of all numerical examples considered in the paper are negligible. The numerical stability of the method is demonstrated in an example involving high-degree polynomials. As a corollary, an upper bound on the size of the support set of the minimally supported optimal designs is also found.