plot.kataegisList: Create a kataegis "rainfall" plot.

Description Usage Arguments Details Value Author(s) Examples

Description

Create a kataegis "rainfall" plot.

Usage

1
2
plot.kataegisList(x, chr = NA, start = NA, end = NA, all = FALSE,
  alpha = 0.5, legend = TRUE, colours = kataegis.colours, ...)

Arguments

x

a kataegisList object, usually created with the kataegis function.

chr

Optionally, only process this chromosome

start

Optionally, only process from this start point (must be specify together with chr and end).

end

Optionally, only process until this point (must be specify together with chr and start).

all

Plot all types of SNVs, default only C>X and T>X are plotted.

alpha

adjust the alpha level of the points

legend

If TRUE display a figure legend

colours

a function specifying the colours to use in the plot, see the default value (kataegis.colours) for an example.

...

passed to plot function

Details

In its simplest form (plot(x)), this function produces a rainfall plot on the current graphical device. By default, the whole genome is displayed. This can be changed by providing either the chr parameter (to show a specific chromosome) or the chr, start and end parameters to show a specific region of a chromosome (note that all three parameters must be provided).

Value

invisibly return data.frame of SNVs plotted.

Author(s)

Jonathan Ellis <jonathan.j.ellis@gmail.com>

Examples

1
2
3
4
5
6
7
## Not run: 
vcf.file <- system.file('extdata', 'example.vcf', package = 'GenomicVis')
x <- kategis(vcf.file)
main <- sprintf("Kataegis: %d", length(x$kataegis))
plot(x, main = main)

## End(Not run)

jjellis/GenomicVis documentation built on May 19, 2019, 11:39 a.m.