-
Notifications
You must be signed in to change notification settings - Fork 0
Open
Description
Hello @dami82,
I am using this package for the first time to analyze mutation signatures, and I would like to express my gratitude as I find it extremely useful.
I have two questions:
I created a plot for 30 COSMIC mutation signatures using the code below and identified the six most frequently occurring mutation signatures across all samples.
Now, I would like to generate a plot only for these six mutation signatures.
How can I do this?
...
# 6. Estimate Exposures to known Mutational Signatures
cosmix <- getCosmicSignatures(forceUseMirror = FALSE, asMutSign = TRUE)
A_matrix <- dcast(merge_vcf, mutType ~ SAMPLEID, fun.aggregate = length, value.var = "SAMPLEID")
View(A_matrix)
A_matrix <- as.data.frame(A_matrix)
rownames(A_matrix) <- A_matrix$mutType
A_matrix$mutType <- NULL
xx <- as.mutation.counts(A_matrix)
xx
blca.expo1 <- resolveMutSignatures(mutCountData = xx,
signFreqData = cosmix)
blca.exp.1x <- blca.expo1$Results$count.result
msigPlot(blca.exp.1x, main = "COSMIC")
Thank you in advance for your help.
Metadata
Metadata
Assignees
Labels
No labels
