plink修改正負鏈(--flip, change the positive and negative stand)

修改正負鏈用到的參數爲--flipcode

假定trial.bim的內容以下:ip

trial.bim
1 rs142578063 0 732746 G A
1 rs144022023 0 732801 G A
1 rs12131618 0 732809 C T
1 rs369225293 0 732909 A G
1 rs373491075 0 732935 C G
1 rs369030935 0 732989 T C
1 rs4951860 0 733013 C T
1 rs367917297 0 733429 T G
1 rs376492866 0 733895 G Atest

如今要把修改rs14257806三、rs367917297這兩個位點的正負鏈;即rs142578063的G變爲C, A變爲T,rs367917297的T 變爲A, G變爲 C, 則須要用到如下命令
/software/plink --bfile trial --flip test.missnp --make-bed --out testfile

其中,test.missnp的內容以下:im

test.missnp
rs142578063
rs367917297make

修改後的結果以下:co

test.bim
1 rs142578063 0 732746 C T
1 rs144022023 0 732801 G A
1 rs12131618 0 732809 C T
1 rs369225293 0 732909 A G
1 rs373491075 0 732935 C G
1 rs369030935 0 732989 T C
1 rs4951860 0 733013 C T
1 rs367917297 0 733429 A C
1 rs376492866 0 733895 G Ablock

相關文章
相關標籤/搜索