R：fRMAを行う

frozen RMA (fRMA)は正規化方法の一つで、サンプル数が少ない場合にRMAよりも優れるとの報告がある。

この計算も、Rのパッケージを用いて行うことができる。

以下は、遺伝子発現情報を抽出するまでのコードである。frmaのパッケージを読みこんでfrma関数を用いれば完了する。

1 2 3 4 5 6

library(affy) library(frma) cel_dat <- ReadAffy() eset <- frma(cel_dat) emat <- exprs(eset)