GWAS -4 VCF格式文件轉(zhuǎn)為Plink文件

1 VCF格式文件為Plink文件

image.png

參考:https://cloud.tencent.com/developer/article/1556166

2 plink轉(zhuǎn)vcf

plink1.9版本支持轉(zhuǎn)化為VCFv4.2格式

plink2.0版本支持轉(zhuǎn)化為VCFv4.3格式

兩個(gè)版本用到的命令不一樣

對(duì)于plink1.9版本啄巧,轉(zhuǎn)化為vcf文件的命令行為:

plink --bfile binary_fileset --recode vcf-iid --out new_vcf

生成的vcf為4.2版本

對(duì)于plink2.0版本悼院,轉(zhuǎn)化為vcf文件的命令行為:

plink --bfile binary_fileset --export vcf --out new_vcf

生成的vcf為4.3版本

參考:https://www.cnblogs.com/chenwenyan/p/8574237.html

3 plink轉(zhuǎn)vcf


3.1 plink --bfile test_1 --recode vcf --out test_vcf  ###這樣就把plink(二進(jìn)制)文件轉(zhuǎn)換成vcf格式了

3.2  plink --file  test_1  --recode vcf --out test_vcf   ### 直接map ped文件轉(zhuǎn)為vcf格式

鏈接:http://www.reibang.com/p/8b4e7b3b7f5e

4 PLNIK 的多種文件格式轉(zhuǎn)換

vcf 轉(zhuǎn)為 ped/map

使用vcftools

vcftools --vcf snp.vcf --plink --out snp

使用plink

plink --vcf snp.vcf --recode --out snp
ped和map文件是Plink的基本格式。

ped文件包含以下幾列:

第一列:Family ID死遭。

第二列:Individual ID偿短。自然群體這列和Family ID是一樣的欣孤。

第三列:Paternal ID。未提供信息的話這列為0昔逗。

第四列:Maternal ID降传。未提供信息的話這列為0。

第五列:Sex勾怒。未提供信息的話這列為0搬瑰。

第六列:Phenotype款票。一般來(lái)說(shuō),直接拿vcf轉(zhuǎn)換的話這列為-9泽论,也就是缺失。

第七列開(kāi)始就是個(gè)體在每個(gè)標(biāo)記位點(diǎn)的基因型卡乾。

map文件包含以下幾列:

第一列:染色體編號(hào)翼悴。

第二列:SNP編號(hào)。

第三列:遺傳距離幔妨。未提供信息的話這列為0鹦赎。

第四列:物理位置。

ped/map 與 tped/tfam 格式互換

ped/map轉(zhuǎn)換為tped/tfam

plink --file snp --recode --transpose --out snp_test

tped/tfam轉(zhuǎn)換為ped/map

plink --tfile snp_test --recode --out snp
ped/map 與 bed/bim/fam互換

ped/map轉(zhuǎn)換為bed/bim/fam

plink --file snp --make-bed --out snp_test

bed/bim/fam轉(zhuǎn)換為ped/map

plink --bfile snp_test --recode --out snp
tped/tfam 與 bed/bim/fam互換

tped/tfam轉(zhuǎn)換為bed/bim/fam

plink --tfile snp --make-bed --out snp_test

bed/bim/fam轉(zhuǎn)換為tped/tfam

plink --bfile snp_test --recode --transpose --out snp
bed/bim/fam 轉(zhuǎn)為 vcf

bed/bim/fam 轉(zhuǎn)為 vcf

plink --bfile snp --export vcf --out snp_test
常用的Plink格式轉(zhuǎn)換就是這些误堡,大家可以根據(jù)自己實(shí)際需要相互轉(zhuǎn)換古话。

染色體的設(shè)置

因?yàn)镻LINK默認(rèn)的設(shè)置是人的染色體, 所以動(dòng)物中锁施,我們應(yīng)該設(shè)置
--chr-set 19 # 豬
已有的選擇:
--cow
--dog
--horse
--mouse
--rice
--sheep
參考:https://zhuanlan.zhihu.com/p/109071456

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末陪踩,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子悉抵,更是在濱河造成了極大的恐慌肩狂,老刑警劉巖,帶你破解...
    沈念sama閱讀 206,126評(píng)論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件姥饰,死亡現(xiàn)場(chǎng)離奇詭異傻谁,居然都是意外死亡,警方通過(guò)查閱死者的電腦和手機(jī)列粪,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,254評(píng)論 2 382
  • 文/潘曉璐 我一進(jìn)店門(mén)审磁,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái),“玉大人岂座,你說(shuō)我怎么就攤上這事态蒂。” “怎么了掺逼?”我有些...
    開(kāi)封第一講書(shū)人閱讀 152,445評(píng)論 0 341
  • 文/不壞的土叔 我叫張陵吃媒,是天一觀的道長(zhǎng)。 經(jīng)常有香客問(wèn)我吕喘,道長(zhǎng)赘那,這世上最難降的妖魔是什么? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 55,185評(píng)論 1 278
  • 正文 為了忘掉前任氯质,我火速辦了婚禮募舟,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘闻察。我一直安慰自己拱礁,他們只是感情好琢锋,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,178評(píng)論 5 371
  • 文/花漫 我一把揭開(kāi)白布。 她就那樣靜靜地躺著呢灶,像睡著了一般吴超。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上鸯乃,一...
    開(kāi)封第一講書(shū)人閱讀 48,970評(píng)論 1 284
  • 那天鲸阻,我揣著相機(jī)與錄音,去河邊找鬼缨睡。 笑死鸟悴,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的奖年。 我是一名探鬼主播细诸,決...
    沈念sama閱讀 38,276評(píng)論 3 399
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼陋守!你這毒婦竟也來(lái)了震贵?” 一聲冷哼從身側(cè)響起,我...
    開(kāi)封第一講書(shū)人閱讀 36,927評(píng)論 0 259
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤嗅义,失蹤者是張志新(化名)和其女友劉穎屏歹,沒(méi)想到半個(gè)月后,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體之碗,經(jīng)...
    沈念sama閱讀 43,400評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡蝙眶,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 35,883評(píng)論 2 323
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了褪那。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片幽纷。...
    茶點(diǎn)故事閱讀 37,997評(píng)論 1 333
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖博敬,靈堂內(nèi)的尸體忽然破棺而出友浸,到底是詐尸還是另有隱情,我是刑警寧澤偏窝,帶...
    沈念sama閱讀 33,646評(píng)論 4 322
  • 正文 年R本政府宣布收恢,位于F島的核電站,受9級(jí)特大地震影響祭往,放射性物質(zhì)發(fā)生泄漏伦意。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,213評(píng)論 3 307
  • 文/蒙蒙 一硼补、第九天 我趴在偏房一處隱蔽的房頂上張望驮肉。 院中可真熱鬧,春花似錦已骇、人聲如沸离钝。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 30,204評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)卵渴。三九已至慧域,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間浪读,已是汗流浹背吊趾。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 31,423評(píng)論 1 260
  • 我被黑心中介騙來(lái)泰國(guó)打工, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留瑟啃,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 45,423評(píng)論 2 352
  • 正文 我出身青樓揩尸,卻偏偏與公主長(zhǎng)得像蛹屿,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子岩榆,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,722評(píng)論 2 345

推薦閱讀更多精彩內(nèi)容