pdat = data.frame(y = rparalogistic(n=3000, 4, 6))
fit = vglm(y ~ 1, paralogistic, pdat, trace=TRUE)
fit = vglm(y ~ 1, paralogistic(init.a=2.3, init.sc=5),
pdat, trace=TRUE, crit="c")
coef(fit, mat=TRUE)
Coef(fit)
summary(fit)
Run the code above in your browser using DataLab