def RunBICTest(chr, beginSeqName, fileName, globalOrLocalTest, organism): """ Runs all BIC test calculations. chr: Current chromosome number. beginSeqName: Beginning of the sequence name of the queried chromosome. (e.g. "Gene") fileName: Output file name. globalOrLocalTest: String identifying whether to use global or local Mt values. organism: Organism's name. """ print ("BIC") if (globalOrLocalTest == "Global"): PalphaGammaHatAsym = printAndParseFiles.parsePalphaGamma("Calculations_" + organism + "\PalphaGammaHat%s\PalphaGammaHatAsymmetric%d.txt" % (globalOrLocalTest,chr), beginSeqName) PalphaGammaHatSym = printAndParseFiles.parsePalphaGamma("Calculations_" + organism + "\PalphaGammaHat%s\PalphaGammaHatSymmetric%d.txt" % (globalOrLocalTest,chr), beginSeqName) LgammaAsym = printAndParseFiles.parseValueFile("Calculations_" + organism + "\Lgamma%s\LgammaAsymmetric%d.txt" % (globalOrLocalTest,chr)) LgammaSym = printAndParseFiles.parseValueFile("Calculations_" + organism + "\Lgamma%s\LgammaSymmetric%d.txt" % (globalOrLocalTest,chr)) BICAsym = BIC.findBIC(PalphaGammaHatAsym, LgammaAsym, "A") BICSym = BIC.findBIC(PalphaGammaHatSym, LgammaSym, "S") printAndParseFiles.compareAndPrintBICs(BICAsym, BICSym, fileName) elif (globalOrLocalTest == "Local"): PalphaGammaHatAsym = printAndParseFiles.parsePalphaGamma("Calculations_" + organism + "\PalphaGammaHat%s\PalphaGammaHatAsymmetric%d.txt" % (globalOrLocalTest,chr), beginSeqName) PalphaGammaHatSym = printAndParseFiles.parsePalphaGamma("Calculations_" + organism + "\PalphaGammaHat%s\PalphaGammaHatSymmetric%d.txt" % (globalOrLocalTest,chr), beginSeqName) LgammaAsym = printAndParseFiles.parseValueFile("Calculations_" + organism + "\Lgamma%s\LgammaAsymmetric%d.txt" % (globalOrLocalTest,chr)) LgammaSym = printAndParseFiles.parseValueFile("Calculations_" + organism + "\Lgamma%s\LgammaSymmetric%d.txt" % (globalOrLocalTest,chr)) BICAsym = BIC.findBIC(PalphaGammaHatAsym, LgammaAsym, "A") BICSym = BIC.findBIC(PalphaGammaHatSym, LgammaSym, "S") printAndParseFiles.compareAndPrintBICs(BICAsym, BICSym, fileName) else: print "Error with Global or Local choice. Please type 'Global' or 'Local'" print ("DONE WITH CHR %d" % chr)