Monte Carlo simulation methods uses repeated random experiments to determine results.

learn more… | top users | synonyms

24
votes
2answers
4k views

How useful is Markov chain Monte Carlo for quantitative finance?

Naively, it seems that Bayesian modeling, structural models particularly, would be quite useful in finance because of their ability to incorporate market idiosyncrasies and produce accurate ...
15
votes
1answer
957 views

Portfolio optimization with monte carlo sampling from predictive distribution

Let's say we have a predictive distribution of expected returns for N assets. The distribution is not normal. We can interpret the dispersion in the distribution as reflection of our uncertainty (or ...
10
votes
5answers
914 views

Monte carlo methods for vanilla european options and Ito's lemma.

I understand that by applying Ito's lemma to the following SDE $$dX=\mu\,X\,dt+\sigma\,X\,dW$$ one obtains a solution to the above SDE which is as follows: $${X}\left( t\right) =\mathrm{X}\left( ...
9
votes
1answer
295 views

Is creating constrained random portfolios a hard problem?

Creating random portfolios with weights $x_i$ can be thought of as sampling from the surface of a simplex given by $$Ex = f$$ and $$Ax \le b$$ Where $E$ and $A$ are constraint matrices for equality ...
8
votes
2answers
2k views

When to use Monte Carlo simulation over analytical methods for options pricing?

I've been using Monte Carlo simulation (MC) for pricing vanilla options with non-lognormal underlyings returns. I'm tempted to start using MC as my primary option-valuating technique as I can get ...
8
votes
3answers
551 views

Reference on Markov chain Monte Carlo method for option pricing?

I have to implement option pricing in c++ using Markov chain Monte Carlo. Is there some paper which describes this in detail so that I can learn from there and implement?
8
votes
1answer
127 views

Simulating property price index

I am trying to write a Monte Carlo simulation to calculate risk associated with some property based products. What is the most reasonable stochastic process to model property price index? Do people ...
7
votes
4answers
1k views

Methods for pricing options

I'm looking at doing some research drawing comparisons between various methods of approaching option pricing. I'm aware of the Monte Carlo simulation for option pricing, Black-Scholes, and that ...
7
votes
2answers
428 views

What are the merits of pseudo random numbers over quasi random numbers in monte-carlo simulation?

I understand that quasi-random numbers have much better convergence, but are there any reasons for me to use pseudo-random numbers instead?
7
votes
1answer
335 views

Simulating the joint dynamics of a stock and an option

I want to know the joint dynamics of a stock and it's option for a finite number of moments between now and $T$ the expiration date of the option for a number of possible paths. Let $r_{\mathrm{s}}$ ...
7
votes
2answers
500 views

Is there any research on applying state-space or dynamic linear models to forecasting equity risk premia?

Is there any research on applying state-space or dynamic linear models to forecasting equity risk premia on a security-by-security basis with a medium term horizon (say 3 month to 12 months horizon)? ...
7
votes
1answer
320 views

How to reduce variance in a Cox-Ingersoll-Ross Monte Carlo simulation?

I am working out a numerical integral for option pricing in which I'm simulating an interest rate process using a Cox-Ingersoll-Ross process. Each step in my Monte Carlo generated path is a ...
7
votes
1answer
529 views

Monte carlo portfolio risk simulation

My objective is to show the distribution of a portfolio's expected utilities via random sampling. The utility function has two random components. The first component is an expected return vector ...
7
votes
2answers
157 views

Estimation of Empirical Expected Shortfall of a heavy tailed distribution

Assume that you have a portfolio for which you have estimated a parametric model to the underlying instruments, but the distribution of the portfolio as a whole is too complicated to compute ...
7
votes
2answers
936 views

How to transform process to risk-neutral measure for Monte Carlo option pricing?

I am trying to price an option using the Monte Carlo method, and I have the price process simulations as an inputs. The underlying is a forward contract, so at all times the mean of the simulations is ...
7
votes
1answer
455 views

Simulating conditional expectations

There is a multidimensional process X defined via its SDE (we can assume that its a diffusion type process), and lets define another process by $g_t = E[G(X_T)|X_t]$ for $t\leq T$. I would like to ...
6
votes
5answers
659 views

portfolio optimization from empirical return distributions

I'd like to do a portfolio optimization of a set of ETF's but want to avoid traditional problems with normality assumptions in returns etc. Are there techniques that let me sample 'draws' from the ...
6
votes
2answers
476 views

Vanilla European options: Monte carlo vs BS formula

I have implemented a monte carlo simulation for a plain vanilla European Option and I am trying to compare it to the analytical result obtained from the BS formula. Assuming my monte carlo pricer is ...
6
votes
2answers
682 views

Principle Component Analysis vs. Cholesky Decomposition for MonteCarlo

Let's assume we have a portfolio containing large number (~500) of risk factors. We want to simulate the portfolio dynamics. PCA based simulation would be faster as we can reduce the dimensionality. ...
6
votes
2answers
526 views

How to minimize the difference between a parametric VaR and a MC-VaR with lognormal assumption?

Given that we want to find the Value at Risk for a portfolio of stocks only, there are two main methods to proceed. In the problem, we also assume that stocks follow a geometric Brownian motion. A ...
6
votes
1answer
138 views

Should we apply practical constraints on the distribution of monte carlo paths?

to limit interest rate paths to a 'reasonable' range (if we could define reasonable). Now we calibrate log-normal skew and mean reversion monthly to robust basket of atm swaptions and in and out ...
6
votes
1answer
231 views

Consistency of economic scenarios in nested stochastics simulation

I am interested in references on research regarding the consistency of economic scenarios in nested stochastics for risk measurement. Background: Pricing by Monte-Carlo: For pricing complex ...
5
votes
4answers
1k views

Stock Price Behavior and GARCH

In my (limited) understanding, the behavior of a stock price can be modeled using Geometric Brownian Motion (GBM). According to the Hull book I'm currently reading, the discrete-time version of this ...
5
votes
4answers
2k views

How to get greeks using Monte-Carlo for arbitrary option?

Let's assume I have an arbitrary option that I can price using Monte-Carlo simulation. What is the general approach (i.e. without relying on specific option type) to calculating the greeks in this ...
5
votes
2answers
611 views

Simulation of GBM

I have a question regarding the simulation of a GBM. I have found similar questions here but nothing which takes reference to my specific problem: Given a GBM of the form $dS(t) = \mu S(t) dt + ...
5
votes
3answers
131 views

Usage of Brownian Bridge?

I was recommended to read something about Brownian Bridge. Could someone familiar with BB give some recommendation? It was mentioned that BB benefits in 2 places BB could reduce the simulation ...
5
votes
1answer
479 views

Sanity check - How to price callables

This question is meant as a sanity check whether i got the workflow right for pricing callable bonds. If anyone finds a mistake, or has a suggestion, please answer. The workflow is: For every call ...
5
votes
1answer
517 views

How to apply quasi-Monte Carlo to path-dependent options?

Following up on my recent question on variance reduction in a Cox-Ingersoll-Ross Monte Carlo simulation, I would like to learn more about using a quasi-random sequence, such as Sobol or Niederreiter, ...
5
votes
1answer
1k views

Longstaff-Schwartz (Least Squares Monte Carlo) applied to American Options

I'm working on an implementation in R of Longstaff & Schwartz method from the this 2001 article. I've managed to build code that replicates their prices in table 1 (p. 127), but only for the ones ...
4
votes
1answer
879 views

Simple model for option premium (for covered call simulation)?

Given a historical distribution of weekly prices and price changes for a stock, how can I estimate the the option premium for a nearly at-the-money (ATM) option, say with an expiration date 3 months ...
4
votes
4answers
349 views

How to deal with extreme cases in normal random numbers generation?

In order to generate normal random numbers, one usually generates random numbers following a uniform distribution $Z \sim \mathcal{U}(0,1)$ and then applies the reverse CDF function on them ...
4
votes
2answers
547 views

Is drift rate the same as interest rate in risk-neutral random walk when using Monte Carlo for option pricing?

When using following risk-neutral random walk $$\delta S = rS \delta t + \sigma S \sqrt{\delta t} \phi$$ where $\phi \sim N(0,1)$. Now when a text mentions drift = 5% does that mean that interest ...
4
votes
1answer
168 views

Estimating early exercise boundary for American put

I am trying to estimate the early exercise boundary for an American put option. I can find the put value through the Longstaff-Schwartz LSM method. How do I obtain the early exercise boundary within ...
4
votes
2answers
228 views

What sort of order submission strategy would result in a random walk of trade prices?

I have written a simulation that matches buy and sell orders, keeps track of an order book and simulates trades. My first pass at order submission was to generate random orders around the bid/ask ...
4
votes
1answer
277 views

Divergence issue with my monte carlo pricer…

I am trying to implement a vanilla European option pricer with Monte Carlo and compare its result to the BS analytical formula's result. I noticed that as I increase (from 1 million to 10 millions) ...
4
votes
2answers
302 views

Generate correlated random variables from Normal and Gamma distributions

I want to generate a random vector $z$ of dimension $k+m$ with some given correlation matrix $\Sigma$, such that the first $k$ elements of the vector are distributed normally and the last $m$ elements ...
3
votes
4answers
587 views

Other means of calibrating Heston models

I understand that the simplest way of calibrating a Heston model for volatility surface is to use Monte-Carlo to simulate the vol and stock price trajectories and then use the observed price to do a ...
3
votes
1answer
2k views

Mersenne twister random number generator in Java for Monte Carlo Sim.

I am using the Mersenne twister random number generator in Java for a Monte Carlo Simulation. I need a uniform distribution of values between -1 and 1. My code is below (I am importing ...
3
votes
2answers
203 views

When to use the real world drift and when the risk neutral one for a Monte-Carlo simulation?

Under what conditions should the drift be real world and when risk neutral when simulating Delta Hedging option pricing trading strategy any other? For 2. it should be risk neutral. For 1., it ...
3
votes
3answers
474 views

Simulating the short rate in the Hull-White model

What is the best way to simulate the short rate $r(t)$ in a simple one factor Hull White process? Suppose I have $$ dr(t) = (\theta(t)-\alpha r(t))dt+\sigma dW_t $$ where $\theta(t)$ is calibrated ...
3
votes
2answers
108 views

How to Calculate a Monte Calo VaR estimation error

I'm performing a Monte Carlo to calculate value at risk (with a 3 dimension risk factor) Now, I would like to calculate the error of the estimation of the VaR with respect to the number of simulations ...
3
votes
1answer
475 views

Monte Carlo simulating Cox-Ingersoll-Ross process

The CIR process is given by the SDE $$ \mathrm dr_t = \theta(\mu-r_t)\mathrm dt + \sigma\sqrt{r_t}\mathrm dW_t $$ where $W_t$ is a Brownian motion. I am interested in finite-difference schemes of ...
3
votes
1answer
385 views

BDT model implementation

I am looking for a nice and readable description of how to implement BDT model: $d log(r(t)) = [\theta(t)-\frac{\sigma'(t)}{\sigma(t)}log(r(t))]dt + \sigma(t) dW$. I assume I already have ...
3
votes
0answers
136 views

How to compute greeks using the adjoint Monte Carlo approach?

Assume I have a stochastic ODE $$dS = a(S)dt + b(S)dW,$$ with Euler approximation $$\hat{S}_{n+1}=F_n(\hat{S}_n)=\hat{S}_n+a(\hat{S}_n)h+b(\hat{S}_n)Z_n\sqrt{h}.$$ This allows me to create sample ...
2
votes
2answers
3k views

How do I estimate convergence in monte carlo methods?

I am experimenting with Monte Carlo methods. I'd like to measure/estimate convergence with a graph/chart. How do I do that? Can anyone please direct me to relevant documentation/links or even give me ...
2
votes
1answer
214 views

Picking from two correlated distributions

Can anyone provide a simple example of picking from two distributions, such that the two generated time series give a specified value of Pearson's correlation coefficient? I would like to do this in ...
2
votes
1answer
105 views

Monte Carlo for MultiFactor Ornstein Uhlenbeck

I'm following loosely the exposition given in "Monte Carlo Methods in Financial Engineering by Glasserman. For a multifactor OU process: $dX(t)=C(b-X(t))dt+DdW(t)$ Where C and D are d*d matrices ...
2
votes
1answer
728 views

Value at Risk Monte-Carlo using Generalized Pareto Distribution(GPD)

I have created a VBA program to calculate VaR by using Monte Carlo, I have simulated Brownian Motion. This method might be ok for 100% equity portfolio, but let's say this portfolio may have fixed ...
2
votes
1answer
73 views

Fitting stochastic variance distributions to index return data

I want to calculate option prices based on a realistic distribution of the underlying. The underlying is a liquid index such as Eurostoxx50. I think of two aproaches, both of them incorporate ...
2
votes
1answer
229 views

Quasi Monte Carlo in Matlab

I want to use Quasi Monte Carlo to try and improve the convergence of a simulation I am running. The random numbers are simply to produce the observation errors for a standard linear regression ...