Mercurial > repos > yufei-luo > differential_expression_analysis_pipeline_for_rnaseq_data

view DiffExpAnal/DESeqTools/barplotTC.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

# barplotTC
# barplot representing total count per sample

# input : counts, target, projectName
# output : barplotTC (png)

# created Feb 7th, 2012
# modified April 30th, 2012 (group instead of target$group)
# Marie-Anges Dillies

barplotTC <- function( counts, group, OUT_barplotTCName, out = TRUE ){

  if (out) png( file=OUT_barplotTCName )

  ylim <- c(0, max(colSums(counts))*1.2)
  barplot( colSums(counts), col=as.integer(group)+1, main = "Total Read Count per Sample",  ylim=ylim )
  legend( "topright", as.character(unique(group)), lty=1,
         col=as.integer(unique(group))+1 )

  if (out) dev.off()
}