Predicting bivariate binary rare events responses using generalised extreme value regression model and copula function