inst/analysis/2_build-lists/051_ALSoD-gene.R

In soderling-lab/geneLists: tools for mapping gene identifiers and a collection of gene lists in a universal format

#!/usr/bin/env Rscript

## Input Parameters:
script <- "051_ALSoD-gene"
short_name <- "alsodgene"
tags <- c("ALS")
ref_url <- "https://alsod.ac.uk/"
data_url <- "https://alsod.ac.uk/"

# Data were  manually cut-and-paste from the website into excel.
data_file <- "ALSod.csv" # in root/downloads

# Load renv.
root <- getrd()
renv::load(root, quiet = TRUE)

# Imports.
suppressPackageStartupMessages({
  library(dplyr)
  library(getPPIs)
  library(data.table)
})

# Load functions in root/R.
suppressMessages({
  devtools::load_all()
})

# Directories.
Rdir <- file.path(root, "R")
gmtdir <- file.path(root, "gmt")
datadir <- file.path(root, "data")
tabsdir <- file.path(root, "tables")
downdir <- file.path(root, "downloads")

# Load the data.
myfile <- file.path(downdir, data_file)
data <- fread(myfile)

# Map human genes to entrez.
data$Entrez <- getIDs(data$"Gene symbol",
  from = "Symbol", to = "Entrez", species = "human"
)

# Map to mouse homologs.
data$msEntrez <- getHomologs(data$Entrez, species = "mouse")

# Drop un-mapped genes.
mygenes <- unlist(data[!is.na(data$msEntrez), "msEntrez"]) %>% unique()
gene_list <- list("ALS" = mygenes)

# Status.
message(paste(
  "Compiled", length(mygenes), "genes associated with",
  "human ALS from ALSoD.ac.uk"
))

# Save as gmt, and then save as rda and generate documentation.
myfile <- file.path(gmtdir, script)
write_gmt(gene_list, ref_url, myfile)
documentDataset(myfile, short_name, Rdir, datadir)