diff --git a/analyses/post_clustering_associations/01-compare-classes.R b/analyses/post_clustering_associations/01-compare-classes.R index b2ede8a..9e8615a 100644 --- a/analyses/post_clustering_associations/01-compare-classes.R +++ b/analyses/post_clustering_associations/01-compare-classes.R @@ -23,6 +23,7 @@ dir.create(output_dir, showWarnings = F, recursive = T) mm_clusters <- read_tsv(file = opt$cluster_file) # combine Multi-modal clusters with RNA-derived molecular subtypes +cat('Comparing multi-omic clusters to known subtypes') anno_file_rna <- read_tsv(file = opt$histology_file) %>% dplyr::select(Kids_First_Biospecimen_ID, molecular_subtype, CNS_region) %>% dplyr::rename("Kids_First_Biospecimen_ID_RNA" = "Kids_First_Biospecimen_ID") %>%