R包 ieugwasr: 如何查找千人基因组对应的EAF
要查找千人基因组对应的EAF,可以使用'ieugwasr'包中的'lookup_1kg_eaf()'函数。该函数需要提供SNP的rs ID和所在染色体的编号。以下是使用该函数的示例代码:
library(ieugwasr)
# 输入SNP的rs ID和所在染色体的编号
rsid <- 'rs123456'
chrom <- 1
# 查找对应的EAF
result <- lookup_1kg_eaf(rsid, chrom)
print(result)
该函数将返回一个数据框,其中包含了对应SNP在千人基因组中的EAF信息。如果找不到对应的EAF信息,则返回一个空的数据框。
原文地址: https://www.cveoy.top/t/topic/mKUb 著作权归作者所有。请勿转载和采集!