[Sashimi-commit] SF.net SVN: sashimi:[5968] trunk/trans_proteomic_pipeline/src/Validation/ DiscriminateFunction/Sequest/SequestDiscrimFunction.cxx

[<ds...@us...>]

Revision: 5968
          http://sashimi.svn.sourceforge.net/sashimi/?rev=5968&view=rev
Author:   dshteyn
Date:     2012-10-31 17:25:29 +0000 (Wed, 31 Oct 2012)
Log Message:
-----------
No expect score for sequest

Modified Paths:
--------------
    trunk/trans_proteomic_pipeline/src/Validation/DiscriminateFunction/Sequest/SequestDiscrimFunction.cxx

Modified: trunk/trans_proteomic_pipeline/src/Validation/DiscriminateFunction/Sequest/SequestDiscrimFunction.cxx
===================================================================
--- trunk/trans_proteomic_pipeline/src/Validation/DiscriminateFunction/Sequest/SequestDiscrimFunction.cxx	2012-10-31 17:20:55 UTC (rev 5967)
+++ trunk/trans_proteomic_pipeline/src/Validation/DiscriminateFunction/Sequest/SequestDiscrimFunction.cxx	2012-10-31 17:25:29 UTC (rev 5968)
@@ -64,14 +64,14 @@
   }
   SequestResult* seqresult = (SequestResult*)(result);
   if (use_expect_) {
-    if (seqresult->expect_ > 0) {
-      return -5-log((double)seqresult->expect_);
-    }
-    else {
-      cerr << "ERROR: Sequest Result for spectrum " << seqresult->spectrum_ 
-	   << " does NOT have a positive expect score." << endl;
-      return -100;
-    }
+    //    if (seqresult->expect_ > 0) {
+    //  return -5-log((double)seqresult->expect_);
+    //}
+    //else {
+    //  cerr << "ERROR: Sequest Result for spectrum " << seqresult->spectrum_ 
+    //	   << " does NOT have a positive expect score." << endl;
+    //  return -100;
+    //}
   }
   double tot = const_;
   tot += xcorr_p_wt_ * getXcorrP(seqresult->xcorr_, getPepLen(seqresult->peptide_));

This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.