miriammiyagi / QuIBL

Quantifying Introgression via Branch Lengths
49 stars 7 forks source link

ValueError: math domain error #32

Open xiaolongtu opened 5 months ago

xiaolongtu commented 5 months ago

@miriammiyagi @nbedelman hi, I meet this problem when running QuIBL.py, could you help me to solve it

/home/tuxl/software/QuIBL-master/QuIBL.py:119: RuntimeWarning: invalid value encountered in double_scalars xq=[l/sum(xq) for l in xq]#xq/sum(xq) /home/tuxl/software/QuIBL-master/QuIBL.py:109: RuntimeWarning: overflow encountered in exp y=np.sinh(x/lmbd)/(lmbd(-1+np.exp(C))) /home/tuxl/software/QuIBL-master/QuIBL.py:107: RuntimeWarning: overflow encountered in exp y=0.5/lmbdnp.exp(-x/lmbd)(1+np.exp(C)) /home/tuxl/software/QuIBL-master/QuIBL.py:107: RuntimeWarning: invalid value encountered in double_scalars y=0.5/lmbdnp.exp(-x/lmbd)*(1+np.exp(C)) Traceback (most recent call last): File "/home/tuxl/software/QuIBL-master/QuIBL.py", line 338, in inputReader(sys.argv[1]) File "/home/tuxl/software/QuIBL-master/QuIBL.py", line 334, in inputReader outputFormatter(outputDict,inputDict) File "/home/tuxl/software/QuIBL-master/QuIBL.py", line 305, in outputFormatter tripletSet=Parallel(n_jobs=num_cores)(delayed(PLexMax)(triple,K,lthresh,numsteps,gAScalar) for triple in trees) File "/home/tuxl/software/anaconda3/envs/R/lib/python2.7/site-packages/joblib/parallel.py", line 1017, in call self.retrieve() File "/home/tuxl/software/anaconda3/envs/R/lib/python2.7/site-packages/joblib/parallel.py", line 938, in retrieve raise exception.unwrap(this_report) joblib.my_exceptions.JoblibValueError: JoblibValueError


........................................................................... /home/tuxl/software/QuIBL-master/QuIBL.py in () 333 outputDict[flag]=config.get('Output',flag) 334 outputFormatter(outputDict,inputDict) 335 336 337 #print exMax(getTripBranches(readin_Newick(sys.argv[1])), 2, 0.01, 5, 0.5)[0].models --> 338 inputReader(sys.argv[1]) 339 340 341 342

...........................................................................