Fix a bug when there is no result.

fcyu · Oct 4, 2017 · d368f5e · d368f5e
1 parent 3325076
commit d368f5e
Showing 1 changed file with 13 additions and 6 deletions.
diff --git a/src/main/java/proteomics/ECL2.java b/src/main/java/proteomics/ECL2.java
@@ -194,12 +194,19 @@ public static void main(String[] args) {
             // save result
             logger.info("Estimating q-value...");
             List<List<FinalResultEntry>> picked_result = pickResult(final_search_results);
-            CalFDR cal_fdr_obj = new CalFDR(picked_result.get(0), cal_evalue);
-            List<FinalResultEntry> intra_result = cal_fdr_obj.includeStats(cal_evalue);
-            intra_result.sort(Collections.reverseOrder());
-            cal_fdr_obj = new CalFDR(picked_result.get(1), cal_evalue);
-            List<FinalResultEntry> inter_result = cal_fdr_obj.includeStats(cal_evalue);
-            inter_result.sort(Collections.reverseOrder());
+            List<FinalResultEntry> intra_result = new LinkedList<>();
+            List<FinalResultEntry> inter_result = new LinkedList<>();
+            if (!picked_result.get(0).isEmpty()) {
+                CalFDR cal_fdr_obj = new CalFDR(picked_result.get(0), cal_evalue);
+                intra_result = cal_fdr_obj.includeStats(cal_evalue);
+                intra_result.sort(Collections.reverseOrder());
+            }
+            if (!picked_result.get(1).isEmpty()) {
+                CalFDR cal_fdr_obj = new CalFDR(picked_result.get(1), cal_evalue);
+                inter_result = cal_fdr_obj.includeStats(cal_evalue);
+                inter_result.sort(Collections.reverseOrder());
+            }
+
             logger.info("Saving results...");
             saveTargetResult(intra_result, build_index_obj.getProAnnotateMap(), spectra_path, true, cal_evalue);
             saveTargetResult(inter_result, build_index_obj.getProAnnotateMap(), spectra_path, false, cal_evalue);