-lmm 9 option leads to bug in mvlmm model

genetics-statistics / GEMMA

Genome-wide Efficient Mixed Model Association

https://github.com/genetics-statistics/GEMMA

GNU General Public License v3.0

333 stars 125 forks source link

-lmm 9 option leads to bug in mvlmm model #255

Open diogro opened 3 years ago

diogro commented 3 years ago

The mvlmm file is missing the case for the -lmm 9 option, leading to an output file with a single line and no significance values.

https://github.com/genetics-statistics/GEMMA/blob/5caf9487f3d69c43ce296345363c52b2bdcf0b7f/src/mvlmm.cpp#L192-L202

I stumbled on this because I'm also experiencing the long run times referenced in issue #243 when using -lmm 1. Option 9 runs, but with the mangled output.

pjotrp commented 3 years ago

Thanks for reporting!