b*****e
发帖数: 223 | 1 我的意思是说，在 sas 里面，我会在一个 variable 里面找出一共有多少个 distinct
group (假设 N 吧），然后对每个 group 去画图，用一个 do i=1 to N 的循环去画。
用 R 怎么做? R 用的不多， so far 只画过单张单张的图，没有这样批量画过。但是.
.....想见将来恐怕有要用的时候。虽然 sas 会做，但有的图 sas 画起来麻烦，想学
会用 R 来画。
R 里面的循环没怎么用过，网上的例子似乎也不多。来求点 hint 啦。谢过。 | j******4
发帖数: 6090 | 2 There is a "par" function that can set up the forms of your output, you can
type in ?par in R to look up those parameters.
Here is an MA plot example(Bioinformatics/Microarray):
par(mfrow=c(8,8), mar=c(0.1, 0.1, 0.1, 0.1), xaxt="n", yaxt="n")
for(i in 1:length(rma)){for(j in 1:length(rma)){
Ai<-((rma[,i])+(rma[,j]))/2
Mj<-((rma[,i])-(rma[,j]))
plot( c(min(Ai), max(Ai)), c(-3,3), type="n", xlab="A", ylab="M")
points(Ai, Mj, pch=".")
lines(c(0, 13), c(0,0), col=2)
}} | g********r
发帖数: 8017 | 3 同一张上就用par分区。不同张可以x11开新图，或者jpeg，tiff之类直接开新图并存在
硬盘上。
distinct
是.
【在 b*****e 的大作中提到】
: 我的意思是说，在 sas 里面，我会在一个 variable 里面找出一共有多少个 distinct
: group (假设 N 吧），然后对每个 group 去画图，用一个 do i=1 to N 的循环去画。
: 用 R 怎么做? R 用的不多， so far 只画过单张单张的图，没有这样批量画过。但是.
: .....想见将来恐怕有要用的时候。虽然 sas 会做，但有的图 sas 画起来麻烦，想学
: 会用 R 来画。
: R 里面的循环没怎么用过，网上的例子似乎也不多。来求点 hint 啦。谢过。
| b*****n
发帖数: 685 | 4 补充一下，开新图的话，windows是windows()，mac是quartz() | g********r
发帖数: 8017 | 5 x11()通用。
【在 b*****n 的大作中提到】
: 补充一下，开新图的话，windows是windows()，mac是quartz()
| b*****e
发帖数: 223 | 6 我的意思是说，做130张类似的图，怎么像 sas 那样让它自动做。大概会用 function,
do loop 之类吧，我都没用过，所以上来想看看你们有没有例子好学习学习。R 里面
，一写长一点绕一点的我就手软了。
回头我可以贴个 sas 里面的 codes。 | b*****e
发帖数: 223 | 7 哦哦哦，仔细看了下，就是差不多这个啦。谢谢
can
【在 j******4 的大作中提到】
: There is a "par" function that can set up the forms of your output, you can
: type in ?par in R to look up those parameters.
: Here is an MA plot example(Bioinformatics/Microarray):
: par(mfrow=c(8,8), mar=c(0.1, 0.1, 0.1, 0.1), xaxt="n", yaxt="n")
: for(i in 1:length(rma)){for(j in 1:length(rma)){
: Ai<-((rma[,i])+(rma[,j]))/2
: Mj<-((rma[,i])-(rma[,j]))
: plot( c(min(Ai), max(Ai)), c(-3,3), type="n", xlab="A", ylab="M")
: points(Ai, Mj, pch=".")
: lines(c(0, 13), c(0,0), col=2)
| t**i
发帖数: 688 | |
|
|
