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

Including -gxe produces error: GSL ERROR: function value is not finite in brent.c at line 58 errno 9 #225

Closed DeWitP closed 4 years ago

DeWitP commented 4 years ago

Hi! Thanks for nice software! I am trying to run GEMMA with the -gxe option. Without the option, the software runs fine, but when I include the -gxe option I get the rror: "GSL ERROR: function value is not finite in brent.c at line 58 errno 9"

I cannot find any information on the formatting of the gxe file, so perhaps my file is wrongly formatted. The GEMMA manual only states: "This gxe file contains a column of environmental variables" so I only include one column of an environmental factor (salinity in my case) in the file, in the same order as the samples are listed in the .fam file (see attached). Do I need to include something more? Thanks!

S_gxe.txt

DeWitP commented 4 years ago

Hi again, nevermind. I added also the sample names in the gxe file and it worked. You might want to clarify this in the manual though. All the best, Pierre

pjotrp commented 4 years ago

You can send a pull request. That would be appreciated!