set.seed(234)
z=runif(10,2,11)# z is independently created
x=sample(1:10)+z/10 #x is partly indep and partly affected by z
y=1+2*x+3*z+rnorm(10)# y depends on x and z not vice versa
mtx=cbind(x,y,z)
parcorVecH2(mtx)
if (FALSE) {
set.seed(34);mtx=matrix(sample(1:600)[1:80],ncol=4)
colnames(mtx)=c('V1', 'v2', 'V3', 'V4')
parcorVecH2(mtx,verbo=TRUE, idep=2)
}
Run the code above in your browser using DataLab