Mercurial > repos > yufei-luo > s_mart

annotate SMART/DiffExpAnal/DESeqTools/boxplotCounts.R @ 33:aa0420172fc6

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Deleted selected files
author m-zytnicki
date Tue, 30 Apr 2013 14:34:53 -0400
parents 94ab73e8a190
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
18
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
1 # boxplotCounts
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
2 # boxplots representing counts distribution per sample
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
3
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
4 # input : counts, target, projectName, type of data (raw or norm)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
5 # output : boxplot (png)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
6
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
7 # created Feb 7th, 2012
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
8 # modified April 30th, 2012
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
9
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
10 boxplotCounts <- function( counts, group, type = c("raw", "norm"), OUT_boxplotCountsName, out = TRUE ){
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
11
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
12 if (out) png( file=OUT_boxplotCountsName )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
13
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
14 boxplot( log2(counts+1), col=as.integer(group)+1, main = paste(type[1], " counts distribution", sep="" ) )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
15 legend( "topright", as.character(unique(group)), lty=1, col=as.integer(unique(group))+1 )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
16
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
17 if (out) dev.off()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
18 }