I am trying to compute the probability density function of the forward rate implied by the SABR formula approximation in order to see how the density implied by the approximation has negative probabilities or explosive behaviour depending on the input parameters (especially for low strikes).

All I have is the stochastic SABR process and the closed form solutions with the first few terms of the asymptotic expansion of the term lognormal volatility.

How would I go about computing the PDF for a range of points? What I want to achieve is something like this:


  • $\begingroup$ compute the undiscounted call $E[(S_T-K)^+]$ using the SABR approximation then obtain the density from $d(K)=\frac{\partial^2}{\partial K^2} E[(S_T-K)^+]$. $\endgroup$ Apr 11 '18 at 6:03
  • $\begingroup$ would that be the density of the forward? $\endgroup$
    – Iliana
    Apr 11 '18 at 12:48

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.