git-svn-id: http://svn.tuebingen.mpg.de/ag-raetsch/projects/QPalma@8604 e1793c9e...
authorraetsch <raetsch@e1793c9e-67f9-0310-80fc-b846ff1f7b36>
Thu, 17 Apr 2008 19:38:50 +0000 (19:38 +0000)
committerraetsch <raetsch@e1793c9e-67f9-0310-80fc-b846ff1f7b36>
Thu, 17 Apr 2008 19:38:50 +0000 (19:38 +0000)
qpalma/computeSpliceAlignWithQuality.py
scripts/PipelineHeuristic.py

index 939a5bb..4734ac4 100644 (file)
@@ -184,7 +184,6 @@ def  computeSpliceAlignScoreWithQuality(dna, exons, est, original_est, quality,
    accSP       = run['numAccSuppPoints']
    mmatrixSP   = run['matchmatrixRows']*run['matchmatrixCols']
    numq        = run['numQualSuppPoints']*5*6
-   print numq
 
    lengthP = weightvector[0:lengthSP]
    donP = weightvector[lengthSP:lengthSP+donSP]
@@ -242,7 +241,6 @@ def  computeSpliceAlignScoreWithQuality(dna, exons, est, original_est, quality,
                    # because we count from 0 in python
                    Lower -= 1
                    Upper = Lower+1 ; # x-werte bleiben fest
-                   #print value,Lower,Upper
                    weightup  = 1.0*(value - cur_plf.limits[Lower]) / (cur_plf.limits[Upper] - cur_plf.limits[Lower])
                    weightlow = 1.0*(cur_plf.limits[Upper] - value) / (cur_plf.limits[Upper] - cur_plf.limits[Lower])
                    
index 98988a9..9d9a6ff 100644 (file)
@@ -437,7 +437,7 @@ class PipelineHeuristic:
          quality, qualityPlifs, run, self.currentPhi)
       
       if score!=0.0:
-          print 'diff %f' % ((trueWeight.T * self.currentPhi)[0,0] - score)
+          print 'diff %f,%f,%f' % ((trueWeight.T * self.currentPhi)[0,0] - score,(trueWeight.T * self.currentPhi)[0,0], score)
 
       # Calculate w'phi(x,y) the total score of the alignment
       return (trueWeight.T * self.currentPhi)[0,0]