Repository 'limma_voom'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/limma_voom

Changeset 18:97e06a4c7c75 (2019-05-05)
Previous changeset 17:86b3df7db58b (2019-02-18) Next changeset 19:c9c6427edfe9 (2019-05-28)
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/limma_voom commit 5f5e2ad38bf004b33fc5da5d9d4e75cf1130f666
modified:
limma_voom.R
limma_voom.xml
b
diff -r 86b3df7db58b -r 97e06a4c7c75 limma_voom.R
--- a/limma_voom.R Mon Feb 18 17:49:24 2019 -0500
+++ b/limma_voom.R Sun May 05 09:02:04 2019 -0400
[
@@ -879,8 +879,15 @@
         geneanno <- y$genes
         geneanno[, 2] <- make.unique(geneanno[, 2])
 
+        # use the logCPMS for the counts
+        if (wantTrend) {
+            cnts <- logCPM
+        } else{
+            cnts <- vData$E
+        }
+
         # MD plot
-        Glimma::glMDPlot(fit, coef=i, counts=y$counts, anno=geneanno, groups=factors[, 1],
+        Glimma::glMDPlot(fit, coef=i, counts=cnts, anno=geneanno, groups=factors[, 1],
              status=status[, i], sample.cols=col.group,
              main=paste("MD Plot:", unmake.names(con)), side.main=colnames(y$genes)[2],
              folder=paste0("glimma_", unmake.names(con)), launch=FALSE)
@@ -889,10 +896,10 @@
         linkData <- rbind(linkData, c(linkName, linkAddr))
 
         # Volcano plot
-        Glimma::glXYPlot(x=fit$coefficients[, i], y=fit$lods[, i], counts=y$counts, anno=geneanno, groups=factors[, 1], 
+        Glimma::glXYPlot(x=fit$coefficients[, i], y=-log10(fit$p.value[, i]), counts=cnts, anno=geneanno, groups=factors[, 1],
             status=status[, i], sample.cols=col.group,
             main=paste("Volcano Plot:", unmake.names(con)), side.main=colnames(y$genes)[2],
-            xlab="logFC", ylab="logodds",
+            xlab="logFC", ylab="-log10(P-value)",
             folder=paste0("glimma_volcano_", unmake.names(con)), launch=FALSE)
         linkName <- paste0("Glimma_VolcanoPlot_", con, ".html")
         linkAddr <- paste0("glimma_volcano_", con, "/XY-Plot.html")
b
diff -r 86b3df7db58b -r 97e06a4c7c75 limma_voom.xml
--- a/limma_voom.xml Mon Feb 18 17:49:24 2019 -0500
+++ b/limma_voom.xml Sun May 05 09:02:04 2019 -0400
b
@@ -1,4 +1,4 @@
-<tool id="limma_voom" name="limma" version="3.38.3+galaxy1">
+<tool id="limma_voom" name="limma" version="3.38.3+galaxy2">
     <description>
         Perform differential expression with limma-voom or limma-trend
     </description>