Simulation¶
For each copula (Gaussian, Archimedean or Student), it is possible to generate samples according to copula’s structure.
Gaussian Copula Sampling¶
We suppose that our copula is a gaussian \(\mathcal{N}(0, \Sigma)\) where \(\Sigma\) is the covariance matrix. Before sampling new variables from quantile functions, we are required to compute the correlation matrix \(P\) as follow :
\[P=\sqrt{\text{Diag}(\Sigma)}^{-1} \Sigma \sqrt{\text{Diag}(\Sigma)}^{-1}\]
According to the process, we extract the Cholesky triangular matrix \(L\) such that \(P=L L^{\top}\). Finally, sampling gaussian copula is now easily possible.
- Sample \(Z \sim \mathcal{N}(0, I_d)\)
- Compute \(V=LZ\)
- Apply quantile function to \(V\)