[BioC] edgeR: plotSmear for multiple detags

Mark Robinson mark.robinson at imls.uzh.ch
Fri Oct 21 21:12:55 CEST 2011


Hi Steve,

The standard plot parameters are passed through.  But, to do want you want, you'll have to construct it by hand.  Here is a worked example:

# concocted example
y <- matrix(rnbinom(10000,mu=5,size=2),ncol=4)
y[1:10,1:2] <- round(y[1:10,1:2]*4)
d <- DGEList(counts=y, group=rep(1:2,each=2), lib.size=colSums(y))
d <- estimateCommonDisp(d)
de <- exactTest(d)

# construct colors, expansion and plot character manually
cols <- rep("grey",nrow(d))
conc <- d$conc$conc.group
k <- conc[,1]==min(conc[,1]) | conc[,2]==min(conc[,2])
cols[k] <- "orange"
cols[de$table$p.value<.05 & de$table$logFC < 0] <- "blue"
cols[de$table$p.value<.05 & de$table$logFC > 0] <- "red"
cexs <- rep(.2, nrow(d))
cexs[k] <- .8
cexs[1:10] <- 2

plotSmear(de, col=cols, cex=cexs)

Of course, don't take my example as a reasonable way to find DE genes (it's not, no mult. test correction), or as a reasonable way to simulate DE (also not!).

Hope that helps,
Mark


On 20.10.2011, at 19:08, steve Shen wrote:

> Hi,
> 
> Is there anyway to plot multiple detags in one figure with plotSmear labeled
> with either different color or shape? Thanks,
> 
> Steve
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list