search

timflutre / eqtlbma

Package to detect eQTLs jointly in multiple subgroups (e.g. tissues) via Bayesian Model Averaging.
http://www.plosgenetics.org/article/info%3Adoi%2F10.1371%2Fjournal.pgen.1003486
GNU General Public License v3.0
22 stars 12 forks source link

ERROR: new_config_prior_[3] is NaN ! #29

Open sunrain1234 opened 1 year ago

sunrain1234 commented 1 year ago

Hi, When I run eqtlbma_hm, I got this error. I wonder what caused this. Is this due to nan being included in the output of the eqtlbma_bf? Here is commond and error. 

zless out_eqtlbma_l10abfs_raw.txt.gz |head
gene    snp     config  l10abf.grid1    l10abf.grid2    l10abf.grid3    l10abf.grid4    l10abf.grid5    l10abf.grid6    l10abf.grid7    l10abf.grid8    l10abf.grid9    l10abf.grid10   l10abf.grid11  l10abf.grid12    l10abf.grid13   l10abf.grid14   l10abf.grid15   l10abf.grid16   l10abf.grid17   l10abf.grid18   l10abf.grid19   l10abf.grid20   l10abf.grid21   l10abf.grid22   l10abf.grid23   l10abf.grid24   l10abf.grid25
7SK     chr1_119228079_G_A      gen     -2.552074e-01   -2.567767e-01   -2.533767e-01   -2.456162e-01   -2.354144e-01   -6.356711e-01   -6.296440e-01   -6.055082e-01   -5.571040e-01   -4.747932e-01  -1.160806e+00    -1.149360e+00   -1.109070e+00   -1.018222e+00   -7.573831e-01   -1.741525e+00   -1.728199e+00   -1.681967e+00   -1.572422e+00   -1.053561e+00   -2.338098e+00   -2.324259e+00   -2.276372e+00   -2.161275e+00   -1.353362e+00
7SK     chr1_119228079_G_A      gen-fix -2.354144e-01   -2.354144e-01   -2.354144e-01   -2.354144e-01   -2.354144e-01   -4.747932e-01   -4.747932e-01   -4.747932e-01   -4.747932e-01   -4.747932e-01  -7.573831e-01    -7.573831e-01   -7.573831e-01   -7.573831e-01   -7.573831e-01   -1.053561e+00   -1.053561e+00   -1.053561e+00   -1.053561e+00   -1.053561e+00   -1.353362e+00   -1.353362e+00   -1.353362e+00   -1.353362e+00   -1.353362e+00
7SK     chr1_119228079_G_A      gen-maxh        -2.552074e-01   -2.552074e-01   -2.552074e-01   -2.552074e-01   -2.552074e-01   -6.356711e-01   -6.356711e-01   -6.356711e-01   -6.356711e-01   -6.356711e-01   -1.160806e+00   -1.160806e+00   -1.160806e+00   -1.160806e+00   -1.160806e+00   -1.741525e+00   -1.741525e+00   -1.741525e+00   -1.741525e+00   -1.741525e+00   -2.338098e+00   -2.338098e+00  -2.338098e+00    -2.338098e+00   -2.338098e+00
7SK     chr1_119228079_G_A      1       -1.811154e-01   -1.811154e-01   -3.965886e-01   -3.965886e-01   -6.696843e-01   -6.696843e-01   -9.631312e-01   -9.631312e-01   -1.262223e+00   -1.262223e+00  nan      nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan
7SK     chr1_119228079_G_A      2       -7.409202e-02   -7.409202e-02   -2.390825e-01   -2.390825e-01   -4.911213e-01   -4.911213e-01   -7.783937e-01   -7.783937e-01   -1.075875e+00   -1.075875e+00  nan      nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan
7SK     chr1_119228079_G_A      1-2     -2.456162e-01   -2.354144e-01   -5.571040e-01   -4.747932e-01   -1.018222e+00   -7.573831e-01   -1.572422e+00   -1.053561e+00   -2.161275e+00   -1.353362e+00  nan      nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan     nan
7SK     chr1_119230040_T_C      gen     -1.723150e-01   -1.643266e-01   -1.528507e-01   -1.376136e-01   -1.194407e-01   -5.136785e-01   -5.026857e-01   -4.744328e-01   -4.214989e-01   -3.291333e-01  -1.022504e+00    -1.009488e+00   -9.678831e-01   -8.753246e-01   -6.017647e-01   -1.598441e+00   -1.584697e+00   -1.538113e+00   -1.428078e+00   -8.952365e-01   -2.193767e+00   -2.179821e+00   -2.131845e+00   -2.016620e+00   -1.194346e+00
7SK     chr1_119230040_T_C      gen-fix -1.194407e-01   -1.194407e-01   -1.194407e-01   -1.194407e-01   -1.194407e-01   -3.291333e-01   -3.291333e-01   -3.291333e-01   -3.291333e-01   -3.291333e-01  -6.017647e-01    -6.017647e-01   -6.017647e-01   -6.017647e-01   -6.017647e-01   -8.952365e-01   -8.952365e-01   -8.952365e-01   -8.952365e-01   -8.952365e-01   -1.194346e+00   -1.194346e+00   -1.194346e+00   -1.194346e+00   -1.194346e+00
7SK     chr1_119230040_T_C      gen-maxh        -1.723150e-01   -1.723150e-01   -1.723150e-01   -1.723150e-01   -1.723150e-01   -5.136785e-01   -5.136785e-01   -5.136785e-01   -5.136785e-01   -5.136785e-01   -1.022504e+00   -1.022504e+00   -1.022504e+00   -1.022504e+00   -1.022504e+00   -1.598441e+00   -1.598441e+00   -1.598441e+00   -1.598441e+00   -1.598441e+00   -2.193767e+00   -2.193767e+00  -2.193767e+00    -2.193767e+00   -2.193767e+00
cmd-line: /data/wjxu/software/eqtlbma-1.3.3/src/.libs/lt-eqtlbma_hm --data out_eqtlbma_l10abfs.txt.gz --nsubgrp 3 --dim 7 --ngrid 10 --out test_eqtlbma_hm.txt.gz
load data ...
nb of input files: 1
==================================================100.00%
finish loading 35252 genes and 174735745 gene-snp pairs (4065.900000 sec, 678283676 kB)
run EM algorithm (classic) ...
iter    0  loglik 18710.126792  pi0 5.0000e-01  configs 1.4286e-01 1.4286e-01 1.4286e-01 1.4286e-01 1.4286e-01 1.4286e-01 1.4286e-01  grid-points 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01 1.0000e-01
ERROR: new_config_prior_[3] is NaN
terminate called after throwing an instance of 'int'
timflutre commented 1 year ago

Hi, It's been a (very) long time since I used this software and it will be hard to find the cause of the error without a small, reproducible example. Moreover, I am too busy right now to do it, sorry. Tim