metaGene: Returns a histogram of the number of reads in each section of...

In coregenomics/groHMM: GRO-seq Analysis Pipeline

Description

Returns a histogram of the number of reads in each section of a moving window centered on a certain feature.

Usage

metaGene(features, reads = NULL, plusCVG = NULL, minusCVG = NULL,
  size = 100L, up = 10000L, down = NULL, BPPARAM = bpparam())

Arguments

features	A GRanges object representing a set of genomic coordinates. The meta-plot will be centered on the transcription start site (TSS)
reads	A GRanges object representing a set of mapped reads. Instead of 'reads', 'plusCVG' and 'minusCVG' can be used Default: NULL
plusCVG	A RangesList object for reads with '+' strand.
minusCVG	A RangesList object for reads with '-' strand.
size	The size of the moving window.
up	Distance upstream of each features to align and histogram. Default: 10 kb.
down	Distance downstream of each features to align and histogram. If NULL, same as up. Default: NULL.
BPPARAM	Registered backend for BiocParallel. Default: BiocParallel::bpparam()

Value

Returns a integer-Rle representing the 'typical' signal centered on a point of interest.

Author(s)

Charles G. Danko and Minho Chae

Examples

library(GenomicRanges)
features <- GRanges("chr7", IRanges(1000, 1000), strand="+")
reads <- GRanges("chr7", IRanges(start=c(1000:1004, 1100),
 width=rep(1, 6)), strand="+")
mg <- metaGene(features, reads, size=4, up=10)

coregenomics/groHMM documentation built on May 7, 2019, 7:57 a.m.