Rev 67901 | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed
R version 3.1.3 RC (2015-03-01 r67921) -- "Smooth Sidewalk"Copyright (C) 2015 The R Foundation for Statistical ComputingPlatform: x86_64-unknown-linux-gnu (64-bit)R is free software and comes with ABSOLUTELY NO WARRANTY.You are welcome to redistribute it under certain conditions.Type 'license()' or 'licence()' for distribution details.R is a collaborative project with many contributors.Type 'contributors()' for more information and'citation()' on how to cite R or R packages in publications.Type 'demo()' for some demos, 'help()' for on-line help, or'help.start()' for an HTML browser interface to help.Type 'q()' to quit R.> ## 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)+ }Loading required package: bootBOOTSTRAP CONFIDENCE INTERVAL CALCULATIONSBased on 1000 bootstrap replicatesCALL :boot.ci(boot.out = cd4.boot, conf = 0.9, type = c("norm", "basic","perc"), h = atanh, hinv = tanh)Intervals :Level Normal Basic Percentile90% ( 0.4705, 0.8589 ) ( 0.4620, 0.8597 ) ( 0.4900, 0.8689 )Calculations on Transformed Scale; Intervals on Original Scale>>> proc.time()user system elapsed0.127 0.020 0.519