#! /usr/bin/env python
from __future__ import print_function
import openturns as ot
ot.TESTPREAMBLE()
ot.RandomGenerator.SetSeed(0)
# Instanciate one distribution object
distribution = ot.MaximumEntropyOrderStatisticsDistribution([
ot.Trapezoidal(-2.0, -1.1, -1.0, 1.0),
ot.LogUniform(1.0, 1.2),
ot.Triangular(3.0, 4.5, 5.0),
ot.Beta(2.5, 3.5, 4.7, 5.2)
])
dim = distribution.getDimension()
print("Distribution ", distribution)
# Is this distribution elliptical ?
print("Elliptical = ", distribution.isElliptical())
# Test for realization of distribution
oneRealization = distribution.getRealization()
print("oneRealization=", repr(oneRealization))
# Test for sampling
size = 10000
oneSample = distribution.getSample(size)
print("oneSample first=", repr(oneSample[0]), " last=",
repr(oneSample[size - 1]))
# In this example we are going to build maximum entropy statistics distribution, which yields ordered realizations: # # .. math:: # X_1 \leq \dots \leq X_n # # %% import openturns as ot import openturns.viewer as viewer from matplotlib import pylab as plt ot.Log.Show(ot.Log.NONE) # %% # create a collection of distribution coll = [ot.Beta(1.5, 1.7, 0.0, 1.0), ot.Beta(2.0, 2.3, 0.5, 1.2)] # %% # create the distribution distribution = ot.MaximumEntropyOrderStatisticsDistribution(coll) print(distribution) # %% # draw a sample (ordered!) distribution.getSample(10) # %% # draw PDF graph = distribution.drawPDF() view = viewer.View(graph) plt.show()
import openturns as ot
from matplotlib import pyplot as plt
from openturns.viewer import View
myDist = [ot.Beta(1.5, 3.2, 0.0, 1.0), ot.Beta(2.0, 4.3, 0.5, 1.2)]
myOrderStatDist = ot.MaximumEntropyOrderStatisticsDistribution(myDist)
myOrderStatDist.setDescription(['$x_1$', '$x_2$'])
graphPDF = myOrderStatDist.drawPDF()
graphCDF = myOrderStatDist.drawCDF()
fig = plt.figure(figsize=(8, 4))
plt.suptitle("Max Entropy Order Statistics Distribution: pdf and cdf")
pdf_axis = fig.add_subplot(121)
cdf_axis = fig.add_subplot(122)
pdf_axis.set_xlim(auto=True)
cdf_axis.set_xlim(auto=True)
View(graphPDF, figure=fig, axes=[pdf_axis], add_legend=True)
View(graphCDF, figure=fig, axes=[cdf_axis], add_legend=True)
