I wonder if someone can reproduce/check if this is a bug.
genome_ssm_data = get_coding_ssm()
ssm_data_capture = get_coding_ssm(this_seq_type = "capture")
all_ssm = bind_rows(ssm_data_capture,genome_ssm_data) #combine all ssm into one large data frame
#Problem:
#this makes a plot with every cohort in the legend even if no samples from the cohort exist in the MAF
prettyOncoplot(maf_df = all_ssm,these_samples_metadata =
all_sample_meta,
minMutationPercent = 5,metadataColumns = c("pathology","cohort"),
sortByColumns = c("pathology","cohort"))
# Expected functionality (requires an explicit filter)
# this only shows the relevant cohorts in the legend (i.e. those in the oncoplot)
prettyOncoplot(maf_df = all_ssm,these_samples_metadata =
filter(all_sample_meta,sample_id %in% all_ssm$Tumor_Sample_Barcode),
minMutationPercent = 5,metadataColumns = c("pathology","cohort"),
sortByColumns = c("pathology","cohort"))```
[Slack Message](https://morinlabsfu.slack.com/archives/C0224H120CU/p1706573886956999?thread_ts=1706573886.956999&cid=C0224H120CU)
I wonder if someone can reproduce/check if this is a bug.