data(beets)
head(beets)
beet0<-lmer(sugpct~block+sow+harvest+(1|block:harvest), data=beets, REML=FALSE)
beet_no.harv <- update(beet0, .~.-harvest)
rr <- PBrefdist(beet0, beet_no.harv, nsim=20)
rr
## Note clearly many more than 10 simulations must be made in practice.
## Computations can be made in parallel using several processors:
## Not run:
# cl <- makeSOCKcluster(rep("localhost", 4))
# clusterEvalQ(cl, library(lme4))
# clusterSetupSPRNG(cl)
# rr <- PBrefdist(beet0, beet_no.harv, nsim=20)
# stopCluster(cl)
# ## End(Not run)
## Above, 4 cpu's are used and 5 simulations are made on each cpu.
Run the code above in your browser using DataCamp Workspace