fraglets-cpp/test2.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
import fraglets
from lifelines import KaplanMeierFitter, NelsonAalenFitter
from multiprocessing import Pool
from matplotlib import pyplot as plt


def getFrag():
    frag = fraglets.fraglets()
    for i in range(500):
        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)

frag = getFrag()

frag.run(1000000,1000,False)