I am trying to process Agilent data by using limma R package. For GSE10469, I used the following code
raw_data <- read.maimages(pdata[,1], source = "agilent") # pdata file is having group information I get the error: Error in readGenericHeader(fullname, columns = columns, sep = sep) : Specified column headings not found in file
When I try
raw_data <- read.maimages(pdata[,1], source = "agilent", green.only = TRUE) Read GSM264878.txt Error in RG[[a]][, i] <- obj[, columns[[a]]] : number of items to replace is not a multiple of replacement length
And also
raw_data <- read.maimages(pdata[,1], source = "agilent", green.only = FALSE) Error in readGenericHeader(fullname, columns = columns, sep = sep) : Specified column headings not found in file
For GSE32006
raw_data <- read.maimages(pdata[,1], source = "agilent") Error in readGenericHeader(fullname, columns = columns, sep = sep) : Specified column headings not found in file
And
raw_data <- read.maimages(pdata[,1], source = "agilent", green.only = TRUE) Read GSM792633.txt Read GSM792634.txt Read GSM792635.txt Read GSM792636.txt Read GSM792637.txt Read GSM792638.txt Error in readGenericHeader(fullname, columns = columns) : Specified column headings not found in file
The files which are read from GSE32006 are gene expression, where as other failed files are from exon array.
So how can I deal with all these issues.
Code snippets are not useful! Unless you show exactly what you did, you are expecting people to guess at what you might have done, and most people are too busy to bother with such things. You need to show a short, self-contained (e.g., anybody can run) bit of code to show exactly what you did and where the error is.