figures/stats: print stats (for paper) when ceating prodigal fig

workflow_figures/scripts/fig_prodigal.R

@@ -35,6 +35,21 @@ for(sname in names(snakemake@config$samples)){
     TABS[[sname]] <- read_prodigal_gcounts(snakemake@input[[sname]])
 }
 
+############################## STATS
+
+print("STATS: Number of total proteins: max vs min:")
+for(sname in names(snakemake@config$samples)){
+    print(sprintf(
+        "STATS: Sample %s: %s vs %s = %d / %d = %.2f",
+        sname,
+        TABS[[sname]]$tool[which.max(TABS[[sname]]$total)],
+        TABS[[sname]]$tool[which.min(TABS[[sname]]$total)],
+        max(TABS[[sname]]$total),
+        min(TABS[[sname]]$total),
+        max(TABS[[sname]]$total) / min(TABS[[sname]]$total)
+    ))
+}
+
 ############################## PLOT
 
 # individual plots