source: BioInfo/ReadRMA.R @ 59

Last change on this file since 59 was 23, checked in by jazz, 17 years ago

Sample R Code from YM. Prof. Chung

File size: 548 bytes
Line 
1######## read RMA  #################
2
3setwd("D:\\Women cancer\\ovary")
4
5rm(list=ls())
6
7memory.limit(size=2000)
8
9memory.size()
10
11library(affy)
12
13fileNames<-dir(".", "[.]CEL")[1:15]
14
15LatinData<-ReadAffy(filenames=fileNames)  ###  Read CEL files
16
17str(LatinData)
18
19k<-LatinData@exprs
20
21LatinData.rma<-expresso(LatinData, bgcorrect.method="rma", normalize.method="quantiles",  pmcorrect.method="pmonly", summary.method="medianpolish")
22## 完成到此即成功 ##
23
24
25write.table(exprs(LatinData.rma), "LatinData_rma.txt", sep="\t")
26
27## total about 15 .CELs
28
Note: See TracBrowser for help on using the repository browser.