Mercurial > repos > yufei-luo > differential_expression_analysis_pipeline_for_rnaseq_data

view deseq/differential_expression_analysis_pipeline_for_rnaseq_data-a03838a6eb54/DiffExpAnal/DESeqTools/barplotNul.R @ 10:6e573fd3c41b draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author yufei-luo
date Mon, 13 May 2013 10:06:30 -0400
parents
children
line wrap: on
line source

# barplotNul
# barplot representing null counts per sample

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

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

barplotNul <- function( counts, group,  OUT_barplotNulName, out = TRUE ){

  if (out) png( file=OUT_barplotNulName )

  N <- apply(counts, 2, function(x){sum(x == 0)})/nrow(counts)
  barplot(N, col=as.integer(group)+1, main = "Proportion of null counts per Sample", ylim = c(0,1))
  legend("topright", as.character(unique(group)), lty=1, col=as.integer(unique(group))+1)

  if (out) dev.off()
}