Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Monte Carlo simulation with inequality constraints
Replies: 11   Last Post: Mar 23, 2012 8:53 AM

 Search Thread: Advanced Search

 Messages: [ Previous | Next ]
 deltaquattro@gmail.com Posts: 77 Registered: 7/21/06
Monte Carlo simulation with inequality constraints
Posted: Mar 15, 2012 6:21 AM
 Plain Text Reply

Hi, all,

I've been given a computer codes with some input variables X1,...Xn and some outputs Y1,...Ym. I want to use a Monte Carlo code to compute the distributions of the Ys, given some distributions for the Xs. The Xs are geometrical parameters of an industrial design, and as a first approximation they are considered independent. Unfortunately the code doesn't run when X1<X2 or X3<X4, so I need to impose the two constraints X1>=X2 and X3>=X4.
How can I do it? I'd prefer a solution which is valid for general distributions, but I can also accept one which works when the Xs are normally distributed.

Some ideas which sprung to my mind:
1. My Monte Carlo code allows to define a linear correlation matrix for the input variables: I don't see how this helps, but maybe you do :)
2. I let the Monte Carlo code to generate freely the sample runs, and whenever a run has X1<X2 or X3<X4, I discard it. However, I'm worried that this "rejection process" may distort the distributions of the Xs. Also, I guess I'll need to perform four times as many Monte Carlo runs as usual, to have the same level of statistical convergence.
3. I change my Monte Carlo variables: instead than X1 and X3, I use Z1 = |X1-X2| and Z2=|X3-X4|. This way, X1 = X2 + Z1 >= X2, and X3 = X4 + Z2 >= X4. However, while I was able to make some reasonable assumptions on the distributions of the original variables, now I have no idea which distributions I should use for Z1 and Z2...Thanks,

Best Regards,

deltaquattro

ps apologies for the double posts, but I'm not sure which ng is more suitable for this thread, and I'm not able to create a crosspost with Google Groups.

Date Subject Author
3/15/12 deltaquattro@gmail.com
3/15/12 Dave Dodson
3/15/12 deltaquattro@gmail.com
3/15/12 Ray Koopman
3/16/12 deltaquattro@gmail.com
3/16/12 Ray Koopman
3/19/12 deltaquattro@gmail.com
3/22/12 Ray Koopman
3/22/12 deltaquattro@gmail.com
3/22/12 Ray Koopman
3/23/12 deltaquattro@gmail.com
3/16/12 AMX

© The Math Forum at NCTM 1994-2018. All Rights Reserved.