DeepAI AI Chat
Log In Sign Up

Weighted Monte Carlo with least squares and randomized extended Kaczmarz for option pricing

by   Damir Filipović, et al.

We propose a methodology for computing single and multi-asset European option prices, and more generally expectations of scalar functions of (multivariate) random variables. This new approach combines the ability of Monte Carlo simulation to handle high-dimensional problems with the efficiency of function approximation. Specifically, we first generalize the recently developed method for multivariate integration in [arXiv:1806.05492] to integration with respect to probability measures. The method is based on the principle "approximate and integrate" in three steps i) sample the integrand at points in the integration domain, ii) approximate the integrand by solving a least-squares problem, iii) integrate the approximate function. In high-dimensional applications we face memory limitations due to large storage requirements in step ii). Combining weighted sampling and the randomized extended Kaczmarz algorithm we obtain a new efficient approach to solve large-scale least-squares problems. Our convergence and cost analysis along with numerical experiments show the effectiveness of the method in both low and high dimensions, and under the assumption of a limited number of available simulations.


page 1

page 2

page 3

page 4


Pricing high-dimensional Bermudan options with hierarchical tensor formats

An efficient compression technique based on hierarchical tensors for pop...

Theory and construction of Quasi-Monte Carlo rules for option pricing and density estimation

In this paper we propose and analyse a method for estimating three quant...

A Note on Monte Carlo Integration in High Dimensions

Monte Carlo integration is a commonly used technique to compute intracta...

Composite Spatial Monte Carlo Integration Based on Generalized Least Squares

Although evaluation of the expectations on the Ising model is essential ...

Exploiting Low Rank Covariance Structures for Computing High-Dimensional Normal and Student-t Probabilities

We present a preconditioned Monte Carlo method for computing high-dimens...

A Generalization of Spatial Monte Carlo Integration

Spatial Monte Carlo integration (SMCI) is an extension of standard Monte...