The R Project SVN R

Rev

Rev 63512 | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed 

## reproducible examples.

library(parallel)

if(require(boot)) {
    run1 <- function(...) {
        library(boot)
        cd4.rg <- function(data, mle) MASS::mvrnorm(nrow(data), mle$m, mle$v)
        cd4.mle <- list(m = colMeans(cd4), v = var(cd4))
        boot(cd4, corr, R = 500, sim = "parametric",
             ran.gen = cd4.rg, mle = cd4.mle)
    }
    cl <- makeCluster(mc <- 2)
    clusterSetRNGStream(cl, 123)
    cd4.boot <- do.call(c, parLapply(cl, seq_len(mc), run1))
    print(boot.ci(cd4.boot,  type = c("norm", "basic", "perc"),
                  conf = 0.9, h = atanh, hinv = tanh))
    stopCluster(cl)
}