SNP的rs号和染色体位置信息相互转换_snp位置转化成snp号

作者：凡人多烦事01 | 2024-04-21 14:40:24

踩

snp位置转化成snp号

BiocManager::install("SNPlocs.Hsapiens.dbSNP144.GRCh37")
library(SNPlocs.Hsapiens.dbSNP144.GRCh37)
snps<-SNPlocs.Hsapiens.dbSNP144.GRCh37
for (i in unique(BMI2$Chr)){
  my_pos<-BMI2$Position[BMI2$Chr==i]
  chr_snps<-snpsBySeqname(snps, as.character(i))
  idx<-match(my_pos,pos(chr_snps))
  rsids<-mcols(chr_snps)$RefSNP_id[idx]
  BMI2$rsid[BMI2$Chr==i]<-rsids
  print(paste(as.character(i),"is ok"))
}
1
2
3
4
5
6
7
8
9
10
11

37/38自主选择

声明：本文内容由网友自发贡献，不代表【wpsshop博客】立场，版权归原作者所有，本站不承担相应法律责任。如您发现有侵权的内容，请联系我们。转载请注明出处：https://www.wpsshop.cn/w/凡人多烦事01/article/detail/463535