@@ -60,7 +60,8 @@ output_path <- here::here('reports')
6060pinvars <-
6161 read_tsv(pinpoints_path,
6262 col_types = cols(.default = col_character())) %>%
63- type_convert(guess_integer = TRUE)
63+ type_convert(guess_integer = TRUE) %>%
64+ mutate(row_label = as.integer(row_label))
6465
6566poolvars <-
6667 read_tsv(all_variants_path,
@@ -71,6 +72,7 @@ poolvars <-
7172``` {r eval=have(params$rescue)}
7273rescue <-
7374 read_tsv(rescue_path) %>%
75+ mutate(row_label = as.integer(row_label)) %>%
7476 select(-sample_id)
7577
7678pinvars <-
@@ -106,43 +108,6 @@ poolvars <-
106108 left_join(annotations, by='varid')
107109```
108110
109-
110- ``` {r eval=F}
111- annotations <-
112- read_tsv(annotations_path,
113- na = '-', col_types = cols(.default = col_character())) %>%
114- type_convert(guess_integer = TRUE) %>%
115- mutate(is_lof = if_else(LoF == 'HC',1,0)) %>%
116- mutate(clinvar_stars = case_when(
117- str_detect(ClinVar_CLNREVSTAT, "practice_guideline") ~ 4,
118- str_detect(ClinVar_CLNREVSTAT, "reviewed_by_expert_panel") ~ 3,
119- str_detect(ClinVar_CLNREVSTAT, "criteria_provided,_multiple_submitters,_no_conflicts") ~ 2,
120- str_detect(ClinVar_CLNREVSTAT, "criteria_provided,_conflicting_interpretations") ~ 1,
121- str_detect(ClinVar_CLNREVSTAT, "criteria_provided,_single_submitter") ~ 1,
122- TRUE ~ NA_real_
123- )) %>%
124- mutate(is_p = if_else(str_detect(ClinVar_CLNSIG, 'Pathogenic') | str_detect(ClinVar_CLNSIG, "Likely_pathogenic"), 1,0)) %>%
125- mutate(is_lofp = if_else(is_lof == 1 | is_p == 1,1,0))
126-
127- rescue <-
128- read_tsv(rescue_path) %>%
129- select(-sample_id)
130-
131- pinvars <-
132- read_tsv(pinpoints_path,
133- col_types = cols(.default = col_character())) %>%
134- type_convert(guess_integer = TRUE) %>%
135- left_join(rescue, by=c('uvarid','row_id','column_id','varid')) %>%
136- left_join(annotations, by='varid')
137-
138- poolvars <-
139- read_tsv(all_variants_path,
140- col_types = cols(.default = col_character())) %>%
141- type_convert(guess_integer = TRUE) %>%
142- left_join(annotations, by='varid')
143-
144- ```
145-
146111``` {r}
147112matrix_ctx <- fromJSON(matrix_context_path)
148113matrix_size <- matrix_ctx$matrix$n_rows
0 commit comments