stat_aggregate {ggbio}R Documentation

Generates summaries on the specified windows

Description

Generates summaries on the specified windows

Usage


## S4 method for signature 'GRanges'
stat_aggregate(data, ..., xlab, ylab, main, by, FUN, start = NULL,
                 end = NULL, width = NULL, y = NULL, frequency = NULL,
                 delta = NULL, simplify = TRUE, window = NULL, facets =
                 NULL, type = c("mean", "median", "max", "min", "sum",
                 "count", "identity"), geom = NULL)

Arguments

data

A GRanges or data.frame object.

xlab

Label for x

ylab

Label for y

main

Title for plot.

by

An object with 'start', 'end', and 'width' methods. Passed to aggreagate.

FUN

The function, found via 'match.fun', to be applied to each window of 'x'. Passed to aggreagate.

start

Start of the window. If 'by' is missing, then must supply two of the 'start', 'end', 'width'. If 'window' is provided then you don't have to specify it.

end

End of the window. If 'by' is missing, then must supply two of the 'start', 'end', 'width'. If 'window' is provided then you don't have to specify it.

width

Width of the window. If 'by' is missing, then must supply two of the 'start', 'end', 'width'. If 'window' is provided then you don't have to specify it.

y

A character indicate the varialbe column for which aggregation is taken on. Notice for geom like 'boxplot', we don't compute or aggregate the variable, we simply want to use the idenitcal y as y axis, in that case please put y in the aes mapping function.

frequency

Optional arguments that specify the sampling frequency within the window.

delta

Optional arguments that specify the sampling increment within the window.

...

Arguments passed to plot function. such as aes() and color.

simplify

A logical value specifying whether or not the result should be simplified to a vector or matrix if possible.

window

Integer value indicate window size.

facets

Faceting formula to use.

type
geom

The geometric object to use display the data.

Value

A 'Layer'.

Author(s)

Tengfei Yin

Examples

library(GenomicRanges)
set.seed(1)
N <- 1000
## ======================================================================
##  simmulated GRanges
## ======================================================================
gr <- GRanges(seqnames = 
              sample(c("chr1", "chr2", "chr3"),
                     size = N, replace = TRUE),
              IRanges(
                      start = sample(1:300, size = N, replace = TRUE),
                      width = sample(70:75, size = N,replace = TRUE)),
              strand = sample(c("+", "-", "*"), size = N, 
                replace = TRUE),
              value = rnorm(N, 10, 3), score = rnorm(N, 100, 30),
              sample = sample(c("Normal", "Tumor"), 
                size = N, replace = TRUE),
              pair = sample(letters, size = N, 
                replace = TRUE))
## 36,37 doesn't work
## ggplot(gr) + stat_aggregate(y = "value",fill = "gray40")
## ggplot(gr) + stat_aggregate(window = 30,  y = "value",fill = "gray40", geom = "histogram")
## ggplot(gr) + stat_aggregate(window = 100, fill = "gray40",y = "value",
##                           type = "max", geom = "histogram")

[Package ggbio version 1.5.20 ]