I have converted a GSE object to an exprSet object "eset2" just in the
vignette but I wanna analyze microarray dataset again for example with:
 
  
 
 usually starting from CEL files i should edit:
 
  
 
 data<-ReadAffy(widget=T)
 
  
 
 pData(data)
 
  
 
 and then:
 
  
 
 data.rma<-rma(data) # for example!
 
  
 
 BUT STARTING FROM THE exprSet object "eset2" HOW CAN PROCEDE IN ORDER 
 TO
 reanalyze microarray
 
if i write:
 
            eset2.rma<-rma(eset2)
Errore in probeNames(object, subset) : no direct or inherited method for
function 'probeNames' for this call
 
 
where 
> eset2
Expression Set (exprSet) with 
      22283 genes
      76 samples
             phenoData object with 1 variables and 76 cases
       varLabels
            : samples
 
how can I do to create an affybatch object in order to use rma function
on my dataset?
 
ps: to obtain eset2 I have followed the instruction on the GEOquery
vignette:
 
> probesets <- Table(GPLList(gse)[[1]])$ID
> data.matrix <- log2(do.call("cbind", lapply(GSMList(gse), function(x)
{
+ tab <- Table(x)
+ mymatch <- match(probesets, tab$ID_REF)
+ return(tab$VALUE[mymatch])
+ })))
> rownames(data.matrix) <- probesets
> colnames(data.matrix) <- names(GSMList(gse))
> pdata <- data.frame(samples = names(GSMList(gse)))
> rownames(pdata) <- names(GSMList(gse))
> pheno <- new("phenoData", pData = pdata, varLabels =
as.list("samples"))
> eset2 <- new("exprSet", exprs = data.matrix, phenoData = pheno)
> eset2
 
 
 
Dr Alberto Goldoni
Medical Genetics Unit
S. Orsola-Malpighi Hospital
Via Massarenti n.9, Pad 11
40100 Bologna, Italy
Mobile Phone: 338-4145970
Fax: 051-636-4004
 <mailto:alberto.goldoni@eurogene.org> alberto.goldoni@eurogene.org
 <http://www.eurogene.org> www.eurogene.org
 <http://www.lagem.it> www.lagem.it
 

	[[alternative HTML version deleted]]

