2018-10-31 GWAS实战(三)plink 进阶之认识常用文件格式

学会了基本的plink之后,就要开始学习了解输入文件的格式了。
详细的格式建议看官网
因为之前以及接触过几次,这里来个详细的归总
主要认识 ped map bed fam bim 这五种格式
其中 :

ped 和 map 是一组的
bed fam bim 是一组的

ped

ped.png

1.ped 包含样本的谱系信息和基因型信息
2.ped 必须与fam 文件一起,前6个字段与fam文件对应
3.使用recode 产生ped文件

ped-eg.png

举个例子
每一行是一个个体,前六列是固定的,从第七列开始后面就是每个snp位点的基因型情况,第七列第八列就是第一个snp位点,第九列第十列就是第二个snp位点,依次类推。。。 其中 0 代表 no calling, 从图中可以看出 第一个个体第一个snp 就是 00 第二个个体第二个snp 就是 AG 。
我们再回过头来看前六列:
第一列 Family ID
第二列 Individual ID
第三列 Paternal ID
第四列 Maternal ID
第五列 Sex (1=male; 2=female; other=unknown)
第六列 Phenotype

这里要主要的是,对于表型


注意.png

map

map

1.map文件和ped文件是一起的,表示每个SNP的信息
2.总共有4列

详细解释

举个例子

第一列代表染色体
第二列代表snp的名称
第三列代表摩尔距离,一般用不上
第四列代表物理距离,这个有用的

第二组开始

bed

bed.png

1.首先这里强调的是bed文件与UCSC Genome Browser's BED format 是完全不一样的。
2.bed 文件与bim fam 文件一起的
3.bed文件是一个二进制文件,所以你是看不来的


一定要和其他两个文件一起

采用二进制的方式保存文件,代表的数据意义和其他两个文件对应

bim

bim

1.从这个说明我们可以看出bim文件是对map文件的拓展
2.总共有六行,包含了snp(variants)的具体信息、
3.这个文件也是要和bed fam 文件结合来使用


eg

我们来看个例子:
1.第一列是染色体信息
2.第二列是snp的名字
3.第三列是摩尔距离,文件中说可以用0,没关系
4.第四列是物理距离
5.第五列是次要等位基因
6.第六列是主要等位基因

这里的信息要和bed文件对应起来

fam

fam

1.fam文件记录了每个样本家系的信息
2.fam文件也是必须和bed bim文件一起使用
3.主要有6 列

fam-eg

1.第一列是Family ID ('FID')
2.第二列是Within-family ID ('IID'; cannot be '0')
3.第三列是Within-family ID of father ('0' if father isn't in dataset)
4.第四列是Within-family ID of mother ('0' if mother isn't in dataset)
5.第五列是性别 ('1' = male, '2' = female, '0' = unknown)
6.第六列是表型 ('1' = control, '2' = case, '-9'/'0'/non-numeric = missing data if case/control)

看到这里我发现其实这个和ped文件格式太像了,是代表信息的提取,这样的格式更加利于观察

以上就会5中常用格式的介绍,中间还有很多的细节,需要去理解。

这两组格式的数据可以相互转化,还可以使用--recode(eg. recodeA recodeAD)来进行编码,这个使用数据更加灵活

最简单的就是入下

ped map 转 bed bim fam

plink --file toy --out toy
ped:map转bed:bim:fam

当然有些时候转要加 --make-bed

现在我再倒过来转回去

bed bim fam 转 ped map

plink --bfile toy --recode --out toy1
反转

OK了

PS:一定要多看官方的网站

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 205,132评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 87,802评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,566评论 0 338
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,858评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,867评论 5 368
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,695评论 1 282
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,064评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,705评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 42,915评论 1 300
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,677评论 2 323
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,796评论 1 333
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,432评论 4 322
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,041评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,992评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,223评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,185评论 2 352
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,535评论 2 343

推荐阅读更多精彩内容

  • 欢迎关注公众号:oddxixplink的主要功能:数据处理,质量控制的基本统计,群体分层分析,单位点的基本关联分析...
    oddxix阅读 12,492评论 0 19
  • introduction Plink算是生物信息最常用的软件之一了,而且很多软件的输入也会使用Plink的文件格式...
    银色麦穗阅读 21,601评论 7 9
  • 在上一篇文章中,分享了有关基因组Phasing的原理,一共有三种,分别是:家系关系分型(Related indiv...
    黄树嘉阅读 6,828评论 3 18
  • 我已经讨厌惨了两个人相处中慢慢展示缺陷后分道扬镳的这个过程了,所以从开始我就要给你看到,我傲慢又自大,冷血又毒舌,...
    kylin进化论阅读 325评论 0 0
  • 今天我们演了绿野仙踪,太开心了我演的狮子,我们的王 睿铭和王译萱忘台词了后来他们想起来了呢托托一直在台上第5幕托托...
    jshs阅读 292评论 0 0