qvalue nan error
Opened this issue · 4 comments
tabix $VCF chr5:150219491-150290291 | rasqual -y $H3K4me3_counts_bin_file -k $H3K4me3_offsets_bin_file -s 150226191 -e 150226801 -j 30657 -l 922 -m 55 -n 14 -z -f "ACC_Plus_K4me3" -x $H3K4me3_covariates_bin_file -t
Qval nan! 0.500000 0.500000
gsl: gamma.c:1180: ERROR: error
Default GSL error handler invoked.
Aborted
Hi, could you upload the input files for this locus to check? It looks the likelihood became NaN at some point in the iterative steps of EM algorithm.
Hi, I am having the same issue for some, but not all, of the features I'm testing. Are there any updates on how to resolve this issue?
Hi, I just want to follow up on this issue. It turns out this issue sometimes happens depending on the dataset. Also, the program gets into an infinite loop without stoping and keeps printing
Qval nan! 0.500000 0.500000
to stdout.
Is there a way that I can share a minimum dataset for testing?
Thanks!
Just to want to help who may have the similar error.
I had the same error and fixed it by removing a column (with the same value in all samples) in -x covariates file.