How can I pre-process illumina raw file
1
0
Entering edit mode
@agaz-hussain-wani-7620
Last seen 6.1 years ago
India

I have some raw `Illumina` data and when processing with `lumi` R package,  gives some error.

Actually I pre-process some other data which works fine by using the following code

x.lumi <- lumiR("filename")

lumiExpr <- lumiExpresso(x.lumi, bg.correct = TRUE, normalize = TRUE, verbose = TRUE)

exprs <- exprs(lumiExpr)

pvalue <- detection(lumiExpr)

But when I try to process the data  https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE52519, which looks something like

ID_REF    Sample 1.2.    Detection Pval    Sample 2.    Detection Pval    Sample 3.    Detection Pval    Sample 4.    Detection Pval    Sample 5.    Detection Pval    Sample 6.    Detection Pval    Sample 7.    Detection Pval    Sample 8.    Detection Pval    Sample 2.2.    Detection Pval    Sample 10.    Detection Pval    Sample 11.    Detection Pval    Sample 12.    Detection Pval
ILMN_1762337    43    0.28063    36.7    0.75758    41.5    0.55336    31.4    0.92754    35.2    0.88406    37.4    0.70619    34.4    0.80632    36.5    0.86561    39    0.89987    36.8    0.83399    44.1    0.36232    42.2    0.29381
ILMN_2055271    41.4    0.4361    51.1    0.01449    46.8    0.13175    35.7    0.5336    44.3    0.13307    43.7    0.13043    38.9    0.34519    44.1    0.16601    46.1    0.278    43.7    0.19499    43.4    0.4361    52    0.01976
ILMN_1736007    43.2    0.26746    41.2    0.29908    46.6    0.14098    38.6    0.23584    44.9    0.10408    38.1    0.639    37    0.54809    42.3    0.28459    52.2    0.04348    41.1    0.4282    46.3    0.20685    41.3    0.37813
ILMN_2383229    46    0.10804    51.6    0.01449    46    0.16996    36.3    0.4585    37.2    0.7365    39.2    0.5112    40.3    0.20553    43.3    0.20158    38.8    0.90777    36.6    0.85771    39.3    0.84058    37.9    0.72596
ILMN_1806310    43    0.28063    62.5    0.00527    47.9    0.09091    36.5    0.42951    43    0.20026    43.2    0.15547    36.2    0.63636    40.6    0.44664    48.6    0.12912    42.1    0.31752    47.7    0.12912    40.4    0.45982

I get the error

Error in gregexpr("\t", dataLine1)[[1]] : subscript out of bounds

Why there is a problem and how can I process such data.

 

illumina microarray • 1.5k views
ADD COMMENT
0
Entering edit mode

What do you mean by 'process the data'? You gave code that worked, but not the code that didn't work! How are people to help if you don't show what you did, and what happened?

ADD REPLY
0
Entering edit mode

'process the data', I mean to generate expression values and p-values from the raw data. The same code is used for both the files.

ADD REPLY
0
Entering edit mode
@james-w-macdonald-5106
Last seen 8 hours ago
United States

You are trying to process already processed data. If you want to re-run using lumi, you need to get the raw data, using e.g., getGEOSuppFiles from the GEOquery package.

ADD COMMENT

Login before adding your answer.

Traffic: 772 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