Question: "Core Dumped" when using GEMMA for GWAS
gravatar for amandinette
2.9 years ago by
amandinette20 wrote:

Dear all,

I am following this tutorial about bacterial GWAS using GEMMA and Bugwas :

When I try correct for population structure, I obtain a segmentation fault.

Rscript bacterialGWAS/kmerGWAS/kmerLMM.R -chisqStat fus300.gwaskmer-out.chisqStat.txt -patternKey fus300.gwaskmer-out.patternKey.txt -patternIndex fus300.gwaskmer-out.patternIndex.txt -signif 5000 -relateMatrix ~/data/fus300_gemma_relmatrixout.cXX.txt -phenotype data/fus300_bugwas.pheno.txt -prefix fus300 -externalSoftware dependencies.txt

Running kmerLMM.R (...) Produce a matrix to fill to create the GEMMA gen file: Create GEMMA snp file: # Find the patterns which each kmer corresponds to out of those used in the LMM: (...) Reading Files ... Segmentation fault (core dumped) Error in file(file, "rt") : impossible d'ouvrir la connexion Calls: read.table -> file De plus : Warning message: In file(file, "rt") : impossible d'ouvrir le fichier './output/fus300_kmer_lmmout.assoc.txt' : Aucun fichier ou dossier de ce type Exécution arrêtée

I can't understand why... Can someone help me ?

Thank you !


gwas bugwas gemma • 1.3k views
ADD COMMENTlink modified 15 months ago by morteza.mahmoudisaber70 • written 2.9 years ago by amandinette20
gravatar for rnnh
2.8 years ago by
rnnh30 wrote:

"Segmentation fault" means that you tried to access memory that you do not have access to. You need more memory to run your job.

Try requesting more memory using qsh, qlogin or qrsh. For example, the following command will request "X" MB:

$ qlogin -l h_vmem=XM
ADD COMMENTlink modified 2.8 years ago • written 2.8 years ago by rnnh30
gravatar for morteza.mahmoudisaber
15 months ago by
morteza.mahmoudisaber70 wrote:

Actually in this case, segmentation fault is not due to memory error. BugWas uses GEMMA which accepts genotype and phenotype in BIMBAM format. For the phenotype it should be just one column without any headers, however, in the provided phenotype it is two columns with header. Delete the ID column and the header and it works fine.

ADD COMMENTlink written 15 months ago by morteza.mahmoudisaber70
Please log in to add an answer.


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