June 23rd, 2003, 12:25 pm
HiI have been trying to implement a Bermudan option on an underlying CDS. I have read Schonbucher'spaper "A Tree Implementation of a Credit Spread Model for Credit Derivatives". However, I mustcode the problem on a symbolic manipulation package for finance similar to Mathematica. Therefore,I had to come up with a more abstract scheme to price the product, rather than just code the tree modelas specified by Schonbucher.The underlying process I am taking is the Hull-White process for the hazard ratedlambda=(k(t)-a*lambda)dt +sigma*dW. Note that I consider interest rates asdeterministic (although with term structure).Following Schonbucher, we know that for given \tau>0, the density of the time of thefirst default as seen from time zero is for T>0Pd(0,t)=Expect Value[\lambda(t) e^{-int_0^t lambda(\tau) d\tau}]so this defines, via the Feynman-Kac theorem a backwards Kolmogorov equation problemfor Pd(0,t) which is\frac{\partial Pd(t,\lambda)}{\partial t}+1/2\sigma^2\frac{\partial^2 Pd(t,\lambda)}{\partial \lambda^2}+[k(t)-a\lambda]{\partial Pd(t,\lambda)}{\partial \lambda}=\lambda Pd(t,lambda)with the terminal data Pd(T,\lambda)=\lambda. Similarly for the survival density,Ps(0,t)=Expect Value[e^{-int_0^t lambda(\tau) d\tau}]So the PDE is the same but the terminal data Ps(T,\lambda)=1.This is similar to the derivation of the Black-Scholes formula, however, I am not totally sure one can apply the Feynman-Kac since the discount term e^{-int_0^t lambda(\tau) d\tau} contains lambda as well, while for Black-Scholes the discount term does not depende on the stock price.So my first question to anyone is, is this right, can I use these two PDE’s to backwards integrate and solve the Pd(0,t) and Ps(0,t) probabilities? If so that makes my life easier, as the symbolic manipulation package for finance I am using allows me to simple write the PDE in the form above and solves it without me having to write any code. Second question is, in order to value the option on a cds, I exercise if the ratio between the sum default payments = (1-Recovery) \sum_{i=1}^{N Periods} DiscountFactor_i * Pd(0,t) and the sum coupon payments = \sum_{i=1}^{N Periods} K DiscountFactor_i * Ps(0,t), where K is the fixed rate coupon, is bigger than K. If the ration is bigger than K, exercise,otherwise do not.Third and last question, how to calibrate k(t), in the Hull-White process for the hazard ratedlambda=(k(t)-a*lambda)dt +sigma*dW? I think “a” and sigma we will get from historical data andgood guesses, but k(t) should be such that the Pd(0,t) at coupon payment dates matches the probability of default coming from stripping cds spread rates in the market, right?Hope someone will be able to clarify these three questions for me. If so, thanks in advance!Eurico