We have data in PGEN format, and as part of our workflow, we initially created temporary PSAM, PVAR, and PGEN files using the following command:
./plink2 --bgen data.bgen ref-first --sample data.sample --set-missing-var-ids @:#:'\$r':'\$a' --new-id-max-allele-len 99 truncate --make-pgen --out data.intermediate
These two steps are essential components of our workflow. However, when attempting to run an association analysis on the resulting BED file using the following command:
./plink2 --bfile cons_imp_data_sep_8th --glm --out test_assoc --covar cons_imp_data_sep_8th_imp_generic.cov --covar-name Sex --pheno constipation_cc.phe --pheno-name constipation_cc
log file information:PLINK v2.00a4.4LM 64-bit Intel (21 Jun 2023)
www.cog-genomics.org/plink/2.0/(C) 2005-2023 Shaun Purcell, Christopher Chang GNU General Public License v3
Logging to test_assoc.log.
Options in effect:
--1
--bfile cons_imp_data_sep_8th
--covar cons_imp_data_sep_8th_imp_generic.cov
--covar-name Sex
--glm
--out test_assoc
--pheno constipation_cc.phe
--pheno-name constipation_cc
Start time: Mon Sep 11 18:21:10 2023
7882 MiB RAM detected, ~7185 available; reserving 3941 MiB for main workspace.
Using up to 4 compute threads.
48352 samples (26548 females, 21804 males; 48352 founders) loaded from
cons_imp_data_sep_8th.fam.
827 variants loaded from cons_imp_data_sep_8th.bim.
Error: Unexpected PLINK 1 .bed file size (expected 9996779 bytes).
End time: Mon Sep 11 18:21:10 2023
We encounter an error message stating, 'Unexpected PLINK 1.bed file size (expected 9996779 bytes).'