[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