def createQQplots(tree, treehome, model1, model2,xlow,xhigh,ylow,yhigh): sturm_version = "SturmMean_130704.jar" analysis_jar = "analysis_140702.jar" model1out = treeout(tree, model1, treehome) model2out = treeout(tree, model2, treehome) #shutil.copyfile(treehome + '/' + model1 + '/' + model1out, treehome + '/' + model1out) #shutil.copyfile(treehome + '/' + model2 + '/' + model2out, treehome + '/' + model2out) t1root = False t2root = False if (model1 == "BEAST"): t1root = True print "model1 is BEAST" if (model2 == "BEAST"): t2root = True print "model2 is BEAST" treeplot = distribcompare(tree + "_" + model1 + "_" + model2, model1out, model2out, treehome, t1root, t2root) treeplot.find_means() treeplot.find_disttomean()
def createQQplots(tree, treehome, model1, model2,xlow=0,xhigh=1,ylow=0,yhigh=1): sturm_version = "SturmMean_130704.jar" analysis_jar = "analysis_140702.jar" shutil.copyfile('c:\\seqgen\\jars\\' + sturm_version, treehome + "\\" + sturm_version) shutil.copyfile('c:\\seqgen\\jars\\' + analysis_jar, treehome + "\\" + analysis_jar) model1out = treeout(tree, model1, treehome) model2out = treeout(tree, model2, treehome) shutil.copyfile(treehome + '/' + model1 + '/' + model1out, treehome + '/' + model1out) shutil.copyfile(treehome + '/' + model2 + '/' + model2out, treehome + '/' + model2out) t1root = False t2root = False if (model1 == "BEAST"): t1root = True print "model1 is BEAST" if (model2 == "BEAST"): t2root = True print "model2 is BEAST" treeplot = distribcompare(tree + "_" + model1 + "_" + model2, model1out, model2out, treehome, t1root, t2root) treeplot.find_means() treeplot.find_disttomean() treeplot.set_dims(xlow,xhigh,ylow,yhigh) treeplot.qq_plot()
''' Created on Jun 24, 2015 @author: alkpongsema ''' from distribcompare import distribcompare if __name__ == '__main__': #treeplot = distribcompare("newbasetreeout.txt", "basenewtreeout.txt") #treeplot = distribcompare("basenew_raxmlout.txt","newbase_raxmlout.txt") #treeplot = distribcompare("newbasetreeout.txt","newbase_raxmlout.txt") treeplot = distribcompare("basenewtreeout.txt","basenew_raxmlout.txt") treeplot.find_means() treeplot.find_disttomean() treeplot.qq_plot()