rqr_nbinom_logs function

compute numerically stable negative binomial rqr