|
Last change
on this file since 153 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 | |
|---|
| 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.
