Last change
on this file since 186 was
23,
checked in by jazz, 16 years ago
|
Sample R Code from YM. Prof. Chung
|
File size:
548 bytes
|
Line | |
---|
1 | ######## read RMA ################# |
---|
2 | |
---|
3 | setwd("D:\\Women cancer\\ovary") |
---|
4 | |
---|
5 | rm(list=ls()) |
---|
6 | |
---|
7 | memory.limit(size=2000) |
---|
8 | |
---|
9 | memory.size() |
---|
10 | |
---|
11 | library(affy) |
---|
12 | |
---|
13 | fileNames<-dir(".", "[.]CEL")[1:15] |
---|
14 | |
---|
15 | LatinData<-ReadAffy(filenames=fileNames) ### Read CEL files |
---|
16 | |
---|
17 | str(LatinData) |
---|
18 | |
---|
19 | k<-LatinData@exprs |
---|
20 | |
---|
21 | LatinData.rma<-expresso(LatinData, bgcorrect.method="rma", normalize.method="quantiles", pmcorrect.method="pmonly", summary.method="medianpolish") |
---|
22 | ## 完成到此即成功 ## |
---|
23 | |
---|
24 | |
---|
25 | write.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.