grm.bin and grm.N.bin have different length when using --make-grm-bin

20 views
Skip to first unread message

Yanyu Liang

unread,
May 14, 2022, 2:01:22 PM5/14/22
to plink2-users
Hi there,

I ran PLINK2 with the following command 

plink2 \
  --pfile mygeno \
  --make-grm-bin \
  --out mygrm

I'd expect the resulting mygrm.grm.bin and mygrm.grm.N.bin have equal size but it turned out that grm.N.bin contains fewer entries. It seems that *.grm.bin contains n*(n+1)/2 entries whereas *.grm.bin contains n*(n-1)/2 entries. From my understanding, the former should be expected (which is also the number of rows given by --make-grm-list). 

I was wondering if there is an issue in --make-grm-bin code and it will become a lot handier if it can be resolved.

Thanks!
Yanyu

Christopher Chang

unread,
May 14, 2022, 3:18:36 PM5/14/22
to plink2-users
Thanks for reporting this; yes, the .grm.N.bin file had too few entries when no missing genotypes were present.  Bugfix has been posted.
Reply all
Reply to author
Forward
0 new messages