hi,
i have some CEL file and i did like below but i have X244901_at while i need AGI code even i used cdf
library(affy)
library(vsn)
library(limma)
library(altcdfenvs)
library(simpleaffy)
# listing the cel files
celFiles <- list.celfiles()
# assigning the cel files to affyraw variable
affyraw=ReadAffy(filenames = celFiles)
# making cdf file
tmp.env=make.cdf.env("ATH1121501_At_TAIRG.cdf")
# performing vsn normalization
vsn.data <- expresso(affyraw, normalize.method="vsn", bg.correct=F, pmcorrect.method="pmonly", summary.method="medianpolish")
# examining the normalization
boxplot(affyraw,col="red")
plot(exprs(affyraw)[,1:2], log = "xy", pch=".",
main="all")
# writing the result
write.table(vsn.data, file = "vsn1.txt", dec = ".", sep = "\t", quote = FALSE)
head(vsn.data[,1:2])
ExpressionSet (storageMode: lockedEnvironment)
assayData: 1 features, 2 samples
element names: exprs, se.exprs
protocolData
sampleNames: Col-0 24h primed.CEL.CEL Col-0 24h unprimed.CEL.CEL
varLabels: ScanDate
varMetadata: labelDescription
phenoData
sampleNames: Col-0 24h primed.CEL.CEL Col-0 24h unprimed.CEL.CEL
varLabels: sample
varMetadata: labelDescription
featureData: none
experimentData: use 'experimentData(object)'
Annotation: ath1121501
>
what is my fault please in the above code??? even i scared maybe i have written an incomplete normalized file

thank you,
my insisting on writing the data in a txt file is because i need the normalized file as an input for another tool for GRN inference anyway i did like below
> setwd("/usr/data/nfs6/izadi/Fereshteh thesis2/Data/Microarray/CEL files") > download.file("http://mbni.org/customcdf/20.0.0/tairg.download/ath1121501attairgcdf_20.0.0.tar.gz", "ath1121501attairgcdf_20.0.0.tar.gz") trying URL 'http://mbni.org/customcdf/20.0.0/tairg.download/ath1121501attairgcdf_20.0.0.tar.gz' Content type 'application/x-gzip' length 1573285 bytes (1.5 MB) ================================================== downloaded 1.5 MB > install.packages("ath1121501attairgcdf_20.0.0.tar.gz", repos=NULL, type="source") Installing package into ‘/usr/people/home/izadi/R/x86_64-redhat-linux-gnu-library/3.2’ (as ‘lib’ is unspecified) * installing *source* package ‘ath1121501attairgcdf’ ... ** R ** data ** inst ** preparing package for lazy loading Creating a generic function for ‘nchar’ from package ‘base’ in package ‘S4Vectors’ ** help *** installing help indices converting help for package ‘ath1121501attairgcdf’ finding HTML links ... done ath1121501attairgcdf html ath1121501attairgdim html geometry html ** building package indices ** testing if installed package can be loaded Creating a generic function for ‘nchar’ from package ‘base’ in package ‘S4Vectors’ * DONE (ath1121501attairgcdf) > library(affy) Loading required package: BiocGenerics Loading required package: parallel Attaching package: ‘BiocGenerics’ The following objects are masked from ‘package:parallel’: clusterApply, clusterApplyLB, clusterCall, clusterEvalQ, clusterExport, clusterMap, parApply, parCapply, parLapply, parLapplyLB, parRapply, parSapply, parSapplyLB The following object is masked from ‘package:stats’: xtabs The following objects are masked from ‘package:base’: anyDuplicated, append, as.data.frame, as.vector, cbind, colnames, do.call, duplicated, eval, evalq, Filter, Find, get, intersect, is.unsorted, lapply, Map, mapply, match, mget, order, paste, pmax, pmax.int, pmin, pmin.int, Position, rank, rbind, Reduce, rep.int, rownames, sapply, setdiff, sort, table, tapply, union, unique, unlist, unsplit Loading required package: Biobase Welcome to Bioconductor Vignettes contain introductory material; view with 'browseVignettes()'. To cite Bioconductor, see 'citation("Biobase")', and for packages 'citation("pkgname")'. > library(vsn) > celFiles <- list.celfiles() > abatch <- ReadAffy(filenames = celFiles, cdfname = "ath1121501attaircdf") > eset <- justvsn(abatch) vsn2: 506944 x 164 matrix (1 stratum). Please use 'meanSdPlot' to verify the fit. > boxplot(eset,col="red") Error in getCdfInfo(object) : Could not obtain CDF environment, problems encountered: Specified environment does not contain ath1121501attaircdf Library - package ath1121501attaircdf not installed Bioconductor - ath1121501attaircdf not available > write.table(eset, file = "eset.txt", dec = ".", sep = "\t", quote = FALSE) Error in as.data.frame.default(x[[i]], optional = TRUE) : cannot coerce class "structure("AffyBatch", package = "affy")" to a data.framei only need a vsn normalized file that the rownames are AGI code not _at
>