[BioC] MAPlot

Mohammad Esad-Djou shahrgol at web.de
Tue Mar 29 17:13:34 CEST 2005


Hello, 

I wrote a small program for MAS5.0:


library(affy)
library(affydata)

data.raw <- ReadAffy(filenames="./R/ME_cel/Expt1_R1.CEL",
			     "./R/ME_cel/Expt2_R1.CEL")

data.exprs <- mas5(data.raw, sc=150)
rm(data.raw)
data.out <- cbind(483, "ME-mas5", "ME_Affy MAS5 Transformation (R)", data.exprs at exprs, 0)
write.table(data.out, file="./R/ME_data/aggregation_1.dat", sep="\t", col.names=NA, quote=FALSE)
rm(data.exprs, data.out)


The result is written in aggregation_1.dat.  Now, I would like to visualize result file with MAPlot, but I don't know, which intermediate steps are necessary.  What can I do?

Thanks,
Mohammad Esad-Djou



More information about the Bioconductor mailing list