library(httr)
library(jsonlite)
library(xml2)
server <- "http://grch37.rest.ensembl.org"
for (m in exposure$SNP) {
path1 = paste0("/variation/human/", m, '?pops=1')
ext <- path1
res <- GET(paste(server, ext, sep = ""), content_type("application/json"))
stop_for_status(res)
res <- httr::content(res)
res <- jsonlite::fromJSON(toJSON(res))$populations
s = res[res$population == '1000GENOMES:phase_3:EUR',]
s_res = s[s$allele == exposure[exposure$SNP == m,"effect_allele.exposure"],]
exposure[exposure$SNP == m, ]$eaf.exposure = s_res$frequency
print(s_res)
}
SNP查找MAF代码
于 2023-09-06 15:46:07 首次发布