← Back to branch summary

~ubuntu-branches/ubuntu/maverick/openturns/maverick

« back to all changes in this revision

Viewing changes to python/test/t_ClaytonCopula_std.py

  • Committer: Bazaar Package Importer
  • Author(s): Fabrice Coutadeur
  • Date: 2010-05-10 17:27:55 UTC
  • mfrom: (1.1.4 upstream) (5.1.5 sid)
  • Revision ID: james.westby@ubuntu.com-20100510172755-cb5ynskknqqi5rhp
Tags: 0.13.2-2ubuntu1
* Merge with Debian testing. No changes left.
* ubuntu_fix-python-2.6.patch: fix detection of python 2.6 libs, to not use
  LOCALMODLIBS. This pulls a dependency on SSL and makes the package FTBFS.

Show diffs side-by-side

added added

removed removed

Lines of Context:
python/test/t_ClaytonCopula_std.py
10
10
    # Instanciate one distribution object
11
11
    dim = 2
12
12
    copula = ClaytonCopula(2.5)
13
 
    print "Copula " , copula 
 
13
    print "Copula " , copula
14
14
    print "Mean " , repr(copula.getMean() )
15
 
    print "Covariance " , copula.getCovariance() 
16
 
 
 
15
    print "Covariance " , repr(copula.getCovariance())
 
16
    
17
17
    # Is this copula an elliptical distribution?
18
 
    print "Elliptical distribution= " , copula.isElliptical() 
19
 
 
 
18
    print "Elliptical distribution= " , copula.isElliptical()
 
19
    
20
20
    # Is this copula elliptical ?
21
21
    print "Elliptical copula= " , copula.hasEllipticalCopula()
22
 
 
 
22
    
23
23
    # Is this copula independent ?
24
 
    print "Independent copula= " , copula.hasIndependentCopula() 
25
 
 
 
24
    print "Independent copula= " , copula.hasIndependentCopula()
 
25
    
26
26
    # Test for realization of distribution
27
27
    oneRealization = copula.getRealization()
28
28
    print "oneRealization=", repr(oneRealization)
29
 
 
 
29
    
30
30
    # Test for sampling
31
31
    size = 10
32
32
    oneSample = copula.getNumericalSample( size )
33
 
    print "oneSample=", oneSample
34
 
 
35
 
   # Test for sampling
 
33
    print "oneSample=", repr(oneSample)
 
34
    
 
35
    # Test for sampling
36
36
    size = 10000
37
37
    anotherSample = copula.getNumericalSample( size )
38
38
    print "anotherSample mean=", repr(anotherSample.computeMean())
39
 
    print "anotherSample covariance=", anotherSample.computeCovariance()
 
39
    print "anotherSample covariance=", repr(anotherSample.computeCovariance())
40
40
    
41
41
    # Define a point
42
42
    point = NumericalPoint(dim, 0.2)
43
 
 
44
 
    # Show PDF and CDF of point 
 
43
    
 
44
    # Show PDF and CDF of point
45
45
    pointPDF = copula.computePDF( point )
46
46
    pointCDF = copula.computeCDF( point )
47
47
    print "Point = ", repr(point), " pdf=%.6f" % pointPDF, " cdf=%.6f" % pointCDF
48
 
 
 
48
    
49
49
    # Get 50% quantile
50
50
    quantile = copula.computeQuantile( 0.5 )
51
51
    print "Quantile=", repr(quantile)
52
52
    print "CDF(quantile)=%.6f" % copula.computeCDF(quantile)
53
 
 
 
53
    
54
54
    # Extract the marginals
55
55
    for i in range(dim) :
56
 
      margin = copula.getMarginal(i)
57
 
      print "margin=", margin 
58
 
      print "margin PDF=%.6f" % margin.computePDF( NumericalPoint(1,0.25))
59
 
      print "margin CDF=%.6f" % margin.computeCDF(NumericalPoint(1,0.25))
60
 
      print "margin quantile=", repr(margin.computeQuantile(0.95)) 
61
 
      print "margin realization=", repr(margin.getRealization()) 
62
 
 
 
56
        margin = copula.getMarginal(i)
 
57
        print "margin=", margin
 
58
        print "margin PDF=%.6f" % margin.computePDF( NumericalPoint(1,0.25))
 
59
        print "margin CDF=%.6f" % margin.computeCDF(NumericalPoint(1,0.25))
 
60
        print "margin quantile=", repr(margin.computeQuantile(0.95))
 
61
        print "margin realization=", repr(margin.getRealization())
 
62
        
63
63
    # Extract a 2-D marginal
64
64
    indices = Indices(2, 0)
65
65
    indices[0] = 1
72
72
    quantile = NumericalPoint(margins.computeQuantile(0.95))
73
73
    print "margins quantile=", repr(quantile)
74
74
    print "margins CDF(qantile)=%.6f" % margins.computeCDF(quantile)
75
 
    print "margins realization=", repr(margins.getRealization()) 
76
 
 
 
75
    print "margins realization=", repr(margins.getRealization())
 
76
    
77
77
except :
78
 
   import sys
79
 
   print "t_NormalCopula_std.py", sys.exc_type, sys.exc_value
 
78
    import sys
 
79
    print "t_NormalCopula_std.py", sys.exc_type, sys.exc_value