import sys
import optparse
-from commoncode import getGeneInfoDict, getConfigParser, getConfigOption
+from commoncode import getGeneInfoDict, getConfigParser, getConfigOption, getConfigIntOption
from cistematic.cisstat.analyzego import calculateGOStats
print "analyzego: version 2.2"
infilename = args[1]
prefix = args[2]
- analyzeGOFromFile(genome, infilename, prefix, options.translateGene, fieldID)
+ analyzeGOFromFile(genome, infilename, prefix, options.translateGene, fieldID, numTests=options.numTests)
def makeParser(usage=""):
help="translate gene")
parser.add_option("--field", type="int", dest="fieldID",
help="column containing gene ID/Name")
+ parser.add_option("--trials", type="int", dest="numTests",
+ help="column containing gene ID/Name")
configParser = getConfigParser()
section = "analyzego"
translateGene = getConfigOption(configParser, section, "translateGene", False)
fieldID = getConfigOption(configParser, section, "fieldID", None)
+ numTests = getConfigIntOption(configParser, section, "numTests", 1)
- parser.set_defaults(translateGene=translateGene, fieldID=fieldID)
+ parser.set_defaults(translateGene=translateGene, fieldID=fieldID, numTests=numTests)
return parser
-def analyzeGOFromFile(genome, infilename, prefix, translateGene=False, fieldID=1):
+def analyzeGOFromFile(genome, infilename, prefix, translateGene=False, fieldID=1, numTests=1):
infile = open(infilename)
- analyzeGO(genome, infile, prefix, translateGene=False, fieldID=1)
+ analyzeGO(genome, infile, prefix, translateGene, fieldID, numTests)
infile.close()
-def analyzeGO(genome, geneInfoList, prefix, translateGene=False, fieldID=1):
+def analyzeGO(genome, geneInfoList, prefix, translateGene=False, fieldID=1, numTests=1):
if translateGene:
symbolToGidDict = getSymbolDict(genome)
locusList.append((genome, gID))
if len(locusList) > 0:
- calculateGOStats(locusList, prefix)
+ calculateGOStats(locusList, prefix, trials=numTests)
def getSymbolDict(genome):