library(randtoolbox) #### outputs #### print(halton(1, 5)) options(digits=15) n <- 100 n <- 5 cbind( 1/get.primes(n), as.vector(halton(1, n) ) ) #### n argument #### try(halton(-1)) #### dim argument #### try(halton(1, 0)) #### init argument #### halton(n) randtoolbox:::.getrandtoolboxEnv(".halton.seed") halton(n, init=TRUE) randtoolbox:::.getrandtoolboxEnv(".halton.seed") halton(n, init=FALSE) try(halton(5, init="a")) #### mixed argument #### try(halton(1, mixed=3)) halton(n, mixed=TRUE) #### usetime argument #### #QMC with time machine start halton(n, usetime=TRUE) halton(n, usetime=TRUE) #hybrid QMC with SFMT : test continuing the sequence (bug reported by Hiroyuki Kawakatsu) setSeed(1234); halton(n, init=TRUE, mixed=TRUE, mexp=607) halton(n, init=FALSE, mixed=TRUE) setSeed(1234); halton(n, init=TRUE, mixed=TRUE, mexp=607) halton(n, init=FALSE, mixed=TRUE) #### method argument #### n <- 5 d <- 4 halton(n, d, method="C") randtoolbox:::.getrandtoolboxEnv(".halton.seed") halton(n, d, method="C", start=5) #### normal argument #### halton(n, normal=TRUE) try(halton(3, normal=1)) #### mexp argument #### halton(n, mexp=607) try(halton(3, mexp=3)) try(halton(3, mexp="3")) #### start argument #### try(halton(3, start="3")) halton(n, start=0) halton(n, start=1)