[BioC] Reading Illumina data from GenomeStudio in limma

Josh [guest] guest at bioconductor.org
Wed Aug 14 04:30:15 CEST 2013

I have data exported from GenomeStudio (1.6) that I can correctly read and pre-process with lumi. However, when I try using read.ilmn in limma, I found that I get an error saying "Error in `rownames<-`(`*tmp*`, value = list(PROBE_ID = c("ILMN_1762337",  : length of 'dimnames' [1] not equal to array extent." I have found that this was due to the inclusion of extra columns in my probe profiles (TargetID, ProbeID, SPECIES, SOURCE, TRANSCRIPT and on). When I removed them, I can correctly read the probe and control profiles with read.ilmn. Is there a better fix around this than manually removing the columns?

 -- output of sessionInfo(): 

R version 3.0.0 (2013-04-03)
Platform: x86_64-apple-darwin10.8.0 (64-bit)

[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] parallel  stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] lumi_2.12.0        Biobase_2.20.0     BiocGenerics_0.6.0 limma_3.16.2

Sent via the guest posting facility at bioconductor.org.

More information about the Bioconductor mailing list