genetics-statistics / GEMMA

Genome-wide Efficient Mixed Model Association
https://github.com/genetics-statistics/GEMMA
GNU General Public License v3.0
330 stars 124 forks source link

GSL ERROR: matrix is singular in lu.c at line 147 errno 1 #249

Open AgatheMaupetit opened 3 years ago

AgatheMaupetit commented 3 years ago

Hi !

Running a multivariate model with gemma 0.98.4: gemma -bfile $DATAWORK/genot -k $DATAWORK/K_Plim.txt -lmm 4 -n 12 41 -o mlmm_2traits I get this wellknown error: GSL ERROR: matrix is singular in lu.c at line 147 errno 1

But I do not get it using only one trait. So there is no issue with the data files.

I have 100 individuals and 98031 SNPs analyzed. Traits have a normal distribution and are not correlated as suggested in other issues.

Here is the output :

image

Do you have any idea to solve this problem ?

Thank you in advance for your help.

Agathe

Check list:

  1. [X] I have found an issue with GEMMA
  2. [ X] I have searched for it on the issue tracker (incl. closed issues)
  3. [X ] I have searched for it on the mailing list
  4. [X ] I have tried the latest release of GEMMA
  5. [X ] I have read and agreed to below code of conduct
  6. [x] If it is a support/install question I have posted it to the mailing list
  7. [x] If it is software development related I have posted a new issue on the issue tracker or added to an existing one
  8. [ X] In the message I have included the output of my GEMMA run
  9. [ ] In the message I have included the relevant .log.txt file in the output directory
  10. [ ] I have made available the data to reproduce the problem (optional)

To find bugs the GEMMA software developers may ask to install a development version of the software. They may also ask you for your data and will treat it confidentially. Please always remember that GEMMA is written and maintained by volunteers with good intentions. Our time is valuable too. By helping us as much as possible we can provide this tool for everyone to use.