francois-a / fastqtl

Enhanced version of the FastQTL QTL mapper
GNU General Public License v3.0
56 stars 22 forks source link

fastQTL no longer normalizes the phenotypes #11

Closed aarbduarte closed 4 years ago

aarbduarte commented 4 years ago

Hello, I've seen that you removed the option to normalize the phenotypes as was doable by using --normal flag. Was it for a particular motive? Is there any work around? Thank you in advance.

francois-a commented 4 years ago

Hi — No, that option is still there: bin/fastQTL.static --help

aarbduarte commented 4 years ago

additionally when running on nominal pass it finishes alright but when doing permutation it throws an error: image

francois-a commented 4 years ago

Can you please run the command (at the bottom) that is throwing the error?

aarbduarte commented 4 years ago

image

francois-a commented 4 years ago

I meant the command that starts with /home/.../fastQTL

aarbduarte commented 4 years ago
Processing gene [AFE_1_+_43531090_43528756_43538198_Hypothetical]
  * Number of variants in cis = 780
  * Best correlation = 0.2902
  * Number of permutations = 1000 / 1203
gsl: beta.c:44: ERROR: domain error
Default GSL error handler invoked.
Aborted (core dumped)
francois-a commented 4 years ago

That error usually arises when the phenotype lacks complexity (0 or same value in most samples). I'm no longer making changes to fastQTL, but you can either try filtering out such phenotypes, or using tensorQTL, which should be more robust to this.

aarbduarte commented 4 years ago

thank you for your help!