10.06.2016 Views

eldo_user

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Monte Carlo Analysis<br />

Importance Sampling Monte Carlo<br />

The two forms of the MCPROB function can be combined to calculate yield, that is, the<br />

probability π L, U that a given response is between bounds y L and y U (with y L ≤ y U ). To calculate<br />

such a probability:<br />

use the following commands:<br />

.EXTRACT MC LABEL=PROB_L MCPROB(OUTPUT, LE, BOUND_L)<br />

.EXTRACT MC LABEL=PROB_U MCPROB(OUTPUT, GE, BOUND_U)<br />

.EXTRACT MC LABEL=PROB_LU '1-(PROB_L+PROB_U)'<br />

where OUTPUT is the name of the output response H(X), and where BOUND_L and<br />

BOUND_U are y L and y U respectively.<br />

Quantile Computation using MCBOUND<br />

The MCBOUND function computes the quantile y π associated with a given probability π.<br />

Formally, the computed value is the solution to the following optimization problem. For 0 < π <<br />

1, the quantity F Y −1 (π) = inf{y | F Y (y) ≥ π} is called the π-quantile of the random variable Y.<br />

The MCBOUND function is inverse to the MCPROB(…, LE, …) function. To calculate a<br />

quantile y π satisfying:<br />

use the following command:<br />

.EXTRACT MC LABEL=BOUND MCBOUND(OUTPUT, PROB)<br />

where OUTPUT is the name of the output response H(X), and where PROB is π.<br />

To determine a yield interval, defined as the interval (y L , y U ) such that the associated coverage<br />

probability π L, U has a given value, use the following commands:<br />

.EXTRACT MC LABEL=BOUND_L MCBOUND(OUTPUT, PROB_L)<br />

.EXTRACT MC LABEL=BOUND_U MCBOUND(OUTPUT, PROB_U)<br />

where OUTPUT is the name of the output response H(X), and where PROB_L and PROB_U<br />

are probability values satisfying PROB_U - PROB_L = π L, U .<br />

Eldo® User's Manual, 15.3 455

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!