annotate_to_gene_models: Annotate genomic regions to gene models

In jokergoo/epik: Integrative Analysis and Visualization of Epigenomic Sequencing Data

Description

Annotate genomic regions to gene models

Usage

annotate_to_gene_models(gr, txdb, gene_model =c("gene", "tx"),
    species = "hg19", promoters_upstream = 2000, promoters_downstream = 200,
    annotation_type = c("percent", "number"),
    annotation_prefix = "overlap_to_")

Arguments

gr	a GRanges object
txdb	a TxDb-class object.
gene_model	type of gene model, by transcript or by genes
species	We need this information to find out proper intergenic regions
promoters_upstream	length of upstream promoter from TSS, pass to promoters
promoters_downstream	length of downstream promoter from TSS, pass o promoters
annotation_type	pass to annotate_to_genomic_features
annotation_prefix	pass to annotate_to_genomic_features

Value

Following columns are attached to gr:

nearest_*_tss

the nearest tss (depending on gene_model)

dist_to_*_tss

distance to the closest tss (depending on gene_model)

nearest_*

the closest gene model (depending on gene_model)

dist_to_*

distance to the closest gene model (depending on gene_model)

'prefix_to'_exon

percent of the region which is covered by exons or number of exons overlapped to the region

'prefix_to'_intron

percent of the region which is covered by introns or number of introns overlapped to the region

'prefix_to'_promoter

percent of the region which is covered by promoters or number of promoters overlapped to the region

'prefix_to'_intergenic

percent of the region which is covered by intergenic regions or number of intergenic regions overlapped to the region

'prefix_to'_fiveUTR

percent of the region which is covered by 5'UTRs or number of 5'UTRs overlapped to the region

'prefix_to'_threeUTR

percent of the region which is covered by 3'UTRs or number of 3'UTRs overlapped to the region

Author(s)

Zuguang Gu <[e-mail]>

Examples

# There is no example
NULL

jokergoo/epik documentation built on Sept. 28, 2019, 9:20 a.m.