Skip to content

LMM output missing something? #237

@pjotrp

Description

@pjotrp

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.

Metadata

Metadata

Assignees

Labels

Type

No type

Projects

No projects

Relationships

None yet

Development

No branches or pull requests

Issue actions