# NOT RUN {
avec <- 5:6 # Alter these (special) values
ivec <- 9:10 # Inflate these (special) values
tvec <- 2:4 # Truncate these (special) values
pobs.a <- c(0.1, 0.2)
pstr.i <- (1:2)/10
size <- 10; prob <- 0.5; xx <- 0:11
y <- rgaitbinom.mlm(1000, size, prob, alter = avec, inflate = ivec,
truncate = tvec, pstr.i = pstr.i, pobs.a = pobs.a, byrow = TRUE)
table(y)
(ii <- dgaitbinom.mlm(xx, size, prob, alter = avec, inflate = ivec,
truncate = tvec, pstr.i = pstr.i, pobs.a = pobs.a, byrow = TRUE))
# }
Run the code above in your browser using DataLab