ThresholdEvent¶

class ThresholdEvent(*args)¶

Random vector defined from a comparison operator and a threshold.

The event occurs when the realization of the underlying random vector exceeds the threshold.

Parameters

antecedentRandomVector of dimension 1: Output variable of interest.
comparisonOperatorComparisonOperator: Comparison operator used to compare antecedent with threshold.
thresholdfloat: threshold we want to compare to antecedent.

See also

ProcessEvent, DomainEvent

Notes

An event is defined as follows:

$\cD_f = \{\vect{X} \in \Rset^n \, / \, g(\vect{X},\vect{d}) \le 0\}$

where $\vect{X}$ denotes a random input vector, representing the sources of uncertainties, $\vect{d}$ is a determinist vector, representing the fixed variables and $g(\vect{X},\vect{d})$ is the limit state function of the model. The probability content of the event $\cD_f$ is $P_f$ :

$P_f = \int_{g(\vect{X},\vect{d})\le 0}f_\vect{X}(\vect{x})\di{\vect{x}}$

Here, the event considered is explicited directly from the limit state function $g(\vect{X}\,,\,\vect{d})$ : this is the classical structural reliability formulation. However, if the event is a threshold exceedance, it is useful to explicite the variable of interest $Z=\tilde{g}(\vect{X}\,,\,\vect{d})$ , evaluated from the model $\tilde{g}(.)$ . In that case, the event considered, associated to the threshold $z_s$ has the formulation:

$\cD_f = \{ \vect{X} \in \Rset^n \, / \, Z=\tilde{g}(\vect{X}\,,\,\vect{d}) > z_s \}$

and the limit state function is:

$g(\vect{X}\,,\,\vect{d}) &= z_s - Z \\ &= z_s - \tilde{g}(\vect{X}\,,\,\vect{d})$

$P_f$ is the threshold exceedance probability, defined as:

$P_f &= P(Z \geq z_s) \\ &= \int_{g(\vect{X}, \vect{d}) \le 0} \pdf\di{\vect{x}}$

Examples

>>> import openturns as ot
>>> dim = 2
>>> X = ot.RandomVector(ot.Normal(dim))
>>> model = ot.SymbolicFunction(['x1', 'x2'], ['x1+x2'])
>>> Y = ot.CompositeRandomVector(model, X)
>>> event = ot.ThresholdEvent(Y, ot.Less(), 1.0)

Methods

`getAntecedent`(self)	Accessor to the antecedent RandomVector in case of a composite RandomVector.
`getClassName`(self)	Accessor to the object’s name.
`getCovariance`(self)	Accessor to the covariance of the RandomVector.
`getDescription`(self)	Accessor to the description of the RandomVector.
`getDimension`(self)	Accessor to the dimension of the RandomVector.
`getDistribution`(self)	Accessor to the distribution of the RandomVector.
`getDomain`(self)	Accessor to the domain of the Event.
`getFunction`(self)	Accessor to the Function in case of a composite RandomVector.
`getId`(self)	Accessor to the object’s id.
`getImplementation`(self)	Accessor to the underlying implementation.
`getMarginal`(self, \*args)	Get the random vector corresponding to the $i^{th}$ marginal component(s).
`getMean`(self)	Accessor to the mean of the RandomVector.
`getName`(self)	Accessor to the object’s name.
`getOperator`(self)	Accessor to the comparaison operator of the Event.
`getParameter`(self)	Accessor to the parameter of the distribution.
`getParameterDescription`(self)	Accessor to the parameter description of the distribution.
`getRealization`(self)	Compute one realization of the RandomVector.
`getSample`(self, size)	Compute realizations of the RandomVector.
`getThreshold`(self)	Accessor to the threshold of the Event.
`intersect`(self, other)	Intersection of two events.
`isComposite`(self)	Accessor to know if the RandomVector is a composite one.
`isEvent`(self)	Whether the random vector is an event.
`join`(self, other)	Union of two events.
`setDescription`(self, description)	Accessor to the description of the RandomVector.
`setName`(self, name)	Accessor to the object’s name.
`setParameter`(self, parameters)	Accessor to the parameter of the distribution.

__init__(self, \*args)¶: Initialize self. See help(type(self)) for accurate signature.

getAntecedent(self)¶

Accessor to the antecedent RandomVector in case of a composite RandomVector.

Returns

antecedentRandomVector: Antecedent RandomVector $\vect{X}$ in case of a CompositeRandomVector such as: $\vect{Y}=f(\vect{X})$ .

getClassName(self)¶

Accessor to the object’s name.

Returns

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

getCovariance(self)¶

Accessor to the covariance of the RandomVector.

Returns

covarianceCovarianceMatrix: Covariance of the considered UsualRandomVector.

Examples

>>> import openturns as ot
>>> distribution = ot.Normal([0.0, 0.5], [1.0, 1.5], ot.CorrelationMatrix(2))
>>> randomVector = ot.RandomVector(distribution)
>>> ot.RandomGenerator.SetSeed(0)
>>> print(randomVector.getCovariance())
[[ 1    0    ]
 [ 0    2.25 ]]

getDescription(self)¶

Accessor to the description of the RandomVector.

Returns

descriptionDescription: Describes the components of the RandomVector.

getDimension(self)¶

Accessor to the dimension of the RandomVector.

Returns

dimensionpositive int: Dimension of the RandomVector.

getDistribution(self)¶

Accessor to the distribution of the RandomVector.

Returns

distributionDistribution: Distribution of the considered UsualRandomVector.

Examples

>>> import openturns as ot
>>> distribution = ot.Normal([0.0, 0.0], [1.0, 1.0], ot.CorrelationMatrix(2))
>>> randomVector = ot.RandomVector(distribution)
>>> ot.RandomGenerator.SetSeed(0)
>>> print(randomVector.getDistribution())
Normal(mu = [0,0], sigma = [1,1], R = [[ 1 0 ]
 [ 0 1 ]])

getDomain(self)¶

Accessor to the domain of the Event.

Returns

domainDomain: Describes the domain of an event.

getFunction(self)¶

Accessor to the Function in case of a composite RandomVector.

Returns

functionFunction: Function used to define a CompositeRandomVector as the image through this function of the antecedent $\vect{X}$ : $\vect{Y}=f(\vect{X})$ .

getId(self)¶

Accessor to the object’s id.

Returns

idint: Internal unique identifier.

getImplementation(self)¶

Accessor to the underlying implementation.

Returns

implImplementation: The implementation class.

getMarginal(self, \*args)¶

Get the random vector corresponding to the $i^{th}$ marginal component(s).

Parameters

iint or list of ints, $0\leq i < dim$: Indicates the component(s) concerned. $dim$ is the dimension of the RandomVector.

Returns

vectorRandomVector: RandomVector restricted to the concerned components.

Notes

Let’s note $\vect{Y}=\Tr{(Y_1,\dots,Y_n)}$ a random vector and $I \in [1,n]$ a set of indices. If $\vect{Y}$ is a UsualRandomVector, the subvector is defined by $\tilde{\vect{Y}}=\Tr{(Y_i)}_{i \in I}$ . If $\vect{Y}$ is a CompositeRandomVector, defined by $\vect{Y}=f(\vect{X})$ with $f=(f_1,\dots,f_n)$ , $f_i$ some scalar functions, the subvector is $\tilde{\vect{Y}}=(f_i(\vect{X}))_{i \in I}$ .

Examples

>>> import openturns as ot
>>> distribution = ot.Normal([0.0, 0.0], [1.0, 1.0], ot.CorrelationMatrix(2))
>>> randomVector = ot.RandomVector(distribution)
>>> ot.RandomGenerator.SetSeed(0)
>>> print(randomVector.getMarginal(1).getRealization())
[0.608202]
>>> print(randomVector.getMarginal(1).getDistribution())
Normal(mu = 0, sigma = 1)

getMean(self)¶

Accessor to the mean of the RandomVector.

Returns

meanPoint: Mean of the considered UsualRandomVector.

Examples

>>> import openturns as ot
>>> distribution = ot.Normal([0.0, 0.5], [1.0, 1.5], ot.CorrelationMatrix(2))
>>> randomVector = ot.RandomVector(distribution)
>>> ot.RandomGenerator.SetSeed(0)
>>> print(randomVector.getMean())
[0,0.5]

getName(self)¶

Accessor to the object’s name.

Returns

namestr: The name of the object.

getOperator(self)¶

Accessor to the comparaison operator of the Event.

Returns

operatorComparisonOperator: Comparaison operator used to define the RandomVector.

getParameter(self)¶

Accessor to the parameter of the distribution.

Returns

parameterPoint: Parameter values.

getParameterDescription(self)¶

Accessor to the parameter description of the distribution.

Returns

descriptionDescription: Parameter names.

getRealization(self)¶

Compute one realization of the RandomVector.

Returns

aRealizationPoint: Sequence of values randomly determined from the RandomVector definition. In the case of an event: one realization of the event (considered as a Bernoulli variable) which is a boolean value (1 for the realization of the event and 0 else).

OpenTURNS

An Open source initiative for the Treatment of Uncertainties, Risks'N Statistics

Previous topic

Next topic

This Page

ThresholdEvent¶