N=11
p=c(0.03667409, 0.09299767, 0.29450436, 0.01733081, 0.04923809, 0.02984361)
si=c(1.69868880, 0.54292079, 0.58489574, 0.92918177, 0.48998484, 0.57004098)
taus=c(-0.52475006, 0.55768873, 0.18454559, 0.02233204, 0.57570506)
tau2par=tau2par.bvn
qcond=qcondbvn
out=rmultinom6dVineCopulaREMADA.norm(N,p,si,taus,qcond,tau2par)
y101=out[,1]
y011=out[,2]
y111=out[,3]
y001=out[,4]
y100=out[,5]
y010=out[,6]
y110=out[,7]
y000=out[,8]
Run the code above in your browser using DataLab