Stats.py
1.7 KB
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
from math import sqrt, ceil
import numpy as np
import matplotlib.pyplot as plt
file = open("/mnt/c/Users/natha/Documents/IBISC/biorseo2/biorseo/data/modules/ISAURE/Motifs_version_initiale/benchmark.dbn", "r")
name = file.readline()
rna = file.readline()
twod = file.readline()
contacts = file.readline()
length = len(rna)
nb_ctc = contacts.count('*')
print("--------------------------------------------------------")
ctc_max = nb_ctc
ctc_min = nb_ctc
np_lgt = []
np_lgt.append(length)
np_ctc = []
np_ctc.append(nb_ctc)
np = []
np.append([length, nb_ctc])
while name:
print(contacts)
print(length)
print(nb_ctc)
print("--------------------------------------------------------")
name = file.readline()
rna = file.readline()
length = len(rna)
if length != 0 :
np_lgt.append(length)
twod = file.readline()
contacts = file.readline()
nb_ctc = contacts.count('*')
if nb_ctc != 0:
np_ctc.append(nb_ctc)
np.append([length, nb_ctc])
if nb_ctc > ctc_max:
ctc_max = nb_ctc
if nb_ctc < ctc_min and nb_ctc != 0:
ctc_min = nb_ctc
file.close()
print(np_lgt)
print(np_ctc)
print(np)
x = np_lgt
y = np_ctc
index = np_ctc.index(ctc_max)
index2 = np_ctc.index(ctc_min)
plt.scatter(x, y, c = 'blue')
plt.annotate("(" + str(np_lgt[index]) + "," + str(ctc_max) + ")", (np_lgt[index], ctc_max),c ='red')
plt.scatter(np_lgt[index], ctc_max,c = 'red')
plt.annotate("(" + str(np_lgt[index2]) + "," + str(ctc_min) + ")", (np_lgt[index2], ctc_min),c ='green')
plt.scatter(np_lgt[index2], ctc_min,c = 'green')
plt.xlabel('longeur de l\'arn')
plt.ylabel('nombre de contacts')
plt.savefig('stats.png')