# NOT RUN {
ymat <- rbifgmcop(n = 1000, apar = rhobit(3, inverse = TRUE))
# }
# NOT RUN {
plot(ymat, col = "blue")
# }
# NOT RUN {
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