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
#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)