TCGAbiolinks GDCprepare produces error in vectbl_as_col_location()
0
0
Entering edit mode
lescano ▴ 10
@76ac7b25
Last seen 2.1 years ago
Canada

Hello all,

Tried downloading raw STAR counts using TCGAbiolinks today using this code snippet


tcga_rna_query <- GDCquery(
  project = "TCGA-LAML",
  data.category = "Transcriptome Profiling",
  data.type = "Gene Expression Quantification",
  workflow.type = "STAR - Counts",
  experimental.strategy = "RNA-Seq"
)
GDCdownload(tcga_rna_query)
tcga_rna_counts <- GDCprepare(tcga_rna_query, summarizedExperiment = FALSE)

Error in `vectbl_as_col_location()`:
! Can't subset columns past the end.
ℹ Locations 2, 3, and 4 don't exist.
ℹ There is only 1 column.
Run `rlang::last_error()` to see where the error occurred.
There were 50 or more warnings (use warnings() to see the first 50)

sessionInfo()
R version 4.1.3 (2022-03-10)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: macOS Big Sur/Monterey 10.16

Matrix products: default
BLAS/LAPACK: /Users/suika-san/opt/miniconda3/lib/libopenblasp-r0.3.20.dylib

locale:
[1] en_GB.UTF-8/en_GB.UTF-8/en_GB.UTF-8/C/en_GB.UTF-8/en_GB.UTF-8

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] forcats_0.5.2       purrr_0.3.4         readr_2.1.2        
 [4] tidyr_1.2.0         tibble_3.1.7        ggplot2_3.3.6      
 [7] tidyverse_1.3.2     edgeR_3.36.0        limma_3.50.1       
[10] stringr_1.4.1       dplyr_1.0.9         TCGAbiolinks_2.25.3

loaded via a namespace (and not attached):
 [1] fs_1.5.2                    bitops_1.0-7               
 [3] matrixStats_0.62.0          lubridate_1.8.0            
 [5] bit64_4.0.5                 filelock_1.0.2             
 [7] progress_1.2.2              httr_1.4.4                 
 [9] GenomeInfoDb_1.30.1         tools_4.1.3                
[11] backports_1.4.1             utf8_1.2.2                 
[13] R6_2.5.1                    DBI_1.1.3                  
[15] BiocGenerics_0.40.0         colorspace_2.0-3           
[17] withr_2.5.0                 tidyselect_1.1.2           
[19] prettyunits_1.1.1           bit_4.0.4                  
[21] curl_4.3.2                  compiler_4.1.3             
[23] cli_3.3.0                   rvest_1.0.3                
[25] Biobase_2.54.0              xml2_1.3.3                 
[27] DelayedArray_0.20.0         scales_1.2.1               
[29] rappdirs_0.3.3              digest_0.6.29              
[31] R.utils_2.12.1              XVector_0.34.0             
[33] pkgconfig_2.0.3             MatrixGenerics_1.6.0       
[35] dbplyr_2.2.1                fastmap_1.1.0              
[37] readxl_1.4.0                rlang_1.0.2                
[39] RSQLite_2.2.8               generics_0.1.3             
[41] jsonlite_1.8.0              vroom_1.5.7                
[43] googlesheets4_1.0.1         R.oo_1.25.0                
[45] RCurl_1.98-1.6              magrittr_2.0.3             
[47] GenomeInfoDbData_1.2.7      Matrix_1.4-1               
[49] Rcpp_1.0.8.3                munsell_0.5.0              
[51] S4Vectors_0.32.3            fansi_1.0.3                
[53] lifecycle_1.0.1             R.methodsS3_1.8.2          
[55] stringi_1.7.6               SummarizedExperiment_1.24.0
[57] zlibbioc_1.40.0             plyr_1.8.7                 
[59] BiocFileCache_2.2.1         grid_4.1.3                 
[61] blob_1.2.3                  parallel_4.1.3             
[63] crayon_1.5.2                lattice_0.20-45            
[65] haven_2.5.0                 Biostrings_2.62.0          
[67] hms_1.1.2                   KEGGREST_1.34.0            
[69] locfit_1.5-9.5              knitr_1.40                 
[71] pillar_1.8.1                GenomicRanges_1.46.1       
[73] TCGAbiolinksGUI.data_1.15.1 biomaRt_2.50.3             
[75] stats4_4.1.3                reprex_2.0.2               
[77] XML_3.99-0.9                glue_1.6.2                 
[79] downloader_0.4              remotes_2.4.2              
[81] modelr_0.1.9                data.table_1.14.2          
[83] BiocManager_1.30.19         png_0.1-7                  
[85] vctrs_0.4.1                 tzdb_0.3.0                 
[87] cellranger_1.1.0            gtable_0.3.1               
[89] assertthat_0.2.1            cachem_1.0.6               
[91] xfun_0.31                   broom_1.0.1                
[93] googledrive_2.0.0           gargle_1.2.1               
[95] AnnotationDbi_1.56.2        memoise_2.0.1              
[97] IRanges_2.28.0              ellipsis_0.3.2

This code has worked for me in the past. Any suggestions for troubleshooting?

Thanks!

TCGAbiolinks • 632 views
ADD COMMENT

Login before adding your answer.

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