Sigma=matrix(c(1,.5,.5,1),ncol=2)
L=t(chol(Sigma))
trunpt=c(0,0,1,1)
above=c(1,1)
# drawn by Halton sequence
ghkvec(L,trunpt,above,r=100)
# use prime number 11 and 13
ghkvec(L,trunpt,above,r=100,HALTON=TRUE,pn=c(11,13))
# drawn by R::runif
ghkvec(L,trunpt,above,r=100,HALTON=FALSE)
Run the code above in your browser using DataLab