fraglets-cpp/test.py at master · kolt-mcb/fraglets-cpp · GitHub

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
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
import fraglets
from lifelines import KaplanMeierFitter, NelsonAalenFitter
from multiprocessing import Pool
from matplotlib import pyplot as plt

def getIter(t):
    frag = getFrag()
    frag.run(1000000,2000,True)
    print(frag.iter)
    return frag.iter


def getFrag():
    frag = fraglets.fraglets()
    for i in range(20):
        frag.parse("[fork nop z match z split match z fork fork fork nop z * split match z fork fork fork nop z * z]")
    frag.parse("[perm z]")


    # alphabet = "abu"

    # permz = "perm z "
    # for symbol in alphabet:
    #     newMol = permz + symbol
    #     frag.parse(newMol)
    #     print(newMol)
    #     newMol = permz + "z " + symbol
    #     frag.parse(newMol)
    #     print(newMol)

    # permz = "perm z "
    # for symbol in frag.unimolTags:
    #     newMol = permz + symbol
    #     frag.parse(newMol)
    #     print(newMol)
    #     newMol = permz + "z " + symbol
    #     frag.parse(newMol)
    #     print(newMol)


    # for symbol in alphabet:
    #     newMol = permz + "match " + symbol
    #     frag.parse(newMol)
    #     print(newMol)
    #     newMol = permz + "z match " + symbol
    #     frag.parse(newMol)
    #     print(newMol)
    #     newMol = permz + "matchp " + symbol
    #     frag.parse(newMol)
    #     print(newMol)
    #     newMol = permz + "z matchp " + symbol
    #     frag.parse(newMol)
    #     print(newMol)

    return frag
    # frag.run(10000,200,True)

# pool = Pool(4)
# data = pool.map(getIter,range(100))

frags = []
data = []

for i in range(200):
    frag=getFrag()
    frag.run(1000,500,True)
    print(frag.iter)
    frags.append(frag)

greatest = frags[0]
for frag in frags:
    if frag.iter > greatest.iter:
        greatest = frag
    data.append(frag.iter)

greatest.drawGraphViz()


print(data)
kmf = KaplanMeierFitter()
kmf.fit(data)
kmf.plot()
plt.show()

naf = NelsonAalenFitter()
naf.fit(data)
naf.plot()
plt.show()