Modelling Frailty for Bivariate Data