Mercurial > repos > yufei-luo > s_mart

annotate SMART/DiffExpAnal/DESeqTools/densityPlot.R @ 31:0ab839023fe4

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author m-zytnicki
date Tue, 30 Apr 2013 14:33:21 -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 # densityPlot
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
2 # density plot of all samples
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
3
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
4 # input : counts, target, projectName
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
5 # output : densplot (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
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
11 densityPlot <- function( counts, group, OUT_densityPlotName, out = TRUE ){
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
12
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
13 if (out) png( file=OUT_densityPlotName )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
14
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
15 couleurs <- as.integer( group ) + 1
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
16 ylim <- c(0, max(density(log2(counts)+1)$y)*1.5)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
17 plot( density(log2(counts[,1])+1), main="Density of counts distribution", col=couleurs[1], ylim = ylim )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
18 for (i in 2:ncol(counts))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
19 lines( density(log2(counts[,i])+1), col=couleurs[i] )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
20 legend( "topright", as.character(unique(group)), lty=1, col=as.integer(unique(group))+1 )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
21
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
22 if (out) dev.off()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
23 }