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

index 028c248..695ef26 100644 (file)
@@ -387,6 +387,13 @@ class PipelineHeuristic:
          trueWeightDon, trueWeightAcc, trueWeightIntron =\
          computeSpliceWeights(d, a, h, trueSpliceAlign, don_supp, acc_supp)
 
+         for i in xrange(0,len(trueWeightDon)):
+             trueWeightDon[0]=0.0
+         for i in xrange(0,len(trueWeightAcc)):
+             trueWeightAcc[0]=0.0
+         for i in xrange(0,len(trueWeightIntron)):
+             trueWeightIntron[0]=0.0
+         
          stop = cpu()
          self.computeSpliceWeightsTime += stop-start
 
@@ -433,9 +440,10 @@ class PipelineHeuristic:
 
          # new score
          DonorScore = calculatePlif(d, [don_score])[0]
-
+         print 'don: %f,%f' % (DonorScore, don_score)
+         
          score = computeSpliceAlignScoreWithQuality(original_est_cut, quality, qualityPlifs, run, self.currentPhi)
-         score += dummyAcceptorScore + IntronScore + DonorScore
+         #score += dummyAcceptorScore + IntronScore + DonorScore
          
          print 'diff %f,%f,%f' % ((trueWeight.T * self.currentPhi)[0,0] - score,(trueWeight.T * self.currentPhi)[0,0], score)