Question: How to extract expression data from RAW data.
gravatar for hkarakurt
2.9 years ago by
hkarakurt90 wrote:

Hello. I downloaded GSE57821 RAW data by getGEOSuppFiles command in R and now I have "filelist.txt" and GSE57821_RAW.tar files.

My command: gse <- getGEO("GSE57821", GSEMatrix = TRUE)

raw <- getGEOSuppFiles("GSE57821")

I want to normalize RAW data. Which way should I follow? How can I retrieve expression data from raw data?

exprs command is not working both in gse and raw files.

I am relatively new at gene expression analysis in R.

Thank you.

raw expression • 1.7k views
ADD COMMENTlink modified 2.9 years ago by h.mon29k • written 2.9 years ago by hkarakurt90
gravatar for h.mon
2.9 years ago by
h.mon29k wrote:

The "raw" data is, in fact, a tar file containing several compressed tables of counts. Maybe someone knows a function which directly reads the "raw" file and creates a ExpressionSet or DGEList - I don't. But getGEOSuppFiles() does not either:

No parsing of the downloaded files is attempted, since the file format is not generally knowable by the computer.

You have to untar GSE57821_RAW.tar and then uncompress the counts files, read one by one with read.table() and finally create the ExpressionSet or DGElist objects yourself.

ADD COMMENTlink written 2.9 years ago by h.mon29k
Please log in to add an answer.


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