首页 > 其他 > 详细

plink软件中--set-missing-var-ids @:#的作用

时间:2020-10-10 23:37:22      阅读:104      评论:0      收藏:0      [点我收藏+]

1、测试数据

[root@linuxprobe test]# ls
test.map  test.ped
[root@linuxprobe test]# cat test.map  ## 测试数据snpID为.
1       .       0       55910
1       .       0       85204
1       .       0       122948
1       .       0       203750
1       .       0       312707
1       .       0       356863
1       .       0       400518
1       .       0       487423
[root@linuxprobe test]# cat test.ped
DOR     1       0       0       0       -9      G A     G G     G G     G C     G G     C C     C C     G G
DOR     2       0       0       0       -9      A A     G G     A G     C C     G G     G C     C C     G G
DOR     3       0       0       0       -9      A A     G G     A G     G C     C G     C C     G G     G G
DOR     4       0       0       0       -9      A A     G G     G G     G G     C G     G G     G G     G G
DOR     5       0       0       0       -9      A A     G G     A G     G C     G G     C C     G G     G G
DOR     6       0       0       0       -9      A A     G G     A A     C C     G G     C C     G G     G G
DOR     7       0       0       0       -9      G G     G G     G G     C C     G G     C C     G G     A G
DOR     9       0       0       0       -9      G G     G G     G G     C C     G G     C C     G G     A G

 

2、

 plink --file test --set-missing-var-ids @:# --recode tab --out a ;rm *.log *.nosex
[root@linuxprobe test]# cat a.map  ## 可见snpID发生变化
1       1:55910 0       55910
1       1:85204 0       85204
1       1:122948        0       122948
1       1:203750        0       203750
1       1:312707        0       312707
1       1:356863        0       356863
1       1:400518        0       400518
1       1:487423        0       487423
[root@linuxprobe test]# md5sum *.ped
09e293ef08104b377871b606816f52a6  a.ped
09e293ef08104b377871b606816f52a6  test.ped

 

 plink --file test --set-missing-var-ids @:#\$1,\$2 --recode tab --out a; rm *.log *.nosex  
[root@linuxprobe test]# cat a.map  ## 加上了snp组合
1       1:55910A,G      0       55910
1       1:852040,G      0       85204
1       1:122948A,G     0       122948
1       1:203750C,G     0       203750
1       1:312707C,G     0       312707
1       1:356863C,G     0       356863
1       1:400518C,G     0       400518
1       1:487423A,G     0       487423
[root@linuxprobe test]# md5sum *.ped
09e293ef08104b377871b606816f52a6  a.ped
09e293ef08104b377871b606816f52a6  test.ped
[root@linuxprobe test]# cat a.ped
DOR     1       0       0       0       -9      G A     G G     G G     G C     G G     C C     C C     G G
DOR     2       0       0       0       -9      A A     G G     A G     C C     G G     G C     C C     G G
DOR     3       0       0       0       -9      A A     G G     A G     G C     C G     C C     G G     G G
DOR     4       0       0       0       -9      A A     G G     G G     G G     C G     G G     G G     G G
DOR     5       0       0       0       -9      A A     G G     A G     G C     G G     C C     G G     G G
DOR     6       0       0       0       -9      A A     G G     A A     C C     G G     C C     G G     G G
DOR     7       0       0       0       -9      G G     G G     G G     C C     G G     C C     G G     A G
DOR     9       0       0       0       -9      G G     G G     G G     C C     G G     C C     G G     A G

 

plink --file test --set-missing-var-ids @:#[b37]\$1,\$2 --recode tab --out a; rm *.log *.nose
[root@linuxprobe test]# md5sum *.ped
09e293ef08104b377871b606816f52a6  a.ped
09e293ef08104b377871b606816f52a6  test.ped
[root@linuxprobe test]# cat a.map ## 同上 ??
1       1:55910[b37]A,G 0       55910
1       1:85204[b37]0,G 0       85204
1       1:122948[b37]A,G        0       122948
1       1:203750[b37]C,G        0       203750
1       1:312707[b37]C,G        0       312707
1       1:356863[b37]C,G        0       356863
1       1:400518[b37]C,G        0       400518
1       1:487423[b37]A,G        0       487423

 

plink软件中--set-missing-var-ids @:#的作用

原文:https://www.cnblogs.com/liujiaxin2018/p/13795556.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!