|
a |
|
b/scripts/ashm_hubs.R |
|
|
1 |
### build track hubs for public SSMs data from aSHM reagions, pathologies BL, |
|
|
2 |
### DLBCL and FL, seq type genome and capture, projection grch37 and hg38 |
|
|
3 |
|
|
|
4 |
library(GAMBLR.utils) |
|
|
5 |
library(GAMBLR.data) |
|
|
6 |
library(dplyr) |
|
|
7 |
|
|
|
8 |
# get metadata |
|
|
9 |
my_meta = get_gambl_metadata(seq_type_filter = c("genome", "capture")) %>% |
|
|
10 |
filter(pathology %in% c("BL", "DLBCL", "FL")) |
|
|
11 |
|
|
|
12 |
# create hub for grch37 |
|
|
13 |
build_browser_hub( |
|
|
14 |
regions_bed = GAMBLR.data::grch37_ashm_regions, |
|
|
15 |
these_samples_metadata = my_meta, |
|
|
16 |
these_seq_types = c("genome", "capture"), |
|
|
17 |
projection = "grch37", |
|
|
18 |
local_web_host_dir = "~/repos/LLMPP", |
|
|
19 |
hub_dir = "hubs/ashm_hubs", |
|
|
20 |
splitColumnName = "pathology", |
|
|
21 |
hub_name = "gamblr_ashm", |
|
|
22 |
shortLabel = "gamblr ashm", |
|
|
23 |
longLabel = "GAMBLR public aSHM mutations from grch37 projection", |
|
|
24 |
contact_email = "rdmorin@sfu.ca" |
|
|
25 |
) |
|
|
26 |
|
|
|
27 |
# create hub for hg38 |
|
|
28 |
build_browser_hub( |
|
|
29 |
regions_bed = GAMBLR.data::hg38_ashm_regions, |
|
|
30 |
these_samples_metadata = my_meta, |
|
|
31 |
these_seq_types = c("genome", "capture"), |
|
|
32 |
projection = "hg38", |
|
|
33 |
local_web_host_dir = "~/repos/LLMPP", |
|
|
34 |
hub_dir = "hubs/ashm_hubs", |
|
|
35 |
splitColumnName = "pathology", |
|
|
36 |
hub_name = "gamblr_ashm", |
|
|
37 |
shortLabel = "gamblr ashm", |
|
|
38 |
longLabel = "GAMBLR public aSHM mutations from hg38 projection", |
|
|
39 |
contact_email = "rdmorin@sfu.ca" |
|
|
40 |
) |
|
|
41 |
|