| ## 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) |
| } |
| |
