Monte Carlo simulation methods are a broad class of computational algorithms that rely on repeated random sampling to obtain numerical results.

learn more… | top users | synonyms

0
votes
1answer
14 views

Interpretation of vega out of BS formula

I am comparing Monte Carlo estimates of VaR (using importance sampling) under both the normal and student distributions. I am also considering risk factors other than log-prices; in particular, ...
1
vote
2answers
81 views

European call down and out option (geometric Brownian motion, Monte Carlo, Euler)

I need to estimate the expected value and the Greeks of an European call down and out option, assuming geometrical Brownian motion of the asset, with Monte Carlo simulation employing Euler ...
1
vote
1answer
361 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 ...
3
votes
1answer
216 views

Do you have a validation set for Libor Market Model implementation?

I'm trying to calibrate a Libor Market Model (LMM) in Matlab with my user-defined function, not their package. I already fitted the market volatilities using SABR but failed to simulate the ...
10
votes
1answer
364 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 ...
2
votes
2answers
123 views

Calculating VaR with Monte Carlo simulation

I would like some help here :) I have a problem calculating VaR with the Monte Carlo Simulation. I have followed then next steps, is this a right way to calculate VaR or I need something more? ...
4
votes
2answers
900 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
91 views

Simulate (imaginary) asset prices using random numbers that follow a Frank Copula

I didn't understand how to simulate asset prices by using non normal random numbers. I am assuming that it would be incorrect to use the standard Geometric Brownian Motion, since it is based solely ...
0
votes
2answers
135 views

What are some examples of non-solvable SDE where Monte Carlo discretization is necessary

Reading Glasserman - "Monte Carlo Methods in Finance" it says in the introduction to Chapter 6 - Discretization Methods, that moste models arising in derivatives pricing can be simulated only ...
0
votes
0answers
48 views

Forecast of ARMA-GARCH model in R

I managed to forecast a GARCH model yesterday and run a Monte Carlo simulation on R. Nevertheless, I can't do the same with an ARMA-GARCH. I tested 4 different method but without achieving an ...
2
votes
1answer
86 views

Evaluation of the semi-closed Heston pricing formula for call options

I'd like to know, how the integral part of the semi-closed Heston pricing formula for call options can be simulated for a given set of model parameters. Monte Carlo simulations shoud work for this ...
0
votes
0answers
81 views

Monte Carlo, convexity and Risk-Neutral ZCB Pricing

I've built a simplistic Excel monte carlo model to price a zero-coupon bond, but it came up with a slightly unepxected result so I wanted to confirm whether my maths is just a little rusty or my model ...
0
votes
1answer
35 views

Should earnings be modelled normally or lognormally?

I am having difficulty deciding whether a company's earnings should be modelled normally or lognormally. If we consider two arguments: (i) The earnings of a company are the returns on the assets of ...
2
votes
0answers
74 views

Simulation of Heston process

I am currently working on implementing Heston model in matlab for option pricing (in this case I am trying to price a European call) and I wanted to compare the results I obtain from using the exact ...
1
vote
1answer
75 views

How to use Halton sequence in monte carlo simulation

Does anybody know how to use the Halton pseudo random technique in monte carlo simulation. I'm able to generate the sequences and I know they are correct. I checked a couple of numbers from different ...
4
votes
1answer
217 views

How to do a Brownian Bridge with quasi-random numbers in the Heston model?

I'm required to use the Euler Monte Carlo method to compute the option price under Heston model settings. I know from some paper that the convergence is volatile for the Heston model with a plain ...
4
votes
2answers
592 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
1answer
41 views

How do you deal with Inflation lag in a MC simulation?

Consider the UK RPI index. This index is published every month around the 15th (give or take a few days). The publication refers to the RPI index of the month before, so there is a lag of a few weeks ...
4
votes
2answers
67 views

Importance Sampling - where to center the sampling distribution?

Consider a Monte Carlo (MC) approximation to a European call with BS parameters $r = 0.05, \sigma = 0.4, T = 10, S_0 = 50$ and $K = 95$. Consider the following results, each using 1M points: plain ...
3
votes
2answers
259 views

Practical implementation of Libor Market Model

I am trying to implement a project about the BGM model, suggested in the book "The Concepts and Practice of mathematical finance" by Mark Joshi. My question is related to the forward volatility ...
2
votes
1answer
113 views

Black-Scholes under stochastic interest rates

I'm trying to implement the Black-Scholes formula to price a call option under stochastic interest rates. Following the book of McLeish (2005), the formula is given by (assuming interest rates are ...
1
vote
0answers
34 views

Initial values for Heston Model calibration

I'm doing a Heston model in Matlab using simple Monte Carlo simulations (5.000 paths and 2 steps per day, simulating 360 days). When I try to calibrate the Heston parameters using fminsearch it takes ...
5
votes
3answers
226 views

Greeks: Why does my Monte Carlo give correct delta but incorrect gamma?

For a vanilla European call, my Monte Carlo method gives the right option price and delta but the wrong gamma. In particular, the value of gamma varies wildly each time I run the method. I estimate ...
0
votes
1answer
67 views

Compute cross-gamma

I am trying to use delta-gamma method with montecarlo simulations to calculate the VAR of a portfolio consisting in options and equities. To use the method I need to compute a gamma matrix, that has ...
7
votes
1answer
171 views

Speeding up computations: when to use Quasi and standard Monte-Carlo in pricing

I am familiar with the theory of Monte-Carlo techniques in the numerical integration, and recently I have started my experiments with these methods applied to derivatives pricing. I am using ...
1
vote
1answer
62 views

Best way to do multithread Monte-Carlo in QuantLib

QuantLib has great facilities for Monte-Carlo pricing engines, classes McSimulation and MonteCarloModel do a lot of work. But they do it in a single thread. What is best way to introduce parallel run ...
3
votes
2answers
2k views

Is there a step-by-step guide for calculating portfolio VaR using monte carlo simulations

I am trying to determine a step-by-step algorithm for calculating a portfolio's VaR using monte carlo simulations. It seems to me that the literature for this is extraordinarily opaque for something ...
3
votes
0answers
40 views

How to price lookback american option when its payment is distributed during its life

I would like to price a floating strike american lookback with a particular feature: I don't want to charge upfront the client, rather I would like to insert a "running fee", some sort of a dividend. ...
4
votes
1answer
85 views

(Re) normalisation of random variable in Monte-Carlo simulations

I have a very simple model (CIR) with a very simple discretisation scheme (Euler) and I use it to do Monte-Carlo Simulations. It is working. Someone insisted that renormalization of my random ...
3
votes
2answers
213 views

How to simulate a CIR process using GPU and Matlab

I am trying to simulate a CIR process using Matlab and my GPU for effeciency. At the moment i run into some implementation problems due to the recursive nature of the discretization. The sheme I ...
1
vote
1answer
103 views

Calculate CVaR for a portfolio

I would like to calculate the Conditional Value at Risk for a portfolio. To be honest, I'm trying for a few days to find an example to calculate for an entire portfolio, not just for one security and ...
1
vote
0answers
92 views

Monte Carlo simulation returns not normal distributed

I am generating 100,000 paths of SPX out to 1 year using Euler discretization. I look at how S is distributed for 100,000 paths at the 1 year point and I find it is lognormally distributed. I look at ...
1
vote
2answers
204 views

What does “convergence” in Monte Carlo simulation mean?

I have read about convergence in terms of MC simulation for derivative pricing, but I am not clear on what it exactly means. Let us suppose I price an option 100,000 paths twice and both result in the ...
8
votes
1answer
311 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) = ...
0
votes
2answers
222 views

How to price exotic options using Monte-Carlo?

I am actually trying to solve some exercise problem using Monte-Carlo and C++ for exotic options. Namely, the exotic options are geometric Asian options and discrete barrier option. It is claimed ...
4
votes
1answer
273 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 ...
0
votes
1answer
78 views

Getting Parameter of Translated Gamma Distribution from Monte Carlo

Spin-off from here. (Edit) Main question: What do I do about a parameter whose suggested values range quite vastly? (Edit) Backstory: I am given data of loss values and the dates that correspond to ...
0
votes
1answer
37 views

Compute moments of aggregate loss using Monte Carlo

Spin-off from here. Richard referred to me an article that tells me how to get parameters of a translated gamma distribution to which I should consider fitting simulated aggregated loss values. The ...
0
votes
1answer
161 views

Get distribution for aggregate loss using Monte Carlo

I am given two data sets containing dates and losses (in some currency). Given a distribution for the amount of losses and an (a,b,0) distribution for frequency of losses, how can I use Monte Carlo ...
3
votes
0answers
307 views

Portfolio Optimization with Monte Carlo Simulation - How to do it with Excel?

If I have three asset classes and their historical weekly returns for five years, how can I construct a minimum variance portfolio and an efficient frontier plot with Excel? To do that do I have to ...
1
vote
1answer
97 views

Median value for geometric brownian motion simulation

I'm trying to simulate stock prices using GBM. I am using the following formula, and MATLAB function, to determine the stock prices: $\nu = \mu - \frac{\sigma^{2}}{2}$; $S = S0*\text{[ones(1,nsims); ...
6
votes
2answers
250 views

How to estimate the greeks with a Monte Carlo simulation?

I am simulating the path of three indices to price a 1 year basket option. All the indices are domestic, so there is no currency component. At each time step I am using the local volatility ...
0
votes
0answers
48 views

How to generate jump times in in Multilevel path simulation for jump-diffusion SDEs?

I am trying to generate jump times in in Multilevel path simulation for jump-diffusion SDEs using the following MATLAB code: I used following Algorithm in Yuan Xia paper: But I have not reached ...
2
votes
0answers
169 views

How to price zero coupon bonds with the Monte Carlo method?

Im trying to calculate monthly ZCB bond prices with a fixed maturity T, over a period of months via Monte Carlo methods. Here is my attempt: For the first month, the price is $P_{t_0}(0,T) = ...
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 ...
3
votes
1answer
202 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 ...
1
vote
2answers
131 views

Importance sampling for barrier option like pricing by Monte carlo

I would like to know some references regarding importance sampling algorithms for variance reduction of Monte Carlo barrier options pricing. Please could someone help me leaving some references? If ...
1
vote
0answers
52 views

Price 3m libor autocap with LMM calibrated on 1y swaption data

I need to calculate a price of an autocap contract which is An autocap is similar to a cap, but at most γ ≤ β caplets can be exercised, and they have to be automatically exercised when in the ...
2
votes
1answer
114 views

Which quantitative tools are actually used for hedging energy price and volume risk?

I'm a finance professor and I am looking for someone with actual trading and risk management knowledge within the energy sector who can tell me about pricing and hedging energy (especially electricity ...
2
votes
0answers
192 views

Law of a geometric brownian motion first hitting time (formula dont match Monte Carlo Simulation)

I posted this question before on MSE I need to use it in a small step in the middle of a simulation and I think I'm not getting correct results to this probabilities and so for my all ...