Mercurial > repos > yufei-luo > s_mart

view SMART/DiffExpAnal/DESeqTools/MAplotDE.R @ 18:94ab73e8a190

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author m-zytnicki
date Mon, 29 Apr 2013 03:20:15 -0400
parents
children
line wrap: on
line source

# MAplotDE
# MAplot of DE genes

# input : res, alpha,OUT_MAplotDEName 
# output : MAplot (png)

MAplotDE <- function( res, alpha, OUT_MAplotDEName, out = TRUE ){

 	if (out) png( file=OUT_MAplotDEName )
  
 	plot( res$baseMean, res$log2FoldChange, pch=".", xlab="Mean expression", ylab="log2FC", main="",
  		log="x", col=ifelse(res$padj < alpha, "red", "black") )
	abline(h=0, col="red")

  	if (out) dev.off()
}