Mercurial > repos > yufei-luo > differential_expression_analysis_pipeline_for_rnaseq_data
view DiffExpAnal/DESeqTools/plotDispEstimates.R @ 10:6e573fd3c41b draft
Uploaded
author | yufei-luo |
---|---|
date | Mon, 13 May 2013 10:06:30 -0400 |
parents | 63799b789162 |
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() }