```
# NOT RUN {
i.mix <- c(5, 10, 15, 20); a.mlm <- 13; mymux <- 2
goffset(mymux, 10, i.mix = i.mix, a.mlm = a.mlm)
# }
# NOT RUN {
org1 <- with(gdata, range(y)) # Original range of the data
vglm(mymux * y ~ 1,
offset = goffset(mymux, nrow(gdata), i.mix = i.mix, a.mlm = a.mlm),
gaitdpoisson(a.mlm = mymux * a.mlm, i.mix = mymux * i.mix,
truncate = Trunc(org1, mymux)),
data = gdata)
# }
```

Run the code above in your browser using DataLab