Abstract
Global sensitivity analysis (GSA) can help practitioners focusing on the inputs whose uncertainties have an impact on the model output, which allows reducing the complexity of the model. Screening, as the qualitative method of GSA, is to identify and exclude non- or less-influential input variables in high-dimensional models. However, for non-parametric problems, there remains the challenging problem of finding an efficient screening procedure, as one needs to properly handle the non-parametric high-order interactions among input variables and keep the size of the screening experiment economically feasible. In this study, we design a novel screening approach based on analysis of variance decomposition of the model. This approach combines the virtues of run-size economy and model independence. The core idea is to choose a low-level complete orthogonal array to derive the sensitivity estimates for all input factors and their interactions with low cost, and then develop a statistical process to screen out the non-influential ones without assuming the effect-sparsity of the model. Simulation studies show that the proposed approach performs well in various settings.
Disclosure statement
No potential conflict of interest was reported by the authors.
Notes
1. A complete OA, also called a complete factorial, is a design consisting of all possible runs in which the first factor takes any level between 1 and
, the second factor takes any level between 1 and
, etc.
2. In the ANOVA, the sum of squares of the th, …,
th factors, denoted as
, is a quadratic form of Y and there exists a unique symmetric matrix
such that
. The matrix
is called the MI of columns
of H.
3. An OA H is called having strength d if each submatrix of H contains all possible
row vectors with the same frequency.
4. A design is called saturated if there are only enough degrees of freedom to estimate the effects specified in the model.