figures/stats: print stats (for paper) when creating barrnap fig

workflow_figures/scripts/fig_barrnap.R

@@ -35,6 +35,24 @@ for(sname in names(snakemake@config$samples)){
     TABS[[sname]] <- proc_barrnap(read_barrnap(snakemake@input[[sname]])) # total counts per tool/kingdom/gene category
 }
 
+############################## STATS
+
+print("STATS: Number of complete genes: max vs min:")
+for(sname in names(snakemake@config$samples)){
+    for(kingdom_name in BARRNAP_KINGDOM_NAMES){
+        tmp_counts <- TABS[[sname]] %>% filter(partial == "complete", kingdom == kingdom_name)
+        print(tmp_counts)
+        print(sprintf(
+            "STATS: Sample %s: %s, %d / %d = %.2f",
+            sname,
+            kingdom_name,
+            max(tmp_counts$x),
+            min(tmp_counts$x),
+            max(tmp_counts$x) / min(tmp_counts$x)
+        ))
+    }
+}
+
 ############################## PLOT
 
 # individual plots