修改正負鏈用到的參數爲--flip
code
假定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 test
file
其中,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