y = rsinmad(n=2000, a=1, 6, 2) # Not genuine data!
fit = vglm(y ~ 1, betaII, trace=TRUE)
fit = vglm(y ~ 1, betaII(init.p=0.7, init.q=0.7), trace=TRUE, crit="c")
coef(fit, mat=TRUE)
Coef(fit)
summary(fit)
Run the code above in your browser using DataLab