Dear all, I am estimating a bivariate probit model using the package VGAM: fit1 = vglm(cbind(y1, y2) ~ x1+x2, binom2.rho, data=dat, trace=TRUE) I would like to estimate this via the method vglm.fit (in an analogous way as lm.fit for lm or glm.fit for glm). Unfortunately my trials did not succeed. Therefore I would like to ask if this is possible at all and - if yes - how one can access this function. Thank you very much for your efforts in advance! Best, Martin [[alternative HTML version deleted]]