从 R biomart getBM 获得 NA 的次要等位基因？

Question

有谁知道为什么下面的代码给出了minor_allele和minor_allele_freq的NA列（通过这种方法提取的其他数据似乎很好）？

library("biomaRt"); snp.db <- useMart(host = "https://www.ensembl.org", "ENSEMBL_MART_SNP", dataset = "hsapiens_snp"); nt.biomart <- getBM(attributes = c("refsnp_id", "minor_allele", "minor_allele_freq"), filters = c("snp_filter"), values = c("rs10788066", "rs765840164"), mart = snp.db, uniqueRows = TRUE); nt.biomart

Ensembl 上绝对提供数据（例如 https://www.ensembl.org/Homo_sapiens/Variation/Explore?db=core;r=10:120308754-120309754;v=rs10788066;vdb=variation;vf =654619804）但是数据提取方法似乎没有提取次要等位基因数据 - 编码有问题吗？

Answer 1

参见https://support.bioconductor.org/p/70935/。 snp.db <- useMart("ENSEMBL_MART_SNP", dataset="hsapiens_snp", host="https://grch37.ensembl.org", path="/biomart/martservice", archive=FALSE). How I'm meant to get this information is anyone's guess. Strangely, replacing grch37 with grch38 doesn't work.