Вы можете сделать это, используя пакет rentrez .
Требуемая функция: entrez_link .
Пример:
library(rentrez)
results <- entrez_link(dbfrom = 'pubmed', id = 19122651, db = 'gds')
results$links$pubmed_gds
[1] "200012783" "200012782" "200012781"
3 результата - это идентификаторы для связанных записей набора данных GEO. Вы можете конвертировать их в образцы GSE, используя entrez_summary
.
Вот несколько уродливый sapply
, который может служить основой для функции:
sapply(results$links$pubmed_gds, function (id) entrez_summary("gds", id)$accession,
USE.NAMES = FALSE)
[1] "GSE12783" "GSE12782" "GSE12781"