This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Gemma linear mixed model analysis with second grm format

Hi

I have been trying to run a univariate linear mixed model using GEMMA. However I was hoping to use a grm I have made else where- it is the second in grm format. The first few lines are as as follows

1       1       8.923316e-01
1       2       -5.544370e-02
1       3       1.903249e-02
1       4       1.880089e-02
1       5       9.801378e-02
1       6       -3.117359e-02
1       7       2.540684e-01
1       8       -1.174230e-01
1       9       1.168132e-01
1       10      1.833550e-01

Then I am trying to run a lmm using the follow command

./gemma.linux -bfile mouse_hs1940 -km 2 -k mouse_plinkgctasorted.txt -lmm 2 -o mouse_plinkgctasorted

However I keep getting the following error. Any advice on how to fix this?

Reading Files ...
## number of total individuals = 1940
## number of analyzed individuals = 1410
## number of covariates = 1
## number of phenotypes = 1
## number of total SNPs = 12226
## number of analyzed SNPs = 10768
Start Eigen-Decomposition...
gsl: newton.c:74: ERROR: derivative is zero
Default GSL error handler invoked.
Aborted

The file will run the analysis when I supply a gemma grm, however the file I intend on supplying will be much larger and would prefer to make the grm elsewhere.

gemma gwas gsl software error

did you end up figuring out what caused this error on yours?

0 answers

No answers yet.

Log in to answer this question.