data(datocc)
str(datocc)
## Not run:
# ## simulation
# n <- 1000
# set.seed(1234)
# x1 <- runif(n, -1, 1)
# x2 <- as.factor(rbinom(n, 1, 0.5))
# x3 <- rnorm(n)
# x4 <- rnorm(n)
# beta <- c(0.6, 0.5)
# theta <- c(0.4, -0.5, 0.3)
# X <- model.matrix(~ x1)
# Z <- model.matrix(~ x1 + x3)
# mu <- drop(X %*% beta)
# nu <- drop(Z %*% theta)
# p.occ <- binomial("cloglog")$linkinv(mu)
# p.det <- binomial("logit")$linkinv(nu)
# Y <- rbinom(n, 1, p.occ)
# W <- rbinom(n, 1, Y * p.det)
# datocc <- data.frame(Y, W, x1, x2, x3, x4, p.occ, p.det)
# ## End(Not run)
Run the code above in your browser using DataLab