B1B2 component of the function pnbd.pmf.General defined in R/pnbd.R has changed. The fix is that the third parameter of the B2 function described by equation (24) here is now correctly set to c <- r + s + x + 1. The previous, incorrect version carried over the i component from the first term by setting c <- w + x + 1.hardie defaults to TRUE everywhere but in bgndb.generalParams, where it defaults to NULL