Jam Sashimi plot

Jam Sashimi plot

Usage

plotSashimi(
  sashimi,
  show = c("coverage", "junction", "junctionLabels"),
  coord_method = c("scale", "coord", "none"),
  exonsGrl = NULL,
  junc_color = jamba::alpha2col("goldenrod2", 0.3),
  junc_fill = jamba::alpha2col("goldenrod2", 0.9),
  junc_alpha = 0.8,
  junc_accuracy = 1,
  junc_nudge_pct = 0.05,
  junc_fontsize = 12,
  fill_scheme = c("sample_id", "exon"),
  color_sub = NULL,
  ylabel = "read depth",
  xlabel = NULL,
  xlabel_ref = TRUE,
  use_jam_themes = TRUE,
  apply_facet = TRUE,
  facet_scales = "free_y",
  ref2c = NULL,
  label_coords = NULL,
  do_highlight = FALSE,
  verbose = FALSE,
  ...
)

Arguments

sashimi

Sashimi data prepared by prepareSashimi() which is a list with covDF coverage data in data.frame format, juncDF junction data in data.frame format, juncLabelDF junction label coordinates in data.frame format, exonLabelDF exon label coordinates per coverage polygon in data.frame format, ref2c list output from make_ref2compressed() to transform genomic coordinates.

show

character vector of Sashimi plot features to include: "coverage" sequence read coverage data; "junction" splice junction read data.

coord_method

character value indicating the type of coordinate scaling to use: "scale" uses ggplot2::scale_x_continuous(); "coord" uses ggplot2::coord_transform(); "none" does not compress genomic coordinates.

exonsGrl

GRangesList object with one or more gene or transcript exon models, where exons are disjoint (not overlapping.)

junc_color, junc_fill

character string with valid R color, used for junction outline, and fill, for the junction arc polygon. Alpha transparency is recommended for junc_fill so overlapping junction arcs are visible.

junc_alpha

numeric value between 0 and 1, to define the alpha transparency used for junction colors, where 0 is fully transparent, and 1 is completely non-transparent.

junc_nudge_pct

numeric value to nudge junction labels by the percent of the maximum y-axis junction label position.

The default junc_nudge_pct=0.05 nudges labels up by 5%, which makes them consistently appear above the top edge of the junction ribbon.
Negative values will place labels just below the top edge of the junction ribbon.
A vector can be supplied to nudge each junction label individually, applied in order the labels appear from left to right.
Note that when the distance from label to the top edge of the ribbon exceeds a threshold, a line segment is drawn from the label to the top edge of the junction ribbon. This threshold is controlled by ggrepel::geom_text_repel(..., min.segment.length=0.5) and is not configurable at this time.

junc_fontsize

numeric default 12, the text font size used when show includes 'junctionLabels' as is default.

fill_scheme

character string for how the exon coverages will be color-filled: "exon" will define colors for each distinct exon, using the GRanges names from flatExonsByGene; "sample_id" to color all exons the same by sample_id.

color_sub

optional character vector of R compatible colors or hex strings, whose names are used to color or fill features in the ggplot object. For example, if fill_sheme="sample_id" the color_sub should have names for each "sample_id" value. If any values are missing, they will be filled in using colorjam::rainbowJam().

ylabel

character string used as the y-axis label, by default "score" reflects the coverage score and junction score, respectively for coverage and junction data. Scores are also adjusted using the scale_factor value for each sample_id as defined in the filesDF. Set to NULL to hide the y-axis label completely.

xlabel

character string used to define the x-axis name, which takes priority over argument xlabel_ref. When xlabel is NULL and xlabel_ref is FALSE, then the x-axis name is "", which displays no x-axis label.

xlabel_ref

logical indicating whether the x-axis name should be determined by the reference (chromosome).

use_jam_themes

logical indicating whether to apply colorjam::theme_jam(), by default for the ggplot theme.

apply_facet

logical indicating whether to apply ggplot2::facet_wrap() with "~sample_id" defining each panel.

facet_scales

character value used as "scales" argument in ggplot2::facet_wrap() when apply_facet=TRUE.

ref2c

optional list output from make_ref2compressed() to compress axis coordinates during junction arc calculations.

label_coords

numeric vector length 2, optional range of genomic coordinates to restrict labels, so labels are not arranged by ggrepel::geom_text_repel() even when coord_cartesian() is used to zoom into a specific x-axis range.

verbose

logical indicating whether to print verbose output.

...

additional arguments are sent to grl2df().

Details

This function uses Sashimi data prepared by prepareSashimi() and creates a ggplot graphical object ready for visualization. As a result, this function provides several arguments to customize the visualization.

Examples

suppressPackageStartupMessages(library(GenomicRanges));
data(test_exon_gr);
data(test_junc_gr);
data(test_cov_gr);
filesDF <- data.frame(url="sample_A",
   type="coverage_gr",
   sample_id="sample_A");
sh1 <- prepareSashimi(
   flatExonsByGene=GRangesList(TestGene1=test_exon_gr),
   filesDF=filesDF,
   gene="TestGene1",
   covGR=test_cov_gr,
   juncGR=test_junc_gr);
plotSashimi(sh1);

Usage

Arguments

Details

See also

Examples