This paper makes three contributions. First, it uses copula functions to obtain a flexible bivariate parametric model for nonnegative integer-valued data (counts). Second, it recovers the distribution of the difference in the two counts from a specifed bivariate count distribution. Third, the...