1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33
|
import openturns as ot
from matplotlib import pyplot as plt
from openturns.viewer import View
distribution = ot.Normal([0.5], [2.], ot.CorrelationMatrix(1))
sample = distribution.getSample(200)
# Construct empirical CDF for the sample
data = ot.NumericalSample(20, 2)
cursor = ot.NumericalPoint(2)
for i in range(20):
cursor[0] = 13. * i / 20 - 6.5
count = 0.
for j in range(200):
if(sample[j, 0] < cursor[0]):
count += 1.
cursor[1] = count / 200
data[i] = cursor
# Create an empty graph
graph = ot.Graph("Some curves", "x1", "x2", True, "topleft")
# Create the staircase
staircase = ot.Staircase(data, "blue", "solid", "s", "Normal CDF")
# Then, draw it
graph.add(staircase)
fig = plt.figure(figsize=(4.5, 4))
plt.suptitle("Staircase example")
axis = fig.add_subplot(111)
axis.set_xlim(auto=True)
View(graph, figure=fig, axes=[axis], add_legend=True)
|