Description Usage Arguments Value See Also Examples
View source: R/plot_compare_dbs.R
Plots two dbs mutation profiles and their difference, reports the residual sum of squares (RSS).
1 2 3 4 5 6 7 | plot_compare_dbs(
profile1,
profile2,
profile_names = c("profile 1", "profile 2"),
profile_ymax = 0.2,
diff_ylim = c(-0.1, 0.1)
)
|
profile1 |
First mutation profile |
profile2 |
Second mutation profile |
profile_names |
Character vector with names of the mutations profiles used for plotting, default = c("profile 1", "profile 2") |
profile_ymax |
Maximum value of y-axis (relative contribution) for profile plotting. This can only be used to increase the y axis. If bars fall outside this limit, the maximum value is automatically increased. default = 0.2. |
diff_ylim |
Y-axis limits for profile difference plot, default = c(-0.1, 0.1) |
A ggplot2 object
plot_compare_profiles
,
plot_compare_indels
,
plot_compare_mbs
Other DBS:
count_dbs_contexts()
,
get_dbs_context()
,
plot_dbs_contexts()
,
plot_main_dbs_contexts()
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 | ## Get the dbs counts
## See 'count_dbs_contexts()' for more info on how to do this.
dbs_counts <- readRDS(system.file("states/blood_dbs_counts.rds",
package = "MutationalPatterns"
))
## Get dbs refit info.
## See 'fit_to_signatures()' for more info on how to do this.
fit_res <- readRDS(system.file("states/dbs_refit.rds",
package = "MutationalPatterns"
))
## Compare the reconstructed profile of sample 1 with the original profile
## The same thing could be done with a reconstructed profile from NMF.
plot_compare_dbs(dbs_counts[, 1], fit_res$reconstructed[, 1])
## You could also compare regular mutation profiles with eachother.
plot_compare_dbs(
dbs_counts[, 1],
dbs_counts[, 2]
)
## Or change the names of the profiles
plot_compare_dbs(dbs_counts[, 1],
dbs_counts[, 2],
profile_names = c("Original", "Reconstructed")
)
## You can also change the y limits.
## This can be done separately for the profiles and the different facets.
plot_compare_dbs(dbs_counts[, 1],
dbs_counts[, 2],
profile_ymax = 0.3,
diff_ylim = c(-0.03, 0.03)
)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.