R包ieugwasr 如何查找千人基因组对应的EAF

要查找千人基因组对应的EAF，可以使用R包“ieugwasr”中的函数“get_eaf()”。该函数需要指定要查询的SNP ID和参考基因组版本。下面是示例代码：

library(ieugwasr)

# 指定要查询的SNP ID和参考基因组版本
snp_id <- "rs123456"
ref_version <- "hg19"

# 调用get_eaf()函数查询EAF
eaf <- get_eaf(snp_id, ref_version)

# 输出查询结果
print(eaf)

其中，snp_id是要查询的SNP ID，ref_version是参考基因组版本。查询结果eaf是一个数据框，包含了该SNP在不同种群中的EAF值。