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()
