The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.

hdbm is a Bayesian inference method that uses continuous shrinkage priors for high-dimensional mediation analysis, developed by Song et al (2018). hdbm provides estimates for the regression coefficients as well as the posterior inclusion probability for ranking mediators.

Install

You can install hdbm via CRAN

install.packages("hdbm")

Or devtools

devtools::install_github("umich-cphds/hdbm", build_opts = c())

If you wish to install the package via devtools, you will need a C++ compiler installed. This can be accomplished by installing Rtools on Windows and Xcode on MacOS.

Example

Taken from the hdbm help file

library(hdbm)

Y <- hdbm.data$y
A <- hdbm.data$a

# grab the mediators from the example data.frame
M <- as.matrix(hdbm.data[, paste0("m", 1:100)], nrow(hdbm.data))

# We just include the intercept term in this example.
C <- matrix(1, nrow(M), 1)
beta.m <- rep(0, 100)
alpha.a <- rep(0, 100)

set.seed(1245)
output <- hdbm(Y, A, M, C, C, beta.m, alpha.a, burnin = 3000, ndraws = 100)

# Which mediators are active?
active <- which(colSums(output$r1 * output$r3) > 50)
colnames(M)[active]

Reference

Yanyi Song, Xiang Zhou et al. Bayesian Shrinkage Estimation of High Dimensional Causal Mediation Effects in Omics Studies. bioRxiv 467399

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.