# NOT RUN {
# following example from https://jenfb.github.io/bkmr/overview.html
set.seed(111)
library(bkmr)
dat <- bkmr::SimData(n = 50, M = 4)
y <- dat$y
Z <- dat$Z
X <- dat$X
set.seed(111)
Sys.setenv(R_FUTURE_SUPPORTSMULTICORE_UNSTABLE="quiet")
future::plan(strategy = future::multiprocess)
# run 4 parallel Markov chains
fitkm.list <- kmbayes_parallel(nchains=4, y = y, Z = Z, X = X, iter = 5000,
verbose = FALSE, varsel = TRUE)
bigkm = comb_bkmrfits(fitkm.list)
ests = ExtractEsts(bigkm)
ExtractPIPs(bigkm)
pred.resp.univar <- PredictorResponseUnivar(fit = bigkm)
risks.overall <- OverallRiskSummaries(fit = bigkm, y = y, Z = Z, X = X,
qs = seq(0.25, 0.75, by = 0.05), q.fixed = 0.5, method = "exact")
# additional objects that are not in a standard bkmrfit object:
summary(bigkm$iters)
table(bigkm$chain)
# }
# NOT RUN {
closeAllConnections()
# }
Run the code above in your browser using DataCamp Workspace