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

learn more… | top users | synonyms

4
votes
1answer
44 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 ...
2
votes
1answer
60 views

American Swaption Pricing with Monte-Carlo method

I want to price an American swaption but I am not sure about what I am doing. Tree methods and PDE discretization seem difficult to adapt to a swaption. I am trying a Monte-Carlo approach. (in ...
1
vote
2answers
68 views

Efficient numerical approaches for pricing American Options with multiple sources of noise

I am looking for efficient numerical approaches for pricing American options when two or more sources of noise are involved (the simplest case coming to mind would be the Heston Model) Eventhough I ...
7
votes
2answers
371 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?
2
votes
1answer
84 views

Does one use the covariance or correlation matrix in cholesky decomposition to generate correlated samples

Can we interchangeably use Cholesky decomposition of covariance and correlation matrix to generate simulations? If not, in which situations do we use one or the other and why? Thanks in advance.
2
votes
1answer
135 views

Distribution of Geometric Brownian Motion

Please let me know where I have been mistaken! Let the SDE satisfied by the GBM $S(t)$ be $$ \frac{dS(t)}{S(t)} = \mu dt + \sigma dW(t). $$ Then, the underlying BM $X(t)$ will satisfy $$ dX(t) = ...
2
votes
0answers
74 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
95 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
70 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
87 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 ...
2
votes
3answers
182 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
141 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 ...
3
votes
0answers
105 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
120 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
34 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 ...
8
votes
2answers
119 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 ...
3
votes
0answers
75 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
186 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 ...
2
votes
1answer
88 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
78 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
155 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 ...
3
votes
0answers
216 views

A question on Monte Carlo method

I am studying Monte Carlo method by reading An Introduction to Financial Option Valuation and my questions come from the last paragraph of the section 15.2 on page 144. The author, Desmond J.Higham, ...
0
votes
1answer
301 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 ...
7
votes
2answers
588 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
417 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 ...
10
votes
1answer
264 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
212 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
370 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
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 ...
1
vote
1answer
120 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
639 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
463 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 ...
7
votes
1answer
212 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
2answers
454 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
514 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 ...
4
votes
2answers
283 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 ...
2
votes
1answer
327 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
540 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
765 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 ...
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 ...
5
votes
1answer
426 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
222 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
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 ...
7
votes
1answer
332 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
495 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 ...
7
votes
1answer
411 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
124 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
479 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, ...