plink格式数据(ped文件第5列、第6列)

 

001、ped文件第5列:表示性别信息(1表示male、2表示female、0或者其他字母表示未知)

root@DESKTOP-1N42TVH:/home/test3# ls
outcome.map  outcome.ped
root@DESKTOP-1N42TVH:/home/test3# cat outcome.ped     ## 测试数据, 第5列
DOR     1       0       0       0       9       G G     C C     G G     G G     A G     A A
DOR     2       0       0       ewr     9       G G     G C     G G     G G     G G     A A
DOR     3       0       0       1       19      G G     C C     G G     G G     G G     A A
DOR     4       0       0       2       39      G G     C C     G G     G G     G G     A A
DOR     5       0       0       2       93      G G     C C     G G     G G     G G     A A
DOR     6       0       0       kku     -9      G G     C C     G G     G G     G G     A A
root@DESKTOP-1N42TVH:/home/test3# plink --file outcome --recode tab --out outcome; rm *.log *.nosex
PLINK v1.90b6.26 64-bit (2 Apr 2022)           www.cog-genomics.org/plink/1.9/
(C) 2005-2022 Shaun Purcell, Christopher Chang   GNU General Public License v3
Logging to outcome.log.
Options in effect:
  --file outcome
  --out outcome
  --recode tab

16007 MB RAM detected; reserving 8003 MB for main workspace.
.ped scan complete (for binary autoconversion).
Performing single-pass .bed write (6 variants, 6 people).
--file: outcome-temporary.bed + outcome-temporary.bim + outcome-temporary.fam
written.
6 variants loaded from .bim file.
6 people (1 male, 2 females, 3 ambiguous) loaded from .fam.   ## 1个male,2个female,0或者其他表示未知
Ambiguous sex IDs written to outcome.nosex .
5 phenotype values loaded from .fam.
Warning: Ignoring phenotypes of missing-sex samples.  If you don't want those
phenotypes to be ignored, use the --allow-no-sex flag.
Using 1 thread (no multithreaded calculations invoked).
Before main variant filters, 6 founders and 0 nonfounders present.
Calculating allele frequencies... done.
Total genotyping rate is exactly 1.
6 variants and 6 people pass filters and QC.
Phenotype data is quantitative.
--recode ped to outcome.ped + outcome.map ... done.
root@DESKTOP-1N42TVH:/home/test3# cat outcome.ped     ## 第5列性别未知的情况下,自动转换为了0
DOR     1       0       0       0       9       G G     C C     G G     G G     A G     A A
DOR     2       0       0       0       9       G G     G C     G G     G G     G G     A A
DOR     3       0       0       1       19      G G     C C     G G     G G     G G     A A
DOR     4       0       0       2       39      G G     C C     G G     G G     G G     A A
DOR     5       0       0       2       93      G G     C C     G G     G G     G G     A A
DOR     6       0       0       0       -9      G G     C C     G G     G G     G G     A A

 

 

 

002、ped文件第6列表示表型,-9或者其他字符表示缺失。

root@DESKTOP-1N42TVH:/home/test3# ls
outcome.map  outcome.ped
root@DESKTOP-1N42TVH:/home/test3# cat outcome.ped     ## 测试数据,第6列表示性别,-9或者其他字符表示缺失
DOR     1       0       0       0       9       G G     C C     G G     G G     A G     A A
DOR     2       0       0       0       9       G G     G C     G G     G G     G G     A A
DOR     3       0       0       1       19      G G     C C     G G     G G     G G     A A
DOR     4       0       0       2       kku     G G     C C     G G     G G     G G     A A
DOR     5       0       0       2       ewr     G G     C C     G G     G G     G G     A A
DOR     6       0       0       0       -9      G G     C C     G G     G G     G G     A A
root@DESKTOP-1N42TVH:/home/test3# plink --file outcome --recode tab --out outcome; rm *.log *.nosex
PLINK v1.90b6.26 64-bit (2 Apr 2022)           www.cog-genomics.org/plink/1.9/
(C) 2005-2022 Shaun Purcell, Christopher Chang   GNU General Public License v3
Logging to outcome.log.
Options in effect:
  --file outcome
  --out outcome
  --recode tab

16007 MB RAM detected; reserving 8003 MB for main workspace.
.ped scan complete (for binary autoconversion).
Performing single-pass .bed write (6 variants, 6 people).
--file: outcome-temporary.bed + outcome-temporary.bim + outcome-temporary.fam
written.
6 variants loaded from .bim file.
6 people (1 male, 2 females, 3 ambiguous) loaded from .fam.
Ambiguous sex IDs written to outcome.nosex .
3 phenotype values loaded from .fam.   ## -9或者其他字符表示缺失
Warning: Ignoring phenotypes of missing-sex samples.  If you don't want those
phenotypes to be ignored, use the --allow-no-sex flag.
Using 1 thread (no multithreaded calculations invoked).
Before main variant filters, 6 founders and 0 nonfounders present.
Calculating allele frequencies... done.
Total genotyping rate is exactly 1.
6 variants and 6 people pass filters and QC.
Phenotype data is quantitative.
--recode ped to outcome.ped + outcome.map ... done.
root@DESKTOP-1N42TVH:/home/test3# cat outcome.ped     ## 表型缺失自动转换为了-9
DOR     1       0       0       0       9       G G     C C     G G     G G     A G     A A
DOR     2       0       0       0       9       G G     G C     G G     G G     G G     A A
DOR     3       0       0       1       19      G G     C C     G G     G G     G G     A A
DOR     4       0       0       2       -9      G G     C C     G G     G G     G G     A A
DOR     5       0       0       2       -9      G G     C C     G G     G G     G G     A A
DOR     6       0       0       0       -9      G G     C C     G G     G G     G G     A A

 

小结:

001、ped文件的第5列表示性别:1表示male; 2表示female; 0或者其他任意字符表示未知,默认缺失为0.

002、ped文件的第6列表示表型:-9或者其他任意字符表示缺失,默认缺失为-9.

 

posted @ 2022-07-08 20:15  小鲨鱼2018  阅读(539)  评论(0编辑  收藏  举报