MAS5で正規化

メモ。こちらを非常に参考にさせて頂く。

CELファイルのあるディレクトリに移動してから

library(affy) 
data <- ReadAffy() 
eset_mas <- mas5(data) 
write.exprs(eset_mas, file="data_mas.txt")

で結果を「data_mas.txt」に出力する。

RMAよりも時間は倍以上かかるし、メモリは食うし（下手するとエラーが出て停止）、結構骨が折れる作業だと実感。