Quantitative Credit Risk · Module 07 量化信用風險 · 模組 07

The Gaussian Copula & the Anatomy of a CDO高斯 Copula 與 CDO 結構解析

Default Correlation Modeling 違約相關性建模
Li (2000)One-factor modelMonte Carlo
A copula separates the marginal behavior of each asset from the dependence structure that binds them. In credit, this is everything — a portfolio of 100 single-B obligors behaves nothing like the same 100 names if their defaults move together.Copula 將每項資產的邊際行為與束縛它們的相依結構分離開來。在信用領域,這至關重要——100 個 B 級債務人的投資組合,若其違約同步發生,行為將與各自獨立時截然不同。 — After Sklar's theorem (1959); operationalized for credit by David X. Li (2000)基於 Sklar 定理(1959);由 David X. Li(2000)應用至信用領域
01 — FOUNDATIONS 基礎概念

From bivariate normal to copula從二元常態分布到 Copula

A bivariate normal mixes the marginals and the dependence into one object. A copula isolates the dependence by feeding it uniform inputs F(x), F(y) ∈ [0,1].二元常態分布將邊際分布與相依結構混合成單一物件。Copula 則透過代入均勻分布輸入 F(x), F(y) ∈ [0,1] 來隔離相依結構。

Compare the surfaces 比較曲面

Slide ρ to see how correlation reshapes the joint behavior. Both surfaces describe the same statistical fact — but the copula is shape-agnostic.拖動 ρ 觀察相關性如何改變聯合行為。兩個曲面描述相同的統計事實——但 Copula 與形狀無關。

Why the change of variable matters 為何變數替換很重要

The bivariate normal density lives on (−∞, +∞)². Its shape changes when you swap normal marginals for t-distributions or lognormals.

The copula density lives on the unit square [0,1]². Its shape depends only on the dependence parameter ρ — not on what the marginals look like.

C(u,v) = Φρ[ Φ⁻¹(u), Φ⁻¹(v) ]

Bivariate normal vs Gaussian copula 二元常態 vs 高斯 Copula

Bivariate Normal Distribution
3-D Surface (rotate · drag)
Contour view
Gaussian Copula
3-D Surface (rotate · drag)
Contour view
3-D surface above the contour for each. Density mode shows the bivariate normal "bell" and the copula's corner spikes; distribution mode shows the smooth CDF climbing to 1.各圖上方為3D曲面。密度模式顯示二元常態「鐘形」及Copula的角落尖峰;分布模式顯示平滑CDF爬升至1。
01b — DERIVATION 推導過程

How the copula equation is built — five conceptual stagesCopula 方程式的建構——五個概念階段

From a single random variable's density to the multivariate dependence object: each step is a coordinate change designed to strip out marginal shape until only the pure dependence remains.從單一隨機變數的密度函數到多元相依物件:每一步都是座標轉換,旨在剔除邊際形狀,直到只剩下純粹的相依結構。
① Density f(x)機率密度
② CDF F(x)累積分布
③ Joint density f(x,y)聯合密度
④ Joint CDF F(x,y)聯合累積分布
⑤ Copula C(u,v)Copula 函數

Stage 1 — Probability density f(x)

Three normal densities. Shaded region under one curve = P(X ≤ x₂). Adjust threshold to see how area changes.

f(x) = (1/σ√2π) · exp[−(x−μ)²/2σ²]

Mathematics 數學推導

Sklar's Theorem (1959) — the foundation Sklar 定理——理論基礎

Every multivariate distribution can be split cleanly into its marginal distributions and a copula that captures the dependence. This is the theoretical license to do everything that follows.每個多元分布都可以被清晰地分解為邊際分布和捕捉相依性的 Copula。這是後續所有操作的理論依據。

Statement

Let H be a joint distribution function with continuous marginals FX and FY. Then there exists a unique function C: [0,1]² → [0,1] — the copula — such that for all (x,y) ∈ ℝ²:

H(x,y) = C( FX(x), FY(y) )

Conversely, if C is any copula and FX, FY are univariate distribution functions, then the function H defined above is a joint distribution function with marginals FX and FY.

Why it matters — the separation principle

Sklar lets you build any joint distribution in two completely independent steps:

Choose marginals for X and Y separately — they don't have to be from the same family. X could be t with df=8, Y could be lognormal, Z could be a generalized extreme value. Whatever fits the data.

Choose a copula to bind them together. The copula governs the dependence structure (linear correlation, tail dependence, asymmetry) without being tied to any particular marginal shape.

This separation is what makes copulas the lingua franca of credit portfolio modeling: every CDS gives you an obligor's marginal default probability, and the copula glues 125 of them into a joint distribution.

Constructing the Gaussian copula explicitly 顯式推導高斯 Copula

Apply Sklar's theorem to the bivariate normal distribution. The marginals are univariate normals; what's left over after stripping them out is the copula.將 Sklar 定理應用於二元常態分布。邊際分布為單變量常態;剝除後剩餘的部分即為 Copula。

Step 1 — Start with a bivariate standard normal

Let (X,Y) be jointly standard normal with correlation ρ. Its joint distribution function is the bivariate standard normal CDF, denoted Φρ:

Φρ(x,y) = P(X ≤ x, Y ≤ y) = ∫−∞x−∞y φρ(s,t) ds dt

The marginals are individually standard normal, FX(x) = Φ(x) and FY(y) = Φ(y).

Step 2 — Apply Sklar's identity

By Sklar, there exists a unique copula CGa such that:

Φρ(x,y) = CGa( Φ(x), Φ(y) )

Now substitute u = Φ(x) and v = Φ(y), which means x = Φ⁻¹(u) and y = Φ⁻¹(v):

Φρ( Φ⁻¹(u), Φ⁻¹(v) ) = CGa(u, v)

Step 3 — Read off the Gaussian copula

This gives the explicit formula for the bivariate Gaussian copula with parameter ρ:

CGaρ(u, v) = Φρ( Φ⁻¹(u), Φ⁻¹(v) )

where u, v ∈ [0,1] are uniform inputs (the marginal CDF values). This is exactly the equation used in the worked example below — Step 1 produces u and v, Step 2 inverts to Φ⁻¹(u) and Φ⁻¹(v), Step 3 evaluates the bivariate normal CDF.

Step 4 — Why this is "the" Gaussian copula

The result is the same for any marginals FX, FY, not just normal ones. If you want a joint distribution with t-marginals tied together by Gaussian dependence:

H(x, y) = CGaρ( Tνx(x), Tνy(y) ) = Φρ( Φ⁻¹[Tνx(x)], Φ⁻¹[Tνy(y)] )

This is what makes the equation so powerful: the same CGaρ dependence object can be glued onto any combination of marginals. The copula carries the correlation; the marginals carry the shape.

Step 5 — The copula density (for completeness)

Differentiate twice to get the copula density. By the change-of-variables formula, with x = Φ⁻¹(u), y = Φ⁻¹(v):

cGaρ(u, v) = ∂²C/∂u∂v = φρ(x, y) / [ φ(x) · φ(y) ]

which simplifies to a clean closed form:

cGaρ(u, v) = (1/√(1−ρ²)) · exp{ [2ρxy − ρ²(x² + y²)] / [2(1−ρ²)] }

The exploding peaks at (0,0) and (1,1) you see in the copula density chart above are the (1−ρ²)−½ factor combined with the corner behavior of Φ⁻¹.

02 — WORKED EXAMPLE 計算範例

Two-asset joint loss via Gaussian copula透過高斯 Copula 計算雙資產聯合損失

t-distributed returns, ρ = 0.7. Three-step recipe: marginal CDFs → normal-equivalent quantiles → bivariate normal CDF.t 分布報酬、ρ = 0.7。三步計算法:邊際 CDF → 標準常態分位數 → 二元常態 CDF。

Inputs — two risky assets 輸入——兩個風險資產

Each asset has t-distributed returns. Scale ≠ standard deviation: σ = s · √[df/(df−2)].每個資產的報酬服從 t 分布。尺度參數 ≠ 標準差:σ = s · √[df/(df−2)]。

Step-by-step calculation 逐步計算過程

Probability that both assets produce a return below the threshold.兩個資產報酬均低於門檻值的機率。

Step 1 — Marginal probabilities under each t-distribution

zx = (k − μx) / sx =
zy = (k − μy) / sy =
Fx(k) = Tνx(zx) =
Fy(k) = Tνy(zy) =

Step 2 — Map to standard normal quantiles

Φ⁻¹[Fx(k)] =
Φ⁻¹[Fy(k)] =

Step 3 — Apply bivariate standard normal CDF

C(u,v) = Φρ[ Φ⁻¹(Fx), Φ⁻¹(Fy) ]
= Φ0.70[ , ]
=
P(X≤k)
P(Y≤k)
P(X≤k ∩ Y≤k)
03 — APPLICATION TO CREDIT 信用應用

The one-factor Gaussian copula for credit baskets信用籃子的單因子高斯 Copula 模型

Vasicek (1987) / Li (2000): each obligor's "asset return" decomposes into a single systematic factor M and an idiosyncratic shock. Default if asset value falls below the default threshold.Vasicek (1987) / Li (2000):每個債務人的「資產報酬」分解為單一系統性因子 M 和特異性衝擊。當資產價值低於違約門檻時即發生違約。

The model 模型架構

Step 1 — Factor decomposition 因子分解

For each obligor i: Ai = √ρ · M + √(1−ρ) · εi, where M and εi are independent standard normals. ρ is the asset correlation — the "copula" parameter governing how defaults cluster.

Step 2 — Default trigger 違約觸發條件

Obligor i defaults over horizon T if Ai ≤ Φ⁻¹(pi), where pi is the marginal default probability (from CDS spreads or rating agency tables).

Step 3 — Conditional independence 條件獨立性

Given a realization of M, defaults are independent. Conditional default probability:

p(M) = Φ( [ Φ⁻¹(p) − √ρ · M ] / √(1−ρ) )

Bad M (negative) → high p(M) → many defaults. This is the contagion mechanism.不利的 M(負值)→ 高 p(M) → 大量違約。這就是傳染機制。

How ρ reshapes the default-count distribution ρ 如何改變違約數量分布

100-name homogeneous basket, 5-year PD = 5%. Watch the tail explode as correlation rises.100 個同質標的的籃子,5 年違約機率 = 5%。觀察隨著相關性上升,尾部如何急劇膨脹。

Distribution of the number of defaults over 5 years, by asset correlation ρ. Higher ρ ⇒ fatter right tail and bimodality emerges at extreme ρ.按資產相關性 ρ 分類的五年違約數量分布。ρ 越高 ⇒ 右尾越厚,極端 ρ 時出現雙峰分布。
04 — MONTE CARLO LABORATORY 蒙地卡羅模擬實驗室

Simulate a credit basket under the copula在 Copula 框架下模擬信用籃子

Set portfolio characteristics, run draws, observe how correlation transforms loss outcomes from binomial-like to wildly skewed.設定投資組合特性,執行抽樣,觀察相關性如何將損失結果從類二項分布轉變為極度偏態分布。

Simulation controls 模擬參數設定

Mean losses 平均損失
Std dev 標準差
95% VaR 95%風險值
99.5% VaR 99.5%風險值

Portfolio loss distribution 投資組合損失分布

Loss in % of notional. Compare independent (ρ=0) vs simulated correlated outcome.以名義本金百分比計算的損失。比較獨立情境(ρ=0)與模擬相關情境的結果。

Histogram of losses from Monte Carlo paths under the one-factor Gaussian copula.
05 — CDO TRANCHE ANALYSIS CDO 分段分析

5-year synthetic CDO tranche pricing & risk五年期合成 CDO 分段定價與風險

A CDO carves the loss distribution into vertical slices. Equity absorbs first losses, mezzanine sits in the middle, senior is "protected" — but only if correlation is well-estimated.CDO 將損失分布切割成垂直切片。股權層吸收初始損失,夾層居中,優先層「受保護」——但前提是相關性估計準確。

Pool & structure 資產池與結構設定

Tranche attachments (% of pool) 分段觸發點(%本金)

Tranche risk profile 分段風險概況

Expected tranche loss, default probability, and implied par spread (bps/yr) under the assumed correlation regime.在假設的相關性情境下,各分段的預期損失、違約機率及隱含平價利差(基點/年)。

Tranche 分段 Attach 觸發點 Detach 脫離點 Width 寬度 P(any loss) 損失機率 E[loss] 預期損失 Spread (bps) 利差
Spread ≈ E[tranche loss] / (T · width) · 10,000 bps. A first-order par approximation, ignoring discounting and timing.利差 ≈ E[分段損失] / (T · 寬度) · 10,000 基點。一階平價近似值,忽略折現和時間因素。

Pool size effect — 10 / 20 / 50 / 100 names 資產池規模效應

Same per-name PD and correlation. Smaller pools have less diversification benefit; the loss distribution stays granular and bumpier.相同的個別違約機率和相關性。較小的資產池分散效益較低,損失分布較粗糙且起伏較大。

Loss distribution over 5 years for varying pool sizes at fixed ρ and PD.

Correlation sensitivity per tranche 各分段的相關性敏感度

The "correlation smile" of structured credit: equity loves low ρ; senior hates high ρ; mezzanine is roughly correlation-neutral somewhere in between.結構性信用的「相關性微笑」:股權層偏好低 ρ;優先層懼怕高 ρ;夾層在兩者之間大致相關性中性。

Expected tranche loss vs ρ (other parameters held fixed).
06 — STRESS SCENARIOS 壓力測試情境

What happens when correlation jumps?相關性驟升時會發生什麼?

2007–2008 lesson: the senior tranche's "AAA-ness" depends critically on the assumed ρ. A jump from 0.20 to 0.50 vaporizes that protection.2007-2008 年的教訓:優先層的「AAA 品質」高度依賴假設的 ρ。ρ 從 0.20 跳升至 0.50,這種保護就煙消雲散。
Benign — ρ = 0.10 溫和情境Senior expected loss 優先層預期損失
Base — ρ = 0.20 基準情境Senior expected loss 優先層預期損失
Stressed — ρ = 0.40 壓力情境Senior expected loss 優先層預期損失
Crisis — ρ = 0.70 危機情境Senior expected loss 優先層預期損失
All four scenarios use 100 names, 5% PD, 40% recovery. The senior tranche (15%–30%) is "untouchable" only under the benign assumption.所有四個情境均使用100個標的、5%違約機率、40%回收率。優先層(15%-30%)僅在溫和假設下才「不受影響」。

Reading the table 閱讀說明

Loss thresholds and the probability the pool exceeds each, by correlation regime.在不同相關性情境下,資產池超過各損失門檻的機率。

Loss threshold 損失門檻 ρ=0.10 ρ=0.20 ρ=0.40 ρ=0.70

Practitioner's note 實務注意事項

What David Li actually proposed David Li 的原始提案

Li (2000) replaced the joint default time distribution with a Gaussian copula on the marginal survival times. Tractable, fast to calibrate, easy to plug into existing pricing systems.

Where it broke 模型失靈之處

Single ρ assumes constant pairwise dependence; the Gaussian copula has zero tail dependence — extreme co-movements are systematically underestimated. In 2007–2008, realized correlations spiked far beyond the 0.20–0.30 calibrated by base correlation models, repricing senior tranches by orders of magnitude.

Modern alternatives 現代替代方案

t-copulas (positive tail dependence), Clayton/Gumbel for asymmetric tails, factor-loading copulas, random-recovery models, and direct Monte Carlo with regime-switching factor M. The plumbing is the same; the dependence kernel is richer.