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]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.