Copula function
- Copula tạo phân phối đồng thời cho nhiều biến bằng cách ánh xạ các phân phối biên chưa biết sang một phân phối đã biết (ví dụ chuẩn đa biến), giữ nguyên phân phối biên và định nghĩa tương quan ở không gian mới.
- C: [0,1]n → [0,1]; C[G1(u1), …, Gn(un)] = Fn[F1−1(G1(u1)), …, Fn−1(Gn(un)); ρF].
- Gi(ui): phân phối biên chưa biết dạng; Fn: phân phối tích lũy đồng thời; Fi−1: hàm ngược dùng để ánh xạ; ρF: ma trận tương quan của Fn.
- Phổ biến quanh năm 2000 (một hàm cho CDO 100+ tài sản), thất sủng khi khủng hoảng 2007–2009 bắt đầu.
Gaussian copula
- Ánh xạ từng phân phối biên sang phân phối chuẩn tắc (trung bình 0, độ lệch chuẩn 1) theo percentile-to-percentile: percentile thứ 5 → −1.645.
- CG[G1(u1), …, Gn(un)] = Mn[N1−1(G1(u1)), …, Nn−1(Gn(un)); ρM].
- Gaussian default time copula: CGD[Q1(t), …, Qn(t)] = Mn[N1−1(Q1(t)), …, Nn−1(Qn(t)); ρM], với Qi(t) là xác suất vỡ nợ tích lũy.
- Hai tài sản: CGD[QB(t), QC(t)] = M2[N−1(QB(t)), N−1(QC(t)); ρ], một hệ số ρ, không cần ma trận.
- Ví dụ B và C, ρ = 0.4: QB(1) = 0.065 → −1.514; QC(1) = 0.238 → −0.713; P(cùng vỡ nợ trong một năm) = 3.4%, so với 1.55% nếu độc lập.
Default time tương quan
- Không có lời giải dạng đóng nên phải mô phỏng.
- Cholesky decomposition của ρM để rút mẫu tương quan Mn(•) ∈ [0,1]; gán Mn(•) = Qi(τi) và giải τi (Excel hoặc Newton-Raphson).
- Lặp lại khoảng 100,000 lần để ước lượng default time kỳ vọng; ví dụ mẫu 25% ứng với τ = 3.5 năm.