ydata <- data.frame(x2 = runif(nn <- 1000))
ydata <- transform(ydata, y = ryules(nn, rho = exp(1.5 - x2)))
with(ydata, table(y))
fit <- vglm(y ~ x2, yulesimon, ydata, trace = TRUE)
coef(fit, matrix = TRUE)
summary(fit)
Run the code above in your browser using DataLab