1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
import openturns as ot
from matplotlib import pyplot as plt
import openturns.viewer as otv
pdf_graph = ot.Graph("PDF graph", "x", "PDF", True, "upper left")
cdf_graph = ot.Graph("CDF graph", "x", "CDF", True, "upper left")
palette = ot.Drawable.BuildDefaultPalette(10)
for i, p in enumerate([-1.0, 0.0, 1.0]):
distribution = ot.GeneralizedExtremeValue(0.0, 1.0, p)
pdf_curve = distribution.drawPDF().getDrawable(0)
cdf_curve = distribution.drawCDF().getDrawable(0)
legend = f"xi={p}"
pdf_curve.setLegend(legend)
cdf_curve.setLegend(legend)
pdf_graph.add(pdf_curve)
cdf_graph.add(cdf_curve)
fig = plt.figure(figsize=(10, 4))
pdf_axis = fig.add_subplot(121)
cdf_axis = fig.add_subplot(122)
otv.View(pdf_graph, figure=fig, axes=[pdf_axis], add_legend=True)
otv.View(cdf_graph, figure=fig, axes=[cdf_axis], add_legend=True)
fig.suptitle("GeneralizedExtremeValue(0,1,xi)")
|
