Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Error when reading microarray idat files from illumina
0
gravatar for salamandra
4 weeks ago by
salamandra • 240
salamandra • 240 wrote:

Hi, I tried to read an idat file with R package illuminaio:

library(illuminaio)
idatfile <- 'C:/Users/Tania/Documents/GSE26677/GSE26677_non-normalized.txt'
idat <- readIDAT(idatfile, what='all')

and also with limma package:

library(limma)
bigxfile <- 'C:/Users/Tania/Documents/GSE26677/FILE/GPL6884_HumanWG-6_V3_0_R0_11282955_A.bgx'
idat <- read.idat(idatfile, bigxfile)

and got this error:

Error in readIDAT(idatfile, what = "all") : 
  Cannot read IDAT file. File format error. Unknown magic: ID_R

Seems something is wrong with idat file, but I tried these commands with several different idat files on GEO (GSE22247, GSE79308, GSE43336, GSE26677) and always got that error.

Does anyone know how to solve this?

The top of idat file is:

ID_REF  NSP Detection Pval  SP  Detection Pval
ILMN_1762337    75.85098    0.2555995   68.39282    0.6166008
ILMN_2055271    92.8892 0.03820817  95.76293    0.02371541
microarray R limma illuminaio illumina • 93 views
ADD COMMENT • link •
modified 4 weeks ago by h.mon ♦ 28k • written 4 weeks ago by salamandra • 240
0
gravatar for h.mon
4 weeks ago by
h.mon ♦ 28k
Brazil
h.mon ♦ 28k wrote:

Use the readBGX() function from the illuminaio package:

library( illuminaio )
bigxfile <- 'C:/Users/Tania/Documents/GSE26677/FILE/GPL6884_HumanWG-6_V3_0_R0_11282955_A.bgx.gz'
bgx <- readBGX( bigxfile )

Note that readBGX() can read compressed and uncompressed files.

You may also want to consider the GEOquery package, see for example Cannot read Illumina data (.bgx and .txt format) in R for microanalysis .
ADD COMMENT • link modified 4 weeks ago • written 4 weeks ago by h.mon ♦ 28k

I can read bgx, but not the txt file. Do you mean the txt file is not a idat file is that it? I tried the link you suggested. Specifically tried:

library(GEOquery)
eset = getGEO('GSE26677')[[1]]

but when doing:

exprs(eset)

gets the data already quantile normalized (I know that cause gives same values as normalized data when we click in specific samples in GEO). How do I get expressionset of raw data of that study?

ADD REPLY • link written 4 weeks ago by salamandra • 240
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1841 users visited in the last hour