Dear sir, Your work on the Zelig package is higly appreciated. It has helped me alot. I'm trying to fit a bivariate probit model using zelig. I'm aware that I can using rhobit transformation on the 3rd coefficient from the fitted model. Please, how can I switch the link function to the identity link directly at the modelling stage just like the one in vglm function (family = binom2.rho(lrho = "identity"). Thank you for your time. Kind regards Segun