I have two tasks:
- Given country's CDS spread draw implied probability of default.
- Given probability of default calculate CDS spread.
If possible, refer to any papers.
I have two tasks:
If possible, refer to any papers.
Risk-neutral default probability implied from CDS is approximately $P=1-e^\frac{-S * t}{1-R}$, where $S$ is the flat CDS spread and $R$ is the recovery rate.
The CDS Spread can be solved using the inverse:
$$S=\ln(1-P) \frac{R-1}{t}$$
Hulls equation is a gross simplification. This equation is not perfect, but is far more accurate and works for all tenor points. It generally works well except when approaching boundary conditions (distressed credits).
I believe the answer can be further improved for all those being directed here by google after 3 years.
A common way to model the default probability is by the hazard rate. As @Bob correctly mentions, a traditional requirement is for it to satisfy (see Option Futures and Other Derivatives section 23.4 in which the author discusses also other more exact approximations): $$\lambda(t)=\frac{S(t)}{1-R}.$$ This is associated with the default probability by (see Poisson Process): $$P(t,t+h)=\lambda(t)h+o(h)\,,$$ with $P(t,t+h)$ the probability of a default occurring between $t$ and $t+h$. Therefore: $$P(0,T)=\int_0^T(1-P(0,t))P(t,t+dt)=\int_0^T\lambda(t)(1-P(0,t))dt\,,$$ where the first term of the integral is "default has not occurred so far" and the second is "default occurs on the next time step". This means that $P$ satisfies: $$\frac{dP(0,t)}{dt}=\lambda(t)(1-P(0,t)).$$ If the CDS is assumed to be constant then $\lambda$ is constant and a solution would be: $$P(0,t)=1-\exp\left(\frac{-St\,\,\,}{1-R}\right).$$ Equivalently solution for the CDS is: $$S=\frac{R-1}{t}\log(1-P(0,t)).$$
The chapter in Hull on Credit Risk gives the same formula as emcor as a first approximation with a justification:
Consider first an approximate calculation. Suppose that a bond yields 200 basis points more than a similar risk-free bond and that the expected recovery rate in the event of a default is 40%. The holder of a corporate bond must be expecting to lose 200 basis points (or 2% per year) from defaults. Given the recovery rate of 40%, this leads to an estimate of the probability of a default per year conditional on no earlier default of $0.02/(1-04)$, or 3.33%. In general $$ \bar{\lambda} = \frac{s}{1-R}$$ where $\bar{\lambda}$ is the average default intensity (hazard rate) per year, $s$ is the spread of the corporate bond yield over the risk-free rate, and $R$ is the expected recovery rate.
This formula can easily be rewritten as $$s = \bar{\lambda} (1 -R)$$ as pointed out by @emcor.
User @lakesh pointed in a deleted question to a blog by Donald van Deventer that analyses this formula and he rejects it.
Both Hull and van Deventer remark that this formula is an imperfect approximation. It makes sense but there are some caveats and a number of improvements can be made and Hull gives one you can readily do yourself. What is best probably depends on the goal of the study.
R
I start to have some questions. Why is it that as R->0 (which should imply more risk?) does lambda become lower. If I assume 0 recovery, i have 100% loss at default. As R->1 lambda->large. Is this not unintuitive? I may be bungling something obvious.
$\endgroup$
From this research from Deutsche Bank, : $$p_{def}=\frac{CDS_{spread}}{1-Rec}$$ $$\Leftrightarrow CDS_{spread}=p_{def}(1-Rec)$$ where $Rec$ is the recovery rate in case of default.