Loading a dataset via acession number
1
0
Entering edit mode
4.4 years ago

I am trying to load this data set: https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA251582

I have installed GEOquerry a while back. When i run the following code I get this result:

library(GEOquery)
eset <- getGEO("PRJNA251582", getGPL = FALSE)


> eset <- getGEO("PRJNA251582", getGPL = FALSE)
Error in file.exists(destfile) : object 'destfile' not found

Also what is the difference betwen calling eset <- getGEO("PRJNA251582", getGPL = FALSE) and eset <- getGEO("PRJNA251582", getGPL = FALSE)[[1]] ?

R sequencing • 841 views
ADD COMMENT
1
Entering edit mode
4.4 years ago
ATpoint 88k

This is RNA-seq, not microarrays which GEOquery assumes, so it is expected that it produces no output. RNA-seq is sequencing data, you therefore will need to download the fastq files, align/quantify to get a matrix of raw counts.

For starters maybe this paper: https://peerj.com/preprints/27283/

Please also search biostars and the web for RNA-seq processing tutorials.

For downloading the files (if you really want to do that), search sra-explorer.info for the dataset to get download links. Beside that, there are plenty of posts and tutorials on how to download and process RNA-seq data, please use the search function.

ADD COMMENT

Login before adding your answer.

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