data(trypanos)
library(npmlreg)
(trypanos.np1 <- alldist(cbind(y, (n - y)) ~ log(x),
random = ~1, data = trypanos, family = binomial,
plot.opt = 0, verbose = FALSE,k=1))
(trypanos.np2 <- update(trypanos.np1,k=2))
Run the code above in your browser using DataLab