Addendum to Rubtsov and Petrov (2016): “A point-in-time–through-the-cycle approach to rating assignment and probability of default calibration”

Torsten Pyttlik; Mark Rubtsov; Alexandre Petrov

In June 2016, The Journal of Risk Model Validation published a paper by Rubtsov and Petrov (2016) called “A point-in-time–through-the-cycle approach to rating assignment and probability of default calibration”. On p. 102 of the paper, the authors solved a system of equations (5.7)–(5.9) numerically; these equations are reproduced below as (1)–(3):

$\displaystyle\mathbb{E}[\varPhi^{-1}(d_{r})]$	$\displaystyle=\frac{\mu_{r}-\sqrt{\rho_{r}}\,\mathbb{E}(\hat{Z})}{\sqrt{1-\rho% _{r}}},$	(1)
$\displaystyle\mathbb{E}[(\varPhi^{-1}(d_{r}))^{2}]$	$\displaystyle=\frac{1}{1-\rho_{r}}[\mathbb{E}(B_{r}^{2})-2\sqrt{\rho_{r}}\,% \mathbb{E}(B_{r}\hat{Z})+\rho_{r}\mathbb{E}(\hat{Z}^{2})]$
	$\displaystyle=\frac{1}{1-\rho_{r}}[(\gamma_{r}+\mu_{r}^{2})-2\sqrt{\rho_{r}}% \mu_{r}\mathbb{E}(\hat{Z})+\rho_{r}\mathbb{E}(\hat{Z}^{2})],$	(2)
$\displaystyle\mathbb{E}[(\varPhi^{-1}(d_{r}))^{3}]$	$\displaystyle=(1-\rho_{r})^{-3/2}[\mathbb{E}(B_{r}^{3})-3\sqrt{\rho_{r}}\,% \mathbb{E}(B_{r}^{2}\hat{Z})+3\rho_{r}\mathbb{E}(B_{r}\hat{Z}^{2})-\rho_{r}^{3% /2}\mathbb{E}(\hat{Z}^{3})]$
	$\displaystyle=(1-\rho_{r})^{-3/2}[(3\mu_{r}\gamma_{r}+\mu_{r}^{3})-3\sqrt{\rho% _{r}}(\mu_{r}^{2}+\gamma_{r})\mathbb{E}(\hat{Z})+3\rho_{r}\mu_{r}\mathbb{E}(% \hat{Z}^{2})-\rho_{r}^{3/2}\mathbb{E}(\hat{Z}^{3})].$	(3)

Torsten Pyttlik has recently proposed an analytical solution to this system, and we present the details of that solution below. We believe it adds substantial extra value to the original material.

Let $Y_{r}:=\varPhi^{-1}(d_{r})$ for brevity. The original equations (5.7)–(5.9) then become

$\displaystyle\mathbb{E}[Y_{r}]$	$\displaystyle=\frac{\mu_{r}-\sqrt{\rho_{r}}\,\mathbb{E}[\hat{Z}]}{\sqrt{1-\rho% _{r}}},$	(4)
$\displaystyle\mathbb{E}[Y_{r}^{2}]$	$\displaystyle=\frac{1}{1-\rho_{r}}[\gamma_{r}+\mu_{r}^{2}-2\sqrt{\rho_{r}}\mu_% {r}\mathbb{E}[\hat{Z}]+\rho_{r}\mathbb{E}[\hat{Z}^{2}]],$	(5)
$\displaystyle\mathbb{E}[Y_{r}^{3}]$	$\displaystyle=(1-\rho_{r})^{-3/2}[\mu_{r}^{3}+3[\mu_{r}\gamma_{r}-\sqrt{\rho_{% r}}(\mu_{r}^{2}+\gamma_{r})\mathbb{E}[\hat{Z}]+\rho_{r}\mu_{r}\mathbb{E}[\hat{% Z}^{2}]]-\rho_{r}^{3/2}\mathbb{E}[\hat{Z}^{3}]].$	(6)

Rearranging (4) gives

\mu_{r}=\sqrt{1-\rho_{r}}\mathbb{E}[Y_{r}]+\sqrt{\rho_{r}}\,\mathbb{E}[\hat{Z}].

(7)

Taking the square of (4) and subtracting that from (5) and then rearranging gives us

	$\displaystyle\mathbb{E}[Y_{r}^{2}]-\mathbb{E}[Y_{r}]^{2}$	$\displaystyle=\frac{1}{1-\rho_{r}}[\gamma_{r}+\rho_{r}(\mathbb{E}[\hat{Z}^{2}]% -\mathbb{E}[\hat{Z}]^{2})],$		(8)
	$\displaystyle\gamma_{r}$	$\displaystyle=(1-\rho_{r})\mathbb{V}[Y_{r}]-\rho_{r}\mathbb{V}[\hat{Z}].$		(9)

Here, we have introduced the variance, defined as

\mathbb{V}[X]:=\mathbb{E}[X^{2}]-\mathbb{E}[X]^{2}.

Note that (9) might result in $\gamma_{r}<0$ if $\rho_{r}>0$ , which is undesirable since $\gamma_{r}$ was defined as a variance when the original system of equations was set up. Negative values of $\rho_{r}$ could therefore be considered, which would require an extensive modification of (1)–(3), using $\smash{\sqrt{-\rho_{r}}}$ and changing signs in several places.

Taking the third power of (4) and subtracting this from (6) gives

\mathbb{E}[Y_{r}^{3}]-\mathbb{E}[Y_{r}]^{3}=(1-\rho_{r})^{-3/2}[3[\mu_{r}% \gamma_{r}-\sqrt{\rho_{r}}\gamma_{r}\mathbb{E}[\hat{Z}]+\rho_{r}\mu_{r}\mathbb% {V}[\hat{Z}]]-\rho_{r}^{3/2}(\mathbb{E}[\hat{Z}^{3}]-\mathbb{E}[\hat{Z}]^{3})].

Inserting (7) and (9) into the inner square brackets on the right-hand side yields, after rearranging, an expression that is solvable for $\rho_{r}$ alone:

	$\displaystyle\mathbb{E}[Y_{r}^{3}]-3\mathbb{E}[Y_{r}]\mathbb{V}[Y_{r}]-\mathbb% {E}[Y_{r}]^{3}$	$\displaystyle=-\bigg(\frac{\rho_{r}}{1-\rho_{r}}\bigg)^{\!3/2}[\mathbb{E}[\hat% {Z}^{3}]-3\mathbb{E}[\hat{Z}]\mathbb{V}[\hat{Z}]-\mathbb{E}[\hat{Z}]^{3}],$
	$\displaystyle\rho_{r}$	$\displaystyle=\frac{1}{[1+(\mathbb{S}[Y_{r}]/\mathbb{S}[\hat{Z}])^{-2/3}]}.$		(10)

Here, we have defined

\mathbb{S}[X]:=\mathbb{E}[X^{3}]-3\mathbb{E}[X]\mathbb{V}[X]-\mathbb{E}[X]^{3},

which is the nonnormalized skewness (to obtain normalized skewness, multiply $\mathbb{S}[X]$ by $\mathbb{V}[X]^{-3/2}$ ).

After evaluating $\rho_{r}$ from (10), use (9) and (7) to obtain values for $\gamma_{r}$ and $\mu_{r}$ , respectively.

Note that if the distribution of $\smash{\hat{Z}}$ is symmetrical, ie, $\mathbb{S}[\smash{\hat{Z}}]=0$ , then (10) has no solution if $\mathbb{S}[Y_{r}]\neq 0$ . There is no unique solution if both $\mathbb{S}[\smash{\hat{Z}}]=0$ and $\mathbb{S}[Y_{r}]=0$ . For the limiting case $\rho_{r}=1$ , the whole system of equations (4)–(6) would be invalid.

References

Rubtsov, M., and Petrov, A. (2016). A point-in-time–through-the-cycle approach to rating assignment and probability of default calibration. The Journal of Risk Model Validation 10(2), 83–112 (http://doi.org/bzcb).