Extract the copula from a distribution

In this example we are going to retrieve the copula from a multidimensional distribution.

from __future__ import print_function
import openturns as ot
import openturns.viewer as viewer
from matplotlib import pylab as plt
ot.Log.Show(ot.Log.NONE)

create a multivariate gaussian distribution

mu = [0.0] * 2
sigma = [1.0]*2
R = ot.CorrelationMatrix(2)
R[0, 1] = 0.25
distribution = ot.Normal(mu, sigma, R)
print(distribution)

Out:

Normal(mu = [0,0], sigma = [1,1], R = [[ 1    0.25 ]
 [ 0.25 1    ]])

extract the copula

copula = distribution.getCopula()
print(copula)

Out:

NormalCopula(R = [[ 1    0.25 ]
 [ 0.25 1    ]])