n <- 2500
x.use = rnorm(n)
w.use = runif(n,-1,1)
eta1.use = c(-2.2, 0.7)
eta2.use = c(-2, 0.8)
beta1.use = c(1,0.5)
beta2.use = c(1,1)
alpha1.use = 7
alpha2.use = 3
tau.use = c(-0.2, .3)
marginals.use <- c("ZINB", "ZIGA")
y.use <- scdeco.sim.cop(marginals=marginals.use, x=x.use,
eta1.true=eta1.use, eta2.true=eta2.use,
beta1.true=beta1.use, beta2.true=beta2.use,
alpha1.true=alpha1.use, alpha2.true=alpha2.use,
tau.true=tau.use, w=w.use)
y.use[1:10,]
Run the code above in your browser using DataLab