Diff of /scripts/ashm_hubs.R [000000] .. [fcb5af]

Switch to unified view

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