对genotype的等位型进行计数,需要用到以下参数:

--freq Allele frequencies
--counts Modifies --freq to report actual allele counts

具体用法如下命令:

/plink-1.07-x86_64/plink --freq --counts --noweb --bfile file  --make-bed --out file

 

.frq.count (basic allele count report)

Produced by '--freq counts'. Valid input for --read-freq.

A text file with a header line, and then one line per variant with the following seven fields:

CHR

Chromosome code

SNP

Variant identifier

A1

Allele 1 (usually minor)

A2

Allele 2 (usually major)

C1

Allele 1 count

C2

Allele 2 count

G0

Missing genotype count (so C1 + C2 + 2 * G0 is constant on autosomal variants)

 

注意:这里用到的是plink的1.07版本,如果是1.9或者2.0版本,则命令可能会有些许的改动,具体参考官网。

参考:http://zzz.bwh.harvard.edu/plink/reference.shtml

相关文章:

  • 2021-11-04
  • 2021-07-10
  • 2021-07-30
  • 2021-11-08
  • 2022-01-01
  • 2022-12-23
  • 2021-08-02
  • 2021-10-14
猜你喜欢
  • 2021-12-21
  • 2021-07-22
  • 2022-01-30
  • 2021-11-11
  • 2021-11-19
  • 2021-06-14
  • 2021-12-08
相关资源
相似解决方案