如何使用 R 搜索 PubMed 或其他数据库

我最近一直在使用优秀的rplos package https://github.com/ropensci/rplos，这使得搜索公共科学图书馆 (PLOS) API 上托管的论文变得非常容易。我遇到了一个障碍，因为 API 本身似乎缺少一些信息 - 其中一个主要问题是至少有 2012 篇关于 API 的论文在“期刊”字段中没有任何信息。我有每篇论文的 DOI，因此很容易通过 Google 搜索 DOI 并显示这些是在真实期刊（通常是 PLoS ONE）上发表的真实论文。显然，重复 2000 次是很愚蠢的。

我想知道是否有人知道如何找到来源期刊，如果我有 DOI 列表？我调查了RISmed 包 http://cran.r-project.org/web/packages/RISmed/index.html，它显然可以从 R 中搜索 PubMed，但我无法弄清楚如何让它提供有用的信息（只是搜索命中的数量，以及一些可能导致我想要的信息的 PubMed ID）。

有人知道如何将 DOI 列表转换为来源期刊名称吗？

编辑：我只是想到了另一个简单的解决方案。 DOI 包含期刊名称的缩写，对于像这样只有少数期刊的情况，只需使用正则表达式来读取 DOI 并选择它们来自哪个期刊。示例：10.1371/期刊。pone.0046711 来自 PLoS ONE。

这是基于 Thomas 建议尝试 rpubmed 的答案。它从有问题的 DOI 列表开始，使用 RISmed 中的 EUtilsSummary 函数查找匹配的 PubMed ID 号，然后使用从rpubmed 的 Github https://github.com/rOpenHealth/rpubmed/blob/master/R/rpubmed_fetch.R并转载如下。抱歉编辑 rpubmed 代码，但第 44 行的对象似乎没有定义或必不可少，所以我把它们删除了。

library(RCurl); library(XML); library(RISmed); library(multicore)

# dummy list of 5 DOIs. I actually have 2012, hence all the multicoring below
dois <- c("10.1371/journal.pone.0046711", "10.1371/journal.pone.0046681", "10.1371/journal.pone.0046643", "10.1371/journal.pone.0041465", "10.1371/journal.pone.0044562")

# Get the PubMed IDs
res <- mclapply(1:length(dois), function(x) EUtilsSummary(dois[x]))
ids<-sapply(res,QueryId)


######## rpubmed functions from https://github.com/rOpenHealth/rpubmed/blob/master/R/rpubmed_fetch.R
fetch_in_chunks <- function(ids, chunk_size = 500, delay = 0, ...){
  Sys.sleep(delay * 3600) # Wait for appropriate time for the server.
  chunks <- chunker(ids, chunk_size)
  Reduce(append, lapply(chunks, function(x) pubmed_fetch(x, ...)))
}

pubmed_fetch <- function(ids, file_format = "xml", as_r_object = TRUE, ...){

  args <- c(id = paste(ids, collapse = ","), db = "pubmed", rettype = file_format, ...)

  url_args <- paste(paste(names(args), args, sep="="), collapse = "&")
  base_url <- "http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?retmode=full"
  url_string <- paste(base_url, url_args, sep = "&")
  records <- getURL(url_string)
  #NCBI limits requests to three per second
  Sys.sleep(0.33)
  if(as_r_object){
    return(xmlToList(xmlTreeParse(records, useInternalNodes = TRUE)))
  } else return(records)
}

chunker <- function(v, chunk_size){
  split(v, ceiling(seq_along(v)/chunk_size))
}
###### End of rpubmed functions

d<-fetch_in_chunks(ids)
j<-character(0)
for(i in 1:2012) j[i]<-as.character(d[[i]][[1]][[5]][[1]][[3]]) # the tortuous path to the journal name