Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Parsing input file '' failed in function AnalyzeBimbamGXE #261

Open
MaxiEstravis opened this issue Feb 21, 2022 · 0 comments
Open

Parsing input file '' failed in function AnalyzeBimbamGXE #261

MaxiEstravis opened this issue Feb 21, 2022 · 0 comments

Comments

@MaxiEstravis
Copy link

MaxiEstravis commented Feb 21, 2022

Hi all, thanks for the great work in GEMMA.
I'm trying to do an analysis with a bim genotype file, a phenotype file, and a GxE file. This is the Linux command line:

~/tools/GEMMA/bin/gemma -g test.bim -p test_pheno.txt -n 1 -a test.anno -k ./output/test_matrix.cXX.txt -gxe test_2021_bio1.data -lmm 4 -o test_BranchQuality_2021_bio1_lmm

This gives the following log with an error:

GEMMA 0.98.5 (2021-08-25) by Xiang Zhou, Pjotr Prins and team (C) 2012-2021
Reading Files ... 
## number of total individuals = 3
## number of analyzed individuals = 3
## number of covariates = 1
## number of phenotypes = 1
## number of total SNPs/var        =  2046023
## number of analyzed SNPs         =  1177875
Start Eigen-Decomposition...
pve estimate =3.06394e-06
se(pve) =6.45794
**** FAILED: Parsing input file '' failed in function AnalyzeBimbamGXE in src/lmm.cpp at line 2329

Because of this error, I suspect that something is wrong with the GxE file, which looks like this:

3.3
4.1
2.9

According to the manual specification (Page 21):

This gxe file contains a column of environmental variables

Besides, this following line without invoking the GxE file works perfectly:

~/tools/GEMMA/bin/gemma -g test.bim -p test_pheno.txt -n 1 -a test.anno -k ./output/test_matrix.cXX.txt -lmm 4 -o test_BranchQuality_2021_bio1_lmm

I noticed there's a related issue that's planning to be incorporated in the new release (Issue #241), but does this mean that I can't incorporate GxE in my analysis for now? Or did I do something wrong?

Cheers,
Maximiliano.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants