-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathSec-Plotting.R
46 lines (31 loc) · 1.3 KB
/
Sec-Plotting.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
## @knitr env, echo=FALSE, massage = FALSE
suppressPackageStartupMessages(library(Biobase))
data(sample.ExpressionSet)
## @knitr plotcode, eval = FALSE, tidy = FALSE
## plot(exprs(sample.ExpressionSet[, 1]),
## exprs(sample.ExpressionSet[, 2]),
## log = "xy",
## xlab = sampleNames(sample.ExpressionSet)[1],
## ylab = sampleNames(sample.ExpressionSet)[2])
## abline(0, 1)
## grid()
## @knitr label=plotfig1, echo=FALSE, fig.width=5, fig.height=4, tidy=FALSE, warning = FALSE
plot(exprs(sample.ExpressionSet[, 1]),
exprs(sample.ExpressionSet[, 2]),
log = "xy",
xlab = sampleNames(sample.ExpressionSet)[1],
ylab = sampleNames(sample.ExpressionSet)[2])
abline(0, 1)
grid()
## @knitr plotcode0, eval = FALSE, tidy = FALSE
## pairs(log2(exprs(sample.ExpressionSet)[, 1:4]),
## pch = 19,
## col = "#0000FF20")
## @knitr label=plotfig2,echo=FALSE,fig.width=3.5,fig.height=3.5,tidy=FALSE, warning = FALSE
pairs(log2(exprs(sample.ExpressionSet)[, 1:3]),
pch = 19,
col = "#0000FF20")
## @knitr label=plotfig3,echo=TRUE,fig.width=5,fig.height=4,tidy=FALSE
boxplot(log2(exprs(sample.ExpressionSet)))
## @knitr label=plotfig4,echo=TRUE,fig.width=5,fig.height=4,tidy=FALSE, warning = FALSE, message = FALSE
smoothScatter(log2(exprs(sample.ExpressionSet)[, 1:2]))