
setwd("C:/statgen")getwd()
library(affy)
library(help=affy)
files<-dir(pattern=".cel")
AffyData<-read.affybatch(filenames=files)
exprsSet.RMA<-rma(AffyData)
exp.RMA<-exprs(exprsSet.RMA)
exprsSet.MAS5<-mas5(AffyData)
write.exprs(exprsSet.MAS5, file="mydata.txt")
这个是别人叫我写的语句,但还是做不出来,不知道是不是cel文件有问题?有问题的话要怎么下载,大四学生在做毕业专题,第一次用R,很多不懂,希望热心的大神能教我一下。谢谢谢谢~