for i in range(1,len(yeastset)):
tree = yeastset[i]
treefiles = [tree + "_bayes_treeout.txt",tree + "_raxml_treeout.txt",tree + "_BEAST_treeout.txt",tree + "_phyml_treeout.txt"]
for m in range(len(models)):
topo_data_file = topo_data_dir + "/yeast_" + models[m] + "_" + str(topoPicked) + "_topomatches.txt"
treehome = basedir + "/" + yeastset[i] + "/" + models[m] + "/"
infile = open(topo_data_file,'r')
sameTop = 1
for line in infile:
temp = line.split()
if (temp[0] == "'" + str(i) + "'"):
if (temp[1] == "TRUE"):
sameTop = 1
break
elif (temp[2] == "TRUE"):
sameTop = 2
break
logplot = LogMapPlotter(tree, treefiles[m], treehome, models[m], rooted=False)
logplot.split_topos()
logplot.make_centres(sameTop)
logplot.make_coords(sameTop)
#os.remove(tree + str(i) + "_" + models[m] + "_allcoords.r")
#os.remove(tree + str(i) + "_" + models[m] + "_allcoords_norm3D.r")
logplot.plot_coords_all()
logplot.plot_coords_quads()
# Next step: figure out edges for a 'progression' in topologies
