Nothing
context("test-plot_rainfall")
# Read data
vcfs <- readRDS(system.file("states/read_vcfs_as_granges_output.rds",
package = "MutationalPatterns"
))
# Specify chromosomes of interest.
chromosomes <- paste0("chr", c(1:22))
# Do a rainfall plot for all chromosomes:
output <- plot_rainfall(vcfs[[1]], title = names(vcfs[1]), chromosomes = chromosomes)
# Plot a single chromosome (chromosome 1):
output_singlechrom <- plot_rainfall(vcfs[[1]], title = names(vcfs[1]), chromosomes = chromosomes[1])
# Plot a subset of the variants
output_subset <- plot_rainfall(vcfs[[1]][1:10], title = names(vcfs[1]), chromosomes = chromosomes)
# Plot an empty gr
output_empty <- plot_rainfall(vcfs[[1]][0], title = names(vcfs[1]), chromosomes = chromosomes)
test_that("Output has correct class", {
expect_true(inherits(output, c("gg")))
expect_true(inherits(output_singlechrom, c("gg")))
expect_true(inherits(output_subset, c("gg")))
expect_true(inherits(output_empty, c("gg")))
})
test_that("Subsetted output contains the correct subset of colours", {
colours_used <- unique(ggplot_build(output_subset)$data[[1]][["colour"]])
expect_equal(colours_used, c("#ADCC54", "#DE1C14", "#2EBAED", "#D4D2D2"))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.