bkmr_wrapper: Wrapper function to fit bkmr model in package bkmr.
In lvhoskovec/mmpack: Methods for multipollutant mixtures analysis

Wrapper function to fit bkmr model in package bkmr.

1	bkmr_wrapper(niter, nburn, Y, X, W, varsel = FALSE, groups = NULL)

`niter`	number of iterations of sampler
`nburn`	number of burn-in iterations
`Y`	n x 1 vector of response data
`X`	n x p matrix of exposure data
`W`	n x q matrix of covariate data
`varsel`	logical; if TRUE variable selection is implemented, default is FALSE
`groups`	optional vector of length p indicating group membership to conduct hierarchical variable selection

a list with components

fit: the original model fit
risk: the posterior mean risk for each subject
pips: component-wise (if groups = NULL) or conditional (if groups is specified) posterior inclusion probabilities for exposures
group.pips: group posterior inclusion probabilities (only if groups is specified))
risk.summary: mean, SD, and 0.95 CI for risk for each subject
preds: predicted values for each subject based on exposures and covariates
risk.distn: posterior distribution of risk for each subject