import sys from PIPRes.util.crude_xml import toPipresXML '''Prints out 8 taxon "Felsenstein zone cherry tree". Args are: TwoEdgeLength ThreeEdgeLength CherryLength A Felsenstein zone tree with TwoEdgeLength and ThreeEdgeLength is created, then each tips is replace with a cherry (each edge of the cherry has length CherryLength).''' twoEdge = float(sys.argv[1]) threeEdge = float(sys.argv[2]) cherryEdge = float(sys.argv[3]) class Generic: pass g = Generic() g.m_newick = '((1:%g,2:%g):%g,(3:%g,4:%g):%g,((5:%g,6:%g):%g,(7:%g,8:%g):%g):%g);' % (cherryEdge, cherryEdge, twoEdge, cherryEdge, cherryEdge, threeEdge, cherryEdge, cherryEdge, threeEdge, cherryEdge, cherryEdge, twoEdge, threeEdge) print toPipresXML([ ['tree', g], ['split', 255L - 3L], ['split', 12L], ['split', 48L], ['split', 192L], ['split', 255L - 195L], ['split', 255L - 15L], ]),