annotate get_data/kegg_glycan/test_findKEGG.py @ 1:0a5e0df17054 draft default tip

Uploaded
author chrisb
date Fri, 06 May 2016 08:05:48 -0400
parents 89592faa2875
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
89592faa2875 Uploaded
chrisb
parents:
diff changeset
1 __license__ = "MIT"
89592faa2875 Uploaded
chrisb
parents:
diff changeset
2
89592faa2875 Uploaded
chrisb
parents:
diff changeset
3 import unittest
89592faa2875 Uploaded
chrisb
parents:
diff changeset
4 import findKEGG as fk
89592faa2875 Uploaded
chrisb
parents:
diff changeset
5
89592faa2875 Uploaded
chrisb
parents:
diff changeset
6
89592faa2875 Uploaded
chrisb
parents:
diff changeset
7 class SimpleUnitTest(unittest.TestCase):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
8 def setUp(self):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
9 import os
89592faa2875 Uploaded
chrisb
parents:
diff changeset
10
89592faa2875 Uploaded
chrisb
parents:
diff changeset
11 os.environ["http_proxy"] = "" # work around for IOError: [Errno url error] invalid proxy for http:
89592faa2875 Uploaded
chrisb
parents:
diff changeset
12 pass
89592faa2875 Uploaded
chrisb
parents:
diff changeset
13
89592faa2875 Uploaded
chrisb
parents:
diff changeset
14 def tearDown(self):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
15 pass
89592faa2875 Uploaded
chrisb
parents:
diff changeset
16
89592faa2875 Uploaded
chrisb
parents:
diff changeset
17 def test_no_db_specified(self):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
18 """should raise error"""
89592faa2875 Uploaded
chrisb
parents:
diff changeset
19 with self.assertRaises(IOError):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
20 m = fk.find_entries_in_KEGG([], " ")
89592faa2875 Uploaded
chrisb
parents:
diff changeset
21 with self.assertRaises(IOError):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
22 m = fk.find_entries_in_KEGG("", " ")
89592faa2875 Uploaded
chrisb
parents:
diff changeset
23 with self.assertRaises(IOError):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
24 m = fk.find_entries_in_KEGG(None, " ")
89592faa2875 Uploaded
chrisb
parents:
diff changeset
25
89592faa2875 Uploaded
chrisb
parents:
diff changeset
26 def test_no_query_specified(self):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
27 """should raise error"""
89592faa2875 Uploaded
chrisb
parents:
diff changeset
28 with self.assertRaises(IOError):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
29 m = fk.find_entries_in_KEGG("glycan", "")
89592faa2875 Uploaded
chrisb
parents:
diff changeset
30 with self.assertRaises(IOError):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
31 m = fk.find_entries_in_KEGG("glycan", [])
89592faa2875 Uploaded
chrisb
parents:
diff changeset
32 with self.assertRaises(IOError):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
33 m = fk.find_entries_in_KEGG("glycan", None)
89592faa2875 Uploaded
chrisb
parents:
diff changeset
34
89592faa2875 Uploaded
chrisb
parents:
diff changeset
35
89592faa2875 Uploaded
chrisb
parents:
diff changeset
36 def test_bad_db(self):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
37 """should raise error"""
89592faa2875 Uploaded
chrisb
parents:
diff changeset
38 import urllib2
89592faa2875 Uploaded
chrisb
parents:
diff changeset
39
89592faa2875 Uploaded
chrisb
parents:
diff changeset
40 with self.assertRaises(urllib2.HTTPError):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
41 m = fk.find_entries_in_KEGG("john", "glucose")
89592faa2875 Uploaded
chrisb
parents:
diff changeset
42
89592faa2875 Uploaded
chrisb
parents:
diff changeset
43 def test_unfindable_entry(self):
1
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
44 """should return an empty string"""
0
89592faa2875 Uploaded
chrisb
parents:
diff changeset
45 m = fk.find_entries_in_KEGG("glycan", "sally")
1
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
46 emptystring=""
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
47 self.assertItemsEqual(m,emptystring,"Expected Empty String for non-existent entry")
0
89592faa2875 Uploaded
chrisb
parents:
diff changeset
48
89592faa2875 Uploaded
chrisb
parents:
diff changeset
49 def test_find_example(self):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
50 m = fk.find_entries_in_KEGG("glycan", "glucose")
89592faa2875 Uploaded
chrisb
parents:
diff changeset
51 self.assertIsNotNone(m)
89592faa2875 Uploaded
chrisb
parents:
diff changeset
52 self.assertIn("GDP-glucose", m)
89592faa2875 Uploaded
chrisb
parents:
diff changeset
53
89592faa2875 Uploaded
chrisb
parents:
diff changeset
54 def test_query_has_newlines(self):
89592faa2875 Uploaded
chrisb
parents:
diff changeset
55 """
1
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
56 Assume a new line query is an AND query.
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
57 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
58 m = fk.find_entries_in_KEGG("glycan", "glucose\nUDP")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
59 self.assertIsNotNone(m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
60 print m
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
61 self.assertIn("UDP-D-glucose", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
62
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
63 def test_query_has_newlines_and_a_space(self):
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
64 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
65 glucose\\n UDP. OR's should pass
0
89592faa2875 Uploaded
chrisb
parents:
diff changeset
66 """
89592faa2875 Uploaded
chrisb
parents:
diff changeset
67 m = fk.find_entries_in_KEGG("glycan", "glucose\n UDP")
89592faa2875 Uploaded
chrisb
parents:
diff changeset
68 self.assertIsNotNone(m)
89592faa2875 Uploaded
chrisb
parents:
diff changeset
69 self.assertIn("UDP-D-glucose", m)
89592faa2875 Uploaded
chrisb
parents:
diff changeset
70
1
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
71 def test_gene_specific_or(self):
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
72 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
73 test "shiga toxin" returns.
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
74 /find/genes/"shiga toxin" for keywords "shiga toxin"
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
75 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
76 m = fk.find_entries_in_KEGG("genes","shiga toxin")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
77 self.assertIsNotNone(m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
78 #. not matching to all entries - that is silly, the db get updated!
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
79 self.assertIn("shiga", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
80 self.assertIn("toxin", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
81 self.assertNotIn("stm:STM0284",m) # should not be found in the OR query
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
82
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
83 def test_gene_specific_and(self):
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
84 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
85 test "shiga+toxin" returns. it should
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
86 /find/genes/shiga+toxin for keywords "shiga" and "toxin"
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
87 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
88 m = fk.find_entries_in_KEGG("genes","shiga+toxin")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
89 self.assertIsNotNone(m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
90 #. not matching to all entries - that is silly, the db get updated!
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
91 self.assertIn("shiga", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
92 self.assertIn("toxin", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
93 self.assertIn("stm:STM0284",m) # should be found in the AND query
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
94
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
95 def test_enzyme_specific_1(self):
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
96 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
97 test 2.4.99.1 returns
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
98 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
99 m = fk.find_entries_in_KEGG("enzyme","2.4.99.1")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
100 self.assertIsNotNone(m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
101 #. not matching to all entries - that is silly, the db get updated!
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
102 self.assertIn("2.4.99.1", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
103
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
104 def test_enzyme_specific_2(self):
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
105 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
106 test 2.4.99.6 returns
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
107 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
108 m = fk.find_entries_in_KEGG("enzyme","2.4.99.6")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
109 self.assertIsNotNone(m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
110 #. not matching to all entries - that is silly, the db get updated!
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
111 self.assertIn("2.4.99.6", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
112
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
113 def test_enzyme_specific_3(self):
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
114 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
115 test ec: 2.4.99.6 with space (tests stripping space functionality)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
116 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
117 m = fk.find_entries_in_KEGG("enzyme","ec: 2.4.99.6")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
118 self.assertIsNotNone(m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
119 #. not matching to all entries - that is silly, the db get updated!
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
120 self.assertIn("2.4.99.6", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
121
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
122 def test_enzyme_specific_4(self):
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
123 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
124 test ec: 2.4.99. , should return 2.4.99.*
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
125 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
126 m = fk.find_entries_in_KEGG("enzyme","ec: 2.4.99.")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
127 self.assertIsNotNone(m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
128 #. not matching to all entries - that is silly, the db get updated!
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
129 self.assertIn("2.4.99.1", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
130 self.assertIn("2.4.99.6", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
131 self.assertIn("2.4.99.11", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
132
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
133 def test_enzyme_specific_5(self):
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
134 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
135 test ec: 2.4.99.1+2.4.99.6 . AND function, so returns nothing here. Do not expect an enzyme to be doubly classified.
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
136 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
137 m = fk.find_entries_in_KEGG("enzyme","2.4.99.1+2.4.99.6")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
138 self.assertIsNotNone(m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
139 emptystring=""
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
140 self.assertItemsEqual(m,emptystring,"Expected Empty String for non-existent entry")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
141
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
142 def test_enzyme_specific_6(self):
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
143 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
144 test ec: 2.4.99.1 2.4.99.6 . OR function, should return something but does not (KEGG FAILURE?).
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
145 note that trying this without "" will return the AND function result which is incorrect.
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
146 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
147 m = fk.find_entries_in_KEGG("enzyme","2.4.99.1 2.4.99.6")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
148 self.assertIsNotNone(m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
149 #self.assertIn("2.4.99.1", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
150 emptystring=""
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
151 print m
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
152 self.assertItemsEqual(m,emptystring,"Expected Empty String for non-existent entry")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
153
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
154
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
155 def test_enzyme_specific_7(self):
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
156 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
157 test deoxy+1.1.1. . AND function, that should return something
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
158 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
159 m = fk.find_entries_in_KEGG("enzyme","deoxy+1.1.1")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
160 self.assertIsNotNone(m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
161 print m
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
162 self.assertIn("deoxy", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
163 self.assertIn("1.1.1", m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
164
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
165 def test_enzyme_specific_8(self):
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
166 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
167 test deoxy 1.1.1. . OR function. should return deoxy or 1.1.1. but does not (KEGG FAILURE?). i
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
168 note that trying this without "" will return the AND function result which is incorrect.
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
169 """
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
170 m = fk.find_entries_in_KEGG("enzyme","deoxy 1.1.1")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
171 self.assertIsNotNone(m)
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
172 emptystring=""
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
173 print m
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
174 self.assertItemsEqual(m,emptystring,"Expected Empty String for non-existent entry")
0a5e0df17054 Uploaded
chrisb
parents: 0
diff changeset
175