annotate deseq/differential_expression_analysis_pipeline_for_rnaseq_data-a03838a6eb54/DiffExpAnal/DESeqTools/boxplotCounts.R @ 10:6e573fd3c41b draft

Uploaded
author yufei-luo
date Mon, 13 May 2013 10:06:30 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
10
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
1 # boxplotCounts
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
2 # boxplots representing counts distribution per sample
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
3
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
4 # input : counts, target, projectName, type of data (raw or norm)
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
5 # output : boxplot (png)
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
6
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
7 # created Feb 7th, 2012
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
8 # modified April 30th, 2012
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
9 # Marie-Anges Dillies
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
10
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
11 boxplotCounts <- function( counts, group, type = c("raw", "norm"), OUT_boxplotCountsName, out = TRUE ){
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
12
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
13 if (out) png( file=OUT_boxplotCountsName )
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
14
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
15 boxplot( log2(counts+1), col=as.integer(group)+1, main = paste(type[1], " counts distribution", sep="" ) )
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
16 legend( "topright", as.character(unique(group)), lty=1, col=as.integer(unique(group))+1 )
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
17
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
18 if (out) dev.off()
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
19 }