Now computes log(lr) when p.in > q.in and returns log(lr)=0 when not