Example #1
0
 def count(self,sequence1,sequence2):
     dec1 = FrequenceAnalysis.FrequenceAnalysis(sequence1)
     dec2 = FrequenceAnalysis.FrequenceAnalysis(sequence2)
     '''
     probability no substitution occurs
     '''
     #sites_count = len(min([dec1,dec2], key = len))
     sites_count = len(dec1)
     different_sites = 0
     for e1,e2 in zip(dec1,dec2):
         different_sites += (1- Model.dotProduct(e1,e2))
     #return (-float(3)/4)* log1p(-(float(4)/3)*different_sites/sites_count)      
     return different_sites/sites_count