git-svn-id: http://svn.tuebingen.mpg.de/ag-raetsch/projects/QPalma@8646 e1793c9e...
authorraetsch <raetsch@e1793c9e-67f9-0310-80fc-b846ff1f7b36>
Thu, 17 Apr 2008 21:58:24 +0000 (21:58 +0000)
committerraetsch <raetsch@e1793c9e-67f9-0310-80fc-b846ff1f7b36>
Thu, 17 Apr 2008 21:58:24 +0000 (21:58 +0000)
scripts/PipelineHeuristic.py

index 8e878e2..8b5e640 100644 (file)
@@ -315,7 +315,9 @@ class PipelineHeuristic:
       dna   = currentDNASeq
 
       alt_don,alt_acc = self.findHighestScoringSpliceSites(currentAcc,currentDon)
-      max_acc=max(currentAcc) ;
+
+      sort_acc=enumerate(currentAcc).sort(lambda x,y: x[1]-y[1]) 
+      (max_acc_pos,max_acc_score)=sort_acc[-1]
       
       #print alt_don
       #print alt_acc
@@ -331,9 +333,9 @@ class PipelineHeuristic:
       # inlined
 
       # compute dummy scores
-      IntronScore = calculatePlif(h, [self.intron_size+1])[0]
+      IntronScore = calculatePlif(h, [max_acc_pos-30])[0]
       #print IntronScore
-      dummyAcceptorScore = calculatePlif(a, [max_acc])[0] 
+      dummyAcceptorScore = calculatePlif(a, [max_acc_score])[0] 
       dummyDonorScore = calculatePlif(d, [0.25])[0]
       
       _start = cpu()