R/get_g_level_annotations.R

In alextidd/tgp: Target Gene Prediction

get_g_level_annotations <- function(vxt_master,
                                    enriched){

  g <- list()

  # gene TPM signal/spec bins
  g <- enriched$expression %>%
    purrr::map( ~ .x %>% dplyr::inner_join(vxt_master %>% dplyr::distinct(ensg), by = "ensg"))
  names(g) <- paste0("expression_", names(g))

  # gene expressed
  g$expressed <- enriched$expressed %>%
    dplyr::inner_join(vxt_master %>% dplyr::distinct(ensg), by = "ensg")

  # return
  names(g) <- paste0("g_", names(g))
  return(g)
}