Credit Risk Models
Why Credit Risk Models?
Credit risk — the possibility that a borrower will fail to meet obligations — is one of the oldest and most significant risks in finance. Unlike market risk, which involves price fluctuations of liquid assets, credit risk involves discrete events (default) with potentially large losses and complex recovery processes.
Modern credit risk modeling became essential after the 1990s, driven by Basel regulations, the growth of credit derivatives, and major defaults like Enron and Lehman Brothers. These models are crucial for loan pricing, portfolio management, regulatory capital, and trading credit derivatives.
Types of Credit Risk
1. Default Risk
The probability that a borrower will fail to make required payments.
2. Migration Risk
The risk that a borrower's creditworthiness deteriorates (rating downgrade) without actual default.
3. Recovery Risk
Uncertainty in the amount recovered if default occurs.
4. Concentration Risk
Risk arising from lack of diversification in the credit portfolio.
5. Counterparty Risk
Credit risk in derivatives and other bilateral contracts.
Structural Models
Merton Model (1974)
Framework : Default occurs when firm value falls below debt level.
Setup :
Firm value: V t = V 0 e ( r − σ 2 2 ) t + σ W t V_t = V_0 e^{(r-\frac{\sigma^2}{2})t + \sigma W_t} V t = V 0 e ( r − 2 σ 2 ) t + σ W t
Debt: Face value D D D maturing at time T T T
Default: Occurs if V T < D V_T < D V T < D
Equity Value :
E 0 = V 0 N ( d 1 ) − D e − r T N ( d 2 ) E_0 = V_0 N(d_1) - De^{-rT} N(d_2) E 0 = V 0 N ( d 1 ) − D e − r T N ( d 2 )
where:
d 1 = ln ( V 0 / D ) + ( r + σ 2 2 ) T σ T , d 2 = d 1 − σ T d_1 = \frac{\ln(V_0/D) + (r + \frac{\sigma^2}{2})T}{\sigma\sqrt{T}}, \quad d_2 = d_1 - \sigma\sqrt{T} d 1 = σ T ln ( V 0 / D ) + ( r + 2 σ 2 ) T , d 2 = d 1 − σ T
Default Probability :
P ( default ) = N ( − d 2 ) \mathbb{P}(\text{default}) = N(-d_2) P ( default ) = N ( − d 2 )
Extensions of Merton Model
1. First-Passage Models
Default can occur
anytime when firm value hits barrier
B < D B < D B < D :
τ = inf { t ≥ 0 : V t ≤ B } \tau = \inf\{t \geq 0 : V_t \leq B\} τ = inf { t ≥ 0 : V t ≤ B }
Default Probability :
P ( τ ≤ T ) = N ( ln ( B / V 0 ) + ( r − σ 2 2 ) T σ T ) + ( B V 0 ) 2 r / σ 2 − 1 N ( ln ( B / V 0 ) − ( r − σ 2 2 ) T σ T ) \mathbb{P}(\tau \leq T) = N\left(\frac{\ln(B/V_0) + (r-\frac{\sigma^2}{2})T}{\sigma\sqrt{T}}\right) + \left(\frac{B}{V_0}\right)^{2r/\sigma^2-1} N\left(\frac{\ln(B/V_0) - (r-\frac{\sigma^2}{2})T}{\sigma\sqrt{T}}\right) P ( τ ≤ T ) = N ( σ T ln ( B / V 0 ) + ( r − 2 σ 2 ) T ) + ( V 0 B ) 2 r / σ 2 − 1 N ( σ T ln ( B / V 0 ) − ( r − 2 σ 2 ) T )
2. Jump-Diffusion Extensions
d V t = r V t d t + σ V t d W t + V t − ( e J − 1 ) d N t dV_t = rV_t dt + \sigma V_t dW_t + V_{t-}(e^J - 1)dN_t d V t = r V t d t + σ V t d W t + V t − ( e J − 1 ) d N t
This captures sudden deteriorations in firm fundamentals.
3. Stochastic Barrier Models
The default barrier itself is stochastic:
d B t = μ B B t d t + σ B B t d W t B dB_t = \mu_B B_t dt + \sigma_B B_t dW_t^B d B t = μ B B t d t + σ B B t d W t B
Calibration Challenges
Unobservable firm value : Must be inferred from equity prices
Volatility estimation : Equity volatility ≠ asset volatility
Parameter stability : Model parameters change over time
Low default frequencies : Difficult to validate empirically
Reduced-Form Models
Hazard Rate Approach
Framework : Model the intensity (hazard rate) of default directly, without modeling firm value.
Default Time :
τ \tau τ has intensity process
λ t \lambda_t λ t
Survival Probability :
P ( τ > t ) = e − ∫ 0 t λ s d s \mathbb{P}(\tau > t) = e^{-\int_0^t \lambda_s ds} P ( τ > t ) = e − ∫ 0 t λ s d s
Default Probability :
P ( τ ≤ t ) = 1 − e − ∫ 0 t λ s d s \mathbb{P}(\tau \leq t) = 1 - e^{-\int_0^t \lambda_s ds} P ( τ ≤ t ) = 1 − e − ∫ 0 t λ s d s
Constant Hazard Rate
If λ t = λ \lambda_t = \lambda λ t = λ (constant):
Survival function : S ( t ) = e − λ t S(t) = e^{-\lambda t} S ( t ) = e − λ t
Default density : f ( t ) = λ e − λ t f(t) = \lambda e^{-\lambda t} f ( t ) = λ e − λ t
Exponential distribution for default times
Stochastic Hazard Rates
Cox-Ingersoll-Ross (CIR) Process
d λ t = κ ( θ − λ t ) d t + σ λ t d W t d\lambda_t = \kappa(\theta - \lambda_t)dt + \sigma\sqrt{\lambda_t} dW_t d λ t = κ ( θ − λ t ) d t + σ λ t d W t
Features :
Mean reversion : High default risk tends to decrease
Non-negative : Square-root diffusion keeps λ t ≥ 0 \lambda_t \geq 0 λ t ≥ 0
Affine structure : Analytical bond pricing formulas
Vasicek Process
d λ t = κ ( θ − λ t ) d t + σ d W t d\lambda_t = \kappa(\theta - \lambda_t)dt + \sigma dW_t d λ t = κ ( θ − λ t ) d t + σ d W t
Features :
Gaussian : Can become negative (less realistic)
Analytical tractability : Closed-form survival probabilities
Correlated Defaults
For portfolio models, hazard rates can be correlated:
d λ i , t = κ i ( θ i − λ i , t ) d t + σ i d W i , t d\lambda_{i,t} = \kappa_i(\theta_i - \lambda_{i,t})dt + \sigma_i dW_{i,t} d λ i , t = κ i ( θ i − λ i , t ) d t + σ i d W i , t
where d W i , t d W j , t = ρ i j d t dW_{i,t} dW_{j,t} = \rho_{ij} dt d W i , t d W j , t = ρ ij d t .
Credit Derivatives Pricing
Credit Default Swaps (CDS)
Structure : Insurance against default of reference entity.
Premium Leg : Periodic payments of spread
s s s times notional
Protection Leg :
( 1 − R ) × Notional (1-R) \times \text{Notional} ( 1 − R ) × Notional if default occurs
Pricing Equation :
Spread × PV of Premium Leg = ( 1 − R ) × Default Probability × Discount Factor \text{Spread} \times \text{PV of Premium Leg} = (1-R) \times \text{Default Probability} \times \text{Discount Factor} Spread × PV of Premium Leg = ( 1 − R ) × Default Probability × Discount Factor
Risky Bond Pricing
Bond Price :
P ( 0 , T ) = E [ e − ∫ 0 τ ∧ T r s d s 1 { τ > T } + e − ∫ 0 τ r s d s R 1 { τ ≤ T } ] P(0,T) = \mathbb{E}\left[e^{-\int_0^{\tau \wedge T} r_s ds} \mathbf{1}_{\{\tau > T\}} + e^{-\int_0^\tau r_s ds} R \mathbf{1}_{\{\tau \leq T\}}\right] P ( 0 , T ) = E [ e − ∫ 0 τ ∧ T r s d s 1 { τ > T } + e − ∫ 0 τ r s d s R 1 { τ ≤ T } ]
Components :
Survival : Full repayment if no default
Default : Recovery R R R if default occurs
Credit Spread
Credit Spread :
s = y risky − y risk-free s = y_{\text{risky}} - y_{\text{risk-free}} s = y risky − y risk-free
Relationship to Default Probability :
s ≈ λ ( 1 − R ) s \approx \lambda (1-R) s ≈ λ ( 1 − R )
for small spreads, where λ \lambda λ is the hazard rate.
Portfolio Credit Risk Models
Gaussian Copula Model
Framework : Joint defaults driven by correlated normal variables.
Setup :
Y i = ρ X + 1 − ρ Z i Y_i = \sqrt{\rho} X + \sqrt{1-\rho} Z_i Y i = ρ X + 1 − ρ Z i
where X , Z i ∼ N ( 0 , 1 ) X, Z_i \sim \mathcal{N}(0,1) X , Z i ∼ N ( 0 , 1 ) independent.
Default : Firm
i i i defaults if
Y i < Φ − 1 ( P D i ) Y_i < \Phi^{-1}(PD_i) Y i < Φ − 1 ( P D i )
Portfolio Loss :
L = ∑ i = 1 N ( 1 − R i ) E A D i 1 { default i } L = \sum_{i=1}^N (1-R_i) EAD_i \mathbf{1}_{\{\text{default}_i\}} L = ∑ i = 1 N ( 1 − R i ) E A D i 1 { default i }
One-Factor Model
Vasicek (1987) : Industry standard for regulatory capital.
Asset Value :
A i = ρ X + 1 − ρ ε i A_i = \sqrt{\rho} X + \sqrt{1-\rho} \varepsilon_i A i = ρ X + 1 − ρ ε i
Default : If
A i < Φ − 1 ( P D i ) A_i < \Phi^{-1}(PD_i) A i < Φ − 1 ( P D i )
Conditional Default Probability :
P ( default i ∣ X = x ) = Φ ( Φ − 1 ( P D i ) − ρ x 1 − ρ ) \mathbb{P}(\text{default}_i | X = x) = \Phi\left(\frac{\Phi^{-1}(PD_i) - \sqrt{\rho} x}{\sqrt{1-\rho}}\right) P ( default i ∣ X = x ) = Φ ( 1 − ρ Φ − 1 ( P D i ) − ρ x )
Multi-Factor Models
Extend to multiple systematic risk factors:
A i = ∑ k = 1 K ρ i , k X k + 1 − ∑ k = 1 K ρ i , k ε i A_i = \sum_{k=1}^K \sqrt{\rho_{i,k}} X_k + \sqrt{1-\sum_{k=1}^K \rho_{i,k}} \varepsilon_i A i = k = 1 ∑ K ρ i , k X k + 1 − k = 1 ∑ K ρ i , k ε i
Applications :
Industry factors : Technology, financial, energy
Geographic factors : US, Europe, Asia
Macroeconomic factors : GDP, interest rates
Copula Models
Definition
A
copula C ( u 1 , … , u n ) C(u_1, \ldots, u_n) C ( u 1 , … , u n ) links univariate marginal distributions to form a multivariate distribution:
F ( x 1 , … , x n ) = C ( F 1 ( x 1 ) , … , F n ( x n ) ) F(x_1, \ldots, x_n) = C(F_1(x_1), \ldots, F_n(x_n)) F ( x 1 , … , x n ) = C ( F 1 ( x 1 ) , … , F n ( x n ))
Popular Credit Copulas
1. Gaussian Copula
C Gauss ( u 1 , … , u n ) = Φ n ( Φ − 1 ( u 1 ) , … , Φ − 1 ( u n ) ; Σ ) C^{\text{Gauss}}(u_1, \ldots, u_n) = \Phi_n(\Phi^{-1}(u_1), \ldots, \Phi^{-1}(u_n); \Sigma) C Gauss ( u 1 , … , u n ) = Φ n ( Φ − 1 ( u 1 ) , … , Φ − 1 ( u n ) ; Σ )
where Φ n \Phi_n Φ n is the multivariate normal CDF with correlation matrix Σ \Sigma Σ .
2. t-Copula
C t ( u 1 , … , u n ) = t ν , Σ ( t ν − 1 ( u 1 ) , … , t ν − 1 ( u n ) ) C^t(u_1, \ldots, u_n) = t_{\nu,\Sigma}(t_\nu^{-1}(u_1), \ldots, t_\nu^{-1}(u_n)) C t ( u 1 , … , u n ) = t ν , Σ ( t ν − 1 ( u 1 ) , … , t ν − 1 ( u n ))
Features :
Heavy tails : More realistic extreme dependence
Tail dependence : Non-zero probability of joint extreme events
3. Archimedean Copulas
Generated by function φ \varphi φ :
C ( u 1 , … , u n ) = φ − 1 ( φ ( u 1 ) + ⋯ + φ ( u n ) ) C(u_1, \ldots, u_n) = \varphi^{-1}(\varphi(u_1) + \cdots + \varphi(u_n)) C ( u 1 , … , u n ) = φ − 1 ( φ ( u 1 ) + ⋯ + φ ( u n ))
Examples :
Clayton : Lower tail dependence
Gumbel : Upper tail dependence
Frank : Symmetric dependence
Calibration
Maximum Likelihood :
max θ ∑ t = 1 T ln c ( u 1 , t , … , u n , t ; θ ) \max_{\theta} \sum_{t=1}^T \ln c(u_{1,t}, \ldots, u_{n,t}; \theta) θ max t = 1 ∑ T ln c ( u 1 , t , … , u n , t ; θ )
where c c c is the copula density.
Credit Risk Metrics
Value at Risk (VaR)
Credit VaR : Maximum expected loss due to credit events at confidence level
α \alpha α :
CVaR α = F L − 1 ( α ) − E [ L ] \text{CVaR}_\alpha = F_L^{-1}(\alpha) - \mathbb{E}[L] CVaR α = F L − 1 ( α ) − E [ L ]
Expected Shortfall (ES)
ES α = E [ L ∣ L ≥ VaR α ] \text{ES}_\alpha = \mathbb{E}[L | L \geq \text{VaR}_\alpha] ES α = E [ L ∣ L ≥ VaR α ]
Economic Capital
Definition : Capital needed to absorb losses at high confidence level (e.g., 99.9%) over one year.
Basel III :
EC = K × EAD \text{EC} = K \times \text{EAD} EC = K × EAD where
K K K is the capital ratio.
Unexpected Loss
UL = Var ( L ) \text{UL} = \sqrt{\text{Var}(L)} UL = Var ( L )
Components :
UL 2 = ∑ i = 1 N UL i 2 + 2 ∑ i < j ρ i j UL i UL j \text{UL}^2 = \sum_{i=1}^N \text{UL}_i^2 + 2\sum_{i<j} \rho_{ij} \text{UL}_i \text{UL}_j UL 2 = i = 1 ∑ N UL i 2 + 2 i < j ∑ ρ ij UL i UL j
Monte Carlo Simulation
Basic Algorithm
Generate systematic factors : X ∼ N ( 0 , 1 ) X \sim \mathcal{N}(0,1) X ∼ N ( 0 , 1 )
Generate idiosyncratic shocks : Z i ∼ N ( 0 , 1 ) Z_i \sim \mathcal{N}(0,1) Z i ∼ N ( 0 , 1 )
Compute asset values : A i = ρ X + 1 − ρ Z i A_i = \sqrt{\rho} X + \sqrt{1-\rho} Z_i A i = ρ X + 1 − ρ Z i
Determine defaults : If A i < Φ − 1 ( P D i ) A_i < \Phi^{-1}(PD_i) A i < Φ − 1 ( P D i )
Calculate portfolio loss : L = ∑ defaults ( 1 − R i ) E A D i L = \sum_{\text{defaults}} (1-R_i) EAD_i L = ∑ defaults ( 1 − R i ) E A D i
Repeat for many scenarios
Variance Reduction
Importance Sampling
Sample from modified distribution to increase default probability in tail scenarios.
Stratified Sampling
Divide systematic factor space into strata and sample within each.
Antithetic Variables
Use pairs ( X , − X ) (X, -X) ( X , − X ) to reduce variance.
Stress Testing
Regulatory Stress Tests
CCAR (US) : Comprehensive Capital Analysis and Review
EBA (EU) : European Banking Authority stress tests
Scenarios :
Baseline : Expected economic conditions
Adverse : Recession scenario
Severely adverse : Financial crisis scenario
Methodologies
Historical scenarios : Replicate past crises
Hypothetical scenarios : Expert judgment on tail risks
Monte Carlo : Simulate many possible scenarios
Machine Learning in Credit Risk
Credit Scoring
Traditional :
Logistic regression , linear discriminant analysis
Modern :
Random forests , gradient boosting, neural networks
Features
Financial ratios : Leverage, profitability, liquidity
Market data : Stock prices, volatility, CDS spreads
Alternative data : Social media, transaction data
Macroeconomic : GDP, unemployment, interest rates
Model Validation
Out-of-sample testing : Performance on unseen data
Backtesting : Historical performance evaluation
Stress testing : Performance under adverse conditions
Model interpretability : Regulatory requirements
Regulatory Framework
Basel III
Key Components :
Minimum capital ratios
Risk-weighted assets calculation
Leverage ratio
Liquidity requirements
Internal Ratings-Based (IRB) Approach
Banks develop own models for:
PD : Probability of default
LGD : Loss given default
EAD : Exposure at default
Capital Requirement :
K = ( L G D × N ( 1 1 − R N − 1 ( P D ) + R 1 − R N − 1 ( 0.999 ) ) − P D × L G D ) × M K = (LGD \times N(\sqrt{\frac{1}{1-R}}N^{-1}(PD) + \sqrt{\frac{R}{1-R}}N^{-1}(0.999)) - PD \times LGD) \times M K = ( L G D × N ( 1 − R 1 N − 1 ( P D ) + 1 − R R N − 1 ( 0.999 )) − P D × L G D ) × M
Connection to Other Topics
Credit risk models integrate multiple quantitative areas:
Built on Stochastic Differential Equations for firm value
Use Jump-Diffusion Processes for sudden deteriorations
Apply Martingale pricing for derivatives
Connect to Interest Rate Models for discounting
Foundation for fixed-income derivatives
Link to survival analysis in statistics
Enable sophisticated risk management and regulatory compliance