deepblue_enrich_regions_overlap: enrich_regions_overlap

Description Usage Arguments Value See Also Examples

View source: R/deepblue.R

Description

Enrich the regions based on regions overlap analysis.

Usage

1
2
deepblue_enrich_regions_overlap(query_id = NULL, background_query_id = NULL,
  datasets = NULL, genome = NULL, user_key = deepblue_options("user_key"))

Arguments

query_id

- A string (Query ID)

background_query_id

- A string (query_id containing the regions that will be used as the background data.)

datasets

- A struct (a map where each key is an identifier and the value is a list containing experiment names or query_ids (you can use both together).)

genome

- A string (the target genome)

user_key

- A string (users token key)

Value

request_id - A string (Request ID - Use it to retrieve the result with info() and get_request_data(). The result is a list containing the datasets that overlap with the query_id regions.)

See Also

Other Enrich the genome regions: deepblue_enrich_regions_fast, deepblue_enrich_regions_go_terms

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
query_id = deepblue_select_experiments(
  experiment_name="S00VEQA1.hypo_meth.bs_call.GRCh38.20150707.bed")

filtered_query_id = deepblue_filter_regions(
  query_id = query_id,
  field = "AVG_METHYL_LEVEL",
  operation = "<",
  value = "0.0025",
  type="number")

rg_10kb_tilling = deepblue_tiling_regions(
    size = 1000,
    genome = "hg19")

# We could have included more Epigenetic Marks here
epigenetic_marks <- c("h3k27ac", "H3K27me3", "H3K4me3")

histones_datasets = c()
for (i in 1:length(epigenetic_marks)) {
  experiments_list <- deepblue_list_experiments(
    epigenetic_mark=epigenetic_marks[[i]],
    type="peaks",
    genome="grch38",
    project="BLUEPRINT Epigenome");

    experiment_names = deepblue_extract_names(experiments_list)
    histones_datasets[[epigenetic_marks[[i]]]] = experiment_names
}

deepblue_enrich_regions_overlap(
  query_id=filtered_query_id,
  background_query=rg_10kb_tilling,
  datasets=histones_datasets,
  genome="grch38")

MPIIComputationalEpigenetics/DeepBlue-R documentation built on Aug. 11, 2021, 3:18 p.m.