This is a test version of Biostars. For the public version, visit https://www.biostars.org.
CPTAC data download

Hi friends

Is there an R code to download CPTAC mutation data? I tried to use the following code but it did not work. I would appreciate your help. Thanks

install.packages("tidyverse")
library(tidyverse)
if (!require("BiocManager"))
    install.packages("BiocManager")
BiocManager::install('GenomicDataCommons')
library(GenomicDataCommons)

qfiles = files() %>% filter( ~ cases.project.project_id == 'CPTAC-3-kirc' &
                            type == 'mutation' &
                            access == "open" #&
                            #analysis.workflow_type == 'STAR - Counts'
                              )
manifest_df = qfiles %>% manifest() # I got error here
nrow(manifest_df)

The error I get is : Error in if (ncol(tmp) < 5) { : argument is of length zero

r cptac
• 1,553 views • link
updated by Ram Ram • written by Rob • Rob
  1. What do you mean "it did not work"? You need to give us as much detail as you can
  2. Your code is incomplete. What is files()? What libraries do you have loaded? Add the output to sessionInfo().
• 0 views • link
written by Ram Ram

Hi thanks for responding. here is the full code with libraries and packages. file here is all the data before filtering for only mutation.

install.packages("tidyverse")
library(tidyverse)
if (!require("BiocManager"))
    install.packages("BiocManager")
BiocManager::install('GenomicDataCommons')
library(GenomicDataCommons)

qfiles = files() %>% filter( ~ cases.project.project_id == 'CPTAC-3-kirc' &
                            type == 'mutation' &
                            access == "open" #&
                            #analysis.workflow_type == 'STAR - Counts'
                              )
manifest_df = qfiles %>% manifest() # I got error here
nrow(manifest_df)

The error I get is : Error in if (ncol(tmp) < 5) { : argument is of length zero

• 0 views • link
updated by Ram Ram • written by Rob • Rob

I've moved the relevant content to your post. You've still not explained where files() comes from and you've also not added the output to sessionInfo() yet.

• 0 views • link
written by Ram Ram

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • cannot install GO.db in windows
    written by Naila • Naila
    #code .libPaths() if (!require("BiocManager", quietly = TRUE)) install.packages("BiocManager") BiocManager::install(c("clusterProfiler", "GO.db", "org.Hs.eg.db"), dependencies = TRUE, force = TRUE) sessionInfo() #output The downloaded source packages are in …
  • Differetial analysis EdgeR
    written by maryak • maryak
    if (!require("BiocManager", quietly = TRUE)) install.packages("BiocManager") BiocManager::install("TCGAbiolinks", force = TRUE) library(edgeR) library(SummarizedExperiment) library(TCGAbiolinks) BiocManager::install("edgeR") #install.packages("DMwR") #library(DMwR) BiocManager::install("minfi") library(minfi) library(limma) library(dplyr) # Load necessary libraries library(edgeR) …
  • DESeq2 Install in Jupyter Notebook
    written by Lolo • Lolo
    I am having trouble installing the DESeq2 package within a Jupyter Notebook. I have tried installing using Conda using a Python kernel, and using BiocManager …
  • ImmuCellAI
    written by sana • sana
    Hello everyone, I'm trying to use the code for ImmuCellAI (the old version) Here is my code: ```r install.packages("tidyverse") library(tidyverse) library(lubridate) install.packages("BiocManager") BiocManager::install("usethis") install.packages("devtools") library(devtools) …
  • Mutation, MAF, CPTAC-3
    written by Rob • Rob
    Hi friends How can I download mutation data from CPTAC for kidney cancer (ccRCC)? Is there any R-code for this?
  • I am trying to install TCGAbiolinks package in both google colab with r runtime and new version of …
    written by Salsabil • Salsabil
    **Code in google colab:** if (!require("BiocManager", quietly = TRUE)) install.packages("BiocManager") BiocManager::install("TCGAbiolinks") **Error message:** Installing package into ‘/usr/local/lib/R/site-library’ (as ‘lib’ is unspecified) 'getOption("repos")' replaces Bioconductor standard …
  • DESeq2 normalization
    written by Rob • Rob
    Hi friends I am using this code to normalize my data. is this code complete for normalizing or the last 3 lines are not needed? …
  • HT-Seq count data coding gene
    written by Rob • Rob
    Hello friends, I want to download HT-seq data from TCGA biolink. How can I download only coding genes? what should I add to my code? …
  • problem installing rtracklayer
    written by Meghan.T • Meghan.T
    I'm pretty new to Rstudio on ubuntu .so I'm trying to install "rtracklayer" package from BiocManager I used this code: if (!requireNamespace("BiocManager", quietly = TRUE)) …
  • heatmap in R
    written by Rob • Rob
    Hello friends I want to make a heatmap for gene data in R. this is the code: geneExp <- read.table("file.txt", header=T, sep="\t") head(geneExp) # Extract …