Question

Download TCGA dataset like a S4 object

0

Entering edit mode

Giuseppe ▴ 20

@giuseppe-16310

Last seen 5.6 years ago

Hi all,

I'm making a function that uses input TCGA datasets. I know that for the purpose of reproducibility, inputs need to be BioC objects, and not text files.

My question is: is there any package that allows me download TCGA datasets like a BioC objects (S4) ?

Thank you so much!

dataset tcga • 1.7k views

ADD COMMENT • link updated 5.7 years ago by mario.zanfardino ▴ 160 • written 5.7 years ago by Giuseppe ▴ 20

score 0 · Answer 1 · 2018-08-01

0

Entering edit mode

mario.zanfardino ▴ 160

@mariozanfardino-15232

Last seen 3.7 years ago

Naples (Italy)

Hi Giuseppe, try curatedTCGA package:

https://bioconductor.org/packages/release/data/experiment/html/curatedTCGAData.html

"This package provides publicly available data from The Cancer Genome Atlas (TCGA) Bioconductor MultiAssayExperiment class objects. These objects integrate multiple assays (e.g. RNA-seq, copy number, mutation, microRNA, protein, and others) with clinical / pathological data. The MultiAssayExperiment class links assay barcodes with patient IDs, enabling harmonized subsetting of rows (features) and columns (patients / samples) across the entire experiment."

ADD COMMENT • link 5.7 years ago mario.zanfardino ▴ 160

1

Entering edit mode

Thank you Mario for your answer. I want to download the following dataset :

I executed the following instruction: curatedTCGAData("SKCM","Methylation",F) to download HumanMethylation450 dataset, but the output is not that I want.

I never used this package, can you help me to resolve this problem?

ADD REPLY • link 5.7 years ago Giuseppe ▴ 20

1

Entering edit mode

This run in my env:

SKCM <- curatedTCGAData(diseaseCode = "SKCM", assays = "Methylation", dry.run = FALSE)

Obviously, the data are in MultiAssayExperiment format (a Bioconductor object-oriented S4 class)(https://bioconductor.org/packages/release/bioc/vignettes/MultiAssayExperiment/inst/doc/MultiAssayExperiment.html).

ADD REPLY • link 5.7 years ago mario.zanfardino ▴ 160

1

Entering edit mode

Thank you very much for your support! I have last question.. I executed this instruction:

> experiments(SKCM)

ExperimentList class object of length 1:

[1] SKCM_Methylation-20160128: SummarizedExperiment with 485577 rows and 475 columns

I want to convert ExperimentList class to a data frame. How can I do it?

ADD REPLY • link 5.7 years ago Giuseppe ▴ 20

2

Entering edit mode

Probably you don't want to convert this to a data frame, but learn to use the SummarizedExperiment class, see for instance here and the package vignette here.

ADD REPLY • link 5.7 years ago Martin Morgan 25k

1

Entering edit mode

df <- as.data.frame(wideFormat(SKCM[1:10,1:10 ,"SKCM_Methylation-20160128"], colDataCols = c(1:10)))

df subsetted for 10 Patient, 10 features and 10 colData

ADD REPLY • link 5.7 years ago mario.zanfardino ▴ 160

0

Entering edit mode

thank you so much! The Only problem is that download Methylation dataset is very slow! And the S4 object occupies so much memory

ADD REPLY • link 5.7 years ago Giuseppe ▴ 20