set.seed(123)
n <- 100
p <- 3
q <- 2
b0 <- rep(1, p)
X <- matrix(rnorm(n * q), n, q)
B <- matrix(rnorm(q * p), q, p)
Sigma <- outer(1:p, 1:p, function(i, j) 0.3^abs(i - j))
probl <- 0.05
probr <- 0.05
probna <- 0.05
Z <- rcggm(n = n, b0 = b0, Sigma = Sigma, probl = probl, probr = probr,
probna = probna)
nobs(Z)
nresp(Z)
npred(Z)
Z <- rcggm(b0 = b0, X = X, B = B, Sigma = Sigma, probl = probl, probr = probr,
probna = probna)
nobs(Z)
nresp(Z)
npred(Z)
Run the code above in your browser using DataLab