|
8 | 8 | from openalea.sequence_analysis import * |
9 | 9 | from openalea.sequence_analysis.estimate import Estimate |
10 | 10 | from openalea.sequence_analysis.compare import Compare |
11 | | -from openalea.sequence_analysis import get_shared_data |
| 11 | +from tools import runTestClass, robust_path as get_shared_data |
12 | 12 |
|
13 | 13 |
|
14 | | -seq20 = Sequences(get_shared_data("belren1.seq")) |
15 | | -seq21 = Sequences(get_shared_data("elstar1.seq")) |
16 | | -seq22 = Sequences(get_shared_data("fuji1.seq")) |
17 | | -seq23 = Sequences(get_shared_data("gala1.seq")) |
18 | | -seq24 = Sequences(get_shared_data("granny1.seq")) |
19 | | -seq25 = Sequences(get_shared_data("reinet1.seq")) |
20 | | -seq26 = Sequences(get_shared_data("wij1.seq")) |
| 14 | +seq20 = Sequences(str(get_shared_data("belren1.seq"))) |
| 15 | +seq21 = Sequences(str(get_shared_data("elstar1.seq"))) |
| 16 | +seq22 = Sequences(str(get_shared_data("fuji1.seq"))) |
| 17 | +seq23 = Sequences(str(get_shared_data("gala1.seq"))) |
| 18 | +seq24 = Sequences(str(get_shared_data("granny1.seq"))) |
| 19 | +seq25 = Sequences(str(get_shared_data("reinet1.seq"))) |
| 20 | +seq26 = Sequences(str(get_shared_data("wij1.seq"))) |
21 | 21 |
|
22 | 22 | Display(seq25, ViewPoint="Data") |
23 | 23 | Plot(seq25, "Intensity") |
|
52 | 52 | matrix30 = Compare(seq32) |
53 | 53 | matrix31 = Compare(seq32, VectorDistance("S")) |
54 | 54 | matrix32 = Compare(seq32, VectorDistance("S"), Transposition=True) |
55 | | -matrix33 = Compare(seq32, VectorDistance(get_shared_data("test_align1.a")), Transposition=True) |
| 55 | +matrix33 = Compare(seq32, VectorDistance(str(get_shared_data("test_align1.a"))), Transposition=True) |
56 | 56 |
|
57 | 57 | Display(Clustering(matrix33, "Partition", 2)) |
58 | 58 | Clustering(matrix33, "Hierarchy", Algorithm="Agglomerative") |
|
76 | 76 | Plot(seq25, 14, 6, "Multinomial", ViewPoint="SegmentProfile", Output="ChangePoint") |
77 | 77 | # hidden semi-Markov chains |
78 | 78 |
|
79 | | -hsmc0 = HiddenSemiMarkov(get_shared_data("belren1.hsc")) |
| 79 | +hsmc0 = HiddenSemiMarkov(str(get_shared_data("belren1.hsc"))) |
80 | 80 | hsmc20 = Estimate(seq20, "HIDDEN_SEMI-MARKOV", hsmc0) |
81 | 81 |
|
82 | | -hsmc0 = HiddenSemiMarkov(get_shared_data("elstar1.hsc")) |
| 82 | +hsmc0 = HiddenSemiMarkov(str(get_shared_data("elstar1.hsc"))) |
83 | 83 | hsmc21 = Estimate(seq21, "HIDDEN_SEMI-MARKOV", hsmc0) |
84 | 84 |
|
85 | | -hsmc0 = HiddenSemiMarkov(get_shared_data("fuji1.hsc")) |
| 85 | +hsmc0 = HiddenSemiMarkov(str(get_shared_data("fuji1.hsc"))) |
86 | 86 | hsmc22 = Estimate(seq22, "HIDDEN_SEMI-MARKOV", hsmc0) |
87 | 87 |
|
88 | | -hsmc0 = HiddenSemiMarkov(get_shared_data("gala1.hsc")) |
| 88 | +hsmc0 = HiddenSemiMarkov(str(get_shared_data("gala1.hsc"))) |
89 | 89 | hsmc23 = Estimate(seq23, "HIDDEN_SEMI-MARKOV", hsmc0) |
90 | 90 |
|
91 | | -hsmc0 = HiddenSemiMarkov(get_shared_data("granny1.hsc")) |
| 91 | +hsmc0 = HiddenSemiMarkov(str(get_shared_data("granny1.hsc"))) |
92 | 92 | hsmc24 = Estimate(seq24, "HIDDEN_SEMI-MARKOV", hsmc0) |
93 | 93 |
|
94 | | -hsmc0 = HiddenSemiMarkov(get_shared_data("reinet1.hsc")) |
| 94 | +hsmc0 = HiddenSemiMarkov(str(get_shared_data("reinet1.hsc"))) |
95 | 95 | hsmc25 = Estimate(seq25, "HIDDEN_SEMI-MARKOV", hsmc0) |
96 | 96 |
|
97 | 97 | Display(hsmc25) |
|
112 | 112 | seq25_1 = ExtractData(hsmc25) |
113 | 113 | Display(seq25_1, ViewPoint="Data", Format="Line") |
114 | 114 |
|
115 | | -hsmc0 = HiddenSemiMarkov(get_shared_data("wij1.hsc")) |
| 115 | +hsmc0 = HiddenSemiMarkov(str(get_shared_data("wij1.hsc"))) |
116 | 116 | hsmc26 = Estimate(seq26, "HIDDEN_SEMI-MARKOV", hsmc0) |
117 | 117 |
|
118 | 118 | # model comparison |
|
0 commit comments