I'm trying to turn a tped & tfam file to bfile format using PLINK 1.9.
My tped is set out like this:
1 AX-XXX 0 1000 A T A A
My tfam something like this:
F1 001 0 0 1 1
I ran:
plink --tfile INPUT --autosome-num 31 --allow-extra-chr --keep-allele-order --debug --make-bed --out OUTPUT
And the output I get is just an OUTPUT-temporary.bed .bim and .fam file. I've run the --debug flag, and the log I get is as so:
Processing .tped file.
OUTPUT-temporary.bed + OUTPUT-temporary.bim + OUTPUT-temporary.fam written.
And everything stops there. I've tried increasing the memory (it reserved 257785 MB in my last run), I've made sure all missing calls are set to 0, but I get no actual error message at all so it's really hard to know whether it's an issue with my files or the commands I'm using or something else.
I also tried --recode and --recode vcd just to see if this issue happened with other commands, and the same thing happens.
Any ideas much appreciated.
Yes:
And that is all I get.
Edit: I have successfully imported a VCF from another dataset to PLINK binary format in the last few days with this same version, so I feel like it's probably an issue with the tfiles, but I had some others from my group look them over and no one can spot an obvious issue with them.
I'm unable to replicate this issue. Is it possible for you to send me a .tped+.tfam which triggers this behavior? If not, I can send you a series of debug builds to try to nail down what's happening, but that process may take much longer.
Let me check whether I have permission to forward the data on - I'll get right back to you.
Yes, I can send you this on. What is the best email? Thanks for your help with this!
[my username] at gmail
Thank you, much appreciated!
This is definitely a plink bug (a segfault occurred), working on a fix now.
Bugfix is now posted. This will now print a regular error message on your data; you need to use --chr-set (which permits chrY) instead of --autosome-num.
That's amazing - thank you very much!