检索biomaRt中的SNP位置，命令挂起

Question

我在试着检索SNP位置的信息。我试着按照这个网站的答案来做，但是这个命令不再起作用了：

library(biomaRt) # biomaRt_2.30.0

snp_mart = useMart("ENSEMBL_MART_SNP", dataset="hsapiens_snp")

snp_ids = c("rs16828074", "rs17232800")
snp_attributes = c("refsnp_id", "chr_name", "chrom_start")

snp_locations = getBM(attributes=snp_attributes, filters="snp_filter", 
                      values=snp_ids, mart=snp_mart)

在等待了很长时间后，我得到了以下错误：

Error in value[[3L]](cond) : 
  Request to BioMart web service failed. Verify if you are still connected to the internet.  Alternatively the BioMart web service is temporarily down.

与上一个版本相比，biomaRt命令有什么变化吗？还是我做错了什么？

Answer 1

这更像是一种变通方法，而不是对这个问题的明确答案。但是现在Ensembl的版本是90。如果我使用以前版本(来自http://may2017.archive.ensembl.org的v89)中的归档主机，则SNP数据集将再次工作。因此，当v90不适用于SNP时，以下是我的临时解决方案：

library("biomaRt")
snp_mart = useMart(biomart = "ENSEMBL_MART_SNP", dataset="hsapiens_snp", host='may2017.archive.ensembl.org')

snp_ids = c("rs16828074", "rs17232800")

snp_attributes = c("refsnp_id", "chr_name", "chrom_start")

snp_locations = getBM(attributes=snp_attributes, filters="snp_filter", 
                       values=snp_ids, mart=snp_mart)

结果应该如下所示：

snp_locations
   refsnp_id chr_name chrom_start
1 rs16828074        2   231454043
2 rs17232800       18    68625022