PyRate icon indicating copy to clipboard operation
PyRate copied to clipboard

[RuntimeWarning: invalid value encountered in scalar divide] error when running RJMCMC

Open isovaline2230 opened this issue 1 year ago • 0 comments

Hi, I've been running PyRate v 3.1.3 using RJMCMC for my dataset. When running the following command,

python PyRate.py ./Data_PyRate.py -j 1 -A 4 -qShift ./epochs.txt -mG -pP 1.5 0 -log_sp_q_rates -edgeShift 144 0 -n 5000000 -s 500

MCMC stops with following warning message after some iterations:

RuntimeWarning: invalid value encountered in scalar divide p1 = (time_i1-t_prime)/(time_i1-time_i2) RuntimeWarning: invalid value encountered in scalar divide p2 = (t_prime-time_i2)/(time_i1-time_i2) RuntimeWarning: divide by zero encountered in log log_q_prob = log(abs(r_time))-prior_sym_beta(u,shape_beta_RJ) # prob latent parameters: Gamma pdf

I'm not sure but I think this issue happens when using -edgeShift option. Hope I could find some solutions. Thanks in advance.

isovaline2230 avatar Aug 31 '24 15:08 isovaline2230