|
Last change
on this file since 158 was
23,
checked in by jazz, 18 years ago
|
Sample R Code from YM. Prof. Chung
|
|
File size:
485 bytes
|
| Line | |
|---|
| 1 | ######### MAS5 ################### |
|---|
| 2 | |
|---|
| 3 | |
|---|
| 4 | setwd("D:\\Women cancer\\ovary") |
|---|
| 5 | |
|---|
| 6 | rm(list=ls()) |
|---|
| 7 | |
|---|
| 8 | memory.limit(size=2000) |
|---|
| 9 | |
|---|
| 10 | memory.size() |
|---|
| 11 | |
|---|
| 12 | library(affy) |
|---|
| 13 | |
|---|
| 14 | fileNames<-dir(".", "[.]CEL")[1:10] |
|---|
| 15 | |
|---|
| 16 | LatinData<-ReadAffy(filenames=fileNames) ### Read CEL files |
|---|
| 17 | |
|---|
| 18 | str(LatinData) |
|---|
| 19 | |
|---|
| 20 | k<-LatinData@exprs |
|---|
| 21 | |
|---|
| 22 | LatinData.mas5<-expresso(LatinData, pmcorrect.method="mas", bgcorrect.method="mas", summary.method="mas") |
|---|
| 23 | |
|---|
| 24 | LatinData.mas5.norm<-affy.scalevalue.exprSet(LatinData.mas5, sc=500) |
|---|
| 25 | ## 完成到此即成功 ## |
|---|
| 26 | |
|---|
Note: See
TracBrowser
for help on using the repository browser.
