######## read RMA  #################

setwd("D:\\Women cancer\\ovary")

rm(list=ls())

memory.limit(size=2000)

memory.size()

library(affy)

fileNames<-dir(".", "[.]CEL")[1:15]

LatinData<-ReadAffy(filenames=fileNames)  ###  Read CEL files

str(LatinData)

k<-LatinData@exprs 

LatinData.rma<-expresso(LatinData, bgcorrect.method="rma", normalize.method="quantiles",  pmcorrect.method="pmonly", summary.method="medianpolish")
## 完成到此即成功 ##


write.table(exprs(LatinData.rma), "LatinData_rma.txt", sep="\t")

## total about 15 .CELs