ymat <- rbifgmcop(1000, apar = rhobitlink(3, inverse = TRUE))
if (FALSE) plot(ymat, col = "blue")
fit <- vglm(ymat ~ 1, fam = bifgmcop, trace = TRUE)
coef(fit, matrix = TRUE)
Coef(fit)
head(fitted(fit))
Run the code above in your browser using DataLab