Search
Question: Error in TCGAquery of TCGAbiolinks Package
1
gravatar for Beginner
11 weeks ago by
Beginner50
Beginner50 wrote:

Dear All,

I'm looking PAM50 molecular Subtypes of TCGA Breast RNA-Seq data. I'm using the following code. But I'm having an error with TCGAquery function. I have used similar code long back and it worked but now it is not working with TCGAquery function.

source("http://bioconductor.org/biocLite.R")

biocLite("TCGAbiolinks")
library(TCGAbiolinks)
cancer <- "BRCA"
PlatformCancer <- "IlluminaHiSeq_RNASeqV2"
dataType <- "rsem.genes.results"
pathCancer <- "TCGAData/miRNA"

datQuery <- TCGAquery(tumor = cancer, platform = PlatformCancer, level = "3") 

Error in TCGAquery(tumor = cancer, platform = PlatformCancer, level = "3") : 
  could not find function "TCGAquery"

lsSample <- TCGAquery_samplesfilter(query = datQuery)

# get subtype information
dataSubt <- TCGAquery_subtype(tumor = cancer)
lumA <- dataSubt[which(dataSubt$PAM50.mRNA == "Luminal A"),1]
allSamples <- lsSample$IlluminaHiSeq_RNASeqV2 #1218 total samples

Any help please !!

 

 

ADD COMMENTlink modified 11 weeks ago by Benjamin Ostendorf50 • written 11 weeks ago by Beginner50
1
gravatar for Benjamin Ostendorf
11 weeks ago by
United States
Benjamin Ostendorf50 wrote:

Hi Beginner,

after the transfer of the TCGA data to GDC the commands changed - try GDCquery instead of TCGAquery. Also make sure you have the latest version of TCGAbiolinks installed and take a look at the current vignette. 

Best,
Benjamin

ADD COMMENTlink written 11 weeks ago by Benjamin Ostendorf50
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.2.0
Traffic: 251 users visited in the last hour