Question: Retrieve TCGA file ID and associated barcode with TCGAbiolinks: "GDC server down"
0
gravatar for user31888
14 months ago by
user3188870
United States
user3188870 wrote:

I am trying to follow the example here to get file names and associated barcodes using TCGAbiolinks.

> library(TCGAbiolinks)
> query <- GDCquery(project = c("TCGA-BRCA"),
              data.category = "Raw Sequencing Data",  
              sample.type = "Primary solid Tumor")
Error in value[[3L]](cond) :
  GDC server down, try to use this package later

I checked the status of the server but apparently it is ok:

> library(GenomicDataCommons)
> status()
$commit
[1] "e9e20d6f97f2bf6dd3b3261e36ead57c56a4c7cc"

$data_release
[1] "Data Release 12.0 - June 13, 2018"

$status
[1] "OK"

$tag
[1] "1.14.1"

$version
[1] 1

Is it possible to similarly get the file names and corresponding barcodes using GenomicDataCommons, which seems to work better?

ADD COMMENTlink modified 14 months ago • written 14 months ago by user3188870
3
gravatar for Kevin Blighe
14 months ago by
Kevin Blighe48k
Kevin Blighe48k wrote:

Hello,

Please install the GitHub version by doing the following:

require(devtools)
install_github("BioinformaticsFMRP/TCGAbiolinks")

Then it should work.

query <- GDCquery(project = c("TCGA-BRCA"),
+               data.category = "Raw Sequencing Data",  
+               sample.type = "Primary solid Tumor")
--------------------------------------
o GDCquery: Searching in GDC database
--------------------------------------
Genome of reference: hg38
--------------------------------------------
oo Accessing GDC. This might take a while...
--------------------------------------------
ooo Project: TCGA-BRCA
--------------------
oo Filtering results
--------------------
ooo By sample.type
----------------
oo Checking data
----------------
ooo Check if there are duplicated cases
ooo Check if there results for the query
-------------------
o Preparing output
-------------------

Kevin

ADD COMMENTlink written 14 months ago by Kevin Blighe48k

Works ! However, it cannot find function datatable in the second part of the code.

> datatable(getResults(query, rows = 1:100,cols = c("file_name","cases")), 
              filter = 'top',
              options = list(scrollX = TRUE, keys = TRUE, pageLength = 5), 
              rownames = FALSE)
Error in datatable(getResults(query, rows = 1:100, cols = c("file_name",  :
  could not find function "datatable"
ADD REPLYlink written 14 months ago by user3188870

My mistake. I just had to get rid of the datatable from the DT package (cannot open browser on the server I use anyway).

ADD REPLYlink modified 14 months ago • written 14 months ago by user3188870

Great! - all good now?

ADD REPLYlink written 14 months ago by Kevin Blighe48k
Please log in to add an answer.

Help
Access

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