diff edgeR_Differential_Gene_Expression.xml @ 109:dfebd3787276 draft

planemo upload for repository https://bitbucket.org/EMCbioinf/galaxy-tool-shed-tools/raw/master/edger_with_design_matrix commit 2700e500a4fb135a20ede7d52221a9d31f1aaa5e-dirty

--- a/edgeR_Differential_Gene_Expression.xml	Tue Sep 01 04:45:24 2015 -0400
+++ b/edgeR_Differential_Gene_Expression.xml	Tue Sep 01 04:53:59 2015 -0400
@@ -236,7 +236,7 @@
     
     
     diff_x <- abs(max(points\$x)-min(points\$x))
-    diff_y &lt;-(max(points\$y)-min(points\$y))
+    diff_y <-(max(points\$y)-min(points\$y))
     plot(c(min(points\$x),max(points\$x) + 0.45 * diff_x), c(min(points\$y) - 0.05 * diff_y,max(points\$y) + 0.05 * diff_y), main="edgeR logFC-MDS Plot on top 500 genes",type="n", xlab="Leading logFC dim 1", ylab="Leading logFC dim 2")
     points(points\$x,points\$y,pch=20)
     text(points\$x, points\$y,rownames(dge\$samples),cex=1.25,col="gray",pos=4)
@@ -329,7 +329,7 @@
       }
       
       with(etable, plot(logCPM, logFC, pch=20, main="edgeR: Fold change vs abundance"))
-      with(subset(etable, FDR &lt; fdr), points(logCPM, logFC, pch=20, col="red"))
+      with(subset(etable, FDR < fdr), points(logCPM, logFC, pch=20, col="red"))
       abline(h=c(-1,1), col="blue")
       dev.off()
     }
@@ -348,8 +348,8 @@
         bitmap(output_PValue_distribution_plot,type="png16m",width=14,height=14)
       }
       
-      expressed_genes <- subset(etable, PValue &lt; 0.99)
-      h <- hist(expressed_genes\$PValue,breaks=nrow(expressed_genes)/15,main="Binned P-Values (&lt; 0.99)")
+      expressed_genes <- subset(etable, PValue < 0.99)
+      h <- hist(expressed_genes\$PValue,breaks=nrow(expressed_genes)/15,main="Binned P-Values (< 0.99)")
       center <- sum(h\$counts) / length(h\$counts)
       lines(c(0,1),c(center,center),lty=2,col="red",lwd=2)
       k <- ksmooth(h\$mid, h\$counts)