GSobolSensitivity¶
- class GSobolSensitivity(a=[0.0, 9.0, 99.0])¶
Class to define the g-Sobol’ sensitivity benchmark problem.
Methods
Returns the first order Sobol' sensitivity indices.
Returns the function.
Returns the input distribution.
getName
()Returns the name of the problem.
Returns the total order Sobol' sensitivity indices.
- __init__(a=[0.0, 9.0, 99.0])¶
Create the g-Sobol sensitivity problem.
The model is:
g(x) = prod_{i=0,…, d-1} g_i(x_i)
where d is the dimension and
x_i = Uniform(0.0, 1.0)
for i = 0, …, d-1. The input random variables are independent.
The default dimension is equal to 3.
- Parameters:
- asequence of floats
The coefficients of the linear sum, with length d + 1.
Notes
The dimension of the problem can be changed. The exact sensitivity indices are computed from the vector a.
The function g has no derivative at X=(1/2,…, 1/2). The function g is symmetric with respect to X=(1/2,…, 1/2).
When a[i] increases, the variable X[i] has a first order indice closer to zero.
The detailed analysis is the following: * if a[i] = 0, then the variable X[i] is ”important”, since 0 ≤ gi (x) ≤ 2. * if a[i] = 9, then the variable X[i] is ”non important”, since 0.90 ≤ gi (x) ≤ 1.10. * if a[i] = 99, then the variable X[i] is ”non significant”, since 0.99 ≤ gi (x) ≤ 1.01.
The model was first introduced in (Saltelli, Sobol’, 1995).
References
Saltelli, A., & Sobol’, I. Y. M. (1994). Sensitivity analysis for nonlinear mathematical models: numerical experience. Matematicheskoe Modelirovanie, 7(11), 16-28.
Saltelli, A., & Sobol’, I. M. (1995). About the use of rank transformation in sensitivity analysis of model output. Reliability Engineering & System Safety, 50(3), 225-239.
Marrel, A., Iooss, B., Van Dorpe, F., & Volkova, E. (2008). An efficient methodology for modeling complex computer codes with Gaussian processes. Computational Statistics & Data Analysis, 52(10), 4731-4744.
Saltelli, A., Chan, K., & Scott, E. M. (Eds.). (2000). Sensitivity analysis (Vol. 134). New York: Wiley.
Examples
>>> import otbenchmark as otb >>> problem = otb.GSobolSensitivity()
- getFirstOrderIndices()¶
Returns the first order Sobol’ sensitivity indices.
- Parameters:
- None.
- Returns:
- firstOrderIndices: ot.Point
The first order sensitivity indices.
- getFunction()¶
Returns the function.
- Parameters:
- None.
- Returns:
- function: ot.Function
The function.
- getInputDistribution()¶
Returns the input distribution.
- Parameters:
- None.
- Returns:
- distribution: ot.Distribution
The distribution.
- getName()¶
Returns the name of the problem.
- Parameters:
- None.
- Returns:
- name: str
The name.
- getTotalOrderIndices()¶
Returns the total order Sobol’ sensitivity indices.
- Parameters:
- None.
- Returns:
- totalOrderIndices: ot.Point
The total order sensitivity indices.