Optimal design to discriminate between rival copula models for a bivariate binary response