set.seed(123)
nn = 1000
y = rbeta(nn, shape1=1, shape2=3)
# parameters are estimated on a log scale
fit = vglm(y ~ 1, betaff, tr=TRUE, crit="c") # intercept-only model
coef(fit, matrix=TRUE) # log scale
Coef(fit) # On the original scale
Run the code above in your browser using DataLab