Apply a transform or inverse transform on your polynomial chaos

Introduction

In this document we present the transformation involved in the creation of a polynomial chaos. Indeed, the polynomial chaos expansion is never directly applied to the input random variable of a model g. Instead, the expansion is expressed based on the associated standard random variable. In this example, we show how to transform a Log-Normal random variable into its standardized variable with the class DistributionTransformation.

Probabilistic transform

Let T be the probabilistic transform which maps the physical random variable X to the standardized random variable \xi:

\xi = T(X).

Let F_X be the Cumulative Distribution Function (CDF) associated with the random variable X and let F_\xi be the CDF associated with the random variable \xi. Therefore, F_\xi is the CDF associated with the orthogonal polynomials involved in the polynomial chaos expansion. In this case, the transform is:

\xi = F_\xi^{-1} \left( F_X(X) \right),

for any X \in \mathbb{R}.

Example

We want to use the Hermite orthogonal polynomials to expand a Log-Normal random variable:

  • let X follow the Log-Normal distribution with the following parameters: Lognormal(\mu=3 \times 10^4, \sigma=9\times 10^3),

  • let Z=\xi follow the Normal distribution with zero mean and unit standard deviation (the letter Z is often used for standard Normal random variables).

Let F_{LN} be the CDF of the Log-Normal distribution associated with X and let \Phi be the CDF of the standard Normal distribution.

Therefore,

Z = T(X) = \Phi^{-1}(F_{LN}(X)),

for any X \in \mathbb{R}.

import openturns as ot

In the first step, we define the LogNormal distribution. Its parameters - mean and standard deviation - have been selected so that there is no ambiguity with the standard Normal distribution. This parametrization can be used thanks to the ParametrizedDistribution class.

Xparam = ot.LogNormalMuSigma(3.0e4, 9.0e3, 15.0e3)  # in N
X = ot.ParametrizedDistribution(Xparam)

Then we generate a 5-point sample.

sampleX = X.getSample(5)
sampleX
X0
024851.57
139644.81
227619.95
328962.72
425575.91


In the second step, we define the random variable Z with standard Normal distribution.

Z = ot.Normal()

We also generate a sample from it.

sampleZ = Z.getSample(5)
sampleZ
X0
0-0.2627547
10.2319856
2-0.3728253
3-0.3093705
4-1.054682


In the third step, we create the transform T which maps X to Z.

T = ot.DistributionTransformation(X, Z)

We apply this transform on the sample generated from X.

T(sampleX)
y0
0-0.4809216
11.172666
2-0.03431651
30.1480268
4-0.3529745


The inverse transform maps Z into X.

Tinverse = T.inverse()
Tinverse(sampleZ)
y0
026118.46
129628.14
225460.13
325834.74
422166.9


Conclusion

The DistributionTransformation class is rarely used directly because the polynomial chaos classes perform the transformation automatically. However, understanding how this transform is done clarifies why the coefficients of the chaos expansion cannot be related to the input random variable X: the expansion is based on the standard variables \xi. Hence, the absolute values of the corresponding coefficients have no straightforward interpretation, even though squaring them yields the part of the global variance associated with each coefficient.

Total running time of the script: ( 0 minutes 0.003 seconds)