Yomi | 1 Jul 10:38

Monte Carlo simulation convergence issue


Hi,

I have created different engiens based on the MC simulation framework.
I am experiencing some convergence issues. If I increase the number of steps
the Max Time Step by year, the price starts to diverge. An example would be
the following graph:

http://www.nabble.com/file/p18211801/MC%2Bissues.jpg 

I have the following questions:

1. have you ever encountered such issues?
2. How to use the error estimate? Whatever the RNG I put (Low
Discrenpancies, Pseudo Random...) I can nto access the error which might be
usefull.

Thanks
--

-- 
View this message in context: http://www.nabble.com/Monte-Carlo-simulation-convergence-issue-tp18211801p18211801.html
Sent from the quantlib-users mailing list archive at Nabble.com.

-------------------------------------------------------------------------
Check out the new SourceForge.net Marketplace.
It's the best place to buy or sell services for
just about anything Open Source.
http://sourceforge.net/services/buy/index.php
Luigi Ballabio | 22 Aug 17:38

Re: Monte Carlo simulation convergence issue

Guillaume,
	first of all, apologies for the delay.

On Tue, 2008-07-01 at 01:40 -0700, Yomi wrote:
> I have created different engiens based on the MC simulation framework.
> I am experiencing some convergence issues. If I increase the number of steps
> the Max Time Step by year, the price starts to diverge. An example would be
> the following graph:
> 
> http://www.nabble.com/file/p18211801/MC%2Bissues.jpg 
> 
> I have the following questions:
> 
> 1. have you ever encountered such issues?

Not that I remember. Did you try doing the same with one of the engines
in the library? Do they have the same behavior?

> 2. How to use the error estimate? Whatever the RNG I put (Low
> Discrenpancies, Pseudo Random...) I can nto access the error which might be
> usefull.

Your engine should access it from the model and write it into the
results. In your engine's calculate() method, you probably have
something like:

this->results_.value = this->mcModel_->sampleAccumulator().mean();

after that, you'll have to add:

(Continue reading)


Gmane