Question: "Core Dumped" when using GEMMA for GWAS
gravatar for amandinette
22 months 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 • 936 views
ADD COMMENTlink modified 12 weeks ago by morteza.mahmoudisaber70 • written 22 months ago by amandinette20
gravatar for r.harrington
21 months ago by
r.harrington30 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 21 months ago • written 21 months ago by r.harrington30
gravatar for morteza.mahmoudisaber
12 weeks 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 12 weeks 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: 852 users visited in the last hour