Mercurial > repos > artbio > repenrich

comparison edgeR_repenrich.R @ 8:5dd4791c7b70 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/repenrich commit ccf4b178d13aeedf7214f2400e7411bcfa3e73b4
author artbio
date Tue, 11 Dec 2018 12:43:07 -0500
parents 308591fd36e5
children 6bba3e33c2e7
comparison
equal deleted inserted replaced
7:308591fd36e5 8:5dd4791c7b70
201 abline(v=0) 201 abline(v=0)
202 # Plot the repeat types 202 # Plot the repeat types
203 types <- with(results, reorder(type, -logFC, median)) 203 types <- with(results, reorder(type, -logFC, median))
204 boxplot(logFC ~ types, data=results, outline=FALSE, horizontal=TRUE, 204 boxplot(logFC ~ types, data=results, outline=FALSE, horizontal=TRUE,
205 las=2, xlab="log2(Fold Change)", main=paste0(allcontrasts, ", by Type"), yaxt="n") 205 las=2, xlab="log2(Fold Change)", main=paste0(allcontrasts, ", by Type"), yaxt="n")
206 axis(2, cex.axis=(1*11)/(length(levels(types))), 206 axis(2, cex.axis=(1*28)/(length(levels(types))),
207 at=seq(from=1, to=length(levels(types))), 207 at=seq(from=1, to=length(levels(types))),
208 labels=levels(types), las=2) 208 labels=levels(types), las=2)
209 abline(v=0) 209 abline(v=0)
210 # volcano plot 210 # volcano plot
211 TEdata = cbind(rownames(results), as.data.frame(results), score=-log(results$FDR, 10)) 211 TEdata = cbind(rownames(results), as.data.frame(results), score=-log(results$FDR, 10))