From 79b91237241512da18b893b6b0fb9a6b6dd06d2f Mon Sep 17 00:00:00 2001
From: Fengchao <fyuab@connect.ust.hk>
Date: Wed, 4 Jul 2018 13:35:18 +0800
Subject: [PATCH] Fall back to the linear FDR estimation approach when
 fuse_count < decoy_count. This change will affect the results a lot in some
 cases.

---
 src/main/java/proteomics/Validation/CalFDR.java | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/main/java/proteomics/Validation/CalFDR.java b/src/main/java/proteomics/Validation/CalFDR.java
index 071f7c9..ff4dc25 100644
--- a/src/main/java/proteomics/Validation/CalFDR.java
+++ b/src/main/java/proteomics/Validation/CalFDR.java
@@ -96,7 +96,7 @@ public static Map<String, Entry> calFDR(String sqlPath, boolean cal_evalue, Stri
                 fdr = 0;
             } else if (fuse_count < decoy_count) {
                 printWarning = true;
-                fdr = 0;
+                fdr = (double) (fuse_count + decoy_count) / (double) target_count;
             } else {
                 fdr = (double) (fuse_count - decoy_count) / (double) target_count;
             }