Mercurial > repos > yufei-luo > s_mart

annotate SMART/DiffExpAnal/DESeqTools/plotDispEstimates.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 # plotDispEstimates
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
2 # scatter plots representing dispersion estimates vs mean expression
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
3
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
4 # input : cds, OUT_plotDispEstimatesName
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
5 # output : scatterplot (png)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
6
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
7 plotDispEstimates <- function( cds, OUT_plotDispEstimatesName, out = TRUE ){
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
8
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
9 if (out) png( file=OUT_plotDispEstimatesName )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
10
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
11 plot(
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
12 rowMeans( counts(cds, normalized=T) ),
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
13 fitInfo(cds)$perGeneDispEsts,
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
14 pch=".", log="xy",
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
15 xlab = "Mean expression strength", ylab = "Dispersion estimate" )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
16
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
17 xg <- 10^seq(-.5, 5, length.out=300)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
18 lines( xg, fitInfo(cds)$dispFun(xg), col="red" )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
19
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
20 if (out) dev.off()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
21 }