def test_search_all_Cnumber_from_label(self): """ test for search_all_Cnumber_from_label """ expected = ([ 'C00001546', 'C00017854', 'C00017855', 'C00027106', 'C00036583', 'C00037061', 'C00045424', 'C00045427' ], 8) label = "N1b-C2c-S2a" actual = tools.search_all_Cnumber_from_label(label) self.assertEqual(actual, expected)
# 属中の種の数とか知りたい # classes.toolsを作りたい。 # In[30]: print(cag.get_Cnumber_from_label("C1a-C1z-C5x")) # In[31]: from tools import search_all_Cnumber_from_label from tools import get_genuses from tools import get_name # In[32]: ac = search_all_Cnumber_from_label("C1a-C1z-C5x") # In[35]: print(len(ac[0])) # In[36]: ac2 = search_all_Cnumber_from_label("C8-C8-C8-C8-C8-C8-C1-C8") # In[75]: print(len(ac2)) # In[37]:
for unit in units[:-1]: temp = unit.split() kosuu[(temp[1], temp[2])] = kosuu.get((temp[1], temp[2]),0) + int(temp[3]) sorted(kosuu.items(), reverse=True, key=lambda x: x[1]) # In[28]: print(len(list(spedict2.items())[0][1])) # 共発現しているのは14って書いてあるところ周辺 # In[29]: all1 = tools.search_all_Cnumber_from_label('C8x-C8x-C8x-C8x-C8y-C8y-C8y-N4x-C8y') print(len(all1)) # In[30]: all2 = tools.search_all_Cnumber_from_label('C8x-C8y-C8y-C8y-C8y-C8y-N4x-C8y-C8y') print(len(all2)) # In[31]: print(all2) # In[32]:
# In[29]: from rdkit.Chem import Draw img = Draw.MolsToGridImage(mol_list, legends=sorted(list(Ery.items())[0][1]), subImgSize=(400, 400)) img.save("Faba_test-2.png") # In[11]: from tools import search_all_Cnumber_from_label # In[12]: ac = search_all_Cnumber_from_label('C1x-N1y-C1z') # In[15]: print(len(ac[0])) # In[16]: from tools import make_kcfs # In[17]: print(Ery) # In[27]:
sorted(kosuu.items(), reverse=True, key=lambda x: x[1]) # In[80]: path = "Poaceae/Oryza/only" kosuu2 = dict() with open(path + "kcfscount.txt") as f: units = f.read().split("\n") for unit in units[:-1]: temp = unit.split() kosuu2[(temp[1], temp[2])] = kosuu2.get( (temp[1], temp[2]), 0) + int(temp[3]) # In[82]: sorted(kosuu2.items(), reverse=True, key=lambda x: x[1]) # In[72]: len(Cnlist) # In[84]: al = search_all_Cnumber_from_label("C1x-C1z-C2b") # In[85]: print(len(al)) # In[ ]:
temp = unit.split() kosuu[(temp[1], temp[2])] = kosuu.get( (temp[1], temp[2]), 0) + int(temp[3]) sorted(kosuu.items(), reverse=True, key=lambda x: x[1]) # In[15]: print(len(list(spedict2.items())[0][1])) # In[17]: import tools # In[21]: all1 = tools.search_all_Cnumber_from_label("C1x-C1x-C1x-C1z-C1y-C1z") print(all1) # 結構複雑な構造なのにこんなにあるのは予想外 # ringのみを抜き出すように変更したい。 # In[60]: all2 = tools.search_all_Cnumber_from_label( 'C1(C1)-C1-C1-C1-C1(C1)(C1)-C1-C1-C1-C1-C1-C2-C2') print(all2) print(len(all2)) # In[61]: all3 = tools.search_all_Cnumber_from_label('C1-C1-C1-C1-C1-N1-C1-C1')