File: Staircase.py

package info (click to toggle)
openturns 1.7-3
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 38,588 kB
  • ctags: 26,495
  • sloc: cpp: 144,032; python: 26,855; ansic: 7,868; sh: 419; makefile: 263; yacc: 123; lex: 44
file content (33 lines) | stat: -rw-r--r-- 898 bytes parent folder | download
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)