The R Project SVN R

Rev

Rev 56621 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 56621 Rev 57296
Line 4... Line 4...
4 
 
 4 
 
5 
library(boot)
 5 
library(boot)
6 
run1 <- function(...) {
 6 
run1 <- function(...) {
7 
   library(boot)
 7 
   library(boot)
8 
   cd4.rg <- function(data, mle) MASS::mvrnorm(nrow(data), mle$m, mle$v)
 8 
   cd4.rg <- function(data, mle) MASS::mvrnorm(nrow(data), mle$m, mle$v)
9 
   cd4.mle <- list(m = mean(cd4), v = var(cd4))
 9 
   cd4.mle <- list(m = colMeans(cd4), v = var(cd4))
10 
   boot(cd4, corr, R = 500, sim = "parametric",
 10 
   boot(cd4, corr, R = 500, sim = "parametric",
11 
        ran.gen = cd4.rg, mle = cd4.mle)
 11 
        ran.gen = cd4.rg, mle = cd4.mle)
12 
}
 12 
}
13 
cl <- makeCluster(mc <- 2)
 13 
cl <- makeCluster(mc <- 2)
14 
clusterSetRNGStream(cl, 123)
 14 
clusterSetRNGStream(cl, 123)