From a14a1a38a6cce1a3dd163574340e16ced959e6ae Mon Sep 17 00:00:00 2001 From: Anh Nguyet Vu Date: Thu, 13 Feb 2025 18:23:59 -0800 Subject: [PATCH] Update source --- R/cbioportal.R | 3 +++ R/cboilerplate.R | 26 ++++++++++++++++++++++++++ 2 files changed, 29 insertions(+) diff --git a/R/cbioportal.R b/R/cbioportal.R index 7300f4d..0b8694f 100644 --- a/R/cbioportal.R +++ b/R/cbioportal.R @@ -188,6 +188,9 @@ cbp_add_maf <- function(maf_data, verbose = TRUE) { if(verbose) checked_message("Making maf meta file") make_meta_maf(cancer_study_identifier, verbose = verbose) + + if(verbose) checked_message("Making required _sequenced case list for mutation data") + make_case_list_maf(cancer_study_identifier) if(verbose) checked_message("Done with adding MAF data") diff --git a/R/cboilerplate.R b/R/cboilerplate.R index 6e326a5..969441d 100644 --- a/R/cboilerplate.R +++ b/R/cboilerplate.R @@ -411,6 +411,31 @@ make_meta_study_generic <- function(cancer_study_identifier, return(meta) } +#--- Generating case list files ------------------------------------------------ # + +#' Case lists for mutation samples +#' +#' https://docs.cbioportal.org/file-formats/#case-lists +#' @keywords internal +make_case_list_maf <- function(cancer_study_identifier, verbose = TRUE) { + + mut <- fread("data_mutations.txt") + mut_samples <- unique(mut$Tumor_Sample_Barcode) + n <- length(mut_samples) + case_list_ids <- paste(mut_samples,collapse = "\t") + meta <- glue::glue("cancer_study_identifier: {cancer_study_identifier}") %>% + append_kv("stable_id", paste0(cancer_study_identifier, "_sequenced")) %>% + append_kv("case_list_name", "Samples with mutation data from sequencing") %>% + append_kv("case_list_description", paste0("Samples with mutation data from sequencing ", "(", n, ")")) %>% + append_kv("case_list_ids", case_list_ids) + + if(!dir.exists("case_lists")) { + if(verbose) checked_message(glue::glue("Creating case_lists study directory")) + dir.create(glue::glue("./case_lists")) + } + + writeLines(meta, "case_lists/case-list.txt") +} # --- Other utils -------------------------------------------------------------- # @@ -443,3 +468,4 @@ use_ref_map <- function(ref_map, as_dt = TRUE) { return(ref_map_ls) } } +