Question: SHAPEIT Segmentation Fault on Input Read-In
0
gravatar for nnn1
13 months ago by
nnn10
nnn10 wrote:

Hello all,

I am attempting to phase an individual's genome using the SHAPEIT software package and am running into errors. I installed the most recent binary (r900) on a linux machine and successfully ran the software on the included example files.

I am running the program with the 1000Genomes Phase 3 reference panel that the manual recommends (link). The input file is in VCF format and was split into chromosomes using tabix. The reference genome used to generate the input file is build 38. As the reference panel is build 38, I used liftOver to lift the input VCF down to hg19/b37.

Upon running the program with the chromosome vcf, reference panel files and genetic map, the following occurs:


Segmented HAPlotype Estimation & Imputation Tool * Authors : Olivier Delaneau, Jared O'Connell, Jean-François Zagury, Jonathan Marchini * Contact : send an email to the OXSTATGEN mail list https://www.jiscmail.ac.uk/cgi-bin/webadmin?A0=OXSTATGEN * Webpage : https://mathgen.stats.ox.ac.uk/shapeit * Version : v2.r904 * Date : 02/07/2018 15:20:55 * LOGfile : [gwas.log]

MODE -phase : PHASING GENOTYPE DATA * Autosome (chr1 ... chr22) * Window-based model (SHAPEIT v2) * Reference panel of haplotypes used * MCMC iteration

Parameters : * Seed : 1530559255 * Parallelisation: 1 threads * Ref allele is NOT aligned on the reference genome * MCMC: 35 iterations [7 B + 1 runs of 8 P + 20 M] * Model: 100 states per window [100 H + 0 PM + 0 R + 0 COV ] / Windows of ~2.0 Mb / Ne = 15000

Reading unphased data in [../../genomes/chr9.vcf] Segmentation fault (core dumped)


I'm not sure how to go about fixing this. Any help is much appreciated!

Thanks!

phasing genome • 770 views
ADD COMMENTlink written 13 months ago by nnn10

Hi,

core dumped seems a kind of memory issue. Have you checked RAM?

ADD REPLYlink written 13 months ago by Mehmet470

Thanks for your reply! How do you suggest going about checking the RAM? I'm not sure what to examine. As I understand it, a core dump is created whenever a program crashes and shows the process state at the time of the crash.

ADD REPLYlink written 13 months ago by nnn10

Can you increase memory that you allocate for this program?

ADD REPLYlink written 13 months ago by Mehmet470
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1120 users visited in the last hour