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

learn more… | top users | synonyms

3
votes
1answer
192 views

Practical implementation of Least Squares Monte Carlo (tweaks and pittfalls)

The Longstaff-Schwartz LSM approach is nowadays ubiquitous(at least in the academic literature) in pricing path dependant derivatives. Up to now I have mostly worked with lattice methods. My ...
3
votes
2answers
158 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 ...
0
votes
1answer
300 views

How do I simulate stock prices for a 10 asset portfolio, over a period of 10 years in MATLAB? [closed]

If I have given vectors for return and volatility (i.e. I have two 1x10 vectors), and I assume at first that their correlation is 0 (meaning my covariance-variance matrix is just diagonal), how do I ...
1
vote
2answers
216 views

Basket Option weight sensitivity calculation

I am looking to find/estimate the "greeks"/option price sensitivities/derivatives for a basket option situation. In specific the change in price of a put option associated with a change in weight of a ...
4
votes
3answers
1k 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 ...
2
votes
1answer
326 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 ...
4
votes
1answer
243 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 ...
1
vote
0answers
295 views

MonteCarlo simulation of stock prices using milstein scheme with dividend yield?

While performing a montecarlo simulation of stock prices using the milstein scheme is it possible to take into account the dividend yield into the simulation itself somehow, if we are given a ...
1
vote
0answers
41 views

Weak convergence of Lookback payoff with correction term

In this article on the Multilevel Monte Carlo method on page 8, http://people.maths.ox.ac.uk/gilesm/files/mcqmc06.pdf, Giles uses a correction term to improve the weak convergence rate of the lookback ...
7
votes
2answers
295 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 ...
2
votes
0answers
85 views

Optimizing stochastic functions numerically

Is there an efficient and commonly used optimization method for "more complex" investment strategies. For instance, say you have a function $f(X_1,...,X_n,c,v)$ where the $X_k$'s are your random ...
1
vote
1answer
435 views

How to explain the path dependency in binomial tree model to price options?

I'm new to quantitative finance, so I'm confused with the so-called path dependency in binomial tree model. Originally I thought the path dependency exists because in binomial tree model, we will ...
1
vote
1answer
183 views

Control variate for Heston model

Does anyone have suggestions for potential control variates for vanillas in a Heston model? I've tried black scholes with implied volatility, average volatility and long term volatility all without ...
1
vote
0answers
95 views

What are the industry standard models for monte carlo simulation of basket options?

I would like to simulate an equity index, a risk free cash account and the yield curve for the purposes of valuing a guarantee on an insurance product that is being backed by both equities and cash. ...
4
votes
1answer
231 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 ...
0
votes
1answer
2k views

Stock prices using a monte carlo simulation with a normal inverse gauss distribution

I am supposed to model daily stock prices with a normal inverse gauss distribution in excel. I feel like I am misssing some basics because I cant transform the information from the academic papers ...
6
votes
2answers
977 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. ...
1
vote
0answers
1k views

Models for simulating FX movements

My goal is to develop a model to simulate long term FX movements. (I am not sure if long term makes any difference, but if it does I am more interested in long term fx movements) These Monte Carlo ...
9
votes
1answer
425 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 ...
1
vote
1answer
396 views

Heston MC Simulations - Speed up in Matlab

At the moment I am running a Quad Core Xeon PC with 12GB of RAM doing crude MC with 10k scenarios and 1000 time steps. And using fminsearch for calibration, and it takes about half an hour to an hour ...
3
votes
1answer
818 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 ...
5
votes
1answer
2k 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 ...
1
vote
1answer
237 views

Greeks of Basket

I am considering a product composed of 10 underlying assets. The maturity is 5 year. Each year if the performance of the equi-weighted portfolio reach a barrier, it pays a coupon. My question concern ...
2
votes
1answer
1k 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 ...
4
votes
2answers
902 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 ...
8
votes
1answer
336 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 ...
7
votes
2answers
994 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 + ...
3
votes
4answers
911 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 ...
5
votes
2answers
421 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 ...
4
votes
2answers
505 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 ...
6
votes
5answers
979 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 ...
8
votes
2answers
1k 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 ...
9
votes
2answers
3k 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 ...
5
votes
1answer
599 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 ...
4
votes
2answers
250 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 ...
5
votes
4answers
2k 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 ...
7
votes
5answers
3k 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 ...
7
votes
1answer
349 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}}$ ...
6
votes
2answers
596 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
139 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 ...
7
votes
1answer
573 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 ...
8
votes
1answer
135 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 ...
5
votes
1answer
614 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, ...
7
votes
1answer
350 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
638 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 ...
8
votes
3answers
670 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?
1
vote
2answers
2k views

Black Scholes and Monte Carlo implementations in Java [duplicate]

Possible Duplicate: Is there an all Java options-pricing library (preferably open source) besides jquantlib? Can anyone recommend a library with an implementation of Black Scholes and Monte ...
7
votes
2answers
577 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)? ...
15
votes
1answer
1k 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 ...
4
votes
1answer
310 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) ...