library(cBioPortalData)
library(AnVIL)
This document serves as a reporting tool for errors that occur when running our utility functions on the cBioPortal datasets.
cBioPortalData()
)Typically, the number of errors encountered via the API are low. There are only a handful of packages that error when we apply the utility functions to provide a MultiAssayExperiment data representation.
First, we load the error Rda
dataset.
api_errs <- system.file(
"extdata", "api", "err_api_info.rda",
package = "cBioPortalData", mustWork = TRUE
)
load(api_errs)
We can now inspect the contents of the data:
class(err_api_info)
## [1] "list"
length(err_api_info)
## [1] 11
lengths(err_api_info)
## 'sampleMap' does not have required columns
## 2
## Bad Request (HTTP 400).
## 4
## Bad Gateway (HTTP 502).
## 11
## Frequency of NA values higher than the cutoff tolerance
## 2
## `n` must be a single number, not an integer `NA`.
## 2
## Barcodes must start with 'TCGA'
## 3
## Stream error in the HTTP/2 framing layer [www.cbioportal.org]: HTTP/2 stream 0 was not closed cleanly: INTERNAL_ERROR (err 2)
## 3
## Internal Server Error (HTTP 500).
## 1
## Inconsistent build numbers found
## 46
## when the length of the supplied 'genome' vector is not 1, then it must be equal to the\n number of sequences
## 1
## Only two build types at a time can be used
## 1
There were about 11 unique errors during the last build run.
names(err_api_info)
## [1] "'sampleMap' does not have required columns"
## [2] "Bad Request (HTTP 400)."
## [3] "Bad Gateway (HTTP 502)."
## [4] "Frequency of NA values higher than the cutoff tolerance"
## [5] "`n` must be a single number, not an integer `NA`."
## [6] "Barcodes must start with 'TCGA'"
## [7] "Stream error in the HTTP/2 framing layer [www.cbioportal.org]: HTTP/2 stream 0 was not closed cleanly: INTERNAL_ERROR (err 2)"
## [8] "Internal Server Error (HTTP 500)."
## [9] "Inconsistent build numbers found"
## [10] "when the length of the supplied 'genome' vector is not 1, then it must be equal to the\n number of sequences"
## [11] "Only two build types at a time can be used"
The most common error was Inconsistent build numbers found
. This is
due to annotations from different build numbers that were not able to
be resolved.
To see what datasets (cancer_study_id
s) have that error we can use:
err_api_info[['Inconsistent build numbers found']]
## [1] "thyroid_gatci_2024" "acc_tcga_gdc" "blca_tcga_gdc"
## [4] "brca_tcga_gdc" "cesc_tcga_gdc" "dlbclnos_tcga_gdc"
## [7] "esca_tcga_gdc" "gbm_tcga_gdc" "hnsc_tcga_gdc"
## [10] "ccrcc_tcga_gdc" "prcc_tcga_gdc" "difg_tcga_gdc"
## [13] "hcc_tcga_gdc" "luad_tcga_gdc" "lusc_tcga_gdc"
## [16] "plmeso_tcga_gdc" "hgsoc_tcga_gdc" "paad_tcga_gdc"
## [19] "mnet_tcga_gdc" "prad_tcga_gdc" "read_tcga_gdc"
## [22] "skcm_tcga_gdc" "stad_tcga_gdc" "nsgct_tcga_gdc"
## [25] "thpa_tcga_gdc" "thym_tcga_gdc" "ucec_tcga_gdc"
## [28] "ucs_tcga_gdc" "um_tcga_gdc" "coad_tcga_gdc"
## [31] "soft_tissue_tcga_gdc" "brain_cptac_gdc" "breast_cptac_gdc"
## [34] "coad_cptac_gdc" "luad_cptac_gdc" "lusc_cptac_gdc"
## [37] "ohnca_cptac_gdc" "ovary_cptac_gdc" "pancreas_cptac_gdc"
## [40] "uec_cptac_gdc" "asclc_msk_2024" "cscc_ranson_2022"
## [43] "ucec_msk_2024" "pancan_mappyacts_2022" "nst_nfosi_ntap"
## [46] "brca_aurora_2023"
We can also have a look at the entirety of the dataset.
err_api_info
## $`'sampleMap' does not have required columns`
## [1] "acyc_fmi_2014" "braf_msk_archer_2024"
##
## $`Bad Request (HTTP 400).`
## [1] "glioma_msk_2018" "mbn_sfu_2023" "pcnsl_msk_2024"
## [4] "pancan_pcawg_2020"
##
## $`Bad Gateway (HTTP 502).`
## [1] "heme_msk_impact_2022" "brca_tcga" "gbm_tcga"
## [4] "difg_glass_2019" "kirc_tcga" "lihc_tcga"
## [7] "lusc_tcga" "mng_utoronto_2021" "paad_tcga"
## [10] "pcpg_tcga" "sarc_tcga"
##
## $`Frequency of NA values higher than the cutoff tolerance`
## [1] "mixed_selpercatinib_2020" "ucec_ccr_msk_2022"
##
## $``n` must be a single number, not an integer `NA`.`
## [1] "msk_met_2021" "msk_chord_2024"
##
## $`Barcodes must start with 'TCGA'`
## [1] "nsclc_tcga_broad_2016" "blca_msk_tcga_2020" "mixed_msk_tcga_2021"
##
## $`Stream error in the HTTP/2 framing layer [www.cbioportal.org]: HTTP/2 stream 0 was not closed cleanly: INTERNAL_ERROR (err 2)`
## [1] "pan_origimed_2020" "coadread_tcga_pub" "meso_tcga"
##
## $`Internal Server Error (HTTP 500).`
## [1] "makeanimpact_ccr_2023"
##
## $`Inconsistent build numbers found`
## [1] "thyroid_gatci_2024" "acc_tcga_gdc" "blca_tcga_gdc"
## [4] "brca_tcga_gdc" "cesc_tcga_gdc" "dlbclnos_tcga_gdc"
## [7] "esca_tcga_gdc" "gbm_tcga_gdc" "hnsc_tcga_gdc"
## [10] "ccrcc_tcga_gdc" "prcc_tcga_gdc" "difg_tcga_gdc"
## [13] "hcc_tcga_gdc" "luad_tcga_gdc" "lusc_tcga_gdc"
## [16] "plmeso_tcga_gdc" "hgsoc_tcga_gdc" "paad_tcga_gdc"
## [19] "mnet_tcga_gdc" "prad_tcga_gdc" "read_tcga_gdc"
## [22] "skcm_tcga_gdc" "stad_tcga_gdc" "nsgct_tcga_gdc"
## [25] "thpa_tcga_gdc" "thym_tcga_gdc" "ucec_tcga_gdc"
## [28] "ucs_tcga_gdc" "um_tcga_gdc" "coad_tcga_gdc"
## [31] "soft_tissue_tcga_gdc" "brain_cptac_gdc" "breast_cptac_gdc"
## [34] "coad_cptac_gdc" "luad_cptac_gdc" "lusc_cptac_gdc"
## [37] "ohnca_cptac_gdc" "ovary_cptac_gdc" "pancreas_cptac_gdc"
## [40] "uec_cptac_gdc" "asclc_msk_2024" "cscc_ranson_2022"
## [43] "ucec_msk_2024" "pancan_mappyacts_2022" "nst_nfosi_ntap"
## [46] "brca_aurora_2023"
##
## $`when the length of the supplied 'genome' vector is not 1, then it must be equal to the\n number of sequences`
## [1] "thca_tcga_pan_can_atlas_2018"
##
## $`Only two build types at a time can be used`
## [1] "pancan_mimsi_msk_2024"
cBioDataPack()
Now let’s look at the errors in the packaged datasets that are used for
cBioDataPack
:
pack_errs <- system.file(
"extdata", "pack", "err_pack_info.rda",
package = "cBioPortalData", mustWork = TRUE
)
load(pack_errs)
We can do the same for this data:
length(err_pack_info)
## [1] 5
lengths(err_pack_info)
## more columns than column names
## 12
## Frequency of NA values higher than the cutoff tolerance
## 5
## invalid class "ExperimentList" object: \n Non-unique names provided
## 2
## non-character argument
## 2
## 'wget' call had nonzero exit status
## 13
We can get a list of all the errors present:
names(err_pack_info)
## [1] "more columns than column names"
## [2] "Frequency of NA values higher than the cutoff tolerance"
## [3] "invalid class \"ExperimentList\" object: \n Non-unique names provided"
## [4] "non-character argument"
## [5] "'wget' call had nonzero exit status"
And finally the full list of errors:
err_pack_info
## $`more columns than column names`
## [1] "ccrcc_utokyo_2013" "gbm_cptac_2021"
## [3] "luad_mskimpact_2021" "mbl_dkfz_2017"
## [5] "pan_origimed_2020" "sarcoma_msk_2022"
## [7] "bowel_colitis_msk_2022" "prad_msk_mdanderson_2023"
## [9] "brca_tcga_pan_can_atlas_2018" "coadread_tcga_pan_can_atlas_2018"
## [11] "ov_tcga_pan_can_atlas_2018" "sarc_tcga_pan_can_atlas_2018"
##
## $`Frequency of NA values higher than the cutoff tolerance`
## [1] "ihch_mskcc_2020" "mixed_selpercatinib_2020"
## [3] "ucec_ccr_msk_2022" "mixed_msk_tcga_2021"
## [5] "ihch_msk_2021"
##
## $`invalid class "ExperimentList" object: \n Non-unique names provided`
## [1] "mpnst_mskcc" "stad_tcga_pub"
##
## $`non-character argument`
## [1] "pcpg_tcga_pub" "mbn_mdacc_2013"
##
## $`'wget' call had nonzero exit status`
## [1] "braf_msk_impact_2024" "braf_msk_archer_2024" "prostate_msk_2024"
## [4] "pcnsl_msk_2024" "pdac_msk_2024" "ucs_msk_2024"
## [7] "asclc_msk_2024" "lms_msk_2024" "crc_orion_2024"
## [10] "brca_aurora_2023" "hcc_msk_2024" "pancreas_msk_2024"
## [13] "pancan_mimsi_msk_2024"
sessionInfo()
## R version 4.4.2 (2024-10-31)
## Platform: x86_64-pc-linux-gnu
## Running under: Ubuntu 24.04.1 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.20-bioc/R/lib/libRblas.so
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.12.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_GB LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## time zone: America/New_York
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] survminer_0.5.0 ggpubr_0.6.0
## [3] ggplot2_3.5.1 survival_3.8-3
## [5] cBioPortalData_2.18.1 MultiAssayExperiment_1.32.0
## [7] SummarizedExperiment_1.36.0 Biobase_2.66.0
## [9] GenomicRanges_1.58.0 GenomeInfoDb_1.42.1
## [11] IRanges_2.40.1 S4Vectors_0.44.0
## [13] BiocGenerics_0.52.0 MatrixGenerics_1.18.0
## [15] matrixStats_1.4.1 AnVIL_1.18.2
## [17] AnVILBase_1.0.0 dplyr_1.1.4
## [19] BiocStyle_2.34.0
##
## loaded via a namespace (and not attached):
## [1] jsonlite_1.8.9 magrittr_2.0.3
## [3] magick_2.8.5 GenomicFeatures_1.58.0
## [5] farver_2.1.2 rmarkdown_2.29
## [7] BiocIO_1.16.0 zlibbioc_1.52.0
## [9] vctrs_0.6.5 memoise_2.0.1
## [11] Rsamtools_2.22.0 RCurl_1.98-1.16
## [13] tinytex_0.54 rstatix_0.7.2
## [15] htmltools_0.5.8.1 S4Arrays_1.6.0
## [17] BiocBaseUtils_1.8.0 lambda.r_1.2.4
## [19] curl_6.0.1 broom_1.0.7
## [21] Formula_1.2-5 SparseArray_1.6.0
## [23] sass_0.4.9 bslib_0.8.0
## [25] htmlwidgets_1.6.4 httr2_1.0.7
## [27] zoo_1.8-12 futile.options_1.0.1
## [29] cachem_1.1.0 commonmark_1.9.2
## [31] GenomicAlignments_1.42.0 mime_0.12
## [33] lifecycle_1.0.4 pkgconfig_2.0.3
## [35] Matrix_1.7-1 R6_2.5.1
## [37] fastmap_1.2.0 GenomeInfoDbData_1.2.13
## [39] shiny_1.10.0 digest_0.6.37
## [41] colorspace_2.1-1 RaggedExperiment_1.30.0
## [43] AnnotationDbi_1.68.0 ps_1.8.1
## [45] RSQLite_2.3.9 labeling_0.4.3
## [47] filelock_1.0.3 RTCGAToolbox_2.36.0
## [49] km.ci_0.5-6 RJSONIO_1.3-1.9
## [51] httr_1.4.7 abind_1.4-8
## [53] compiler_4.4.2 bit64_4.5.2
## [55] withr_3.0.2 backports_1.5.0
## [57] BiocParallel_1.40.0 carData_3.0-5
## [59] DBI_1.2.3 ggsignif_0.6.4
## [61] rappdirs_0.3.3 DelayedArray_0.32.0
## [63] rjson_0.2.23 tools_4.4.2
## [65] chromote_0.3.1 httpuv_1.6.15
## [67] glue_1.8.0 restfulr_0.0.15
## [69] promises_1.3.2 gridtext_0.1.5
## [71] grid_4.4.2 generics_0.1.3
## [73] gtable_0.3.6 KMsurv_0.1-5
## [75] tzdb_0.4.0 tidyr_1.3.1
## [77] websocket_1.4.2 data.table_1.16.4
## [79] hms_1.1.3 car_3.1-3
## [81] xml2_1.3.6 utf8_1.2.4
## [83] XVector_0.46.0 markdown_1.13
## [85] pillar_1.10.0 stringr_1.5.1
## [87] later_1.4.1 splines_4.4.2
## [89] ggtext_0.1.2 BiocFileCache_2.14.0
## [91] lattice_0.22-6 rtracklayer_1.66.0
## [93] bit_4.5.0.1 tidyselect_1.2.1
## [95] Biostrings_2.74.1 miniUI_0.1.1.1
## [97] knitr_1.49 gridExtra_2.3
## [99] bookdown_0.41 futile.logger_1.4.3
## [101] xfun_0.49 DT_0.33
## [103] stringi_1.8.4 UCSC.utils_1.2.0
## [105] yaml_2.3.10 evaluate_1.0.1
## [107] codetools_0.2-20 tibble_3.2.1
## [109] BiocManager_1.30.25 cli_3.6.3
## [111] xtable_1.8-4 munsell_0.5.1
## [113] processx_3.8.4 jquerylib_0.1.4
## [115] survMisc_0.5.6 Rcpp_1.0.13-1
## [117] GenomicDataCommons_1.30.0 dbplyr_2.5.0
## [119] png_0.1-8 XML_3.99-0.17
## [121] rapiclient_0.1.8 parallel_4.4.2
## [123] TCGAutils_1.26.0 readr_2.1.5
## [125] blob_1.2.4 bitops_1.0-9
## [127] scales_1.3.0 purrr_1.0.2
## [129] crayon_1.5.3 rlang_1.1.4
## [131] KEGGREST_1.46.0 rvest_1.0.4
## [133] formatR_1.14