[BioC] ExpressionSet created by procset method in ArrayExpress

Md.Mamunur Rashid mamunur.rashid at kcl.ac.uk
Tue Jul 28 15:34:18 CEST 2009


Hi everyone,

I have recently posted a query names "How to use ExpressionSet for 
further 
Analysis"(http://article.gmane.org/gmane.science.biology.informatics.conductor/24343/match=use+expressionset+further+analysis)

I have created an ExpressionSet object using procset method in 
ArrayExpress package.
But the problem I found out that in the ExpressionSet the numeric values 
are in double quote.  (i.e. "123.02".  ) Now when I pass this expression 
set in to the following R code

 > library(lumi)
 > b <- lumiB("p1")    // p1 is the expression set I have created using 
procset
 > MAplot(b)

it shows me an error that Error in log(c("7413.841", "347.639", 
"233.5174", "7120.314", "35.45117",  :  Non-numeric argument to 
mathematical function

can anyone suggest me why the numeric data is converted to 
string/characters while creating ExpressionSet. I will really appriciate 
any help as I am really stuck in this point. Thanks in advance.

regards,
Mamun


below is the code I have used:


> library(ArrayExpress)
 > etabm421 <- list(path = ".",rawadata = NULL, rawfiles = NULL, 
procdata = "E-TABM-421.processed.zip",procfile = 
"E-TABM-421-processed-data-1583036605.txt",sdrf = "E-TABM-421.sdrf.txt", 
idf = "E-TABM-421.idf.txt", adf = "A-MEXP-931.adf.txt")
 > cn = getcolproc(ETABM421)
 > cn
[1] "Reporter REF"              "BeadStudio:AVG_Signal"   
[3] "BeadStudio:Detection Pval"
 > proset = procset(ETABM421,cn[2])



More information about the Bioconductor mailing list