# NOT RUN {
set.seed(0)
n <- 100 # sample size
p <- 10 # dimensionality
rho <- 0.5 # the correlation between pairwise predictors.
Sigma <- matrix(rho, p, p)
diag(Sigma) <- 1
require(MASS)
x <- mvrnorm(n,rep(0,p), Sigma)
y <- x[,2] + x[,3] + rnorm(n)
z <- x[,1]
cdcsis(x,y,z,2)
# }
Run the code above in your browser using DataLab