diff test-data/out.rscript @ 10:99ace6c1ff57 draft default tip

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/volcanoplot commit 086dfdddb6fbf8acb98fdcc3654b0afee6eff442-dirty
author iuc
date Wed, 11 Dec 2024 20:36:56 +0000
parents ab01e379d29e
children
line wrap: on
line diff
--- a/test-data/out.rscript	Mon Dec 09 16:04:24 2024 +0000
+++ b/test-data/out.rscript	Wed Dec 11 20:36:56 2024 +0000
@@ -20,20 +20,16 @@
 
 # Import data  ------------------------------------------------------------
 
-results <- read.delim('/tmp/tmpl4o1f_bf/files/5/2/5/dataset_52538741-d085-42da-817b-263bd4f7cf98.dat', header = TRUE)
+results <- read.delim('/private/var/folders/zn/m_qvr9zd7tq0wdtsbq255f8xypj_zg/T/tmprh4qip75/files/d/2/2/dataset_d2255b46-f0f6-4900-8b9e-bd352e34f303.dat', header = TRUE)
 
 
 # Format data  ------------------------------------------------------------
 
-# Create columns from the column numbers specified and use the existing category_symbol column for shapes
-results <- results %>% mutate(
-    fdr = .[[4]],
-    pvalue = .[[3]],
-    logfc = .[[2]],
-    labels = .[[1]],
-)
-
-# Check if shape_col is provided 
+# Create columns from the column numbers specified
+results <- results %>% mutate(fdr = .[[4]],
+                              pvalue = .[[3]],
+                              logfc = .[[2]],
+                              labels = .[[1]])
 
 # Get names for legend
 down <- unlist(strsplit('Down,Not Sig,Up', split = ","))[1]
@@ -53,7 +49,7 @@
 # Specify genes to label --------------------------------------------------
 
 # Import file with genes of interest
-labelfile <- read.delim('/tmp/tmpl4o1f_bf/files/5/d/4/dataset_5d401b02-f6af-4ed9-b853-992fd4a4d044.dat', header = TRUE)
+labelfile <- read.delim('/private/var/folders/zn/m_qvr9zd7tq0wdtsbq255f8xypj_zg/T/tmprh4qip75/files/5/e/5/dataset_5e5b8fb0-bf65-438e-9b5b-03a540d9aa5d.dat', header = TRUE)
 
 # Label the genes of interest in results table
 results <- mutate(results, labels = ifelse(labels %in% labelfile[, 1], labels, ""))
@@ -65,17 +61,15 @@
 # Open file to save plot as PDF
 pdf("volcano_plot.pdf")
 
-# Set up base plot with faceting by category_symbol instead of shapes
+# Set up base plot
 p <- ggplot(data = results, aes(x = logfc, y = -log10(pvalue))) +
+    geom_point(aes(colour = sig)) +
     scale_color_manual(values = colours) +
     theme(panel.grid.major = element_blank(),
-          panel.grid.minor = element_blank(),
-          panel.background = element_blank(),
-          axis.line = element_line(colour = "black"),
-          legend.key = element_blank()) 
-
-# Conditional logic to use either shape or facet based on user selection
-p <- p + geom_point(aes(colour = sig)) #only add color
+        panel.grid.minor = element_blank(),
+        panel.background = element_blank(),
+        axis.line = element_line(colour = "black"),
+        legend.key = element_blank())
 
 # Add gene labels
 p <- p + geom_text_repel(data = filter(results, labels != ""), aes(label = labels),