TCGA, TCGAbiolink, gene expression
0
0
Entering edit mode
14 months ago
Rob ▴ 170

Hi friends, I m using the following code to download STAR count of KIRP. How should I download just the primary tumor & just coding genes (not all genes)? Also, how should I download the corresponding clinical data? Thanks

query <- GDCquery(project = "TCGA-KIRP",
                  data.category = "Transcriptome Profiling",
                  data.type = "Gene Expression Quantification",
                  experimental.strategy = "RNA-Seq",
                  workflow.type = "STAR - Counts")


# Download data using api
GDCdownload(query, method = "api")

#Prepare data
data <- GDCprepare(query)
?assay()
#generate count matrix
rna_TCGA_STAR_KIRP <- as.data.frame(SummarizedExperiment::assay(data))
#rna_STAR_test <- as.data.frame(SummarizedExperiment::assay(data_TCGA_STAR_KIRP))

write.table(rna_TCGA_STAR_KIRP, "rna_TCGA_STAR_KIRP.txt", sep='\t')
data RNA-seq count STAR • 433 views
ADD COMMENT

Login before adding your answer.

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