Plink 1.9 "File read failure" error when merging with --bmerge
0
2
Entering edit mode
2.5 years ago
svendula ▴ 20

When trying to merge two bfiles, I keep getting the "File read failure" error just after the merging starts, even though the files were obviously read and recognised. Any suggestions what I might be doing wrong? See below.

PLINK v1.90b6.6 64-bit (10 Oct 2018)           www.cog-genomics.org/plink/1.9/
(C) 2005-2018 Shaun Purcell, Christopher Chang   GNU General Public License v3

Logging to file1_file2_merged.log.
Options in effect:

  --allow-no-sex
  --bfile file1
  --bmerge file2
  --keep-allele-order
  --make-bed
  --out file1_file2_merged

128823 MB RAM detected; reserving 64411 MB for main workspace.
399 people loaded from file1.fam.
1239 people to be merged from file2.fam.
Of these, 1239 are new, while 0 are present in the base dataset.
580127 markers loaded from file1.bim.
562582 markers to be merged from file2.bim.
Of these, 155 are new, while 562427 are present in the base dataset.
Performing single-pass merge (1638 people, 580282 variants).

Pass 1: fileset #1 complete.
Error: File read failure.
plink gwas error • 1.8k views
ADD COMMENT
0
Entering edit mode

Please use the formatting bar (especially the code option) to present your post better. You can use backticks for inline code (`text` becomes text), or select a chunk of text and use the highlighted button to format it as a code block. I've done it for you this time.
code_formatting

ADD REPLY
0
Entering edit mode

Any news on this, @svendula? I'm facing the same issue.

ADD REPLY
0
Entering edit mode

same problem here

ADD REPLY
0
Entering edit mode

Can you post a set of files I can use to replicate what you’re seeing?

ADD REPLY
0
Entering edit mode

I'm trying to merge chr10-12 to chr1. and Here is my log file looks like:

    PLINK v1.90b3.42 64-bit (20 Sep 2016)
Options in effect:
  --bed ukb_cal_chr1_v2.bed
  --bim ukb_snp_chr1_v2.bim
  --fam ukb41999_cal_chr1_v2_s488264.fam
  --make-bed
  --merge-list list_beds.txt
  --out ukb_cal_allChrs_sample

Hostname: nodea1401
Working directory: /gpfs/group/home/yucheng/genotype_beds
Start time: Fri Jun 18 22:20:39 2021

Random number seed: 1624080039
96602 MB RAM detected; reserving 48301 MB for main workspace.
Performing single-pass merge (488377 people, 179921 variants).

Error: File read failure.

and my list_beds.txt looks like:

ukb_cal_chr10_v2.bed ukb_snp_chr10_v2.bim ukb41999_cal_chr10_v2_s488264.fam
ukb_cal_chr11_v2.bed ukb_snp_chr11_v2.bim ukb41999_cal_chr11_v2_s488264.fam
ukb_cal_chr12_v2.bed ukb_snp_chr12_v2.bim ukb41999_cal_chr12_v2_s488264.fam
ADD REPLY
1
Entering edit mode

You should be able to use plink 2.0 --pmerge-list for this. If it fails, it should provide a more useful error message.

ADD REPLY
0
Entering edit mode

Yeah now I'm able to run it with plink 2.0! Thank you so much!!

ADD REPLY

Login before adding your answer.

Traffic: 2345 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6