Error in inputting multiple illumina microarray idat files to limma
1
0
Entering edit mode
2.4 years ago
maria2019 ▴ 210

Hi,

I have illumina microarray idat files and I want to read them to limma. Below is my code:

library(limma)

idatfiles = dir(path="File_path", pattern = ".idat")
bgxfile = dir(path = "File_path", pattern = ".bgx")
data = readIDAT(idatfiles, bgxfile)

However, I get the following error:

**Error in readIDAT(idatfiles, bgxfile) : 
  Unable to find file 7963862004_A_Grn.idat7963862004_B_Grn.idat7963862004_C_Grn.idat7963862004_D_Grn.idat7963862004_E_Grn.idat7963862004_F_Grn.idat7963862004_G_Grn.idat7963862004_H_Grn.idat7963862004_I_Grn.idat7963862004_J_Grn.idat7963862004_K_Grn.idat7963862004_L_Grn.idat
In addition: Warning message:
In if (!file.exists(file)) { :
  the condition has length > 1 and only the first element will be used
Execution halted**

Any comment on how I can fix this?

I also tried this:

idatfiles <- list.files(path="File_path", pattern ="\\.idat")

and it did not work either

R idat limma microarray • 689 views
ADD COMMENT
1
Entering edit mode

What are the contents of idatfiles and bgxfile? Also, form where did you obtain the data? I processed a handful of Illumina chips relatively recently.

ADD REPLY
3
Entering edit mode
2.4 years ago
Gordon Smyth ★ 4.4k

First problem. The limma function is called read.idat, not readIDAT.

Second problem. You haven't told read.idat (or readIDAT) which directory the files are in, so naturally the function can't find them.

You need

x <- read.idat(file.path("File_path",idatfiles), file.path("File_path",bgxfile"))
ADD COMMENT

Login before adding your answer.

Traffic: 1799 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6