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 openturns.viewer import View
ot.RandomGenerator.SetSeed(0)
dimension = 2
R = ot.CorrelationMatrix(dimension)
R[0, 1] = 0.8
distribution = ot.Normal([3.0] * dimension, [2.0] * dimension, R)
size = 100
sample = distribution.getSample(size)
firstSample = ot.Sample(size, 1)
secondSample = ot.Sample(size, 1)
for i in range(size):
firstSample[i] = ot.Point(1, sample[i, 0])
secondSample[i] = ot.Point(1, sample[i, 1])
lmtest = ot.LinearModelAlgorithm(firstSample, secondSample).getResult()
drawLinearModelResidual = ot.VisualTest.DrawLinearModelResidual(lmtest)
View(drawLinearModelResidual, figure_kw={"figsize": (5, 5)})
|