Plink 命令汇总

  • 格式:docx
  • 大小:31.01 KB
  • 文档页数:4

分型数据结果出来后,首先要做的是按照染色体将SNP位点分类(在位点位于不同染色体的情况下),
1)如上图所示,按照不同染色体分好位点数;
2)首先分析的是control数据里的LD,如果存在LD,即D’接近1,r2大小似乎关系不大,就可以判定有LD,这个在Haploview软件里面分析。

分析时的数据格式整理成Linkage format。

可以把一个数据弄成ped,一个是map的。

(奇怪的是我用plink format,系统总说我少SNP column,给ped加head,genotype那里换成rs号或SNP1,SNP2之类的都不行,map文件是标准的plink格式,headless),有时间再找出问题在哪。

Map的文件:只有2列,一列是rs号,一列是染色体上的起始位置,至于这个位置,我用的是Hapmap数据库上的。

不同数据库之间的相差不大,起始需要的是每2个SNP之间的相对距离,估计软件自己会计算。

Ped文件:
前6列是固定的,一次为FID,IID,father ID,mother ID,sex(1 male; 2 female),phenotype(1 control;2 case),后面就是基因型的了。

按照SNP号排列好,但是不要表头。

注意在map文件中,SNP的顺序要和这个对上。

3)结果可以看到LD,右键点击,可以看到D’,r2;点Haplotype,不知道怎么有的没有结果出来,可能是没有单倍型对于有单倍型的,比如,对照中有这样的结果,我再计算case中的单倍型结果。

得出这些数据后,我是再把数据导入到SPSS,运用卡方分析case control中数据是否有差异。

可惜今天分析的结果都没有差异,奇怪为什么直接用plink的haps命令结果有一个是单倍型显著的呢?(P=0.045)。