Merge pull request #3 from LiChenPU/dev

Dev

Author: LiChenPU

FDR_example/Target-Decoy Generator.R

@@ -209,7 +209,7 @@ FDR_Calculator <- function(decoy_address, output_address,
     anns_elements <- sapply(anns, function(x) unlist(strsplit(x, "\\d+")))
     decoys <- sapply(anns_elements, function(x) sum(unlist(x) %in% imp_elements))
     FP <- length(which(decoys > 0)) * (f - 1)
-    FDR <- FP / (length(anns) - FP)
+    FDR <- FP / (length(anns))
 }
 
 NetID <- function(work_dir, file_name, ion_mode, decoy_lib) {

FDR_example/plot.R

@@ -51,7 +51,7 @@ load("FDR.RData")
   ## Plot ####
   {
     p_all_precision = ggplot(data = gt_all, mapping = aes(x = method, y = precision, color = method)) + 
-      geom_boxplot() + #geom_point(data = stan, mapping = aes(x = method, y = precision, fill = method), shape = 22, color = "black") +
+      geom_boxplot(outlier.shape = NA) + #geom_point(data = stan, mapping = aes(x = method, y = precision, fill = method), shape = 22, color = "black") +
       geom_point(alpha = 1 / 5, position = "jitter") +
       # ggtitle("Precision by Ground Truth") + 
       theme_classic(base_size = 12 # edit font size for all non-data text
@@ -76,7 +76,7 @@ load("FDR.RData")
       facet_wrap(~ library, nrow = 2) #+coord_flip()
 
     p1_precision = ggplot(data = gt_all %>% filter(library=="HMDB"), mapping = aes(x = method, y = precision, color = method)) + 
-      geom_boxplot() + #geom_point(data = stan, mapping = aes(x = method, y = precision, fill = method), shape = 22, color = "black") +
+      geom_boxplot(outlier.shape = NA) + #geom_point(data = stan, mapping = aes(x = method, y = precision, fill = method), shape = 22, color = "black") +
       geom_point(alpha = 1 / 5, position = "jitter") +
       # ggtitle("Precision by Ground Truth") +
       # guides(color = F) +
@@ -147,7 +147,7 @@ load("FDR.RData")
   ## Plot ####
   {
     p_all_FDR = ggplot(data = fdr_all, mapping = aes(x = method, y = FDR, color = method)) + 
-      geom_boxplot() + #geom_point(data = stan, mapping = aes(x = method, y = precision, fill = method), shape = 22, color = "black") +
+      geom_boxplot(outlier.shape = NA) + #geom_point(data = stan, mapping = aes(x = method, y = precision, fill = method), shape = 22, color = "black") +
       geom_point(alpha = 1 / 5, position = "jitter") +
       # ggtitle("Precision by Ground Truth") + 
       theme_classic(base_size = 12 # edit font size for all non-data text
@@ -172,7 +172,7 @@ load("FDR.RData")
       facet_wrap(~ library, nrow = 2) #+coord_flip()
 
     p1_FDR = ggplot(data = fdr_all %>% filter(library=="HMDB"), mapping = aes(x = method, y = FDR, color = method)) + 
-      geom_boxplot() + #geom_point(data = stan, mapping = aes(x = method, y = precision, fill = method), shape = 22, color = "black") +
+      geom_boxplot(outlier.shape = NA) + #geom_point(data = stan, mapping = aes(x = method, y = precision, fill = method), shape = 22, color = "black") +
       geom_point(alpha = 1 / 5, position = "jitter") +
       # ggtitle("Precision by Ground Truth") + 
       guides(color = guide_legend(