Download T1D summary statistics

dl_T1D_sumstats.R

urls <- gsubfn::strapply(
  readLines("https://datadryad.org//resource/doi:10.5061/dryad.ns8q3"),
  "<a href=\"(/bitstream/handle/10255/dryad\\.[0-9]+/meta_chr_[0-9]+\\?sequence=1)\">",
  simplify = 'c')

sumstats <- purrr::map_dfr(urls, ~ {
  download.file(paste0("https://datadryad.org", .x),
                destfile = (tmp <- tempfile(fileext = ".txt")))
  sumstats <- bigreadr::fread2(
    tmp, select = c("chromosome", "position", "a0", "a1",
                    "beta.meta", "p.meta", "info_score.A", "info_score.I"))
  na.omit(sumstats)
})