[BioC] Affymetrix preprocessing, MAS5 vs. RMA

Dozmorov, Mikhail G. (HSC) Mikhail-Dozmorov at ouhsc.edu
Thu Jan 22 21:47:58 CET 2009


Hello,

I'm learning basics of importing Affy's .cel files to R, and do it with affy library:
	R> Data <- ReadAffy()
When I do preprocessing
	R> eset <- rma(Data) OR eset <- mas5(Data)
I'm getting very different expression values.
	From MAS5:
ID	21.CEL	33.CEL	36.CEL
1415670_at	2494.54	2576.29	2240.92
1415671_at	6292.17	5013.13	3075.89
1415672_at	6727.66	5858.31	5690.70
1415673_at	567.38	718.02	875.25
	From RMA, Antilogged from Log2
ID	21.CEL	33.CEL	36.CEL
1415670_at	74.16	69.30	59.08
1415671_at	191.52	171.38	94.30
1415672_at	194.18	199.57	175.13
1415673_at	14.93	13.16	21.50

Can't find any explanation why it could happen. And what to use? Please, advice! Thank you. Mikhail.



More information about the Bioconductor mailing list