[RuntimeWarning: invalid value encountered in scalar divide] error when running RJMCMC
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.