-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy pathThreshold.py
105 lines (101 loc) · 3.84 KB
/
Threshold.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
#!/usr/bin/env python
# -*- coding:utf-8-*-
# author: Liukewei time:2020/1/21 QQ:422209303 e-mail:Liukeweiaway@hotmail.com
# ----------------------------------------------------------------------------
def choiceThreshold(species, modification, threshold):
threshold = str(threshold)
if species == 'Human':
if modification == 'm1A':
if threshold == 'low': # low
M_threshold = 0.0237
elif threshold == 'normal': # normal
M_threshold = 0.0368
elif threshold == 'high': # high
M_threshold = 0.1071
elif modification == 'm6A':
if threshold == 'low':
M_threshold = 0.2282
elif threshold == 'normal':
M_threshold = 0.2596
elif threshold == 'high':
M_threshold = 0.2882
elif modification == 'm5C':
if threshold == 'low':
M_threshold = 0.7057
elif threshold == 'normal':
M_threshold = 0.8438
elif threshold == 'high':
M_threshold = 0.9052
elif modification == 'pseudouridine':
if threshold == 'low':
M_threshold = 0.7314
elif threshold == 'normal':
M_threshold = 0.7683
elif threshold == 'high':
M_threshold = 0.8137
elif modification == 'A-to-I':
if threshold == 'low':
M_threshold = 0.4484
elif threshold == 'normal':
M_threshold = 0.5286
elif threshold == 'high':
M_threshold = 0.6273
elif species == 'Yeast':
if modification == 'm1A':
if threshold == 'low':
M_threshold = 0.8039
elif threshold == 'normal':
M_threshold = 0.9419
elif threshold == 'high':
M_threshold = 0.9419
elif modification == 'm6A':
if threshold == 'low':
M_threshold = 0.7314
elif threshold == 'normal':
M_threshold = 0.8130
elif threshold == 'high':
M_threshold = 0.8839
elif modification == 'm5C':
if threshold == 'low':
M_threshold = 0.1429
elif threshold == 'normal':
M_threshold = 0.1838
elif threshold == 'high':
M_threshold = 0.7069
elif modification == 'pseudouridine':
if threshold == 'low':
M_threshold = 0.6633
elif threshold == 'normal':
M_threshold = 0.7364
elif threshold == 'high':
M_threshold = 0.7706
elif species == 'Mouse':
if modification == 'm1A':
if threshold == 'low':
M_threshold = 0.0667
elif threshold == 'normal':
M_threshold = 0.0979
elif threshold == 'high':
M_threshold = 0.1225
elif modification == 'm6A':
if threshold == 'low':
M_threshold = 0.5854
elif threshold == 'normal':
M_threshold = 0.6650
elif threshold == 'high':
M_threshold = 0.9196
elif modification == 'm5C':
if threshold == 'low':
M_threshold = 0.2834
elif threshold == 'normal':
M_threshold = 0.9792
elif threshold == 'high':
M_threshold = 0.9792
elif modification == 'pseudouridine':
if threshold == 'low':
M_threshold = 0.6723
elif threshold == 'normal':
M_threshold = 0.7058
elif threshold == 'high':
M_threshold = 0.7256
return M_threshold