plotCoverageFromEnsembldb: Plot read coverage directly from ensembldb object.

Description Usage Arguments Value Examples

Description

A wrapper around the plotCoverage function. See the documentation for (plotCoverage) for more information.

Usage

1
2
plotCoverageFromEnsembldb(ensembldb, gene_names, transcript_ids = NULL,
  ...)

Arguments

ensembldb

ensembldb object.

gene_names

List of gene names to be plotted.

transcript_ids

Optional list of transcript ids to be plotted.

...

Additional parameters to be passed to plotCoverage.

Value

ggplot2 object

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
require("EnsDb.Hsapiens.v86")
require("dplyr")
require("GenomicRanges")
sample_data = dplyr::data_frame(sample_id = c("aipt_A", "aipt_C", "bima_A", "bima_C"), 
 condition = factor(c("Naive", "LPS", "Naive", "LPS"), levels = c("Naive", "LPS")), 
 scaling_factor = 1) %>%
 dplyr::mutate(bigWig = system.file("extdata",  paste0(sample_id, ".str2.bw"), package = "wiggleplotr"))
 
track_data = dplyr::mutate(sample_data, track_id = condition, colour_group = condition)
## Not run: 
plotCoverageFromEnsembldb(EnsDb.Hsapiens.v86, "NCOA7", transcript_ids = c("ENST00000438495", "ENST00000392477"), 
track_data, heights = c(2,1), fill_palette = getGenotypePalette())

## End(Not run)

kauralasoo/wiggleplotr documentation built on July 21, 2021, 4:15 a.m.