-
Notifications
You must be signed in to change notification settings - Fork 0
/
plotLymphNode.py
118 lines (94 loc) · 2.9 KB
/
plotLymphNode.py
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
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
import matplotlib.pyplot as plt
import numpy as np
import sys
# sns.set()
T_final = int(sys.argv[1])
h_t = float(sys.argv[2])
populationTitle = {
"thp": "Linfonodo - T $CD4^+$",
"adc": "Linfonodo - Células Dendríticas Ativadas",
"bcl": "Linfonodo - Células B",
"pcl": "Linfonodo - Plasmócitos",
"tke": "Linfonodo - T $CD8^+$",
"ant": "Linfonodo - Anticorpos IgG"
}
t = np.linspace(0, T_final, int(T_final/h_t))
TL_c_vetor = np.zeros(len(t))
TL_h_vetor = np.zeros(len(t))
B_vetor = np.zeros(len(t))
FL_vetor = np.zeros(len(t))
PL_vetor = np.zeros(len(t))
DL_vetor = np.zeros(len(t))
with open("./result/tCyto.txt", 'r') as f:
lines = f.readlines()
TL_c = [line.rstrip() for line in lines]
with open("./result/tHelper.txt", 'r') as f:
lines = f.readlines()
TL_h = [line.rstrip() for line in lines]
with open("./result/bCell.txt", 'r') as f:
lines = f.readlines()
BV = [line.rstrip() for line in lines]
with open("./result/antibody.txt", 'r') as f:
lines = f.readlines()
FL = [line.rstrip() for line in lines]
with open("./result/plasmaCell.txt", 'r') as f:
lines = f.readlines()
PL = [line.rstrip() for line in lines]
with open("./result/dendritic.txt", 'r') as f:
lines = f.readlines()
DL = [line.rstrip() for line in lines]
for i in range(0,len(TL_c)):
TL_c_vetor[i] = TL_c[i]
TL_h_vetor[i] = TL_h[i]
B_vetor[i] = BV[i]
FL_vetor[i] = FL[i]
PL_vetor[i] = PL[i]
DL_vetor[i] = DL[i]
ax = plt.gca()
ax.set_xticks([0,7,14,21,28])
plt.plot(t,TL_c_vetor, "-r")
plt.title(populationTitle["tke"])
plt.xlabel("Tempo (dias)")
plt.ylabel("Concentração (Células/$mm^2$)")
plt.savefig('result/t_cito_linfonodo.png', dpi = 300)
plt.clf()
ax = plt.gca()
ax.set_xticks([0,7,14,21,28])
plt.plot(t,TL_h_vetor, "-r")
plt.title(populationTitle["thp"])
plt.xlabel("Tempo (dias)")
plt.ylabel("Concentração (Células/$mm^2$)")
plt.savefig('result/t_helper_linfonodo.png', dpi = 300)
plt.clf()
ax = plt.gca()
ax.set_xticks([0,7,14,21,28])
plt.plot(t,B_vetor, "-r")
plt.title(populationTitle["bcl"])
plt.xlabel("Tempo (dias)")
plt.ylabel("Concentração (Células/$mm^2$)")
plt.savefig('result/b_cell_linfonodo.png', dpi = 300)
plt.clf()
ax = plt.gca()
ax.set_xticks([0,7,14,21,28])
plt.plot(t,FL_vetor, "-r")
plt.title(populationTitle["ant"])
plt.xlabel("Tempo (dias)")
plt.ylabel("Concentração (Células/$mm^2$)")
plt.savefig('result/anticorpo_linfonodo.png', dpi = 300)
plt.clf()
ax = plt.gca()
ax.set_xticks([0,7,14,21,28])
plt.plot(t,PL_vetor, "-r")
plt.title(populationTitle["pcl"])
plt.xlabel("Tempo (dias)")
plt.ylabel("Concentração (Células/$mm^2$)")
plt.savefig('result/pl_cell_linfonodo.png', dpi = 300)
plt.clf()
ax = plt.gca()
ax.set_xticks([0,7,14,21,28])
plt.plot(t,DL_vetor, "-r")
plt.title(populationTitle["adc"])
plt.xlabel("Tempo (dias)")
plt.ylabel("Concentração (Células/$mm^2$)")
plt.savefig('result/dc_linfonodo.png', dpi = 300)
plt.clf()