Skip to content
Prev 2932 / 15274 Next

Cumulative Multivariate Normal Distribution

Hi Ben,

Have you looked at pmvnorm  in package mvtnorm.? This is Genz's method, 
I am working with Jan-Dash on a new approximation
for cumulative normal and I have some semi-tested code of that, we are 
currently testing it for Worst-of and other rainbow type options.
http://arxiv.org/pdf/cs.CE/0611061

The alternative to price these would be to use monte-carlo to do the 
integration and just price them using that, see for example
 the attached R code which generates multi-variate paths. You can then 
use these paths to then compute the discounted payoff function.

Hope this helps,

Best
Krishna

ps:  the drift is simply the risk-neutral drift however for currencies 
you have to factor in the quanto adjustment.
Also this piece of code was written about 4 years ago so pardon the 
coding standards.

 >assetPath<-multiassetpath(c(2,3),c(-0.1,-0.2),2,matrix(c(1,0.2,0.2,1),2,2),1,100,2)
Chiquoine, Ben wrote:
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: worstof.r
URL: <https://stat.ethz.ch/pipermail/r-sig-finance/attachments/20080919/ef51ae38/attachment.pl>