cr_enriched_heatmap_at_genomic_features: Visualizing enrichment for epigenomic signals at TSS-CGIs

Description Usage Arguments Details Value Author(s) Examples

View source: R/correlated_regions_enrichedheatmap.R

Description

Visualizing enrichment for epigenomic signals at TSS-CGIs

Usage

1
2
3
cr_enriched_heatmap_at_genomic_features(cr, txdb, expr, gf,
    fdr_cutoff = 0.05, meth_diff_cutoff = 0.1, marks = NULL, type = "neg", extend = 5000,
    min_reduce = 1, min_width = 1000, nearest_by = "tss", expr_ha)

Arguments

cr

correalted regions

txdb

transcriptome annotation which was used in correlated_regions

expr

expression matrix which was used in correlated_regions

gf

genomic features, a GRanges object

fdr_cutoff

cutoff for fdr

meth_diff_cutoff

cutoff for methylation difference. If there are no subgroup information or only one subgroup, meth_IQR column is used for filtering. If there are more than one subgroups, meth_diameter column is used for filtering.

marks

names of histone marks, should be supported in chipseq_hooks

type

visualize negative correlated regions or positive correlated regions

extend

base pairs extended to upstream and downstream

min_reduce

base pairs for merging neighbouring regions

min_width

minimal width of regions

nearest_by

"tss" or "gene", how to connect genomic features to genes

expr_ha

a HeatmapAnnotation class object.It is used for the expression heatmap

Details

There are several heatmaps visualize various signals enriched at genomic features. After annotate to genes, in the extended regions, each region can only have one gene.

There are following heatmaps:

If there are more than 12 heatmaps, they will be put into two pages.

Heatmaps are split into two sub-clusters by k-means clustering on the mean methylation matrix. If there are two subgroups in all samples, each subcluster are split by high expression/low expression in subgroup 1. In each high expression/low expression, rows are split by the k-means clusters calculated in cr_enriched_heatmap. Finally, rows are clustered by considering closeness of signals in the extended gf regions.

Value

no value is returned

Author(s)

Zuguang Gu <z.gu@dkfz.de>

Examples

1
2
# There is no example
NULL

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