“…Hence, its applicability is limited in situations involving sparse and noisy training data (Vohra et al, 2020). The parametric method mainly include the high-order response surface method (HRSM) (Gavin and Yau, 2008), polynomial chaos expansion (PCE) (Ghanem and Spanos, 1992;Xiu and Karniadakis, 2002;Zhang and Xu, 2021;SalehiA et al, 2018;Xu et al, 2017;Wang and Matthies, 2018), polynomial dimensional decomposition (PDD) (Rahman, 2008(Rahman, , 2009(Rahman, , 2011(Rahman, , 2015Yadav and Rahman, 2013;Ren et al, 2016;Tang et al, 2016Tang et al, , 2019Lu, 2018), etc. In the HRSM, the Chebyshev polynomials are used to determine the polynomial degree of each variable, and the constitution of the response surface is obtained based on an assumed criterion, and then the coefficients in the response surface is evaluated by uniformly distributed random sample points.…”