Hello,
I am trying to do background correction/ quantile normalization on a .txt file of expression values (.CEL data is unavailable to me). GSE_data is a 47,000 x 73 dataframe with headers of expression values. The data was gathered on an Illumina microarray.
I'm trying to use lumi like so:
library(lumi)
testMat <- as.matrix(GSE_data)
testMat <- testMat[,-1] # to get rid of probe IDs in first column
eset<-new("ExpressionSet",exprs=testMat)
# background correction/RMA normalization
lumiObj <- lumiB(eset, method='bgAdjust',verbose=TRUE)
error in evaluating the argument 'object' in selecting a method for function 'exprs': Error in bgAdjust(x.lumi, ...) : The object should be class "LumiBatch"!
Can anyone see what I'm doing wrong? What's the best way to do background correction/quantile normalization on .txt data?
Thank you so much!
EDIT: When I rewrite the eset definition, I get an error:
eset<-new('LumiBatch', exprs = test)
Error in `featureNames<-`(`*tmp*`, value = c("1", "2", "3", "4", "5", :
'value' length (47308) must equal feature number in AssayData (47308)
EDIT 2: Also, using lumiR doesn't seem to work - inputting this (importme
= file in my path that is a .txt file of the expression values only)
try.lumi <- lumiR('importme')
Error in gregexpr("\t", dataLine1)[[1]] : subscript out of bounds