annotate rest_tool_comp_for_assay.py @ 4:54358dfa62c0 draft

Uploaded
author bernhardlutz
date Fri, 28 Mar 2014 13:41:22 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
4
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
1 #!/usr/bin/env python
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
2 # Aufruf convert_graph.py --aid list of ids --aid-from-file file
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
3
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
4 import sys, os
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
5 import networkx as nx
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
6 import argparse
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
7 import urllib2, urllib, httplib
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
8
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
9 import readfile
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
10 #supported graph_types
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
11 #output_types = ["tsv", "csv", "png", "json", "txt", "xml", "sdf", "asnt", "asnb", "jsonp"]
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
12
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
13
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
14 #get the cids for bioassay aid
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
15 def getCompoundList(aidlist):
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
16 aidliststring=",".join(aidlist)
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
17 url="http://pubchem.ncbi.nlm.nih.gov/rest/pug/assay/aid/"+aidliststring+"/cids/txt"
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
18 data=readfile.getresult(url)
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
19 return data
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
20
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
21 def main(args):
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
22 if args.aidfile is None:
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
23 aidlist=args.aid.split(",")
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
24 else:
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
25 aidlist=readfile.getListFromFile(args.aidfile)
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
26 cids=getCompoundList(aidlist)
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
27 args.outfile.write(cids)
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
28 args.outfile.close()
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
29
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
30
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
31 if __name__ == "__main__":
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
32 parser = argparse.ArgumentParser()
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
33 parser.add_argument('--aid', type=str,
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
34 help="AIDs of the BioAssay")
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
35 parser.add_argument('--aidfile', type=argparse.FileType('r'),
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
36 help="Specify a file with a list of aids, one per line")
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
37 parser.add_argument('--outfile', type=argparse.FileType('w'),
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
38 help="Specify output file")
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
39 if len(sys.argv) < 2:
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
40 print "Too few arguments..."
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
41 parser.print_help()
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
42 exit(1)
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
43 args = parser.parse_args()
54358dfa62c0 Uploaded
bernhardlutz
parents:
diff changeset
44 main( args )