Actuarium

Dependence and Copulas

Advanced
13 min read·Risk & Reinsurance
On this page
Key formulas
Sklar's theorem
F(x,y)=C(FX(x),FY(y))F(x,y) = C(F_X(x), F_Y(y))
Gaussian copula
C(u,v)=Φρ(Φ1(u),Φ1(v))C(u,v) = \Phi_\rho\big(\Phi^{-1}(u), \Phi^{-1}(v)\big)
Clayton copula
C(u,v)=(uθ+vθ1)1/θC(u,v) = \big(u^{-\theta}+v^{-\theta}-1\big)^{-1/\theta}
Kendall's tau (Clayton)
τ=θθ+2\tau = \frac{\theta}{\theta+2}
Upper tail dependence
λU=limu1P(V>uU>u)\lambda_U = \lim_{u\to1^-} P\big(V>u \mid U>u\big)

Aggregating risk across lines of business, perils, or geographies requires modeling not just each risk's marginal distribution but how they move together — and correlation alone is a poor, sometimes misleading, summary of dependence, especially in the tails that drive capital. Copulas separate the marginal behavior of each risk from the dependence structure linking them.

Sklar's theorem

Sklar's theorem (1959) states that any joint distribution FF with continuous marginals FX,FYF_X, F_Y can be written as

F(x,y)=C(FX(x),FY(y))F(x,y) = C\big(F_X(x), F_Y(y)\big)

for a unique function C:[0,1]2[0,1]C:[0,1]^2\to[0,1] called a copula — itself a joint CDF of two Uniform(0,1) random variables. This decomposition is powerful precisely because it is a decomposition: any marginal distributions can be paired with any copula, letting an actuary model, say, a heavy-tailed lognormal severity and a light-tailed frequency with an arbitrary dependence structure between two lines of business, entirely independent of the choice of marginals.

Major copula families

Gaussian copula: C(u,v)=Φρ(Φ1(u),Φ1(v))C(u,v) = \Phi_\rho(\Phi^{-1}(u), \Phi^{-1}(v)), the copula implicit in a bivariate normal with correlation ρ\rho. It has zero tail dependence — extreme joint events become vanishingly unlikely relative to what a fat-tailed dependence structure would imply — a well-documented shortcoming when used for credit and catastrophe risk aggregation (implicated in the 2008 financial crisis mispricing of CDO tranches).

t-copula: the copula of a multivariate Student-t distribution, with parameters ρ\rho and degrees of freedom ν\nu. Unlike the Gaussian, it has symmetric, positive tail dependence in both tails, controlled by ν\nu (lower ν\nu → fatter tails → stronger tail dependence), making it a common choice when joint extreme moves are plausible but no directional asymmetry is expected.

Clayton copula: C(u,v)=(uθ+vθ1)1/θC(u,v) = (u^{-\theta}+v^{-\theta}-1)^{-1/\theta}, θ>0\theta>0. Exhibits lower tail dependence only — variables tend to crash together but do not necessarily boom together — a natural fit for joint default/loss events (e.g., correlated reserve deterioration across lines in a soft market, or correlated large losses in a recession).

Gumbel copula: C(u,v)=exp{[(lnu)θ+(lnv)θ]1/θ}C(u,v) = \exp\{-[(-\ln u)^\theta+(-\ln v)^\theta]^{1/\theta}\}, θ1\theta\ge1. Exhibits upper tail dependence only, appropriate for perils that tend to produce joint extreme large losses (e.g., regional catastrophe losses across correlated property exposures).

Kendall's tau and calibration

Because copula parameters are not directly comparable to Pearson correlation, actuaries typically calibrate via Kendall's tau τ\tau, a rank-based (and hence marginal-distribution-free) dependence measure, which has closed forms for common families:

  • Gaussian/t: τ=2πarcsin(ρ)\tau = \dfrac{2}{\pi}\arcsin(\rho)
  • Clayton: τ=θθ+2\tau = \dfrac{\theta}{\theta+2}
  • Gumbel: τ=11θ\tau = 1 - \dfrac{1}{\theta}

For example, an empirical Kendall's tau of τ=0.4\tau = 0.4 between two lines' annual loss ratios implies a Clayton θ\theta solving 0.4=θ/(θ+2)θ=0.4(θ+2)θ0.4θ=0.80.6θ=0.8θ=1.3330.4 = \theta/(\theta+2) \Rightarrow \theta = 0.4(\theta+2) \Rightarrow \theta - 0.4\theta = 0.8 \Rightarrow 0.6\theta = 0.8 \Rightarrow \theta = 1.333, versus a Gumbel θ\theta solving 0.4=11/θθ=1/0.6=1.6670.4 = 1-1/\theta \Rightarrow \theta = 1/0.6 = 1.667. The same τ\tau implies materially different tail behavior depending on which family is chosen — a reminder that matching a single rank correlation does not pin down the dependence structure.

Tail dependence

The upper tail dependence coefficient

λU=limu1P(V>uU>u)\lambda_U = \lim_{u\to1^-} P(V>u \mid U>u)

measures the probability that one variable is extreme given the other is, in the limit. λU=0\lambda_U=0 for Gaussian (asymptotic independence in the tail, regardless of ρ<1\rho<1); λU=221/θ\lambda_U = 2-2^{1/\theta} for Gumbel; Clayton has λU=0\lambda_U=0 but λL=21/θ>0\lambda_L = 2^{-1/\theta}>0. This asymmetry — a copula can look similarly correlated in the body of the distribution yet behave completely differently in the tail — is exactly why tail dependence, not correlation, should drive the choice of copula for capital and reinsurance-aggregation purposes.

Simulating from a copula

A standard algorithm to simulate a pair (X,Y)(X,Y) with copula CC and given marginals:

  1. Simulate (U,V)(U,V) from the copula CC (e.g., for Gaussian/t, simulate correlated normals/t-variates and apply Φ\Phi/t-CDF to get uniforms; for Archimedean copulas such as Clayton/Gumbel, use the conditional-distribution or Marshall–Olkin/Laplace-transform algorithms specific to the family).
  2. Apply the inverse marginal CDFs: X=FX1(U)X = F_X^{-1}(U), Y=FY1(V)Y = F_Y^{-1}(V).

This two-step separation is the practical payoff of Sklar's theorem: simulate dependence and marginals independently, then combine.

Stress-testing dependence

Because correlation/copula parameters are themselves highly uncertain (data windows are short relative to how infrequently true tail-correlation events occur), sound practice is to stress the dependence assumption itself, not just the marginals:

  • Refit or override ρ\rho/θ\theta at levels seen in the worst historical periods (e.g., financial crisis correlations, which spike well above long-run averages — "correlation breakdown").
  • Swap families entirely (e.g., replace an assumed Gaussian copula with a Gumbel or t-copula at matched Kendall's tau) and re-measure capital/TVaR to quantify model risk, not just parameter risk.
  • Run a "perfect correlation" (ρ=1\rho=1 or comonotonic) scenario as an upper bound and independence as a lower bound, bracketing the plausible range of aggregate capital.

A capital model that only stresses marginal severities while holding the copula fixed is testing half the problem; dependence assumptions typically drive more of the diversification credit — and its potential erosion in a real stress — than any single marginal parameter.

Pitfalls

  • "Correlation" without specifying which measure (Pearson vs. Kendall vs. Spearman) or which copula it is meant to calibrate.
  • Assuming Gaussian dependence for tail-risk aggregation (e.g., aggregating cat and casualty capital) when the true joint tail behavior is fatter.
  • Fitting a copula to a short or benign historical window and extrapolating that same dependence into a stress scenario.
  • Ignoring that the same Kendall's tau can come from very different tail dependence depending on copula family, as shown above.

Exam relevance

Copulas and dependence measures are tested on CAS Exam 8 and SOA's ERM/QFI Advanced tracks.

Further reading

  • Sklar (1959), Fonctions de répartition à n dimensions et leurs marges
  • Nelsen, An Introduction to Copulas
  • McNeil, Frey & Embrechts, Quantitative Risk Management: Concepts, Techniques and Tools

Related

References

  • Sklar (1959), Fonctions de répartition à n dimensions
  • Nelsen, An Introduction to Copulas
  • McNeil, Frey & Embrechts, Quantitative Risk Management

Want a deeper conversation on this topic?

Ask the agent about this
Ask the tutor