Entering edit mode
7.4 years ago
dp0b
▴
80
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.
did you end up figuring out what caused this error on yours?