-
Notifications
You must be signed in to change notification settings - Fork 135
Closed
Description
Describe the bug
Running -lmm [1-3] gives the following output:
time ./bin/gemma-0.98.3-linux-static -lmm 1 -bfile example/mouse_hs1940 -k output/mouse_hs1940_data.cXX.txt -o lmm1
head -3 lmm*assoc*
==> lmm1.assoc.txt <==
chr rs ps n_miss allele1 allele0 af beta se logl_H1 l_remle p_wald
1 rs3683945 3197400 0 A G 0.443 -7.965330e-02 6.185099e-02 -1.582126e+03 4.318993e+00 1.980182e-01
1 rs3707673 3407393 0 G A 0.443 -6.654282e-02 6.210234e-02 -1.582377e+03 4.316144e+00 2.841271e-01
==> lmm2.assoc.txt <==
chr rs ps n_miss allele1 allele0 af logl_H1 l_mle p_lrt
1 rs3683945 3197400 0 A G 0.443 -1.583892e+03 4.313353e+00 1.977836e-01
1 rs3707673 3407393 0 G A 0.443 -1.584148e+03 4.310280e+00 2.840289e-01
==> lmm3.assoc.txt <==
chr rs ps n_miss allele1 allele0 af beta se logl_H1 p_score
1 rs3683945 3197400 0 A G 0.443 -7.965744e-02 6.191532e-02 0.000000e+00 1.984062e-01
1 rs3707673 3407393 0 G A 0.443 -6.650240e-02 6.217649e-02 0.000000e+00 2.848479e-01
At least one expects lmm3 to show all 3 calculations! Probably a regression on my end.
Reactions are currently unavailable