# HG changeset patch # User nanette # Date 1377002207 14400 # Node ID 900854ec76e6627bf674535bf29a22adf72a1cbe # Parent a0cbbd0dfcc00e3ccd04381f2a7275eb611b4eaa Uploaded diff -r a0cbbd0dfcc0 -r 900854ec76e6 test_script_path.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test_script_path.py Tue Aug 20 08:36:47 2013 -0400 @@ -0,0 +1,49 @@ +""" +@summary: GO enrichment analysis (hotspots) +@author: nanette.coetzer@gmail.com +@version 5 + +""" +import optparse, sys +import subprocess +import tempfile +import os, re + +def stop_err( msg ): + sys.stderr.write( "%s\n" % msg ) + sys.exit() + +def __main__(): + #Parse Command Line + parser = optparse.OptionParser() + parser.add_option("-i", "--input1", default=None, dest="input1", + help="genes") + parser.add_option("-o", "--output1", default=None, dest="output1", + help="star genes") + + (options, args) = parser.parse_args() + + try: + open(options.input1, "r").close() + except TypeError, e: + stop_err("You need to supply the Gene Universe file:\n" + str(e)) + except IOError, e: + stop_err("Can not open the Gene Universe file:\n" + str(e)) + + + ########################################################## + + infile = open(options.input1,"r") + inlist = [] + for line in infile: + inlist.append(line.strip()) + infile.close() + outfile = open(options.output1,"w") + for l in inlist: + outfile.write("* "+str(l)+"\n") + outfile.close() + + ############################################## + +if __name__=="__main__": + __main__()