Gemma linear mixed model analysis with second grm format
0
0
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.

Gemma gwas gsl software error • 2.1k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 1802 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6