This is a test version of Biostars. For the public version, visit https://www.biostars.org.
PLINK BED format: how interpret 0, 1, 2

Hello,

I'm working with RAD-seq data in Stacks. And I obtained my final dataset in .PED format. Then I converted the plink PED file (.ped and .map) into a .BED format file using the following command:

plink -file input.plink --set-missing-var-ids "@_#" --make-bed --out output.file --allow-extra-chr --mind 0.3 --geno 0

And I was wondering how I should interpret the coding of the alleles: 0, 1 , 2. Given that I did not specify any '--recode' command I'm not sure what would be the correct interpretation:

Is it correct to interpret as:

0 homozygous for the major allele 1 heterozygous 2 homozygous for minor alleles

snp plink
• 2,245 views • link
updated by chrchang523 chrchang523 • written by CaffeSospeso • CaffeSospeso

1 answer

0 values correspond to the allele in the 5th column of the .bim file, and 2 values correspond to the allele in the 6th column. The 6th column value actually defaults to major, not minor.

• 0 views • link
written by chrchang523 chrchang523

Log in to answer this question.

More posts like this