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

index 770b77b..e38b5b8 100644 (file)
@@ -387,12 +387,12 @@ class PipelineHeuristic:
          trueWeightDon, trueWeightAcc, trueWeightIntron =\
          computeSpliceWeights(d, a, h, trueSpliceAlign, don_supp, acc_supp)
 
-         for i in xrange(0,len(trueWeightDon)):
-             trueWeightDon[i]=0.0
-         for i in xrange(0,len(trueWeightAcc)):
-             trueWeightAcc[i]=0.0
-         for i in xrange(0,len(trueWeightIntron)):
-             trueWeightIntron[i]=0.0
+         #for i in xrange(0,len(trueWeightDon)):
+         #    trueWeightDon[i]=0.0
+         #for i in xrange(0,len(trueWeightAcc)):
+         #    trueWeightAcc[i]=0.0
+         #for i in xrange(0,len(trueWeightIntron)):
+         #    trueWeightIntron[i]=0.0
          
          stop = cpu()
          self.computeSpliceWeightsTime += stop-start
@@ -445,7 +445,7 @@ class PipelineHeuristic:
          #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)