added R script for exercise 3

ahryhorzhevska · Apr 11, 2022 · 031f10a · 031f10a
1 parent dad4a63
commit 031f10a
Showing 1 changed file with 15 additions and 0 deletions.
diff --git a/plot.R b/plot.R
@@ -0,0 +1,15 @@
+library(ggplot2)
+
+gene <- read.csv2(snakemake@input[[1]], dec = ",")
+colnames(gene) <- "id"
+gene$id <- as.numeric(gene$id)
+
+print(paste0("I'm plotting the file ", snakemake@wildcards[["file"]]))
+
+pdf(file = snakemake@output[[1]])
+
+ggplot(gene, aes(x = id)) +
+    geom_density(alpha = .2, fill = "#FF6666") +
+    ggtitle(paste0("Plot "), snakemake@wildcards[["file"]])
+
+dev.off()