stuff/dip_N_do.R

In mmaechler/diptest: Hartigan's Dip Test Statistic for Unimodality - Corrected

### This is sourced from a file that sets  N , n.sim  , .Random.seed  ...
###
stopifnot(exists("N"),     is.numeric(N),     length(N) == 1, N == round(N),
          exists("n.sim"), is.numeric(n.sim), length(n.sim) == 1,
          n.sim == round(n.sim))
if(exists(".Random.seed")) {
    cat("using \n")
    dump(".Random.seed", "")
} else
    set.seed(12)
stopifnot(require("diptest"))
setwd("/u/maechler/R/Pkgs/diptest/stuff")

dd <- numeric(n.sim)
.pt <- proc.time()

for(i in 1:n.sim) {
    dd[i] <- dip(runif(N))
    if(i %% 100 == 0) {
        cat(".")
        if(i %% 10000 == 0) {
            cat("",i,"\n")
            .opt <- .pt ; .pt <- proc.time(); .pt - .opt
        }
    }
}

nam <- paste("dip", floor(N/1000),"k", sep='')
assign(nam, dd)
(outf <- paste(nam, "rda", sep="."))
save(list = c(nam, "N"), file = outf)

sessionInfo()