Here is the story: GTEx V6 data is not in VCF format. When analyzing it, GSL SVD crashed on me with a complaint message (gsl: svd.c:286: ERROR: svd of MxN matrix, M<N, is not implemented) so I created this branch and tried to create a toy data to trouble-shoot. Yet the program crashed with segfault on the toy data which led to investigation and thus this patch. Now it works on my toy set but I am not sure if this gets rid of the GSL error (investigating). Creating this pull request because I think this is a critical bug that should best be fixed.
Please read the commit message here:
https://github.com/gaow/eqtlbma/commit/2c080acc2ea19c2070d67785caeaa53eb8a36aca
Here is the story: GTEx V6 data is not in VCF format. When analyzing it, GSL SVD crashed on me with a complaint message (gsl: svd.c:286: ERROR: svd of MxN matrix, M<N, is not implemented) so I created this branch and tried to create a toy data to trouble-shoot. Yet the program crashed with segfault on the toy data which led to investigation and thus this patch. Now it works on my toy set but I am not sure if this gets rid of the GSL error (investigating). Creating this pull request because I think this is a critical bug that should best be fixed.