Last change
on this file since 185 was
23,
checked in by jazz, 16 years ago
|
Sample R Code from YM. Prof. Chung
|
File size:
548 bytes
|
Rev | Line | |
---|
[23] | 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.