Mercurial > repos > yufei-luo > differential_expression_analysis_pipeline_for_rnaseq_data

view DiffExpAnal/DESeqTools/plotDispEstimates.R @ 0:63799b789162 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author yufei-luo
date Tue, 22 Jan 2013 10:07:03 -0500
parents
children
line wrap: on
line source

# Marie-Anges Dillies
# plotDispEstimates
# scatter plots representing dispersion estimates vs mean expression

# input : cds, OUT_plotDispEstimatesName
# output : scatterplot (png)

plotDispEstimates <- function( cds, OUT_plotDispEstimatesName, out = TRUE ){

  if (out) png( file=OUT_plotDispEstimatesName )
  
  plot(
  	rowMeans( counts(cds, normalized=T) ),
  	fitInfo(cds)$perGeneDispEsts,
  	pch=".", log="xy",
  	xlab = "Mean expression strength", ylab = "Dispersion estimate" )
  	
  xg <- 10^seq(-.5, 5, length.out=300)
  lines( xg, fitInfo(cds)$dispFun(xg), col="red" )

  if (out) dev.off()
}