一行命令批量修改染色体和位置为RS号

鉴于太多人问我怎么批量根据chr:pos查找RS号，在这里出一个教程。

注意以下教程展示的是修改hg19基因组版本的RS号，如果你的数据是其他版本的，请修改为对应版本的数据。

假定数据是test.txt，内容如下：

现在希望根据第一列chr:pos找到对应的RS号，实现以下的效果：

则可以用dplyr::left_join参数，具体实现过程如下所示：

install.packages(dplyr)
library(dplyr)
tes = read.table("test.txt",header=T,check.names=F,sep="\t") #注意这里我设置的是制表符分隔符，如果你的文件不是制表符的话，需要修改成对应的分隔符
match = read.table("snp150_hg19.txt",header=T,check.names=F,sep="\t")
need=dplyr::left_join(tes,match,by="chromosome:start") #如果snp150_hg19.txt文件中有对应的RS号，则比对到test.txt文件中，如果没有的话，就变为NA
write.table(need, file ="clean.txt", sep ="\t", row.names =FALSE, col.names =TRUE, quote =FALSE) #保存文件

上述命令需要用到snp150_hg19.txt文件，其内容如下所示：

snp150_hg19.txt文件可以从https://hgdownload.soe.ucsc.edu/downloads.html网站下载，也可以通过后台回复RS号修改获取，下载后请解压后再使用;

这个文件很大，总共有234104111个位点，解压后5G左右，只能通过网盘传输；

另外，我也提供了hg38的基因组版本snp151_hg38.txt.gz，需要的同样后台回复RS号修改。

由于本人只做人类，所以只提供了人类基因组版本的教程，做其他物种的小伙伴们需要自行从https://hgdownload.soe.ucsc.edu/downloads.html网站下载。

致谢橙子牛奶糖（陈文燕），请用参考模版：We thank the blogger (orange_milk_sugar, Wenyan Chen) for XXX

感谢小可爱们多年来的陪伴，我与你们一起成长~

一行命令批量修改染色体和位置为RS号相关推荐

批量修改linux换行格式,linux中sed命令批量修改
sed命令下批量替换文件内容格式: sed -i "s/查找字段/替换字段/g" `grep 查找字段 -rl 路径` 文件名 -i 表示inplace edit,就地修改文件 ...
使用命令批量修改文件名
目录需求思路使用awk命令将修改前.修改后的文件名打出来使用xargs命令将修改前.修改后的文件名传入mv命令作为参数需求从某地下载了一些文件,这些文件名字都是aaaa-xxssxxx的格 ...
woo语言一行代码批量修改文件名称
woo语言一行代码批量修改文件名称 1.windows平台 – 下方代码在需要修改的文件夹内使用cmd执行即可,把wo后缀修改为woo后缀: woo.exe -r "a = _args() ...
powershell一行代码批量修改文件名（附命令详解）
文章目录写在前面具体步骤命令详解注意事项简单示例总结写在前面在日常办公中,批量修改照片/视频的文件名是一件让很多人头疼却又不得不面对的事情.我就经常会遇到下载的视频文件名中带有广告等信 ...
使用chpasswd命令批量修改系统用户密码
chpasswd命令工作原理: 从系统的标准输入读入用户的名称和口令,并利用这些信息来更新系统上已存在的用户的口令! 语法: 1:# echo 用户名:密码 | chpasswd 2:# chpass ...
linux chpasswd命令批量修改用户密码
chpasswd 批量修改用户密码工作原理: 从系统的标准输入读入用户的名称和口令,并利用这些信息来更新系统上已存在的用户的口令语法: 1:# echo 用户名:密码 | chpasswd 2: ...
一行命令同时修改maven项目中多个mudule的版本号
这是Hollis的第 215 篇原创分享作者 l Hollis 来源 l Hollis(ID:hollischuang) Maven,是一个Java开发比较常用的项目管理工具,可以对 Java 项目 ...
【BAT-表姐御用05ren命令批量修改后缀】
一.简析从网上下载图片的时候,后缀类型那是五花八门的,不利于后续开展工作,故,这里演示一下如何批量修改后缀,统一文件格式. 二.命令表达式这里修改后缀是属于修改文件名类的,所以我们用ren 命令 ...
使用ren命令批量修改照片名
dos命令ren修改文件名的格式 ren [需要修改的文件的绝对路径包含后缀名] [修改后的名称,包含后缀名] 如果文件名有空格,需要使用英文双引号包裹. 举个栗子: 批量修改照片名,以身份证命名.输 ...

一行命令批量修改染色体和位置为RS号

一行命令批量修改染色体和位置为RS号相关推荐

最新文章

热门文章