Comments (4)
- Check whether the latest stable build also segfaults here.
- If it does, check whether the --make-bed component fails in isolation, with --bfile [plink-temporary file] --allow-extra-chr --make-bed after the first crash; this should indicate whether the bug is in the --bcf or the --make-bed code.
from plink-ng.
Latest version has the same result:
$ time ~/src/plink1.90b2t/plink --make-bed --bcf ALL.phase3_shapeit2_mvncall_integrated_v5_extra_anno.20130502.genotypes.bcf --out ALL.phase3_shapeit2_mvncall_integrated_v5_extra_anno.20130502.genotypes.bcf.plink --allow-extra-chr
PLINK v1.90b2t 64-bit (20 Dec 2014) https://www.cog-genomics.org/plink2
(C) 2005-2014 Shaun Purcell, Christopher Chang GNU General Public License v3
Logging to ALL.phase3_shapeit2_mvncall_integrated_v5_extra_anno.20130502.genotypes.bcf.plink.log.
32054 MB RAM detected; reserving 16027 MB for main workspace.
--bcf: 84739k variants complete.
ALL.phase3_shapeit2_mvncall_integrated_v5_extra_anno.20130502.genotypes.bcf.plink-temporary.bed
+
ALL.phase3_shapeit2_mvncall_integrated_v5_extra_anno.20130502.genotypes.bcf.plink-temporary.bim
+
ALL.phase3_shapeit2_mvncall_integrated_v5_extra_anno.20130502.genotypes.bcf.plink-temporary.fam
written.
84739846 variants loaded from .bim file.
2504 people (0 males, 0 females, 2504 ambiguous) loaded from .fam.
Ambiguous sex IDs written to
ALL.phase3_shapeit2_mvncall_integrated_v5_extra_anno.20130502.genotypes.bcf.plink.nosex
.
Using 1 thread (no multithreaded calculations invoked).
Before main variant filters, 2504 founders and 0 nonfounders present.
Calculating allele frequencies... done.
Total genotyping rate is 0.980324.
84739846 variants and 2504 people pass filters and QC.
Note: No phenotypes present.
--make-bed to
ALL.phase3_shapeit2_mvncall_integrated_v5_extra_anno.20130502.genotypes.bcf.plink.bed
+
ALL.phase3_shapeit2_mvncall_integrated_v5_extra_anno.20130502.genotypes.bcf.plink.bim
+
ALL.phase3_shapeit2_mvncall_integrated_v5_extra_anno.20130502.genotypes.bcf.plink.fam
Segmentation fault (core dumped)
real 43m27.661s
user 37m54.756s
sys 1m6.058s
When I try your suggestion to run --make-bed
using --bfile
on the tmp file I get another seg fault:
$ time ~/src/plink1.90b2t/plink --make-bed --bfile ALL.phase3_shapeit2_mvncall_integrated_v5_extra_anno.20130502.genotypes.bcf.plink-temporary --out plink.test --allow-extra-chr
PLINK v1.90b2t 64-bit (20 Dec 2014) https://www.cog-genomics.org/plink2
(C) 2005-2014 Shaun Purcell, Christopher Chang GNU General Public License v3
Logging to plink.test.log.
32054 MB RAM detected; reserving 16027 MB for main workspace.
84739846 variants loaded from .bim file.
2504 people (0 males, 0 females, 2504 ambiguous) loaded from .fam.
Ambiguous sex IDs written to plink.test.nosex .
Using 1 thread (no multithreaded calculations invoked).
Before main variant filters, 2504 founders and 0 nonfounders present.
Calculating allele frequencies... done.
Total genotyping rate is 0.980324.
84739846 variants and 2504 people pass filters and QC.
Note: No phenotypes present.
Segmentation fault (core dumped)
real 3m3.203s
user 0m40.631s
sys 0m13.533s
from plink-ng.
Thanks, I will try to reproduce and fix this tonight, and let you know if I need any more information to do so.
from plink-ng.
The January 8 development build should fix this; let me know if you still have any problems. (This was actually supposed to crash, but with a "sorting of files too large to fit in RAM has not been implemented yet" error message, rather than a segfault. The multipass sorting code was already 90% written, though, so I went ahead and added the last 10%.)
from plink-ng.
Related Issues (20)
- Feature Request: Upstream Dockerfile
- CentOs7.9 install plink1.9 error HOT 3
- how to install in arm64 ? thks HOT 4
- Request: plink2 --extract/exclude-if-info [key(s)...] HOT 1
- Compilation failure HOT 1
- Can't run anny PLINK command on Ubuntu 22.04 HOT 1
- File read error when converting .vcf to PLINK .ped and .map file HOT 1
- Request for pgen files to include a version number in the header HOT 8
- can plink handle vcf files with * notation? HOT 1
- Pgenlib 0.90.1 is not installable from PyPI with Python 3.12.2 on Mac aarch64 HOT 1
- Logistic regression returns NAs output---FIRTH_CONVERGE_FAIL HOT 1
- plink2 --update-name issue HOT 1
- Obtaining different results from the same data in different formats HOT 2
- PLINK v1.9 `sprintf()` compilation warnings: Does a way to silence/fix them exist? HOT 3
- Persistent download URLs for built binaries? HOT 2
- pgenlibr: ReadList question HOT 7
- Problem with combining different samples into a single plink 1 file for later input HOT 1
- --homozyg-group question HOT 1
- clumping
- Error: Duplicate phenotype/covariate ID HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from plink-ng.