GSobolSensitivity

class GSobolSensitivity(a=[0.0, 9.0, 99.0])

Class to define the g-Sobol’ sensitivity benchmark problem.

Methods

getFirstOrderIndices()

Returns the first order Sobol' sensitivity indices.

getFunction()

Returns the function.

getInputDistribution()

Returns the input distribution.

getName()

Returns the name of the problem.

getTotalOrderIndices()

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

g_i(x_i) = (|4 * x_i - 2.0| + a_i)/ (1 + a_i)

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.