|
14 | 14 | # |
15 | 15 | ######################################################################### |
16 | 16 | """ |
| 17 | + |
17 | 18 | __revision__ = "$Id: test_exploratory.py 8676 2010-04-20 15:28:46Z cokelaer $" |
18 | 19 |
|
| 20 | +from os.path import join as pj |
| 21 | + |
19 | 22 | from openalea.sequence_analysis import * |
20 | | -from openalea.sequence_analysis.estimate import Estimate |
21 | | -from openalea.sequence_analysis.data import path |
22 | 23 | from openalea.sequence_analysis import get_shared_data as path |
23 | | -from os.path import join as pj |
24 | | -seq0 = Sequences(pj(path ,"chene_sessile_15pa.seq")) |
| 24 | +from openalea.sequence_analysis.data import path |
| 25 | +from openalea.sequence_analysis.estimate import Estimate |
25 | 26 |
|
26 | | -#Plot(seq0, ViewPoint="Data") |
| 27 | +seq0 = Sequences(pj(path("chene_sessile_15pa.seq"))) |
| 28 | + |
| 29 | +# Plot(seq0, ViewPoint="Data") |
27 | 30 |
|
28 | 31 |
|
29 | 32 | # change of unit for the variable diameter of the annual shoot |
30 | 33 |
|
31 | 34 | marginal3 = ExtractHistogram(seq0, "Value", 3) |
32 | | -#Plot(Cluster(marginal3, "Information", 0.75)) |
33 | | -#Plot(Cluster(marginal3, "Information", 0.61)) |
34 | | -#Plot(Cluster(marginal3, "Step", 10)) |
| 35 | +# Plot(Cluster(marginal3, "Information", 0.75)) |
| 36 | +# Plot(Cluster(marginal3, "Information", 0.61)) |
| 37 | +# Plot(Cluster(marginal3, "Step", 10)) |
35 | 38 |
|
36 | 39 | vec10 = Vectors(seq0) |
37 | 40 |
|
38 | 41 | # plot of the average sequence |
39 | | -#Plot(Regression(vec10, "MovingAverage", 1, 2, [1])) |
| 42 | +# Plot(Regression(vec10, "MovingAverage", 1, 2, [1])) |
40 | 43 |
|
41 | 44 | vec95 = ValueSelect(vec10, 1, 95) |
42 | 45 | vec96 = ValueSelect(vec10, 1, 96) |
|
45 | 48 | VarianceAnalysis(vec10, 1, 2, "N") |
46 | 49 |
|
47 | 50 |
|
48 | | -print type(ExtractHistogram(vec95, 2)) |
| 51 | +print(type(ExtractHistogram(vec95, 2))) |
49 | 52 |
|
50 | | -Compare(ExtractHistogram(vec95, 2), ExtractHistogram(vec96, 2), ExtractHistogram(vec97, 2), "N") |
51 | | -#Plot(ExtractHistogram(vec95, 2), ExtractHistogram(vec96, 2), ExtractHistogram(vec97, 2)) |
| 53 | +Compare( |
| 54 | + ExtractHistogram(vec95, 2), |
| 55 | + ExtractHistogram(vec96, 2), |
| 56 | + ExtractHistogram(vec97, 2), |
| 57 | + "N", |
| 58 | +) |
| 59 | +# Plot(ExtractHistogram(vec95, 2), ExtractHistogram(vec96, 2), ExtractHistogram(vec97, 2)) |
52 | 60 |
|
53 | 61 | ContingencyTable(vec10, 1, 4) |
54 | 62 |
|
55 | 63 | # one-way variance analysis based on ranks |
56 | 64 |
|
57 | 65 | VarianceAnalysis(vec10, 1, 4, "O") |
58 | | -Compare(ExtractHistogram(vec95, 4), ExtractHistogram(vec96, 4), ExtractHistogram(vec97, 4), "O") |
59 | | -#Plot(ExtractHistogram(vec95, 4), ExtractHistogram(vec96, 4), ExtractHistogram(vec97, 4)) |
| 66 | +Compare( |
| 67 | + ExtractHistogram(vec95, 4), |
| 68 | + ExtractHistogram(vec96, 4), |
| 69 | + ExtractHistogram(vec97, 4), |
| 70 | + "O", |
| 71 | +) |
| 72 | +# Plot(ExtractHistogram(vec95, 4), ExtractHistogram(vec96, 4), ExtractHistogram(vec97, 4)) |
60 | 73 |
|
61 | | -#Plot(ExtractHistogram(vec95, 5), ExtractHistogram(vec96, 5), ExtractHistogram(vec97, 5)) |
62 | | -#Plot(ExtractHistogram(vec95, 6), ExtractHistogram(vec96, 6), ExtractHistogram(vec97, 6)) |
| 74 | +# Plot(ExtractHistogram(vec95, 5), ExtractHistogram(vec96, 5), ExtractHistogram(vec97, 5)) |
| 75 | +# Plot(ExtractHistogram(vec95, 6), ExtractHistogram(vec96, 6), ExtractHistogram(vec97, 6)) |
63 | 76 |
|
64 | 77 | vec11 = ValueSelect(vec10, 4, 1) |
65 | 78 | vec12 = ValueSelect(vec10, 4, 2) |
66 | 79 | vec13 = ValueSelect(vec10, 4, 3, 4) |
67 | 80 |
|
68 | | -#Plot(ExtractHistogram(vec11, 2), ExtractHistogram(vec12, 2), ExtractHistogram(vec13, 2)) |
69 | | -#Plot(ExtractHistogram(vec11, 5), ExtractHistogram(vec12, 5), ExtractHistogram(vec13, 5)) |
70 | | - |
71 | | -mixt20 = Estimate(ExtractHistogram(vec10, 2), "MIXTURE", "NB", "NB", "NB", "NB", NbComponent="Estimated") |
| 81 | +# Plot(ExtractHistogram(vec11, 2), ExtractHistogram(vec12, 2), ExtractHistogram(vec13, 2)) |
| 82 | +# Plot(ExtractHistogram(vec11, 5), ExtractHistogram(vec12, 5), ExtractHistogram(vec13, 5)) |
| 83 | + |
| 84 | +mixt20 = Estimate( |
| 85 | + ExtractHistogram(vec10, 2), |
| 86 | + "MIXTURE", |
| 87 | + "NB", |
| 88 | + "NB", |
| 89 | + "NB", |
| 90 | + "NB", |
| 91 | + NbComponent="Estimated", |
| 92 | +) |
72 | 93 | Display(mixt20) |
73 | 94 | Plot(mixt20) |
74 | 95 | Plot(ExtractDistribution(mixt20, "Mixture")) |
75 | 96 |
|
76 | | -mixt21 = Estimate(ExtractHistogram(vec10, 5), "MIXTURE", "NB", "NB", "NB", "NB", NbComponent="Estimated") |
77 | | - |
78 | | - |
79 | | - |
80 | | - |
| 97 | +mixt21 = Estimate( |
| 98 | + ExtractHistogram(vec10, 5), |
| 99 | + "MIXTURE", |
| 100 | + "NB", |
| 101 | + "NB", |
| 102 | + "NB", |
| 103 | + "NB", |
| 104 | + NbComponent="Estimated", |
| 105 | +) |
0 commit comments