ProbabilitySimulationResult

class ProbabilitySimulationResult(*args)

Probability simulation result.

Methods

drawImportanceFactors()

Draw the importance factors.

getBlockSize()

Accessor to the block size.

getClassName()

Accessor to the object's name.

getCoefficientOfVariation()

Accessor to the coefficient of variation.

getConfidenceLength(*args)

Accessor to the confidence length.

getEvent()

Accessor to the event.

getImportanceFactors()

Accessor to the importance factors.

getMeanPointInEventDomain()

Accessor to the mean point conditioned to the event realization.

getName()

Accessor to the object's name.

getOuterSampling()

Accessor to the outer sampling.

getProbabilityDistribution()

Accessor to the asymptotic probability distribution.

getProbabilityEstimate()

Accessor to the probability estimate.

getStandardDeviation()

Accessor to the standard deviation.

getTimeDuration()

Accessor to the elapsed time.

getVarianceEstimate()

Accessor to the variance estimate.

hasName()

Test if the object is named.

setBlockSize(blockSize)

Accessor to the block size.

setEvent(event)

Accessor to the event.

setName(name)

Accessor to the object's name.

setOuterSampling(outerSampling)

Accessor to the outer sampling.

setProbabilityEstimate(probabilityEstimate)

Accessor to the probability estimate.

setTimeDuration(time)

Accessor to the elapsed time.

setVarianceEstimate(varianceEstimate)

Accessor to the variance estimate.

Notes

Structure created by the method run() of a EventSimulation, and obtained thanks to the method getResult().

Examples

>>> import openturns as ot
>>> ot.RandomGenerator.SetSeed(0)
>>> limitState = ot.SymbolicFunction(['E', 'F', 'L', 'I'], ['-F*L^3/(3.*E*I)'])
>>> # Enable the history mechanism in order to use the getImportanceFactors method
>>> limitState = ot.MemoizeFunction(limitState)
>>> distribution = ot.Normal([50.0, 1.0, 10.0, 5.0], [1.0]*4, ot.IdentityMatrix(4))
>>> output = ot.CompositeRandomVector(limitState, ot.RandomVector(distribution))
>>> event = ot.ThresholdEvent(output, ot.Less(), -3.0)
>>> experiment = ot.MonteCarloExperiment()
>>> algo = ot.ProbabilitySimulationAlgorithm(event, experiment)
>>> algo.run()
>>> result = algo.getResult()
>>> importanceFactors = result.getImportanceFactors()
__init__(*args)
drawImportanceFactors()

Draw the importance factors.

Warning

It is necessary to enable the history of the model to perform this analysis (see MemoizeFunction).

Returns:
graphGraph

Importance factor graph.

getBlockSize()

Accessor to the block size.

Returns:
blockSizeint

Number of terms in the probability simulation estimator grouped together.

getClassName()

Accessor to the object’s name.

Returns:
class_namestr

The object class name (object.__class__.__name__).

getCoefficientOfVariation()

Accessor to the coefficient of variation.

Returns:
coefficientfloat

Coefficient of variation of the simulated sample which is equal to \sqrt{Var_e} / P_e with Var_e the variance estimate and P_e the probability estimate.

getConfidenceLength(*args)

Accessor to the confidence length.

Parameters:
levelfloat, level \in ]0, 1[

Confidence level. By default, it is 0.95.

Returns:
confidenceLengthfloat

Length of the confidence interval at the confidence level level.

getEvent()

Accessor to the event.

Returns:
eventRandomVector

Event we want to evaluate the probability.

getImportanceFactors()

Accessor to the importance factors.

Returns:
importanceFactorsPointWithDescription

Sequence containing the importance factors with a description for each component.

Notes

The importance factors \alpha_i are evaluated from the coordinates of the mean point of event domain \vect{X}^*_{event}, mapped into the standard space as follows:

\alpha_i = \displaystyle \frac{\left(U_{i}^*\right)^2}{||\vect{U}^*||^2}

where \vect{U}^* = T(\vect{X}^*_{event}) with T the iso-probabilistic transformation and the mean point \vect{X}^*_{event} = \displaystyle \frac{1}{n} \sum_{i=1}^{n} \vect{X}_i 1_{event}(\vect{X}_i).

Warning

This notion is only available if the history mechanism of the model is activated (see MemoizeFunction).

getMeanPointInEventDomain()

Accessor to the mean point conditioned to the event realization.

Returns:
meanPointPoint

Mean point in the physical space of all the simulations generated by the EventSimulation algorithm that failed into the event domain.

Notes

Warning

This notion is only available if the history mechanism of the model is activated (see MemoizeFunction).

getName()

Accessor to the object’s name.

Returns:
namestr

The name of the object.

getOuterSampling()

Accessor to the outer sampling.

Returns:
outerSamplingint

Number of groups of terms in the probability simulation estimator.

getProbabilityDistribution()

Accessor to the asymptotic probability distribution.

Returns:
probaDistributionNormal

Asymptotic normal distribution of the event probability estimate.

getProbabilityEstimate()

Accessor to the probability estimate.

Returns:
probaEstimatefloat

Estimate of the event probability.

getStandardDeviation()

Accessor to the standard deviation.

Returns:
sigmafloat

Standard deviation of the estimator at the end of the simulation.

getTimeDuration()

Accessor to the elapsed time.

Returns:
timefloat

Simulation duration in seconds

getVarianceEstimate()

Accessor to the variance estimate.

Returns:
varianceEstimatefloat

Variance estimate.

hasName()

Test if the object is named.

Returns:
hasNamebool

True if the name is not empty.

setBlockSize(blockSize)

Accessor to the block size.

Parameters:
blockSizeint, blockSize \geq 0

Number of terms in the probability simulation estimator grouped together.

setEvent(event)

Accessor to the event.

Parameters:
eventRandomVector

Event we want to evaluate the probability.

setName(name)

Accessor to the object’s name.

Parameters:
namestr

The name of the object.

setOuterSampling(outerSampling)

Accessor to the outer sampling.

Parameters:
outerSamplingint, outerSampling \geq 0

Number of groups of terms in the probability simulation estimator.

setProbabilityEstimate(probabilityEstimate)

Accessor to the probability estimate.

Parameters:
probaEstimatefloat, 0 \leq P_e \leq 1

Estimate of the event probability.

setTimeDuration(time)

Accessor to the elapsed time.

Parameters:
timefloat

Simulation duration in seconds

setVarianceEstimate(varianceEstimate)

Accessor to the variance estimate.

Parameters:
varianceEstimatefloat, Var_e \geq 0

Variance estimate.

Examples using the class

Create a process from random vectors and processes

Create a process from random vectors and processes

Estimate a probability with Monte Carlo

Estimate a probability with Monte Carlo

Use a randomized QMC algorithm

Use a randomized QMC algorithm

Use the Adaptive Directional Stratification Algorithm

Use the Adaptive Directional Stratification Algorithm

Use the Directional Sampling Algorithm

Use the Directional Sampling Algorithm

Specify a simulation algorithm

Specify a simulation algorithm

Use the Importance Sampling algorithm

Use the Importance Sampling algorithm

Estimate a probability with Monte-Carlo on axial stressed beam: a quick start guide to reliability

Estimate a probability with Monte-Carlo on axial stressed beam: a quick start guide to reliability

Estimate a buckling probability

Estimate a buckling probability

Exploitation of simulation algorithm results

Exploitation of simulation algorithm results

Subset Sampling

Subset Sampling

Axial stressed beam : comparing different methods to estimate a probability

Axial stressed beam : comparing different methods to estimate a probability

Using the FORM - SORM algorithms on a nonlinear function

Using the FORM - SORM algorithms on a nonlinear function

Control algorithm termination

Control algorithm termination