create_enh_region: Create enhancer regions from gene annotation data.

Description Usage Arguments Value Author(s) Examples

Description

create_enh_region creates enhancer region from gene annotation data.

Usage

1
2
create_enh_region(annot_data, chrom_size = NULL, upstream = -15000,
  downstream = 15000)

Arguments

annot_data

A GRanges object containing the enhancer annotation data.

chrom_size

Optional data.table containing chromosome sizes, e.g. using the read_chrom_size function.

upstream

Integer defining the length of bp upstream of enhancer centre.

downstream

Integer defining the length of bp downstream of enhancer cetre.

Value

A GRanges object containing the enhancer regions data.

Author(s)

C.A.Kapourani C.A.Kapourani@ed.ac.uk

Examples

1
2
3
4
5
enh_beatson_file <- system.file("extdata", "enhancers_beatson.bed", package = "processHTS")
gene_beatson_file <- system.file("extdata", "rna_beatson2.bed", package = "processHTS")
enh_beatson_data <- read_enh_beatson(enh_beatson_file, gene_beatson_file)
enh_beatson_data <- enh_beatson_data$enhancer_data
enh_regions <- create_enh_region(enh_beatson_data)

andreaskapou/processHTS documentation built on May 12, 2019, 3:33 a.m.