This paper proposes and implements simulated maximum likelihood estimation of bivariate count models with unrestricted correlation pattern of unobserved heterogeneity. The implementation incorporates both antithetic acceleration and adjustment for first-order simulation bias. Both the Monte...