Is it possible to add additional gene information? about gwasrapidd HOT 4 CLOSED

library(dplyr, warn.conflicts = FALSE)
library(glue, warn.conflicts = FALSE)
library(httr, warn.conflicts = FALSE)
library(tidyjson, warn.conflicts = FALSE)
library(purrr)
library(gwasrapidd, warn.conflicts = FALSE)

# For retrieving info about one gene only.
get_gene <- function(gene_name) {
  
  gene_name <- URLencode(gene_name)
  url <- glue::glue("https://www.ebi.ac.uk/gwas/api/search?q=title%3A+%22{gene_name}%22+AND+resourcename%3Agene")
  response <- httr::GET(url)
  json <- httr::content(response, 'text', encoding = 'UTF-8')
  
  tbl <-
    json %>%
    tidyjson::enter_object('response', 'docs') %>%
    tidyjson::gather_array() %>%
    tidyjson::spread_values(
      gene_name = tidyjson::jstring('title'),
      description = tidyjson::jstring('description'),
      chromosome_name = tidyjson::jstring('chromosomeName'),
      chromosome_start = tidyjson::jnumber('chromosomeStart'),
      chromosome_end = tidyjson::jnumber('chromosomeEnd'),
      cytogenetic_band = tidyjson::jstring('cytobands'),
      biotype = tidyjson::jstring('biotype')
    ) %>%
    tidyjson::as_tibble() %>%
    dplyr::select(-c('document.id', 'array.index'))
  
  return(tbl)
}

# Maps over genes
get_genes <- function(gene_names) {
  purrr::map_dfr(gene_names, get_gene)
}

# As an example, let us get some studies about rheumatoid arthritis.
my_studies <- get_studies(efo_trait = 'rheumatoid arthritis')

# Going after those associations related to the first study found only (for the
# sake of simplicity)
one_study <- my_studies[1]
my_associations <- get_associations(study_id = one_study@studies$study_id)
gene_info <- get_genes(my_associations@genes$gene_name)

my_associations@genes <- dplyr::left_join(my_associations@genes, gene_info, by = 'gene_name')

# The slot `genes` of the `associations` object is now augmented to include
# those extra details about each gene.
my_associations@genes
#> # A tibble: 5 × 9
#>   association_id locus_id gene_name description chromosome_name chromosome_start
#>   <chr>             <int> <chr>     <chr>       <chr>                      <dbl>
#> 1 17729                 1 BLK       BLK proto-… 8                       11486894
#> 2 17731                 1 ARHGEF3   Rho guanin… 3                       56727418
#> 3 17883                 1 HLA-DRB1  major hist… 6                       32578769
#> 4 17730                 1 TRHDE     thyrotropi… 12                      72087266
#> 5 40454                 1 PADI4     peptidyl a… 1                       17308195
#> # … with 3 more variables: chromosome_end <dbl>, cytogenetic_band <chr>,
#> #   biotype <chr>