R/get_d_level_annotations.R

In alextidd/tgp: Target Gene Prediction

get_d_level_annotations <- function(variants,
                                    DHSs){
  cat("Annotating DHSs...\n")

  d <- list()

  d$inv_n_variants <- variants %>%
    bed_intersect_left(DHSs, keepBcoords = F) %>%
    dplyr::group_by(DHS) %>%
    dplyr::summarise(value = 1/dplyr::n_distinct(variant))

  # return
  names(d) <- paste0("d_", names(d))
  return(d)
}