This paper makes three contributions. Firstly, it uses copula functions to obtain a flexible bivariate parametric model for non-negative integer-valued data (counts). Secondly, it recovers the distribution of the difference in the two counts from a specified bivariate count distribution....