In [1]:

    

import dendropy
from dendropy.utility.fileutils import find_files


    

In [3]:

    

pb = dendropy.Tree.get_from_path('./Trees/garli_opt/ranked/2598501_pb', 'nexus')


    

In [6]:

    

flist = find_files(top='Trees/garli_opt/ranked', filename_filter='*')
del flist[9]


    

In [22]:

    

sqtrees = [dendropy.Tree.get_from_path(filename,"nexus") for filename in flist]


    

In [ ]:

    




    

In [23]:

    

oneb = sqtrees[0]
one_list = [tree.symmetric_difference(oneb) for tree in sqtrees]


    

In [24]:

    

one_list


    

    
Out[24]:





[0, 736, 648, 732, 0, 734, 709, 728, 654, 728, 694]

In [28]:

    

twob = sqtrees[1]
two_list = [tree.symmetric_difference(twob) for tree in sqtrees]


    

In [29]:

    

two_list


    

    
Out[29]:





[736, 0, 716, 72, 736, 708, 653, 644, 720, 70, 724]

In [30]:

    

threeb = sqtrees[2]
three_list = [tree.symmetric_difference(threeb) for tree in sqtrees]
three_list


    

    
Out[30]:





[648, 716, 0, 728, 648, 764, 705, 726, 10, 730, 732]

In [31]:

    

fourb = sqtrees[3]
four_list = [tree.symmetric_difference(fourb) for tree in sqtrees]
four_list


    

    
Out[31]:





[732, 72, 728, 0, 732, 728, 667, 652, 732, 8, 732]

In [34]:

    

fiveb = sqtrees[4]
five_list = [tree.symmetric_difference(fiveb) for tree in sqtrees]
five_list


    

    
Out[34]:





[0, 736, 648, 732, 0, 734, 709, 728, 654, 728, 694]

In [36]:

    

sixb = sqtrees[5]
six_list = [tree.symmetric_difference(sixb) for tree in sqtrees]
six_list


    

    
Out[36]:





[734, 708, 764, 728, 734, 0, 729, 706, 772, 728, 724]

In [18]:

    

pb = dendropy.Tree.get_from_path('Trees/garli_opt/pb', 'nexus')
mle = dendropy.Tree.get_from_path('Trees/garli_opt/2b', 'nexus')


    

In [19]:

    

pb.symmetric_difference(mle)


    

    
Out[19]:





732

In [ ]: