# NOT RUN {
set.seed(1)
# load the distance correlation for comparison
require(energy)
## independent case
x <- rnorm(100)
y <- rnorm(100)
z <- rnorm(100)
cdcor.ada(x,y,z)
dcor(x,y)
## conditional dependent case
x <- rnorm(100)
y <- x + 0.1*rnorm(100)
z <- rnorm(100)
cdcor.ada(x,y,z)
dcor(x,y)
## conditional independent case: x and z are correlated
require(MASS)
data <- mvrnorm(100,rep(0,2),matrix(c(1,0.8,0.8,1),2,2))
x <- data[,1]
z <- data[,2]
y <- z + 0.1*rnorm(100)
cdcor.ada(x,y,z)
dcor(x,y)
# }
Run the code above in your browser using DataLab